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

Garfield

先将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

感动啊！

mofangbao

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

易小凯

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

iloveplmmok

学习下们嘿嘿

wode_q_x

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

灏海明丽

学习学习。。

QQtangxian

正在学习这个

quentin011099

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

zzfrain

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