24小时热门版块排行榜

返回列表

当前主题已经存档。

既傻又憨

至尊木虫 (知名作家)

山野渔夫

应助: 19 (小学生)
贵宾: 1.514
金币: 19784.1
散金: 1173
红花: 45
沙发: 21
帖子: 5312
在线: 1948.4小时
虫号: 697583
注册: 2009-02-06
性别: GG
专业: 环境污染化学

[交流] 【求助】grads用的好的帮看一下这个程序

reinit
yrstrt=1985
yrend=1996
root_nc='e:\ncdata\hgt\'
fn='hgt'
yy=yrstrt
while(yy<=yrend)
'sdfopen '%root_nc%fn%'.'%yy%'.nc'
say 'sdfopen '%root_nc%fn%'.'%yy'.nc'
'set lev 500'
'set lat 40 80'
'set t 10 33'
'set lon 75'
'define hgt1=ave(hgt/10,lon=75,lon=85)'
'define hgtpj=hgt1'
n=1986
while(n<=1996)
'define hgtpj=hgtpj+hgt1'
n=n+1
endwhile
yy=yy+1
endwhile
'define hgtave=hgtpj/12'
'set gxout shaded'
'set grid off'
'set grads off'
'set cint 4'
'set ylint 10'
'set xlint 4'
'set xyrev on'
'd hgtave'
'cbar 0.6 0'
'draw string 5.8 0.2 (a)'
为什么只能读出1985年的位势高度呢？不能求出85-96年的平均值呢，我觉得可能是在
n=1986
while(n<=1996)
'define hgtpj=hgtpj+hgt1'
这里有问题，怎么改一下就好了？

回复此楼

» 猜你喜欢

hopeeverythinggoeswell. http://db.tt/WTNw1wGM

1楼 2009-11-07 15:19:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

digua88

至尊木虫 (著名写手)

EPI: 1
应助: 10 (幼儿园)
金币: 10883
散金: 1214
红花: 11
帖子: 1646
在线: 1945.9小时
虫号: 819370
注册: 2009-07-31
性别: GG
专业: 数值预报与数值模拟

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xinmeng8803(金币+2,VIP+0):谢谢！ 11-10 08:44

看了一下，貌似求个平均每你这么麻烦，
n=1986
while(n<=1996)
'define hgtpj=hgtpj+hgt1'
n=n+1
endwhile
这个循环中套循环不妥，
'set t 10 33'，这个本身定义了，那你循环的目的呢？

赞一下(11人)

回复此楼

2楼2009-11-09 14:29:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

既傻又憨

至尊木虫 (知名作家)

山野渔夫

应助: 19 (小学生)
贵宾: 1.514
金币: 19784.1
散金: 1173
红花: 45
沙发: 21
帖子: 5312
在线: 1948.4小时
虫号: 697583
注册: 2009-02-06
性别: GG
专业: 环境污染化学

引用回帖:

Originally posted by digua88 at 2009-11-9 14:29:
看了一下，貌似求个平均每你这么麻烦，
n=1986
while(n<=1996)
'define hgtpj=hgtpj+hgt1'
n=n+1
endwhile
这个循环中套循环不妥，
'set t 10 33'，这个本身定义了，那你循环的目的呢？

我想求的是1985-1996年t 10 33这个时间段的逐日平均。我用的资料是逐年的nc文件。

赞一下

回复此楼

hopeeverythinggoeswell. http://db.tt/WTNw1wGM

3楼2009-11-09 19:28:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主既傻又憨的主题更新

返回列表