利用NCEP资料画剖面图的一点儿经验

wdd1992

学习了，帮助很大

wdd1992

火速mark。。。。

白令海

老师，您好。我的gs如下，还是报错了。

'reinit'

'sdfopen F:\Data\GrADS\2016\slp.2016.nc'
*'set dfile 1'
'sdfopen F:\Data\GrADS\2016\hgt.2016.nc'
*'set dfile 2'
'sdfopen F:\Data\GrADS\2016\air.2016.nc'
*'set dfile 3'
'sdfopen F:\Data\GrADS\2016\uwnd.2016.nc'
*'set dfile 4'
'sdfopen F:\Data\GrADS\2016\vwnd.2016.nc'
*'set dfile 5'
'sdfopen F:\Data\GrADS\2016\shum.2016.nc'
*'set dfile 6'
'sdfopen F:\Data\GrADS\2016\rhum.2016.nc'
*'set dfile 7'
'sdfopen F:\Data\GrADS\2016\omega.2016.nc'
*'set dfile 8'

'set grads off'
'set grid off'
*'set parea 1.0 8.0  1.3 9.8'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.25'
*'set xlabs 25`3.|26`3.|27`3.|28`3.|29`3.|30`3.`1N'
'set ylevs 800 700 600 500 400 300'
'set xlab on'
'set ylab on'
'set zlog on'

*************************计算水汽通量散度*************************
'set lon 80 120'
'set lat 10 50'
'set z 1000 300'
'set t 561 672'
'define qdivg=(hdivg(shum.6*uwnd.4,shum.6*vwnd.5))/9.8'
******************************************************************

'set lon 102.5'
'set lat 27.5'
'd qdivg'

*'draw ylab Level/hPa'

'set string 1 c 9 0'
'set strsiz 1.0'
'draw string 3.1  9.5  `4  (b)'

'enable print F:\Data\GrADS\2016\test.gmf'
'print'
'disable print'

浮生未歇

辛苦了楼主，十分感谢你的分享

firebee

{:eb502:}

river

白令海 发表于 2018-5-8 16:57
老师，您好。我的gs如下，还是报错了。

'reinit'

'set z 1000 300' 这是什么意思？

白令海

river 发表于 2018-5-11 13:02
'set z 1000 300' 这是什么意思？

写错了……应该是 'set lev 1000 300'

白令海

river 发表于 2018-5-11 13:02
'set z 1000 300' 这是什么意思？

我改正过来了，但还是报错。

river

白令海 发表于 2018-5-14 19:33
我改正过来了，但还是报错。

你这个问题好像不是脚本的问题了，你试试水汽通量的算法用我帖子里的行不行······

白令海

river 发表于 2018-5-14 21:18
你这个问题好像不是脚本的问题了，你试试水汽通量的算法用我帖子里的行不行······

是电脑的问题。换了一台电脑，用同样版本的 openGrADS 就能出图了。多谢 river老师。