[gradsusr] [EXT] Re: (no subject)

Luthiene Alves Dalanhese de Almeida ludalanhese at usu.edu
Mon May 10 14:48:55 EDT 2021


Hi Evigeny

Try to set your temperature values before this looping.
‘t=t-273.15’ *from Kelvin to Celsius

If you do this, will be less confusing after, and you will get what Is wrong faster.

Sent from my iPhone

On May 10, 2021, at 12:31 PM, James Potemra <jimp at hawaii.edu> wrote:



Evgeny,

Again you are changing your script in between questions.  Please just fix one thing at a time.  You have repeated the confusion of the variable "t" in this line:

'define dewp850mb=(t-273.15)-((14.55+0.114*(t-273.15))*(1-0.01*rh)+pow((2.5+0.007*(t-273.15))*(1-0.01*rh),3)+(15.9+0.117*(t-273.15))*pow((1-0.01*rh),14))'

I presume by "t" you mean "temperature" and not "time".  This won't work since "t" at this point is not defined.

Second, you need to put all the code within the while loop.  For example, you have "printim" presumably to make images for each time-step, but since this is outside the while loop, you'll only get one.

Finally, there does not seem to be "rhprs" in this data set, or at least it's not being read reliably via OPeNDAP.

Comment out the lines defining rh and dewp850mb, move all the stuff below "endwhile" into the "while" loop, and it works, at least for me.

Jim

On 5/10/21 7:49 AM, Evgeny Gavrilov wrote:
but why? can you explain the reasons?

'reinit'

'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210510/gfs_0p25_12z'

'set lon 30 50'
'set lat 40 50'
'set t 1 80'

'define temp=tmpprs(lev=850)'
'define rh=rhprs(lev=850)'
'define dewp850mb=(t-273.15)-((14.55+0.114*(t-273.15))*(1-0.01*rh)+pow((2.5+0.007*(t-273.15))*(1-0.01*rh),3)+(15.9+0.117*(t-273.15))*pow((1-0.01*rh),14))'
'define vapr850mb=6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'
'define e850mb=vapr850mb*1.001+(850-100)/900*0.0034'
'define w850mb=0.62197*(e850mb/(850-e850mb))'
'define te850mb=(t+(2260000*w850mb/1004))'
'define te=(te850mb*pow((1000/850),(287/1004)))-273.15'

maps = 20
'set display color white'
i = 1
while ( i<maps )
  'set t ' i
  'set gxout shaded'
  'd te - te(t='i-1')'
  'cbarn'
  'clear'
  i=i+1
endwhile

'q dims'
times = sublin(result,5)
mse = subwrd(times,6)

'set strsiz 0.17'
'set string 1 r 6 0' ; 'draw string 4.85 7.98 GFS forecast: 'mse
'set string 1 r 6 0' ; 'draw string 10.1 7.98 Run: 'met

'printim C:\Users\123\for\te_'i'.png x800 y650'

'clear'

return

On Mon, May 10, 2021 at 8:14 PM Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:
Do you see it? That very first error message points to the fact that your data file either was not opened successfully or  you are not able to access the data within it.

First thing you should do is confirm you can plot basic fields from the data file.

Jeff

On Mon, May 10, 2021 at 10:05 AM Evgeny Gavrilov <eugen420 at gmail.com<mailto:eugen420 at gmail.com>> wrote:
<778.png>


On Mon, May 10, 2021 at 6:31 PM Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:
at least one of your variables or function calls above is wrong based on the error messages. However, because GrADS throws out more messages than the window you showed can contain, I can't see the first error message, which is likely the most important one. Show us that message.

Jeff Duda

On Mon, May 10, 2021 at 3:13 AM Evgeny Gavrilov <eugen420 at gmail.com<mailto:eugen420 at gmail.com>> wrote:
  the map is not drawn again

https://www.youtube.com/watch?v=s8JG6BUVh8c


On Mon, May 10, 2021 at 4:57 AM James Potemra <jimp at hawaii.edu<mailto:jimp at hawaii.edu>> wrote:

Evgeny,

You have a lot of curiosities in your script.  First, it's never a good to mix GrADS variables (e.g., "t") with your own:

'set t 1 20'
'define t=tmpprs(lev=850)'

instead of 'define t=' maybe make this 'define temp=' or something like that.  Second, you have

