| 查看: 253 | 回复: 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' 这里有问题，怎么改一下就好了？ |
回复此楼
digua88
至尊木虫 (著名写手)
- EPI: 1
- 应助: 10 (幼儿园)
- 金币: 10912.5
- 散金: 1214
- 红花: 11
- 帖子: 1646
- 在线: 1961小时
- 虫号: 819370
- 注册: 2009-07-31
- 性别: GG
- 专业: 数值预报与数值模拟
2楼2009-11-09 14:29:57
既傻又憨
至尊木虫 (知名作家)
山野渔夫
- 应助: 19 (小学生)
- 贵宾: 1.514
- 金币: 19784.1
- 散金: 1173
- 红花: 45
- 沙发: 21
- 帖子: 5312
- 在线: 1948.4小时
- 虫号: 697583
- 注册: 2009-02-06
- 性别: GG
- 专业: 环境污染化学
3楼2009-11-09 19:28:19
投票:
