分享个作站点降水的EOF程序，另附160站点资料（在国家气候中心也可下载）

我是小猪

最近再做毕业论文，顺便分享一下，该程序可分解站点降水（我在气候中心下载的降水资料）EOF主要的四个模态，分解出的vec1-vec4再用站点转格点的方式在grads上画图（详情请搜索）

我是小猪

洪hong 发表于 2016-5-7 23:09
您好，您的模态图画得非常好看，全中国都有阴影，而我插值的图青藏高原西侧就没有阴影，不知道您能不能发画 ...

图表主要参考清风版主的“GrADS站点文件作图详细解决方案”一贴，详情可搜索，具体gs如下：
'reinit'
'open C:\EOF\cREOF\grid.ctl'
'open C:\EOF\cREOF\sta.ctl'
'set lon 70 140'
'set lat 15 55'
'set grads off'
'set mpdset hires'
'define rgrid=oacres(g,rain.2)'
'set gxout shaded'
'd rgrid'
'set gxout contour'
'cnbasemap rgrid'
*加入色标
'cbarn'
'printim C:\EOF\cREOF\vec3.png white'
;

这是综合两个CTL的gs，但是具体的站点ctl即sta.ctl和格点ctl即grid.ctl还要参考上述帖子，其中，sta.ctl如下：
dset C:\EOF\cREOF\vec3.grd
dtype station
stnmap C:\EOF\cREOF\vec3.map
undef -32766
title rain data sample
tdef 1 linear jul1951 1yr
vars 1
rain 0 99 rain data
endvars

grid.ctl如下：
dset C:\EOF\cREOF\grid.grd
title grid data sample
undef -9.99E33
xdef 29 linear 70.0 2.5
ydef 17 linear 15.0 2.5
zdef 1 levels 1000
tdef 1 linear jul1951 1yr
vars 1
g 0 99 grid data
endvars
，但是要生成grid.grd需要先写入fFortran生成

星期六小妹妹

求分享图作为参考

我是小猪

传上去了，时间序列的红线可忽略

我是小猪

传上去了，时间序列的红线可忽略

贫道敬孔

多谢多谢{:5_267:}

洪hong

您好，您的模态图画得非常好看，全中国都有阴影，而我插值的图青藏高原西侧就没有阴影，不知道您能不能发画这个模态的grid.grd和gs给我看看？我也可以发我的grid.grd和gs给您看。

我是小猪

星期六小妹妹 发表于 2016-5-7 16:40
求分享图作为参考

谢谢提醒，忘记了，嘿嘿

洪hong

如果您需要看我的gs，我也可以共享给您哦。

洪hong

我傻傻的gs如下：
'reinit'
'open D:****.ctl'
'open D:****.ctl'
'enable print D:****.gmf'
'set lon 75 135.5'
'set lat 17 55'
'set mpdset cn cnriver'
'set mpdset cnworld'
'set mpdset cnscs'
i=1
while(i<=3)
'define a1=oacres(g,r'i'.2)'
'define b1=maskout(a1,g-0.5)'
'define bb1=smth9(b1)'
'set xlopts 1 10 0.18'
'set ylopts 1 10 0.18'
'set gxout shaded'
'cnbasemap bb1'
'set gxout contour'
'cnbasemap bb1'
'draw title r '
'cbarn.gs'
'set cthick 8'
'set clopts 1 6 0.1'
'print'
i=i+1
'c'
endwhile
'disable print'
;