- 积分
- 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输出,没有报错信息,但就是画出来的图每张都一样。
|
|