17个站点EOF分解话特征向量的图出现问题

乌云珠

这是我的步骤
1 eof.f90 作17个站点降水的eof分解，
2 将特征向量整理成grads识别的格式 sta.grd
3 为sta.grd写ctl文件 sta.ctl
4 生成map文件
5 写格点文件grid.grd
6 为格点文件写ctl文件 grd.ctl
7 sta.gs 画特征向量图 sta.gs  

其中我的sta.ctl文件为
dset d:\my31\grid.grd
title grid data sample
undef -9.99e33
xdef 11 linear 90.0 2.5
ydef 7 linear 30.0 2.5
zdef 1 levels 1000
tdef 17 linear jan1951 1mon
vars 1
g 0 99 grid data
endvars


grd.ctl文件为
dset d:\my31\grid.grd
title grid data sample
undef -9.99e33
xdef 11 linear 90.0 2.5
ydef 7 linear 30.0 2.5
zdef 1 levels 1000
tdef 17 linear jan1951 1mon
vars 1
g 0 99 grid data
endvars


sta.gs文件为（选了一个模态画图）
'open d:\my31\grd.ctl'
'open d:\my31\sta.ctl'
'set lon 95 115'
'set lat 32.5 42.5'
'set t 1'
'set grads off'
'set mpdset cnworld'
'd oacres(g,tzxl.2)'
结过运行sta.gs时出现了上图的错误信息，实在找不出原因呀？我的站点数据是从160个站中选取的17各站（区域范围是32.5-41N,100-112.5E）谢谢帮忙呀！

river

姐，您不识别英文啊，首先插值函数报错说少于两个站点，您回头看看您第二个打开的sta.ctl里面有tzxl这个变量么。在后面那个错误告诉你在那个路径下面打开不了cnworld，您那个路径下有cnworld这个文件么。

站点绘图论坛里的帖子已经很多了，您可以在论坛里搜索看看学习一下吧，照着做一遍也就会了

葡萄儿

为什么不用surfer画图呢？简单很多的

ayzqs

站点太少，EOF意义不大吧

乌云珠

river 发表于 2013-5-7 07:21

                               
登录/注册后可看大图

姐，您不识别英文啊，首先插值函数报错说少于两个站点，您回头看看您第二个打开的sta.ctl里面有tzxl这个变量 ...

谢谢，是经纬度写反了

river

乌云珠 发表于 2013-5-8 09:54

                               
登录/注册后可看大图

谢谢，是经纬度写反了

对哦，最近经常有人出这个问题呢

zl.87625

站点数据做eof，请问数据应该怎么读入？是某一年，然后所有站的数据，然后再接来下一年？

还是某一个站，所有年份的数据，然后再下一个站？

w1215lin

zl.87625 发表于 2013-7-20 08:31

                               
登录/注册后可看大图

站点数据做eof，请问数据应该怎么读入？是某一年，然后所有站的数据，然后再接来下一年？

还是某一个站， ...

看看这个帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=13319

hxyj

请问能否提供一个站点ＥＯＦ的程序？

流水2012

对哦，还是没弄懂，Eof程序帖子多，但是要是能详细说下就好了