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

PIZI

很好很好{:5_213:}

PIZI

很好很好啊

PIZI

赞一个{:eb302:}

PIZI

好。。。。。。。。。。。

343625321

7iong 发表于 2013-3-19 14:31
非常感谢，学习了，请问我要是想要画剖面涡度，'define vor=hcurl(coll2gr(1,-u),coll2gr(2,-u))'！！！！ ...

请问你最后画出来了吗？我用同样的方法老报错呀~~求给个脚本！ 谢谢！！

343625321

楼主你好，我按照你的方法用fnl资料画散度任意剖面图，但是出现错误，麻烦帮我看一下！万分感谢！
Error from GR2STN:  3rd argument invalid.
Operation Error:  Error from gr2stn function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = gr2stn(hdivg(ugrdprs,vgrdprs)*100000,118.0,lat)
Error from GR2STN:  3rd argument invalid.
Operation Error:  Error from gr2stn function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = gr2stn(hdivg(ugrdprs,vgrdprs)*100000,119,lat)
Error from GR2STN:  3rd argument invalid.
Operation Error:  Error from gr2stn function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = gr2stn(hdivg(ugrdprs,vgrdprs)*100000,120,lat)
Error from GR2STN:  3rd argument invalid.
Operation Error:  Error from gr2stn function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = gr2stn(hdivg(ugrdprs,vgrdprs)*100000,121,lat)
Error from COLL2GR:  Collection 1 empty
Operation Error:  Error from coll2gr function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = coll2gr(1,-u)




gs如下：
'reinit'
'open E:/grads/fnl_2015.ctl'
'enable print E:\grads\slice\slice.gmf'
'set mpdset cnhimap'
'set grads off'
'set grid off'
'set zlog on'
lon1=118.0
lon2=121.0
lat1=24.0
lat2=27.0
'define alfa=atan2('lat2-lat1','lon2-lon1')'

'set x 1'
'set y 1'
'set lev 1000 100'
'set t 16'
'set t 6'
'q dim'
tt=sublin(result,5)
tt=subwrd(tt,6)

lon=lon1
'collect 1 free'

while(lon<=lon2)

'collect 1 gr2stn(hdivg(UGRDprs,VGRDprs)*100000,'lon','lat')'

lon=lon+1
endwhile
'set lon 100 130'
'set lev 1000 100'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'
'set xlabs 24N,118E|25N,119E|26N,120E|27N,121E'
'set gxout shaded'
'set csmooth on'
'set clab on'
'd coll2gr(1,-u)'
'draw title 'tt'-slice-div'
'printim E:\grads\slice\'tt'-slice-div.jpg white'
'print'
'disable print'
   

麻烦帮我看看 万分感谢！

wj2608150

好，用用试试

wj2608150

很详细 很好用

wj2608150

请问怎么用wrf结果做这种图吗？地形不会做

river

343625321 发表于 2015-9-29 20:01
楼主你好，我按照你的方法用fnl资料画散度任意剖面图，但是出现错误，麻烦帮我看一下！万分感谢！
Error f ...

你的涡度和散度不是原始资料有的，而是需要通过计算得到。而要得到涡度和散度是需要插值的，一个点怎么插值？