- 积分
- 187
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-9
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2013-3-21 23:18:51
|
显示全部楼层
这里感谢river.下面是我的gs文档。gs 文档中有问题大家一起来交流
#计算5月到九月的时间-纬度的水汽通量散度###
'reinit'
'sdfopen f:/data/201205-09/uwnd.2012.nc'
'sdfopen f:/data/201205-09/vwnd.2012.nc'
'sdfopen f:/data/201205-09/air.2012.nc'
'sdfopen f:/data/201205-09/pres.sfc.2012.nc'
'sdfopen f:/data/201205-09/shum.2012.nc'
'set lon 50 145'
'set lat -30 70'
'set t 122 274'
#当初做不出图来是因为时间放错了位置,现在我个人觉得现在的图还让自己满意#
'set map 1 1 20'
'set annot 1 50'
'set font 4'
'set mpdset cnworld'
'define usq=vint(pres.4(z=1),shum.5*uwnd.1,300)'
'define vsq=vint(pres.4(z=1),shum.5*vwnd.2,300)'
'define qdiv=hdivg(usq,vsq)'
'set lon 1'
'define w1=ave(qdiv,lon=90,lon=120)'
'define w2=ave(uwnd.1,lon=90,lon=120)'
'define w3=ave(vwnd.2,lon=90,lon=120)'
'set lat -10 60'
'set grads off'
'set grid off'
'set xyrev on'
'set timelab off'
'set gxout shaded'
'd w1*1000'
#不知道这里*1000量级对不对?#
'cbarn'
'set xyrev on'
'set gxout vector'
'd w2;w3'
'printim f:\17.gif x4000 y3000 white'
|
|