Function to add or substract the temperature effect to data in a tomodir

def main():
    """Function to add or substract the temperature effect to data in a tomodir
    """
    options = handle_options()

    # read in temperature and resistivity data
    tempdata = readin_temp(options.temp_file)
    magdata = readin_rho(options.filename,
                         options.rhofile,
                         aniso=options.aniso)
    # calculate corrected data
    mag_corr = calc_correction(temp=tempdata,
                               mag=magdata,
                               add=options.add,
                               T_std=options.T_std,
                               m=options.m,)
    # save data
    save_mag_to_file(mag_corr,
                     options.output,
                     options.rhofile)