[gradsusr] (no subject)

Jeff Duda jeffduda319 at gmail.com
Mon May 10 10:25:41 EDT 2021


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> 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> 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> 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> 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> 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>
>>>>> 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>
>>>>>> 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>
>>>>>>> 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>
>>>>>>>> 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> 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
>>>>>>>>>> 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
>>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> gradsusr mailing listgradsusr at gradsusr.orghttp://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>> gradsusr mailing list
>>>>>> 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
>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>
>>>>
>>>> _______________________________________________
>>>> gradsusr mailing listgradsusr at gradsusr.orghttp://gradsusr.org/mailman/listinfo/gradsusr
>>>>
>>>> _______________________________________________
>>>> gradsusr mailing list
>>>> 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 listgradsusr at gradsusr.orghttp://gradsusr.org/mailman/listinfo/gradsusr
>>
>> _______________________________________________
> gradsusr mailing list
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20210510/95a0f0d5/attachment-0001.html>


More information about the gradsusr mailing list