[gradsusr] (no subject)

Evgeny Gavrilov eugen420 at gmail.com
Sat May 8 15:47:10 EDT 2021


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


More information about the gradsusr mailing list