[gradsusr] Number of Full Domain Grid Values Available in HRRR Domain

[Perry, Aaron]

Hi Jeff,


The HRRR files are GRIB2 format and don't seem to be able to be opened by ncdump.


Aaron D. Perry
Class of 2016
Lyndon State College
Lyndonville, VT 05851
Mobile: (617) 780-4312
Twitter: @arnpry<https://twitter.com/arnpry>
FB: fb.com/arnpry<https://www.facebook.com/arnpry>
LinkedIn: arnpry<https://www.linkedin.com/in/arnpry><https://www.linkedin.com/in/upstream>
Web: shearanalytics.net<http://shearanalytics.net/>
Email: Aaron.Perry at lyndonstate.edu<http://Aaron.Perry@lyndonstate.edu/>
________________________________
From: gradsusr-bounces at gradsusr.org <gradsusr-bounces at gradsusr.org> on behalf of Jeff Duda <jeffduda319 at gmail.com>
Sent: Wednesday, November 1, 2017 12:37:46 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Number of Full Domain Grid Values Available in HRRR Domain

First note - I wouldn't prefer to use GrADS if I had other software available, as using GrADS isn't guaranteed to give a 100% accurate answer. If you have access to the NCO utilities (ncdump) and the raw HRRR output, I would just do an 'ncdump -h' on the output file.

If you really need to use GrADS, one trick you can try is this

'd asum(const(refdtotal,1),lon=276,lon=295,lat=37,lat=48.5)'

This function creates a new field such that all valid data points in the field refdtotal (which should match the HRRR domain) are set to 1, and then the sum of all those 1s is taken over the domain. Thus, it should spit out the number of grid points in the domain. You can use any field in the data set and any domain bounds so long as they fully encompass the valid data, as the asum() function won't consider missing data points (i.e., where the field is invalid outside of the model domain).

Jeff Duda

On Wed, Nov 1, 2017 at 11:12 AM, Perry, Aaron <Aaron.Perry at lsc.vsc.edu<mailto:Aaron.Perry at lsc.vsc.edu>> wrote:

I am trying to figure out how to output the exact number of grid values across my domain in GrADS using the HRRR model.


GrADS Code:


*Domain

'set lat 37 48.5'
'set lon 276.0 295.0'

*Generate a text file
'set gxout print'
'set prnopts %0.0f 219 1 u' <--- What should my numl value be? How do I determine to print out a value for every grid cell in the domain?
'd refdtotal'
dummy=write('/home/weather/hrrr.txt',result)


Best regards,

Aaron D. Perry
Class of 2016
Lyndon State College
Lyndonville, VT 05851
Mobile: (617) 780-4312
Twitter: @arnpry<https://twitter.com/arnpry>
FB: fb.com/arnpry<https://www.facebook.com/arnpry>
LinkedIn: arnpry<https://www.linkedin.com/in/arnpry><https://www.linkedin.com/in/upstream>
Web: shearanalytics.net<http://shearanalytics.net/>
Email: Aaron.Perry at lyndonstate.edu<http://Aaron.Perry@lyndonstate.edu/>

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr




--
Jeff Duda
Post-doctoral research fellow
University of Oklahoma School of Meteorology
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171102/cb5df1cd/attachment-0001.html