'define a = thetae'

but thetae is not defined anywhere.  Third, no need set time range every time you define the variables; set it once.  Finally, as Jeff suggested a while back, you'll need to add quotes to get the time index as a number:

'd a - a(t-1)' --> 'd a - a(t='i-1')

---------------------------

* set URL for data
'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210509/gfs_0p25_12z'

* define lat/lon/time ranges
'set lon 30 50'
'set lat 40 50'
'set t 1 81'

* define variables
'define temp=tmpprs(lev=850)'
'define rh=rhprs(lev=850)'
'define dewp850mb=(t-273.15)-((14.55+0.114*(t-273.15))*(1-0.01*rh)+pow((2.5+0.007*(t-273.15))*(1-0.01*rh),3)+(15.9+0.117*(t-273.15))*pow((1-0.01*rh),14))'
'define vapr850mb=6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'
'define e850mb=vapr850mb*1.001+(850-100)/900*0.0034'
'define w850mb=0.62197*(e850mb/(850-e850mb))'
'define te850mb=(t+(2260000*w850mb/1004))'
'define epa=(te850mb*pow((1000/850),(287/1004)))-273.16'


* draw maps
maps = 20
'set display color white'
i = 1
while ( i<maps )
  'set t ' i
  'set gxout shaded'
* I have not idea what you want to do, but this will display the difference between the current timestep and the one prior
  'd temp - temp(t='i-1')'
  'cbarn'
  'clear'
  i=i+1
endwhile

-------------------------


On 5/9/21 9:02 AM, Evgeny Gavrilov wrote:
All the same. The maps are not drawn

You can run my script:

'reinit'

'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210509/gfs_0p25_12z'

'set lon 30 50'
'set lat 40 50'

'set display color white'

maps = 20
  i = 1
  while ( i<maps )
'set t ' i

'set t 1 20'
'define t=tmpprs(lev=850)'
'set t 1 20'
'define rh=rhprs(lev=850)'
'set t 1 20'
'define dewp850mb=(t-273.15)-((14.55+0.114*(t-273.15))*(1-0.01*rh)+pow((2.5+0.007*(t-273.15))*(1-0.01*rh),3)+(15.9+0.117*(t-273.15))*pow((1-0.01*rh),14))'
'set t 1 20'
'define vapr850mb=6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'
'set t 1 20'
'define e850mb=vapr850mb*1.001+(850-100)/900*0.0034'
'set t 1 20'
'define w850mb=0.62197*(e850mb/(850-e850mb))'
'set t 1 20'
'define te850mb=(t+(2260000*w850mb/1004))'
'set t 1 20'
'define epa=(te850mb*pow((1000/850),(287/1004)))-273.16'

'set gxout shaded'
'define a = thetae'
'd a - a(t-1)'
'cbarm'

'printim C:\Users\123\for\temp_'i'.png x800 y600'

'clear'

i=i+1
endwhile

return

On Sun, May 9, 2021 at 8:54 PM Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:
Oh, yeah...I missed that. Good catch, James!

But the display command I suggested to him (d a - a(t-1)) certainly is valid for plotting the 1-time-step change in thetae, which I presume is what he wanted originally. It depends on the final graphic type he is trying to produce.

On Sun, May 9, 2021 at 10:41 AM James Potemra <jimp at hawaii.edu<mailto:jimp at hawaii.edu>> wrote:

You are not defining the time correctly.  In this new script you define "thetae", and this is still only defined for a single time-step.  You need to put the "set t 1 81" before any variable you want to be a function of time.  Also, your statement "d a-a(t-1)" doesn't make sense since "a" (if defined properly) will be a function of time/lat/lon.  If you want to plot 81 contour plots, each one as a different between current time and time before, you'll have to script that like you did in your first example.

Jim

On 5/8/21 12:52 PM, Evgeny Gavrilov wrote:
this command: 'set t 1 81' before 'define' does not solve this problem

On Sun, May 9, 2021 at 1:11 AM Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:
What do you mean by "doesn't help"? We need specifics to be able to fully help you.

Jeff Duda

On Sat, May 8, 2021 at 1:53 PM Evgeny Gavrilov <eugen420 at gmail.com<mailto:eugen420 at gmail.com>> wrote:
  still doesn't help 'set t 1 81'

'reinit'

'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z'

