24小时热门版块排行榜    

查看: 257  |  回复: 2
当前主题已经存档。

既傻又憨

至尊木虫 (知名作家)

山野渔夫

优秀版主

[交流] 【求助】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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digua88

至尊木虫 (著名写手)

★ ★ ★
小木虫(金币+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',这个本身定义了,那你循环的目的呢?
2楼2009-11-09 14:29:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

既傻又憨

至尊木虫 (知名作家)

山野渔夫

优秀版主

引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 既傻又憨 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中,共勉 +8 eulota 2026-08-26 8/400 2026-08-27 12:34 by lice_1987
[基金申请] 基金未中,这种答复是模板吗? +4 zhaosm1982 2026-08-27 5/250 2026-08-27 12:17 by zhaosm1982
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 11:29 by xxniao123
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
信息提示
请填处理意见