用GrADS绘制气候态候平均风速

千里草

       用NCEP日再分析资料，绘制气候态候平均风速，就是想用GrADS绘制10年（1994-2013年）5月第二候850hPA平均风速，gs文件如下，但画不出图（如果只打开一个文件，比如13年的U、V风速，就能画出图来。）。请大家指点，无论是在数据处理方面还是方法上面，非常感谢！

'reinit'


'sdfopen E:\zuoye\uwnd.2004.nc'
'sdfopen E:\zuoye\uwnd.2005.nc'
'sdfopen E:\zuoye\uwnd.2006.nc'
'sdfopen E:\zuoye\uwnd.2007.nc'
'sdfopen E:\zuoye\uwnd.2008.nc'
'sdfopen E:\zuoye\uwnd.2009.nc'
'sdfopen E:\zuoye\uwnd.2010.nc'
'sdfopen E:\zuoye\uwnd.2011.nc'
'sdfopen E:\zuoye\uwnd.2012.nc'
'sdfopen E:\zuoye\uwnd.2013.nc'
'sdfopen E:\zuoye\vwnd.2004.nc'
'sdfopen E:\zuoye\vwnd.2005.nc'
'sdfopen E:\zuoye\vwnd.2006.nc'
'sdfopen E:\zuoye\vwnd.2007.nc'
'sdfopen E:\zuoye\vwnd.2008.nc'
'sdfopen E:\zuoye\vwnd.2009.nc'
'sdfopen E:\zuoye\vwnd.2010.nc'
'sdfopen E:\zuoye\vwnd.2011.nc'
'sdfopen E:\zuoye\vwnd.2012.nc'
'sdfopen E:\zuoye\vwnd.2013.nc'


'set lon 30 150'
'set lat -10 45'
'set lev 850'
'set gxout shaded'
'set grads off'
'set grid off'
'set arrlab off'
'set arrscl 0.5 12'
'set cmin 8'
'define umean1=ave(uwnd,t=157,t=161)'
'define umean2=ave(uwnd.2,t=156,t=160)'
'define umean3=ave(uwnd.3,t=156,t=160)'
'define umean4=ave(uwnd.4,t=156,t=160)'
'define umean5=ave(uwnd.5,t=157,t=161)'
'define umean6=ave(uwnd.6,t=156,t=160)'
'define umean7=ave(uwnd.7,t=156,t=160)'
'define umean8=ave(uwnd.8,t=156,t=160)'
'define umean9=ave(uwnd.9,t=157,t=161)'
'define umean10=ave(uwnd.10,t=156,t=160)'
'define umean=(umean1+umean2+umean3+umean4+umean5+umean6+umean7+umean8+umean9+umean10)/10.0'
'define vmean1=ave(vwnd.11,t=157,t=161)'
'define vmean2=ave(vwnd.12,t=156,t=160)'
'define vmean3=ave(vwnd.13,t=156,t=160)'
'define vmean4=ave(vwnd.14,t=156,t=160)'
'define vmean5=ave(vwnd.15,t=157,t=161)'
'define vmean6=ave(vwnd.16,t=156,t=160)'
'define vmean7=ave(vwnd.17,t=156,t=160)'
'define vmean8=ave(vwnd.18,t=156,t=160)'
'define vmean9=ave(vwnd.19,t=157,t=161)'
'define vmean10=ave(vwnd.20,t=156,t=160)'
'define vmean=(vmean1+vmean2+vmean3+vmean4+vmean5+vmean6+vmean7+vmean8+vmean9+vmean10)/10.0'
'define fs=mag(umean,vmean)'
'd fs'
'set gxout vector'
'd umean;vmean'
'draw title 850hPa Winds May06-10,Climatology'
'draw line 5.5 1 6 1'
'draw line 5.95 0.95 6 1'
'draw line 5.95 1.05 6 1'
'draw string 5.75 0.8 10m/s'
'printim E:\zuoye\tu\Winds_May06_10.png x1000y800 white'
;

yx88229

用fortran调用处理数据吧

守望之月/;-D

不知道grads里最多可以打开几个文件？可能是打开文件太多了~~

同意LS，可以考虑先用FORTRAN计算以后再画图。

千里草

yx88229 发表于 2014-5-5 12:41
用fortran调用处理数据吧

谢谢，其实就是个作业，FORTRAN当然可以，就是每个文件都要处理觉得麻烦，但还是谢谢了。

千里草

守望之月/;-D 发表于 2014-5-5 12:47
不知道grads里最多可以打开几个文件？可能是打开文件太多了~~

同意LS，可以考虑先用FORTRAN计算以后再画 ...

恩，谢了！

琼小爷

可以先写个ctl把nc文件一次读取呀，就不用这么麻烦了

兰蓝岛

不知道问题解决了没

xuchen

推荐