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

wj2608150

好，用用试试

wj2608150

很详细 很好用

wj2608150

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

river

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

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

343625321

river 发表于 2015-9-29 21:38
你的涡度和散度不是原始资料有的，而是需要通过计算得到。而要得到涡度和散度是需要插值的，一个点怎么插 ...

你的意思是要先计算涡度  然后再用gr2stn插值 涡度的值，再coll2gr进行画图吗？

river

343625321 发表于 2015-9-30 10:31
你的意思是要先计算涡度  然后再用gr2stn插值 涡度的值，再coll2gr进行画图吗？

对啊，做剖面图的思路肯定都是得先有这个物理量，然后再做剖面。尤其是这类需要计算得到的，必须是这样的。你自己也可以想象，做剖面你至少得保证物理量在你要的剖面附近从东到西，从南到北，从高到低都有数值才行啊；就好比切西瓜，必须要等西瓜长出来了，你才能一切两半或者更多。如果连西瓜都没有，你直接切是不会有结果的

343625321

river 发表于 2015-9-30 11:53
对啊，做剖面图的思路肯定都是得先有这个物理量，然后再做剖面。尤其是这类需要计算得到的，必须是这样的 ...

楼主给的解释很到位。但是我改了之后还是报错。。麻烦再看一下啦
Error from HDIVG:  Invalid dimension environm
  Horizontal environment (X, Y Varying) is re
Operation Error:  Error from hdivg function
  Error ocurred at column 1
DEFINE error:  Invalid expression.
Error from GR2STN:  3rd argument invalid.
Operation Error:  Error from gr2stn function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = gr2stn(div,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(div,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(div,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(div,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)




脚本：
'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)
'define div=hdivg(UGRDprs,VGRDprs)*100000'
lon=lon1
'collect 1 free'

while(lon<=lon2)

'collect 1 gr2stn(div,'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'

river

343625321 发表于 2015-10-9 09:21
楼主给的解释很到位。但是我改了之后还是报错。。麻烦再看一下啦
Error from HDIVG:  Invalid  ...

用一个点怎么计算散度？

jiaohe1990

haohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohaohao

开卷远行

值得仔细研究。