search

CEMPD / VERDI

This is the repo for the VERDI project, written in java.
GNU General Public License v3.0
17 stars 13 forks source link

Need to improve Obs Overlay example #245

Open lizadams opened 3 years ago

lizadams commented 3 years ago

Current VERDI documentation example for Observational Overlays is available here: https://github.com/CEMPD/VERDI/blob/master/doc/User_Manual/VERDI_ch11.md#observational-data-overlays

To load the model dataset ./verdi.sh -f $cwd/data/model/CCTM46_P16.baseO2a.36k.O3MAX -s "O3[1]" -g tile -f $cwd/data/obs/AQS_overlay_2002.2hr.2002152.ncf

The model dataset that we are currently using for this example is 36km. It would be better to have a 4km dataset to compare with the obs, so that we can ensure that the observational data is overlaying correctly on the grid cells and on the underlying map.

The observational data needs to be loaded interactively, as there isn't a way to load obs data using a command-line script. It would be great if we could create a command-line script that would load the observational data overlay to the tile plot:

./verdi.sh -f $cwd/data/model/CCTM46_P16.baseO2a.36k.O3MAX -s "O3[1]" -f $cwd/obs/AQS_overlay_2002.2hr.2002152.ncf -obs "O3[2]"

There are also only 2 hours in the observational dataset that we provide AQS_overlay_2002.2hr.2002152.ncf. It would be better to provide 24 hours, so that the following error doesn't occur.

Screen Shot 2021-04-01 at 5 49 11 PM
dschwede commented 3 years ago

I put a file containing 2 days of ozone data for a 4 km domain over Denver on atmos in /work/MOD3DEV/sgf/data/verdi along with the AQS overlay for all of 2014. When I plotted the overlay on the tile plot, I did not get any errors.

image