求助用grads画经过EOF的160站降水图

Garfield · 发表于 2012-3-14 14:24:52

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

先将160站降水进行了EOF分析，得到一个egvt.txt文件，然后有根据站点的写入格式用Fortran写成了egvt.grd形式的资料，也生成了相应的map文件。然后描述文件如下：（有20个模态）
dset g:\shixi2\egvt.grd
dtype station
stnmap g:\shixi2\r160.map
undef -999.0
title rain
tdef 1 linear jun1951 1yr
vars 20
c1 0 99 rainfall data
c2 0 99 rainfall data
c3 0 99 rainfall data
c4 0 99 rainfall data
c5 0 99 rainfall data
c6 0 99 rainfall data
c7 0 99 rainfall data
c8 0 99 rainfall data
c9 0 99 rainfall data
c10 0 99 rainfall data
c11 0 99 rainfall data
c12 0 99 rainfall data
c13 0 99 rainfall data
c14 0 99 rainfall data
c15 0 99 rainfall data
c16 0 99 rainfall data
c17 0 99 rainfall data
c18 0 99 rainfall data
c19 0 99 rainfall data
c20 0 99 rainfall data
endvars

.gs如下：
'reinit'
'open g:\shixi2\map.ctl'
'open g:\shixi2\egvt.ctl'
'enable print g:\shixi2\egvt1.gmf'
'set mpdset lowres'
'set poli off'
'set lat 18 55'
'set lon 73 135'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'define aa=oacres(g.1,c1.2)'
'define rgg=maskout(aa,g.1-0.5)'
'define bb=smth9(rgg)'
'set csmooth on'
'set gxout contour'
'set cint 0.1'
'set cthick 3'
'd bb'
'print'
'disable print'
;

但是画不出图来，图上没有格点值，这是为什么？

jqh99999 · 发表于 2012-3-14 14:32:48

感动啊！

mofangbao · 发表于 2012-3-14 15:43:35

楼主没有自己先分析一下吗？

易小凯 · 发表于 2012-3-15 10:50:51

要写一个grid文件,然后再插值到格点上才行.
具体大可以参照<grads 站点作图详细解决方案>

iloveplmmok · 发表于 2012-3-19 16:24:57

学习下们嘿嘿

wode_q_x · 发表于 2012-3-17 18:23:22

楼主的ctl文件t应该设为20吧，变量是一个，给你我的ctl，gs文件，希望对你有所帮助。另外论坛里关于画站点图的资料可以看一下，大侠们整理的非常好！

灏海明丽 · 发表于 2012-4-28 11:14:16

学习学习。。

QQtangxian · 发表于 2012-5-18 16:29:42

正在学习这个

quentin011099 · 发表于 2012-7-16 09:25:24

请问楼主可以把你的egvt的txt转换成grd的格式的程序分享一下吗，谢谢了···

zzfrain · 发表于 2012-11-26 20:51:52

同求egvt的txt转换成grd的格式的程序，谢谢

		自动登录	找回密码
密码			立即注册

[求助] 求助用grads画经过EOF的160站降水图

登录后查看更多精彩内容~