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

[Perry, Aaron]

Thanks for the code Jennifer! It returned a 425 for the grid size, which seems about right. Also, "refdtotal" is just a name I chose for reflectivity.


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 Jennifer M Adams <jadams21 at gmu.edu>
Sent: Wednesday, November 1, 2017 1:56:31 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Number of Full Domain Grid Values Available in HRRR Domain


> On Nov 1, 2017, at 12:37 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:
> 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.
Jeff,
Why would you begin your answer with such an irrelevant, vague, and negative statement? It taints everything else you might have to say. If you don’t like GrADS, unsubscribe or lurk quietly, but don’t answer questions with mean-spirited misinformation.


Aaron,
> I am trying to figure out how to output the exact number of grid values across my domain in GrADS using the HRRR model.
> 'set lat 37 48.5'
> 'set lon 276.0 295.0’

At this point you can get the number of X grid points like this:
‘set gxout stat'
‘d refdtotal’
* The following assumes the first line of the stat output starts with “Notice: Automatic Grid Interpolation…"
* If you are not using a PDEF , then the idims will be in line 3 of the stat output
idims=sublin(result,4)
x1=subwrd(idims,4); x2=subwrd(idims,6)
xsize=x2-x1+1

>
> *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?

Now you can use the ‘xsize’ script variable in your ‘set prnopts’ statement. Your format statment %0.0f needs to be adjusted — the number to the left of the decimal indicates how many total places (columns) your number will occupy. The number to the right of the decimal indicates how many digits will appear to the right of the decimal in your formatted number.  I suggest something like this:
   ‘set prnopts %5.1f ‘xsize’ 1 u’
But it depends a lot on what variable you are displaying. I don’t know what refdtotal is, so this is just a guess.
—Jennifer

> '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
> FB: fb.com/arnpry
> LinkedIn: arnpry
> Web: shearanalytics.net
> Email: Aaron.Perry at lyndonstate.edu
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
>
> --
> Jeff Duda
> Post-doctoral research fellow
> University of Oklahoma School of Meteorology
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer Miletta Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
George Mason University




_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171102/1422b669/attachment.html