GrADS in LINUX for WRF

Got one problem today while scripting in GrADS, in hand-on training. we are animating and printing streamline and rainfall information.

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

WRF training

Having training on January's weekend. Saturday and Sunday morning. WRF(Weather Research and Forecasting) and MM5 (Mesoscale Modelling).. Nice one. On weekend, because the trainers are from Malaysian Meteorological Dept, and they are working on weekdays.. anyhow, i'm looking forward for the end result of this training.. then i can start doing modelling! wow...can't wait (even though i feel lousy right now..aarrgghh!)

and currently doing training, accessing linux server system using windoes, xmanager. cool~