- 积分
 - 846
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2020-3-25
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
GrADS
| 系统平台: | 
GrADs  | 
 
| 问题截图: | 
- | 
 
| 问题概况: | 
运行的时候出现file number out of range  | 
 
| 我看过提问的智慧: | 
没看过(看完再来提问)  | 
 
| 自己思考时长(天): | 
1  | 
 
 
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
x
 
'reinit' 
'sdfopen E:\biyelunwen\ziliao\water_vapour.nc' 
'open E:\biyelunwen\ziliao2\topo10_2007.ctl' 
'enable print E:\biyelunwen\ziliao\water_vapour.gmf X3000 Y3000' 
 
'set annot 1 1' 
'set xlopts 1 4 0.12' 
'set ylopts 1 4 0.12' 
'set xlint 3' 
'set ylint 2' 
 
'set lon 50 120' 
'set lat 20 50'  
'set lev 0 0' 
 
 
 
 
********************************** 
********************************** 
i=1 
while(i<=12) 
'set vpage 2 9 0.5 8.2' 
*'set parea 1 8.5 5 8.5' 
'set grads off' 
'set grid off' 
 
'set t 1 12' 
'define aven=ave(p72.162.1,t=0,t=492,12)' 
'define avee=ave(p71.162.1,t=0,t=492,12)' 
 
'set t 'i'' 
'set map 9 1 3' 
'draw map' 
 
'set gxout contour' 
'set ccolor 4' 
'd aven;avee' 
 
 
*********绘制高原 
'set z 1' 
'set t 1' 
********** 
*'define h3000=maskout(tp.2,tp.2-3000)' 
*'define hconst=const(h3000,3000)' 
*'set gxout fgrid' 
*'set fgvals 3000 15' 
*'d hconst' 
*********** 
'set gxout contour' 
'set clab off' 
'set clevs 3000' 
'set ccolor 4' 
'd tp.2' 
 
 
 
 
 
*********绘制图标 
'set annot 1 1.5' 
'draw xlab longitude' 
'draw ylab latitude' 
'set strsiz 0.16 0.16' 
*'draw string 3.8 6.7 10m wind 'i' Month Height' 
'draw title 'i' Month Water Vapour Flux' 
'print' 
'c' 
i=i+1 
endwhile 
'disable print' 
; 
 
 
 |   
- 
运行gs就是这个样子了 
 
 
 
- 
 
 
 
 
 
 
 
 |