[gradsusr] (no subject)

[Jeff Duda]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20210508/299cfefb/attachment-0001.html>