- 积分
- 2231
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-18
- 最后登录
- 1970-1-1
|
GrADS
| 系统平台: |
|
| 问题截图: |
- |
| 问题概况: |
最近敝人要做中小尺度的天气诊断，我想计算某层的变温，主要内容如下，目的是算出24小时变温。ctl文件就是正常的FNL资料多时次的CTL，共56时次。
'open E:\work\date\all.ctl'
'enable print E:\work\date\bianwen850.gmf'
'set grads off'
'set mpdset hires'
'set mpdset cnworld cnriver'
'set lon 106 125'
'set lat 14 35'
nt=5
while (nt<=56)
'set t 'nt''
'set lev 850'
'set cint 1'
'set gxout contour'
'define w=t-4'
bw=tmpprs(t)-tmpprs(w)
'd bw'
'q dim'
a=sublin(result,5)
a=subwrd(a,6)
'draw title BW'a
cbar
'print'
'clear'
nt=nt+1
endwhile
'reinit'
;
主要问题在于不知道如何实现两个时次的值做差，后来改动也试了tmpprs(t=nt)-tmpprs(t=w),但是依旧显示nt not d variable or function name。才疏学浅，还望各位大大不吝赐教,如何才能实现变温的计算。谢谢。 |
| 我看过提问的智慧: |
看过 |
| 自己思考时长(天): |
1 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
|
|
窥视卡
雷达卡
发表于 2015-1-28 15:36:37
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
楼主
谢谢楼主！我正在愁这个问题！