- 积分
- 202
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-7-6
- 最后登录
- 1970-1-1
|
GrADS
| 系统平台: |
windows |
| 问题截图: |
- |
| 问题概况: |
nc文件用grads插值成1.0*1.0后，怎么求模型间的相关系数和标准差？我看贴后只会用re函数插值，但是不会生成新的nc文件，（不太会写循环，就用了最笨的方法求得26年春季平均）我是想把两个分辨率不同的模式的数据都插值成1.0*1.0后再求相关系数等，然后用这个数据来做泰勒图。。求大神指点下！！
'reinit'
'sdfopen C:\air.mon.mean.nc'
'set gxout fwrite'
'set fwrite E:\gcm\air.mon.mean.grd'
'set lon 70 120'
'set lat 0 40'
'set lev 1000'
'set gxout shaded'
have1=0
have2=0
have3=0
have4=0
have5=0
have6=0
have7=0
have8=0
have9=0
have10=0
have11=0
have12=0
have13=0
have14=0
have15=0
have16=0
have17=0
have18=0
have19=0
have20=0
have21=0
have22=0
have23=0
have24=0
have25=0
have26=0
'have1=ave(air,t=387,t=389)'
'have2=ave(air,t=399,t=401)'
'have3=ave(air,t=411,t=413)'
'have4=ave(air,t=423,t=425)'
'have5=ave(air,t=435,t=437)'
'have6=ave(air,t=447,t=449)'
'have7=ave(air,t=459,t=461)'
'have8=ave(air,t=471,t=473)'
'have9=ave(air,t=483,t=485)'
'have10=ave(air,t=495,t=497)'
'have11=ave(air,t=507,t=509)'
'have12=ave(air,t=519,t=521)'
'have13=ave(air,t=531,t=533)'
'have14=ave(air,t=543,t=545)'
'have15=ave(air,t=555,t=557)'
'have16=ave(air,t=567,t=569)'
'have17=ave(air,t=579,t=581)'
'have18=ave(air,t=591,t=593)'
'have19=ave(air,t=603,t=605)'
'have20=ave(air,t=615,t=617)'
'have21=ave(air,t=627,t=629)'
'have22=ave(air,t=639,t=641)'
'have23=ave(air,t=651,t=653)'
'have24=ave(air,t=663,t=665)'
'have25=ave(air,t=675,t=677)'
'have26=ave(air,t=687,t=689)'
'set lat 5 30'
'set lon 89 110'
'd re((have1+have2+have3+have4+have5+have6+have7+have8+have9+have10+have11+have12+have13+have14+have15+have16+have17+have18+have19+have20+have21+have22+have23+have24+have25+have26)/26,50,linear,70,1,40,linear,0,1,bs)'
'disable fwrite'
; |
| 我看过提问的智慧: |
看过 |
| 自己思考时长(天): |
2 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
|
|
窥视卡
雷达卡
发表于 2017-4-1 09:56:50
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
楼主
