docs: 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)