任意两点间斜剖面做法总结【利用GRADS自带函数】

红儿老子

请问楼主，有1x1的地形文件吗？

river

红儿老子 发表于 2014-3-1 13:42

                               
登录/注册后可看大图

请问楼主，有1x1的地形文件吗？

没有，不过你可以利用这个2.5*2.5的资料插值成1*1的

sun_shine_Xia

river 发表于 2014-3-1 21:25

                               
登录/注册后可看大图

没有，不过你可以利用这个2.5*2.5的资料插值成1*1的

LZ，我也用了你的这个方法去试着画了下，用的脚本就是你贴出来的画流场的那个，为什么他一直提示这个错误啊~~

脚本是下面这个：
'reinit'
'open F:\Data\EC_Data\W\wNDJ_Dtrend.ctl'
'open F:\Data\EC_Data\U\uNDJ_Dtrend.ctl'
'open F:\Data\EC_Data\V\vNDJ_Dtrend.ctl'
'open F:\Data\Topography\orog.ctl'
'set grads off'
'set grid off'
'set zlog on'
lon1=70.0
lon2=130.0
lat1=35.0
lat2=45.0
'define alfa=atan2('lat2-lat1','lon2-lon1')'
'set x 1'
'set y 1'
'set lev 1000 100'
'set t 1'
lon=lon1
'collect 1 free'
'collect 2 free'
'collect 3 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(-w.1*100,'lon','lat')'
'collect 2 gr2stn(u.2*cos(alfa)+v.3*sin(alfa),'lon','lat')'
'collect 3 gr2stn(orog.4(t=1),'lon','lat')'
lon=lon+1
endwhile
'set lon 70 130'
'set lev 1000 100'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'
'set xlabs 35N,70E|37N,82E|39N,94E|41N,106E|43N,118E|45N,130E'
'set gxout shaded'
'set csmooth on'
'set clab off'
'set clevs 0'
'set ccols 0 5 0'
'd coll2gr(3,-u)'
'set gxout stream'
'set clab on'
'd maskout(coll2gr(2,-u),-coll2gr(3,-u));maskout(coll2gr(1,-u),-coll2gr(3,-u));mag(coll2gr(2,-u),coll2gr(1,-u))'

'draw title along (35N,70E) to (45N,130E)'
'printim F:\Cop_Zonal-vertical\anypoumian2.png white'
say 'ok'
;

river

sun_shine_Xia 发表于 2014-3-4 21:36

                               
登录/注册后可看大图

LZ，我也用了你的这个方法去试着画了下，用的脚本就是你贴出来的画流场的那个，为什么他一直提示这个错误 ...

脚本的话看着没有问题，是不是资料有问题呢？介绍下你用的资料

sun_shine_Xia

river 发表于 2014-3-5 10:36

                               
登录/注册后可看大图

脚本的话看着没有问题，是不是资料有问题呢？介绍下你用的资料

我的资料使用NCL处理过然后写出来的~配上ctl之后用grads也没啥问题
这个是垂直速度的ctl文件~
dset F:\Data\EC_Data\W\wNDJ_Dtrend.grd
title NDJ mean  Dtrend w wind
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 14 levels 1000 925 850 700 600 500 400 300 250 200 150 100 50 10
tdef 35 linear 00Z01Feb1979 1yr
vars 1
w    14 99  Dtrend w wind NDJ-mean year
endvars
但是如果我不画流线图，画的只是单要素的等值线就不会有这个问题，就能出图，资料还是这个资料~

river

sun_shine_Xia 发表于 2014-3-5 11:21

                               
登录/注册后可看大图

我的资料使用NCL处理过然后写出来的~配上ctl之后用grads也没啥问题
这个是垂直速度的ctl文件~
dset F:\ ...

U和V也是这样的数据结构？那不应该的

sun_shine_Xia

river 发表于 2014-3-5 15:15

                               
登录/注册后可看大图

U和V也是这样的数据结构？那不应该的

对，也是这样的~而且不管是单独画u,v还是w，都能出图，就是想画流线的时候就不行了~

river

sun_shine_Xia 发表于 2014-3-5 16:06

                               
登录/注册后可看大图

对，也是这样的~而且不管是单独画u,v还是w，都能出图，就是想画流线的时候就不行了~

脚本没有问题，一时半会儿解决不了了

sun_shine_Xia

river 发表于 2014-3-5 15:15

                               
登录/注册后可看大图

U和V也是这样的数据结构？那不应该的

我现在知道问题应该就在'collect 2 gr2stn((u.2*cos(alfa)+v.3*sin(alfa)),'lon','lat')'这句话上~但是不知道该怎么弄啊！

river

sun_shine_Xia 发表于 2014-3-5 17:21

                               
登录/注册后可看大图

我现在知道问题应该就在'collect 2 gr2stn((u.2*cos(alfa)+v.3*sin(alfa)),'lon','lat')'这句话上~但是不 ...

但是这句写的没有问题啊，我也不知道怎么改。你换成nc资料试一下，看看到底是不是脚本的问题。