def geo_correction(date_pair, dem_header, dem, incidence_angle_or_map): """ using geo coordinates to remove APS """ aps1 = pa.PyAPSPyRateGeo(os.path.join(ECMWF_DIR, ECMWF_PRE + date_pair[0] + ECMWF_EXT), dem_header=dem_header, dem=dem, grib=ECMWF, verb=True) aps2 = pa.PyAPSPyRateGeo(os.path.join(ECMWF_DIR, ECMWF_PRE + date_pair[1] + ECMWF_EXT), dem_header=dem_header, dem=dem, grib=ECMWF, verb=True) phs1 = np.zeros((aps1.ny, aps1.nx)) phs2 = np.zeros((aps2.ny, aps2.nx)) print('Without Lat Lon files') aps1.getdelay_pyrate(phs1, dem, inc=incidence_angle_or_map) aps2.getdelay_pyrate(phs2, dem, inc=incidence_angle_or_map) aps_delay = phs2 - phs1 # delay in meters as we don't provide wavelength return aps_delay