the script is like this (shorter version):
********************************************************************************
'reinit'
'open test1.ctl'
tt=2
while (tt<74)
'c'
'set t'tt
'set gxout shaded'
'define rain = rainc+rainnc'
'set lev 1000'
'd rain'tt
'set lev 850'
'set gxout stream'
'd u;v'
'printim test'tt'.png white'
tt=tt+1
endwhile
******************************************************************
note that rain and rainnc data are available in level 1000 only.
and u;v data available in several atmospheric level (19 to be exact)
The output data works for t=2, but then the rain data became 'undefined' (for the rest of t)...
After a few time of figuring out how to solve the problem, we came with moving 'set lev 1000'
to above 'define rain=rainc+rainnc' ,
My theory is the 'define rain' data need also to be set to 'set lev 1000' for the 2nd looping, because the second execution bringing the value of level=850 (that being used for wind output) into the 'define rain' calculation..so, the defined 'rain' is now in value at 850 hPa, instead of level 1000...
so, we adjusted >>
********************************************************************************
'reinit'
'open test1.ctl'
tt=2
while (tt<74)
'c'
'set t'tt
'set lev 1000'
'set gxout shaded'
'define rain = rainc+rainnc'
!****!
'd rain'tt
'set lev 850'
'set gxout stream'
'd u;v'
'printim test'tt'.png white'
tt=tt+1
endwhile
******************************************************************
me: GrADS help is so limited especially in simple-easy problem, i'll try to write something more complicated whenever i'm in intermediate-advanced level..hihihi
