*7月平均の地表風と気圧場を描く

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

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

define meanp = ave(p,time=00z1jul2004,time=12z31jul2004)
set cthick 5
set clopts 1 1 0.1
set cint 4
d meanp/100


define meanus = ave(us,time=00z1jul2004,time=12z31jul2004)
define meanvs = ave(vs,time=00z1jul2004,time=12z31jul2004)

set arrowhead 0.1
set arrscl 0.5 15
d skip(meanus,6,4);meanvs

draw title July monthly-mean pressure [hPa] & wind [m/s] at surface

enable print meanslp.gm
print
disable print


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

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