<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>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.</p>
<p>Jim<span style="color:rgb(80,0,80)"><br>
</span></p>
<div class="moz-cite-prefix">On 5/8/21 12:52 PM, Evgeny Gavrilov
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAKh9uO8WABZANOeSDz9zCAbteJ5VYA0q2YHDHYtwwv4U6pppNw@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr">this command: 'set t 1 81' before 'define' does not
solve this problem<br>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sun, May 9, 2021 at 1:11 AM
Jeff Duda <<a href="mailto:jeffduda319@gmail.com"
moz-do-not-send="true">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>What do you mean by "doesn't help"? We need specifics
to be able to fully help you.</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
1:53 PM Evgeny Gavrilov <<a
href="mailto:eugen420@gmail.com" target="_blank"
moz-do-not-send="true">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><span 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" moz-do-not-send="true">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 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" target="_blank"
moz-do-not-send="true">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" moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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"
moz-do-not-send="true">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"
moz-do-not-send="true">gradsusr@gradsusr.org</a><br>
<a
href="http://gradsusr.org/mailman/listinfo/gradsusr"
rel="noreferrer" target="_blank"
moz-do-not-send="true">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" moz-do-not-send="true">gradsusr@gradsusr.org</a><br>
<a
href="http://gradsusr.org/mailman/listinfo/gradsusr"
rel="noreferrer" target="_blank"
moz-do-not-send="true">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" moz-do-not-send="true">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank" moz-do-not-send="true">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</div>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank"
moz-do-not-send="true">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr"
rel="noreferrer" target="_blank" moz-do-not-send="true">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"
moz-do-not-send="true">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr"
rel="noreferrer" target="_blank" moz-do-not-send="true">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
</body>
</html>