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

守望之月/;-D

已经在用了，画了垂直速度的剖面，赞~

D:\omega section.png

leswell

谢谢楼主分享~很有用

守望之月/;-D

任意两点间斜剖面，垂直速度和带地形的流场，终于都试验成功啦~~

撒花~~

风浪

不错，但是任意剖面的画法怎么没有讲！！

river

风浪 发表于 2014-3-28 10:26
不错，但是任意剖面的画法怎么没有讲！！

看不懂那就说明缺少基础，还是去看看基础知识吧

逐风

画出来了，效果还行吧，多谢

逐风

这个
set x 1
set y 1
是什么意思啊？不太懂唉。。。

river

逐风 发表于 2014-3-28 22:20
这个
set x 1
set y 1

自己去看collect的用法

hubery

river老师您好，小弟想画温度的剖面图，然后叠加上地形，但是老是出错，请river老师指教，gs如下：
'reinit'
'sdfopen D:\20140308paper\outputair.nc'
'open D:\2010\nceporo\orog.ctl'

'set grads off'
'set zlog on'
'set x 1'
'set y 1'
'set lev 1000 300'
'set parea 1 9.25 1 7.75'

lon1 = 78.0
lon2 = 87.0
lat1 = 27.0
lat2 = 35.0
lon = lon1
'collect 1 free'
'collect 2 free'

while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(t_gds0_isbl*pow((1000./lev),0.286),'lon','lat')'
  'collect 2 gr2stn(orog.2,'lon','lat')'
  
lon = lon + 1.0
endwhile



'set lon 78 87'
'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(2,-u)'


'set gxout contour'
'set clab on'
'd maskout(coll2gr(1,-u),-coll2gr(2,-u))'

;

hubery

错误显示如下