- 积分
- 1164
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-16
- 最后登录
- 1970-1-1

|

楼主 |
发表于 2012-8-20 16:10:17
|
显示全部楼层
做个霸气的木头 发表于 2012-8-20 12:35
grads存图有两种格式啦,一种是printim的,一种是enable print的,你两种都混着来的,这是基础知识,你选 ...
惭愧啊, ,我还有个问题想请教哈你,还望斑竹不吝赐教啊。
关于求64年(1948-2011)年平均水汽通量,我初步思路是:将求出的1-12月的一共12个月的月平均水气通量求和然后再除以12,自己改了一下脚本命令,但最后求出的图有问题,还得你帮我看看是哪儿出了问题??? 是红色标注的命令有问题吗???
'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen f:\data\monthly\shum.mon.mean.nc'
'sdfopen f:\data\monthly\uwnd.mon.mean.nc'
'sdfopen f:\data\monthly\vwnd.mon.mean.nc'
'sdfopen f:\data\monthly\pres.mon.mean.nc'
'set lat 4 53'
'set lon 73 135'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld hires'
'set map 1 1 1'
'draw map'
'set ylopts 1 4 0.12'
'set xlopts 1 4 0.12'
'set cthick 4'
'set t 1 768'
'set z 1 8'
'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'set z 1'
'set t 1 768'
'define uqs=vint(pres.4(z=1),uq,300)/1000'
'define vqs=vint(pres.4(z=1),vq,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 1 12'
'define w1=ave(sq,t+0,t=768,12)'
'define w2=ave(uqs,t+0,t=768,12)'
'define w3=ave(vqs,t+0,t=768,12)'
'define m1=sum(w1,t=1,t=12)'
'define m2=sum(w2,t=1,t=12)'
'define m3=sum(w3,t=1,t=12)'
'set gxout shaded'
'd m1'
'set grads off'
'set gxout contour'
'd m1'
'set grads off'
'set gxout vector'
'd m2;m3'
'run cbarn'
'draw title 1948-2011 year average vopor flux china'
'printim g:\year.tif x1000 y800 white'
;
|
|