code: def main():
      
      options = handle_options()
  
      # read in observed and synthetic data
      elecs, d_obs = readin_volt(options.d_obs)
      elecs, d_est = readin_volt(options.d_est)
      elecs, d_estTC = readin_volt(options.d_estTC)
      # calculate corrected data
      volt_corr = calc_correction(d_obs,
                                  d_est,
                                  d_estTC,
                                  )
      # save data
      save_volt(elecs,
                volt_corr,
                options.output,
                )
docs: Function to remove temperature effect from field data