关于GrADS中插值问题oacres

just_a_bird

各位大侠，我在用oacres进行插值时出问题了，具体情况见图。另附本实验的数据（附件），希望能得到各位指点，不胜感激！

lqouc

数据没有用的，也没人会浪费金币下载，完全可以把你的ctl和gs贴上来。而且报错很明显了，就是你的空间网格点变量g有问题。

Aires

你生成的格点变量名检查一下是g吗

mofangbao

遇到这种情况都要问的就说明自己都没有把提示的英文翻译下，翻译过来你就知道要去gs里面找错误了，请帖代码，下次请将帖子发到编程作图求助专区

just_a_bird

mofangbao 发表于 2013-4-2 12:48

                               
登录/注册后可看大图

遇到这种情况都要问的就说明自己都没有把提示的英文翻译下，翻译过来你就知道要去gs里面找错误了，请帖代码 ...

原始数据的CTL：
DSET   D:\ENVISAT\hjj.grd
DTYPE  station  
STNMAP D:\ENVISAT\hjj.map
UNDEF  -32766
TITLE  Global Hs ENVISAT
TDEF 1 linear 01may2011 1mo
VARS 1
hs 0 99 SWH Data
ENDVARS
格网的CTL：
DSET  D:\ENVISAT\hjj0.grd
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 360 LINEAR 0 1.0
YDEF 181 LINEAR 0 1.0
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR 01may2011 1mo
VARS 1
g 0 99 Grid Data
ENDVARS
GS里面内容：
'reinit'
'open D:\ENVISAT\hjj.ctl'
'open D:\ENVISAT\hjj0.ctl'
'set lon -180 180'
'set lat -90 90'
'set mpdset hires'
'd oacres(g,hs.2)'
好像问题出在最后一句话。

lqouc

just_a_bird 发表于 2013-4-2 14:42

                               
登录/注册后可看大图

原始数据的CTL：
DSET   D:\ENVISAT\hjj.grd
DTYPE  station  

你打开的顺序对么........

Aires

g.2  hs.1

river

just_a_bird 发表于 2013-4-2 14:42

                               
登录/注册后可看大图

原始数据的CTL：
DSET   D:\ENVISAT\hjj.grd
DTYPE  station  

g才是你打开的第二文件里的变量啊
'd oacres(g.2,hs)'

just_a_bird

river 发表于 2013-4-2 19:27

                               
登录/注册后可看大图

g才是你打开的第二文件里的变量啊
'd oacres(g.2,hs)'

你说的对，已经能画出来图了，但是我还有有两个疑问，第一，上面的命令我是按照（http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903）做出来的。我是先按照里面重做一遍（能得到最终结果），然后把原始站点数据换成自己的站点数据；第二，我生成的图形是这样的，不知道怎么才能拟合在一起，而不是条带形式的。谢谢，各位大侠指点。

river

just_a_bird 发表于 2013-4-4 16:25

                               
登录/注册后可看大图

你说的对，已经能画出来图了，但是我还有有两个疑问，第一，上面的命令我是按照（http://bbs.06climate.c ...

您这数据是全球啊，怎么生成的二进制资料，把过程和原始资料都讲讲吧