- 积分
 - 336
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 2015-4-2
 
	- 最后登录
 - 1970-1-1
 
 
 
 
 
 
 | 
	
 
GrADS
| 系统平台: | 
  | 
 
| 问题截图: | 
- | 
 
| 问题概况: | 
原来用fortran画的,但好像只能输出一个层次,否则输出的文件是0k,后面想检验一下前面的一层次,顺便把其它层次的季节平均用grads计算一下。 
但是输出来的图都是一样的233333  | 
 
| 我看过提问的智慧: | 
看过  | 
 
| 自己思考时长(天): | 
1  | 
 
 
 
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册 
 
 
 
x
 
利用原来的nc,打算计算每个季节每个高度的平均值 
'reinit' 
'set gxout fwrite' 
'set fwrite i:\data\wind\uvvector\seasontry.grd' 
'sdfopen i:\data\uwnd.mon.mean.nc' 
itt=1 
while(itt<=56) 
'set x 1 144' 
'set y 1 73' 
iz=1 
while(iz<=12) 
'set z 'iz'' 
'set t 'itt'' 
'define tt='itt'+12*('itt'-1)' 
it=0 
while(it<=9) 
'define u31=ave(uwnd,t=155+tt+'it',t=155+2+tt+'it')'   *有的人说需要t='155+tt+it',但我这样会报错 
*'define u32=ave(uwnd,t=155+3+tt,t=155+5+tt)' 
*'define u33=ave(uwnd,t=155+6+tt,t=155+8+tt)' 
*'define u34=ave(uwnd,t=155+9+tt,t=155+11+tt)' 
'd u31' 
it=it+3 
endwhile 
*'d u32' *原本是不一样的变量,但由于一样的图,我就怕是覆盖了,但变成一个变量还是一样。。。 
*'d u33' 
*'d u34' 
iz=iz+1 
endwhile 
itt=itt+1 
endwhile 
'disable fwrite' 
; 
 
用上面的gs输出,没有报错信息,但就是画出来的图每张都一样。 
 
 
 
 
 |   
 
 
 
 |