I am trying to create a L3 gridded daily AIRS ozone profile product from ~20 L2 chunks per day over Europe.
I am using coda.fetch to read in the variables, e.g.
latitude = coda.fetch(pf, ‘L2_Standard_atmospheric_surface_product/Geolocation_Fields/Latitude’)
latitude is a in (GeoTrack=45, GeoXTrack=30) dimensions, so I collapse them onto a time dimension so that I can assign it to a harp.Variable, as follows:
product.latitude = harp.Variable(latitude.flatten(), [‘time’])
Per day I have a a list of ‘harp.Product()’ instances which I then wish to bin_spatial and output into a netcdf file, as follows:
average = harp.execute_operations(productlist,
operations = ‘keep(latitude,longitude,tropopause_temp, tropopause_press, tropopause);
The final ncdf is not geo-referenced, as you can well imagine, since there are no latitude/longitude dimensions there, latitude/longitude are not kept in the operations. The lat/lon_bounds are output, of course. Is there a way to add the latitude/longitude dimension somehow in the post-operations? maybe I should not collapse the original geotrack dimensions onto “time”?
Sorry if this is a bit too confusing, I can add the whole python script is preferable.