deeplycloudy / glmtools

GOES-R Geostationary Lightning Mapper Tools
BSD 3-Clause "New" or "Revised" License
62 stars 33 forks source link

Infer domain and time from ABI file? #74

Open gerritholl opened 4 years ago

gerritholl commented 4 years ago

This is either a question or a feature proposal, depending on whether it's already possible. I imagine a rather common use case to be that the user wishes to aggregate GLM flashes over the domain and time coverage corresponding to the spatial coverage and interval between two ABI files. For example, consider:

OR_ABI-L1b-RadM1-M6C10_G16_s20201041800242_e20201041800311_c20201041800343.nc
OR_ABI-L1b-RadM1-M6C10_G16_s20201041801213_e20201041801282_c20201041801312.nc
OR_ABI-L1b-RadM1-M6C10_G16_s20201041802213_e20201041802282_c20201041802328.nc
OR_ABI-L1b-RadM1-M6C10_G16_s20201041803213_e20201041803282_c20201041803323.nc
OR_ABI-L1b-RadM1-M6C10_G16_s20201041804213_e20201041804282_c20201041804320.nc

then a user processing these alongside GLM L2 data may desire to:

Is either of that functionality currently possible with glmtools, or otherwise desirable?

deeplycloudy commented 4 years ago

Yes, it's possible to pass an arbitrary custom start and end time to the gridding system (--start and --end in make_glm_grids.py), and if you have the x,y fixed grid coords from ABI data then it should be straightforward to use those to drive the gridding process. make_glm_grids.py also already manages the generation of a conus and full disk grids that match ABI, as well as ABI mesoscale domains given a specified center point. When you say --domain=conus it generates 1D arrays of x,y fixed grid coordinates that are used as the target grid spec. If you already have explicit x,y from an ABI file those could be used directly, too. So, at first glance, this looks like a straightforward and welcome modification to make_glm_grids.py.