'set mpdset hires'
'set lon 10 30'
'set lat 40 50'

'set display color white'

'set gxout shaded'
'define tc = tmpprs(lev=850)'
'define rh = rhprs(lev=850)'
'define dewp850mb = (tc-273.15)-((14.55+0.114*(tc-273.15))*(1-0.01*rh)+pow((2.5+0.007*(tc-273.15))*(1-0.01*rh),3)+(15.9+0.117*(tc-273.15))*pow((1-0.01*rh),14))'
'define vapr850mb = 6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'
'define e850mb = vapr850mb*1.001+(850-100)/900*0.0034'
'define w850mb = 0.62197*(e850mb/(850-e850mb))'
'define te850mb = (tc+(2260000*w850mb/1004))'
'define thetae = (te850mb*pow((1000/850),(287/1004)))-273.16'

'set t 1 81'
'define a = thetae'
'd a-a(t-1)'
'cbarm'

'printim C:\Users\for\temp.png x800 y600'

'clear'

return

On Sat, May 8, 2021 at 11:15 PM Evgeny Gavrilov <eugen420 at gmail.com<mailto:eugen420 at gmail.com>> wrote:
still doesn't help 'set t 1 81'

'reinit'

'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z'

'set mpdset hires'
'set lon 10 30'
'set lat 40 50'

'set display color white'

'set gxout shaded'
'define tc = tmpprs(lev=850)'
'define rh = rhprs(lev=850)'
'define dewp850mb = (tc-273.15)-((14.55+0.114*(tc-273.15))*(1-0.01*rh)+pow((2.5+0.007*(tc-273.15))*(1-0.01*rh),3)+(15.9+0.117*(tc-273.15))*pow((1-0.01*rh),14))'
'define vapr850mb = 6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'
'define e850mb = vapr850mb*1.001+(850-100)/900*0.0034'
'define w850mb = 0.62197*(e850mb/(850-e850mb))'
'define te850mb = (tc+(2260000*w850mb/1004))'
'define thetae = (te850mb*pow((1000/850),(287/1004)))-273.16'

'set t 1 81'
'define a = thetae'
'd a-a(t-1)'
'cbarm'

'printim C:\Users\for\temp.png x800 y600'

'clear'

return

On Sat, May 8, 2021 at 10:54 PM James T. Potemra <jimp at hawaii.edu<mailto:jimp at hawaii.edu>> wrote:

You will have to set the time range before defining variable 'a', otherwise it will only be defined for a single time step.  Try add 'set t 1 81' before defining 'a'.

Jim

On 5/8/21 8:42 AM, Evgeny Gavrilov wrote:
"i" from while.

maps = 20
  i = 1
  while ( i<maps )
'set t ' i

2) a-a(t=1) also has a report about constant field (value=0)

'reinit'

'sdfopen https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z'

'set mpdset hires'
'set lon 10 30'
'set lat 40 50'

'set display color white'

'set gxout shaded'
'define a = tmp2m'
'd a-a(t-1)'
'cbarm'

'printim C:\Users\for\temp.png x800 y600'

'clear'

return

On Sat, May 8, 2021 at 10:11 PM Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:
1a) How is 'i' set?
1b) I don't think the grads scripting language supports math in a command line command like that. You'd have to do something like
i = 1
i2 = i-1
'd a(t='i') - a(t='i2')'

2) This command is equivalent and does not require a script variable (provided your dimension environment is set appropriately): 'd a - a(t-1)'

Jeff Duda

On Sat, May 8, 2021 at 12:01 PM Evgeny Gavrilov <eugen420 at gmail.com<mailto:eugen420 at gmail.com>> wrote:
why this method does not working?

'set gxout shaded'
'define a=tmp2m'
'd a(t='i') - a(t='i-1')'

report: constant field. value = 0
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr


--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr



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


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


--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr



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


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


--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO



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


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


--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr


--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO



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


CAUTION: This email originated from outside of USU. If this appears to be a USU employee, beware of impersonators. Do not click links, reply, download images, or open attachments unless you verify the sender’s identity and know the content is safe.


_______________________________________________
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/20210510/5a588f0e/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 778.png
Type: image/png
Size: 27757 bytes
Desc: 778.png
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20210510/5a588f0e/attachment-0001.png>


More information about the gradsusr mailing list