[gradsusr] (no subject)

Sat May 8 18:52:43 EDT 2021

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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20210509/93845d7b/attachment.html>