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

楼主 |
发表于 2012-10-20 18:49:22
|
显示全部楼层
river 发表于 2012-10-20 18:34
如果是NC资料的话是不应该出现这种问题的。你把gs贴上来吧。我怀疑是你计算的时候设置的经纬度比你要画的 ...
'reinit'
'sdfopen e:\data\ncep\uwnd.2012.nc'
'sdfopen e:\data\ncep\vwnd.2012.nc'
'sdfopen e:\data\ncep\shum.2012.nc'
'set grads off'
'set grid off'
'set mpdset cnworld cnriver'
'set map 15 1 2'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set parea 1.5 10 1 8'
nt=809
while(nt<=812)
'set t 'nt''
'set lat 20 45'
'set lon 100 135'
'set lev 850'
'set ylint 5'
'set xlint 5'
'set csmooth on'
*'run define_color.gs'
*'set ccols 4 11 5 13 3 10 12 8 2 27 68'
*'set clevs 2 4 6 8 10 12 14 16 18 20'
'set gxout shaded'
'd hdivg(shum.3*uwnd.1/9.8,shum.3*vwnd.2/9.8)'
'cbarn 0.85 0 5.75 0.35'
'q time'
rc=sublin(result,1)
tt=subwrd(rc,3)
'draw title 850hPa Divergence of Moisture Flux '
'printim e:\IAP\721\850hPa-divergence-of-moisture-flux'tt'.gif white'
'enable print e:\IAP\721\850hPa-divergence-of-moisture-flux'tt'.gmf'
'print'
'disable print'
'c'
'set grads off'
nt=nt+1
endwhile
;
经纬度就设置了一次,不太理解你所说的。。。。请指教 |
|