边界层高度grads脚本求。急急急

river

tangjie199365 发表于 2017-10-15 06:30
是z 1 10 那里的问题 。因为ctl文件里面只有一层的数据 画不出来 1到10 层的。谢谢你把纵坐标设 ...

PBL HEIGHT (m) 你资料里这个变量不就是高度吗？单位是米，那出来的图纵坐标就是高度啊。你这个高度有没有负值，控制一下纵坐标数值范围可以不画负值的部分。具体命令自己去官网学习吧http://cola.gmu.edu/grads/gadoc/gadocindex.html

tangjie199365

river 发表于 2017-10-15 15:14
PBL HEIGHT (m) 你资料里这个变量不就是高度吗？单位是米，那出来的图纵坐标就是高度啊。你这个高度有没 ...

好的好的      谢谢

tangjie199365

river 发表于 2017-10-15 15:14
PBL HEIGHT (m) 你资料里这个变量不就是高度吗？单位是米，那出来的图纵坐标就是高度啊。你这个高度有没 ...

把图画出来之后 显示这样的错误 Data Request Error: Invalid grid coordinates
   World coordinates convert to non-integer  grid coordinates
     Variable = pblh.3  Dimension = 0
   Error ocurred at column 1
DISPLAY error:  Invalid expression
   Expression = pblh.3
Data Request Error: Invalid grid coordinates
   World coordinates convert to non-integer  grid coordinates
     Variable = pblh.4  Dimension = 0
   Error ocurred at column 1
DISPLAY error:  Invalid expression
   Expression = pblh.4
这是我脚本 除了我在图中标记出来的3和4  不知道为什么 画不出来 求大神 解决其余的都画出来了
'reinit'
'open 20091109/1/test02z.ctl'
'open 20091109/2/test02z.ctl'
'open 20091109/3/test02z.ctl'
'open 20091109/4/test02z.ctl'
'open 20091109/5/test02z.ctl'
'open 20091109/6/test02z.ctl'
'open 20091109/7/test02z.ctl'
'set grads off'
'set grid off'
'set lon 119.2'
'set lat 30'
'set z 1'
'set t 1 73'

'set font 4'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set csmooth on'
*'set ylevs 0 100 200 300 400 500 600 700 800 900 '
'set ccolor 2'
'set cstyle 1'
'set cmark 2'
*'set digsiz 0.12'
'set cthick 6'
'd PBLH.1'
'set ccolor 3'
'set cstyle 1'
'set cmark 5'
'set cthick 6'
'd PBLH.2'
'set ccolor 9'
'set cstyle 1'
'set cmark 7'
'set cthick 6'
'd PBLH.3'
'set ccolor 1'
'set cstyle 1'
'set cmark 8'
'set cthick 6'
'd PBLH.4'
'set ccolor 4'
'set cstyle 1'
'set cmark 4'
'set cthick 6'
'd PBLH.5'
'set ccolor 5'
'set cstyle 1'
'set cmark 3'
'set cthick 6'
'd PBLH.6'
'set ccolor 7'
'set cstyle 1'
'set cmark 1'
'set cthick 6'
'd PBLH.7'

'cbar_line -x 9.1 -y 7.7 -c 2 3 9 1 4 5 7 -m 2 5 7 8 4 3 1 -t "1" "2" "3" "4" "5" "6" "7"'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
'draw title 20091109-PBLH'

'printim /nuist/scratch/tj199365/wrfout/20091109/PBLH/pblhd2-'6'.png white'

tangjie199365

请大神帮忙解决 求教导{:eb303:}{:eb303:}

river

tangjie199365 发表于 2017-10-18 09:28
请大神帮忙解决 求教导

我不知道怎么帮你了，反正脚本上基本没问题。你自己看错误提示啊，不是说格点坐标有问题吗，你自己看看你画不出来的两个资料和其他几个有没有不一样的地方

陈陈yu

tangjie199365 发表于 2017-10-12 06:34
ncl也可以

这种用ncl可以这么处理呀？请问大佬有用ncl处理.grib边界层高度的脚本吗？

clarawho

                               
登录/注册后可看大图

                               
登录/注册后可看大图

来自苹果APP客户端