*7月1日から31日まで(月平均場)の平均的な風と高度場を描く

reinit
open ../ctl/200407.ctl
set grads off

set lon 0 180
set lat 10 60
set time 1jul2004

*高度を250hPa気圧面に設定
set lev 250

*--高度場のための準備--

*13日から20日までの期間平均をとり、meanz に代入

define meanz = ave(z,time=00z1jul2004,time=12z31jul2004)

set cthick 5
set clopts 1 1 0.1
set cint 50
d meanz

*--風速場のための準備--

define meanu = ave(u,time=00z1jul2004,time=12z31jul2004)
define meanv = ave(v,time=00z1jul2004,time=12z31jul2004)

set arrowhead 0.1
set arrscl 0.5 40
d skip(meanu,6,4);meanv

draw title July monthly-mean height [m] & wind [m/s] at 250hPa

enable print meanUVZ250.gm
print
disable print


!gxps -c -i meanUVZ250.gm -o meanUVZ250.ps
!gxeps -2 -c -i meanUVZ250.gm -o meanUVZ250.eps

!convert -rotate 90 -geometry 600x400  meanUVZ250.eps meanUVZ250.gif