【已解决】不同资料时间序列图

wuyejincao

river 发表于 2012-11-7 20:12

                               
登录/注册后可看大图

我画了一张，美化啊什么的自己搞吧。
我就是多用一个1.gs计算输出资料，又变了一个1.ctl。

谢啦~~~

river

wuyejincao 发表于 2012-11-7 23:18

                               
登录/注册后可看大图

谢啦~~~

不客气。今天无意间看到grads里有一个close的命令，突然灵光一闪发现它可以比较简单的解决你的这个问题。只需要把时间序列较短的那个放在前面计算出来，然后close掉，再计算其它的画图（也可以计算完一起画图），最后加上这条线的display就出来了。
具体gs如下（就是在你的基础上稍微改了一下）：
'reinit'
'open f:\help\spring_sm_16yr_marked_ave.ctl'
'define mean6=ave(OBS_spring,t=1,t=16)'
'define sd6=sqrt(sum(pow(OBS_spring-mean6,2),t=1,t=16)/15.)'
'set t 1 16'
'define line6=(OBS_spring-mean6)/sd6'
'close 1'

'open f:\help\CFSR_springchina_ave.ctl'
'open f:\help\Inspringchina_ave.ctl'
'open f:\help\MERRAspringchina_ave.ctl'
'open f:\help\JRA_spring_china_30yr_ave.ctl'
'open f:\help\NCEP_china_spring_ave.ctl'

'set grid off'
'set grads off'
*'set vrange -3.5 3.5'
*'set yaxis -3.5 3.5 0.5'
'set xaxis 1981 2010 2'
'enable print f:\help\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)'
'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)'
'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)'
'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)'
'd (NCEP_spring.5-mean5)/sd5'



'set cstyle 1'
'set cmark 6'
'set ccolor 2'
'set cthick 6'
'd line6'

'define y=0'
'set cmark 0'
'set ccolor 22'
'set cstyle 3'
'd y'

'printim f:\help\domain_china.gif white'
'print'
'disable print'
;
出来的图和之前我画的那张是一模一样滴！（突然感觉这个close命令好棒啊，以后可以省掉好多提前提取资料的麻烦啦）

wuyejincao

river 发表于 2012-11-12 21:56

                               
登录/注册后可看大图

不客气。今天无意间看到grads里有一个close的命令，突然灵光一闪发现它可以比较简单的解决你的这个问题。 ...

好，原来这么亮！

chenggao

这是什么？看不太懂

river

chenggao 发表于 2012-11-13 10:48

                               
登录/注册后可看大图

这是什么？看不太懂

gs文件啊，画折线的。这有什么不懂的

subtropical

river 发表于 2012-11-12 21:56
不客气。今天无意间看到grads里有一个close的命令，突然灵光一闪发现它可以比较简单的解决你的这个问题。 ...

大神不仅学术好，人也很好，谢谢，我也在计算两个不同时间段的ctl画在同一张图上，学习了

river

subtropical 发表于 2016-1-19 11:22
大神不仅学术好，人也很好，谢谢，我也在计算两个不同时间段的ctl画在同一张图上，学习了

我都不记得我当年有这么大劲头帮人解决问题呢···

lasa123

bangbangde

huangyao

楼主，请问你这个问题解决了吗，我也遇到同样的

qxtlyf

好厉害,学习了