- 积分
- 6723
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-3
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 wuyejincao 于 2012-11-10 21:32 编辑
前5套资料时次是1981~2010年(1~30个时次),最后一套资料是1993~2008年(16个时次),出图后发现,就最后一套资料,就是红线不对,它对应的应该是中间一段时间~ 请高人指导!
'reinit'
'open D:\model\Domain\CFSR_springchina_ave.ctl'
'open D:\model\Domain\Inspringchina_ave.ctl'
'open D:\model\Domain\MERRAspringchina_ave.ctl'
'open D:\model\Domain\JRA_spring_china_30yr_ave.ctl'
'open D:\model\Domain\NCEP_china_spring_ave.ctl'
'open D:\model\Domain\spring_sm_16yr_marked_ave.ctl'
'set grid off'
'set grads off'
*'set xaxis 1981 2010 2'
'enable print D:\model\Domain\domain_china.gmf'
'set cstyle 1'
'set cmark 3'
'set ccolor 1'
'set cthick 6'
'define mean=ave(CFSR_spring,t=1,t=30)'
'define sd=sqrt((sum(pow(CFSR_spring-mean,2),t=1,t=30))/29)'
'set t 1 30'
'd (CFSR_spring-mean)/sd'
'cbar_line -x 3.5 -y 0.3 -c 1 -m 3 -l 1 -t "CFSR"'
'set cstyle 1'
'set cmark 2'
'set ccolor 1'
'set cthick 6'
'define mean2=ave(ERA_spring.2,t=1,t=30)'
'define sd2=sqrt((sum(pow(ERA_spring.2-mean2,2),t=1,t=30))/29)'
'set t 1 30'
'd (ERA_spring.2-mean2)/sd2'
'cbar_line -x 5 -y 0.3 -c 1 -m 3 -l 1 -t "ERA"'
'set cstyle 1'
'set cmark 5'
'set ccolor 1'
'set cthick 6'
'define mean3=ave(MERRA_spring.3,t=1,t=30)'
'define sd3=sqrt((sum(pow(MERRA_spring.3-mean3,2),t=1,t=30))/29)'
'set t 1 30'
'd (MERRA_spring.3-mean3)/sd3'
'set cstyle 1'
'set cmark 4'
'set ccolor 1'
'set cthick 6'
'define mean4=ave(JRA_spring.4,t=1,t=30)'
'define sd4=sqrt((sum(pow(JRA_spring.4-mean4,2),t=1,t=30))/29)'
'set t 1 30'
'd (JRA_spring.4-mean4)/sd4'
'set cstyle 1'
'set cmark 8'
'set ccolor 1'
'set cthick 6'
'define mean5=ave(NCEP_spring.5,t=1,t=30)'
'define sd5=sqrt((sum(pow(NCEP_spring.5-mean5,2),t=1,t=30))/29)'
'set t 1 30'
'd (NCEP_spring.5-mean5)/sd5'
'set cstyle 1'
'set cmark 6'
'set ccolor 2'
'set cthick 6'
'define mean6=ave(OBS_spring.6,t=1,t=16)'
'define sd6=sqrt((sum(pow(OBS_spring.6-mean6,2),t=1,t=16))/15)'
'set t 1 16'
'd (OBS_spring.6-mean6)/sd6'
'print'
'disable print'
;
|
|