【求助】如何将两个变量分别插入格点场的横纵坐标

shangbn305

我想做站点数据和wrfout的风场对比。

首先wrfout的风场直接用 set gxout vector ; d skip(u,10,10);v; 画出来了。

对于站点数据，我使用最大风速及最大风速方向，使用fortran生成sta.grd文件，主要语句是u=big*sin(21.5*(bigd-1)*3.14/180.0)； v=big*cos(21.5*(bigd-1)*3.14/180.0)，其中bigd是风向，u、v是转换之后的。

站点的描述文件如下：其实这里有个问题，u的UNDEF=3.789826 v的UNDEF=1.810310，在sta.ctl里不知道怎么写两个

sta.ctl

DSET C:/GrADS20/win32/sta.grd
DTYPE station
STNMAP C:/GrADS20/win32/sta.map
UNDEF 3.789826
TITLE Wind Data Sample
TDEF 1 linear 00jun2012 1dy
VARS 2
u 0 99 Wind Data
v 0 99 Wind Data
ENDVARS

然后使用grads作图，使用如下语句：

'reinit'
'open C:\GrADS20\win32\sta.ctl'
'set t 1'
'set z 1'
'set lon 104.54355 148.08408996'
'set lat 32.28431 60.41944476'
'd u;v'#'d u'#'d v'
;
这是 d u;v的结果

d u

d v

原本以为d u;v自动生成和wrf一样的风场，也就是根号下(u平方+v平方)，但是我想错了。仔细想想应该是把u,v分别插值到grid格点中的横纵坐标轴，wrf的是已经插值了的，站点这个还未进行插值。

现在的问题是如何把u,b插值到横纵坐标中，大家又遇到这种问题的么？可否告诉我一下，万分感谢！

lqouc

用oacres这个函数吧，很简单的，但是要自己做一个格点场。

shangbn305

lqouc 发表于 2013-3-12 19:48

                               
登录/注册后可看大图

用oacres这个函数吧，很简单的，但是要自己做一个格点场。

格点场我做完了，下边是我的格点场的描述文件
DSET C:/GrADS20/win32/grid.grd
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 537 LINEAR 104.54355 0.08108108
YDEF 347 LINEAR 32.28431 0.08108108
ZDEF  27 LINEAR 1 1
TDEF  4 linear 00z01jun2012 480mn
VARS  1
g     0 99 Grid Data
ENDVARS

我使用如下语句进行插入操作的：
'reinit'
'open C:\GrADS20\win32\grid.ctl'
'open C:\GrADS20\win32\sta.ctl'
'set t 1'
'set z 1'
'define a=oacres(g,u.2)'
'define b=oacres(g,v.2)'
'd u;v'
;
出错了

lqouc

shangbn305 发表于 2013-3-13 12:51

                               
登录/注册后可看大图

格点场我做完了，下边是我的格点场的描述文件
DSET C:/GrADS20/win32/grid.grd
TITLE Grid Data Sample ...

.......这个，好奇怪啊，报错说你的站点都是缺省啊，格点场的时间最好和数据的时间长度匹配，话说你生成map文件了吧？

shangbn305

lqouc 发表于 2013-3-13 12:59

                               
登录/注册后可看大图

.......这个，好奇怪啊，报错说你的站点都是缺省啊，格点场的时间最好和数据的时间长度匹配，话说你生成m ...

生成map文件了，可能是我的思路不太对，我在论坛里搜搜帖子，看大家都是怎么做风场的