GREMLIN CONUS3 Dataset for 2022

Geostationary Operational Environmental Satellite (GOES) Radar Estimation via Machine Learning to Inform NWP (GREMLIN) is a machine learning model that produces composite radar reflectivity using data from the Advanced Baseline Imager (ABI) and Geostationary Lightning Mapper (GLM). GREMLIN is useful for observing severe weather and providing information during convective initialization especially over regions without ground-based radars. Previous research found good skill compared to ground-based radar products, however, the analysis was done over a dataset with similar climatic and precipitation characteristics as the training dataset: warm season Eastern CONUS in 2019. This study expands the analysis to the entire contiguous United States, during all seasons, and covering the period 2020-2022. Several validation metrics including root-mean-square difference (RMSD), probability of detection (POD), and false alarm ratio (FAR) are plotted over CONUS by season, day-of-year, and time-of-da..., The methodology is described in detail by Hilburn et al. (2021). The ABI, GLM, and MRMS data sets were resampled to a common 3 km grid. A cloud height of 10 km was used for removing parallax displacements. Satellite and radar samples were matched in time with a maximum time difference of 2.5 minutes. GLM lightning groups were accumulated over 15-minute time periods., Code for reading the data: The data can be read using the provided Python code (\"read_conus3_file.py\" and \"test_read.py\") or the provided Fortran code (\"gzmodule.f90\", \"test_read.f90\", \"compile.sh\"). Note that this code reads the data without having to unzip the files.
The “test_read.py” and “test_read.f90” use one sample file (ABI_C13_202001010000.bin.gz) to verify the code is reading correctly by checking against seven points within the image and against the minimum and maximum over the full image. This file has been included in the software.tar package.
To use the Python code, at the command line simply invoke:
 python test_read.py
The function that actually reads a data file, in “read_conus3_file.py”, makes use of the gzip module, which is part of the Python Standard Library. If your code is reading correctly, you should see this output:
testfile= ABI_C13_202001010000.bin.gzmin data, expected, isclose= -1e+30 -1e+30 Truemax data, expected, isclose= 296.77463 296.77463 Truedata[   0,...