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

李文帅

river 发表于 2017-7-10 22:10
因为插值的时候是每次把一个值插值在一个固定的点上，就是为了表示这个意思。而且应该和GrADS软件本身也 ...

好的，谢谢啦

不能用真名924

666.赶快存下来学习。

flying_mouse

这个好，正找呢！

孙亦

好厉害，我要快点升级

liuyw

大神，我使用的是0.25*0.25的资料，然后根据楼主给的脚本得出0.25*0.25的地形资料，套用脚本运行，结果只有地形图像，没有流线图像。请教大神指导！

                               
登录/注册后可看大图

'reinit'
'open D:\data\20170707.ctl'
'open E:\gs\dixing\1X1\orog1.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(-VVELprs*10,'lon','lat')'
'collect 2 gr2stn(UGRDprs*cos(alfa)+VGRDprs*sin(alfa),'lon','lat')'
'collect 3 gr2stn(orog.2(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 coll2gr(2,-u);coll2gr(1,-u);mag(coll2gr(2,-u),coll2gr(1,-u))'
****mag求平方根
'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 E:\gs\xiepoumian\anypoumian_liuxian.gif white'
say 'ok'
;

river

没有报错？那是不是经度增量太大造成的lon=lon+1，你把增量改小试试

xibao

怎么就是不行呢，唉。。。画不出来啊

xibao

有人画出来能教教我吗。

xibao

毕业论文急需要帮忙呀呀呀呀呀！

river

liuyw 发表于 2017-9-6 22:43
大神，我使用的是0.25*0.25的资料，然后根据楼主给的脚本得出0.25*0.25的地形资料，套用脚本运行，结果只有 ...

没有报错？那是不是经度增量太大造成的lon=lon+1，你把增量改小试试