<div dir="ltr"> still doesn't help 'set t 1 81' <div><br></div><div><span class="gmail-im" style="color:rgb(80,0,80)">'reinit'<br><br>'sdfopen <a href="https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z" target="_blank">https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z</a>'<br><br>'set mpdset hires'<br>'set lon 10 30'<br>'set lat 40 50'<br><br>'set display color white'<br><br>'set gxout shaded'<br></span>'define tc = tmpprs(lev=850)'<br>'define rh = rhprs(lev=850)'<br>'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))'<br>'define vapr850mb = 6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'<br>'define e850mb = vapr850mb*1.001+(850-100)/900*0.0034'<br>'define w850mb = 0.62197*(e850mb/(850-e850mb))'<br>'define te850mb = (tc+(2260000*w850mb/1004))'<br>'define thetae = (te850mb*pow((1000/850),(287/1004)))-273.16'<br><br>'set t 1 81'<br>'define a = thetae'<span class="gmail-im" style="color:rgb(80,0,80)"><br>'d a-a(t-1)'<br>'cbarm'<br><br>'printim C:\Users\for\temp.png x800 y600'<br><br>'clear'<br><br>return</span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 8, 2021 at 11:15 PM Evgeny Gavrilov <<a href="mailto:eugen420@gmail.com">eugen420@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">still doesn't help 'set t 1 81' <div><br></div><div>'reinit'<br><br>'sdfopen <a href="https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z" target="_blank">https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z</a>'<br><br>'set mpdset hires'<br>'set lon 10 30'<br>'set lat 40 50'<br><br>'set display color white'<br><br>'set gxout shaded'<br>'define tc = tmpprs(lev=850)'<br>'define rh = rhprs(lev=850)'<br>'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))'<br>'define vapr850mb = 6.112*exp((17.67*dewp850mb)/(dewp850mb+243.5))'<br>'define e850mb = vapr850mb*1.001+(850-100)/900*0.0034'<br>'define w850mb = 0.62197*(e850mb/(850-e850mb))'<br>'define te850mb = (tc+(2260000*w850mb/1004))'<br>'define thetae = (te850mb*pow((1000/850),(287/1004)))-273.16'<br><br>'set t 1 81'<br>'define a = thetae'<br>'d a-a(t-1)'<br>'cbarm'<br><br>'printim C:\Users\for\temp.png x800 y600'<br><br>'clear'<br><br>return<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 8, 2021 at 10:54 PM James T. Potemra <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<p>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'.</p>
<p>Jim<br>
</p>
<div>On 5/8/21 8:42 AM, Evgeny Gavrilov
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">"i" from while.
<div><br>
</div>
<div>maps = 20<br>
i = 1<br>
while ( i<maps )<br>
'set t ' i <br>
</div>
<div><br>
</div>
<div>2) a-a(t=1) also has a report about constant field
(value=0) </div>
<div><br>
</div>
<div>'reinit'<br>
<br>
'sdfopen <a href="https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z" target="_blank">https://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs20210508/gfs_0p25_12z</a>'<br>
<br>
'set mpdset hires'<br>
'set lon 10 30'<br>
'set lat 40 50'<br>
<br>
'set display color white'<br>
<br>
'set gxout shaded'<br>
'define a = tmp2m'<br>
'd a-a(t-1)'<br>
'cbarm'<br>
<br>
'printim C:\Users\for\temp.png x800 y600'<br>
<br>
'clear'<br>
<br>
return<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sat, May 8, 2021 at 10:11
PM Jeff Duda <<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div>1a) How is 'i' set?</div>
<div>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</div>
<div>i = 1</div>
<div>i2 = i-1</div>
<div>'d a(t='i') - a(t='i2')'</div>
<div><br>
</div>
<div>2) This command is equivalent and does not require a
script variable (provided your dimension environment is
set appropriately): 'd a - a(t-1)'</div>
<div><br>
</div>
<div>Jeff Duda<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sat, May 8, 2021 at
12:01 PM Evgeny Gavrilov <<a href="mailto:eugen420@gmail.com" target="_blank">eugen420@gmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">why this method does not working? <br>
<div><br>
</div>
<div>'set gxout shaded'<br>
'define a=tmp2m'<br>
'd a(t='i') - a(t='i-1')' <br>
</div>
<div><br>
</div>
<div>report: constant field. value = 0 </div>
</div>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote>
</div>
<br clear="all">
<br>
-- <br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><font size="2"><span>Jeff
Duda, Research Scientist</span></font></div>
<div dir="ltr"><font size="2"><span></span></font>University
of Colorado Boulder</div>
<div dir="ltr"><font size="2"><span></span></font>Cooperative
Institute for Research in Environmental
Sciences</div>
<div dir="ltr">NOAA/OAR/ESRL/Global Systems
Laboratory<br>
<font size="2"><span>
<span>Boulder, CO<br>
</span></span></font>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
</div>
</blockquote></div>
</blockquote></div>