爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4489|回复: 2

grads循环输出的数据层次不变

[复制链接]

新浪微博达人勋

发表于 2017-4-22 14:14:12 | 显示全部楼层 |阅读模式
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输出,没有报错信息,但就是画出来的图每张都一样。




密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-4-22 14:25:09 | 显示全部楼层
fortran大概抽搐好了。。我又用fortran算了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-1 21:11:07 | 显示全部楼层
set t 之后来个空格
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表