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

canghaidie

沿任意线画垂直剖面，学习中

canghaidie

非常棒，已完美作图，谢谢！

--惊鸿-

小白感谢river大神的分享，但想请教如下问题：
运行grads，报错如下，我没搞懂它报错的意思，请求指点一下，谢谢了。
【grads报错代码如下】
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Cannot use log scaling when coordinates <= 0
Linear scaling used
Cannot contour grid - all undefined values
No hardcopy metafile open

【gs如下】
'reinit'
'open G:\b\pou\CF_3D_season_clim_WNP.ctl'
'set grads off'
'set grid off'
'set zlog on'
'set csmooth on'
'set x 1'
'set y 1'
'set z 0 40'
'set t 1'
lon1=110.0
lon2=125.0
lat1=5.0
lat2=20.0
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(cf,'lon','lat')'
lon=lon+1
endwhile
'set lon 110 125'
'set z 0 40'
'set ylabs 0 2 4 6 8 10 12 14 16 18'
'set xlabs 5N,110E|20N,125E'
'set gxout shaded'
*'set clevs 26 29 32 35 38 41 44 47 50 53 56 59 62 65 68 71 74 77 80 83 86 89 92 95 98 101'
*'set ccols 59 58 57 56 49 48 47 46 45 44 13 38 37 36 35 34 33 32 21 22 23 24 25 26 27 28 29'
'd coll2gr(1,-u)'
'draw title cloud_fraction along (20N,110E) to (5N,125E)'
'printim G:\b\pou\xie.png x1600 y1200 white'
'disable print'

river

--惊鸿- 发表于 2020-2-25 22:34
小白感谢river大神的分享，但想请教如下问题：
运行grads，报错如下，我没搞懂它报错的意思，请求指点一下 ...

超出数据范围，或者全是缺测值。而且还提示了小于0的时候不能使用log函数，不清楚你这个怎么会有这种错误提示。你好歹说一下你用的资料，把描述文件贴出来看看吧

--惊鸿-

river 发表于 2020-2-26 21:33
超出数据范围，或者全是缺测值。而且还提示了小于0的时候不能使用log函数，不清楚你这个怎么会有这种错误 ...

谢谢river对grads报错的解答。
我用的数据是三维的云数据（cloud_fraction）
之前批量的月资料预处理成4个季节的数据
【ctl如下】
dset G:\b\pou\CF_3D_season_clim_WNP.dat
undef -9999.0
title   historical 30
xdef  81  linear  100.5 1
ydef  61  linear  0.5   1
zdef  40  linear  0  0.48
tdef  4   linear jan2007  6mo
vars  1
cf    40  99   
endvars

river

--惊鸿- 发表于 2020-2-29 09:25
谢谢river对grads报错的解答。
我用的数据是三维的云数据（cloud_fraction）
之前批量的月资料预处理成 ...

你这个资料垂直分辨率这么高啊，而且是均匀层次，不是气压的单位啊。那你得试试把一个参数改改
'd coll2gr(1,40)'
-u这个参数只适用于垂直坐标是气压的情况。
你试试看吧，实在不行你把资料发给我（weather-man@foxmail.com），我试试。

--惊鸿-

river 发表于 2020-2-29 20:44
你这个资料垂直分辨率这么高啊，而且是均匀层次，不是气压的单位啊。那你得试试把一个参数改改
'd coll2 ...

啊啊啊啊把-u改成40一下子就出图了！可以了！
太感谢您了，这个斜剖面图困扰了我好久，真的太谢谢你的指点了。{:eb511:}

river

--惊鸿- 发表于 2020-3-1 18:47
啊啊啊啊把-u改成40一下子就出图了！可以了！
太感谢您了，这个斜剖面图困扰了我好久，真的太谢谢你的指 ...

没事儿多去GrADS官网逛逛，了解了解命令和和函数的用法，说不准就不要困扰喇嘛久了

lemonade

大神请问，我按您的方法画水汽通量qv的垂直剖面图，出来的图地形阴影是有的，但qv等值线有问题，不知道是什么原因。

lemonade

lemonade 发表于 2020-3-26 12:20
大神请问，我按您的方法画水汽通量qv的垂直剖面图，出来的图地形阴影是有的，但qv等值线有问题，不知道是什 ...

补，这是我的ctl和gs文件：（多谢大神喽~）
out.ctl文件：
dset F:\pou\out.grd
title 25Nmay
undef -9.96921e+36
xdef 144 linear 0 2.5
ydef 1 linear 25 2.5
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 1 linear 00Z13May2007 360mn
vars 1
qva 8 99 poumian
endvars
gs文件：
'reinit'
'open F:\pou\out.ctl'
'open F:\NCEPdx\orog.ctl'
'set zlog on'
lon1=80.0
lon2=110.0
'set x 1'
'set y 1'
'set z 1 8'
'set t 1'
lon=lon1
'collect 1 free'
'collect 2 free'
while(lon<=lon2)
'collect 1 gr2stn(qva,'lon',25)'
'collect 2 gr2stn(orog.2(t=1),'lon',25)'
lon=lon+1
endwhile
'set lon 80 110'
'set z 1 8'
'set xlint 10'
'set ylint 100'
'set gxout shaded'
'set csmooth on'
'set clab off'
'set clevs 0'
'set ccols 0 1 0'
'd coll2gr(2,-u)'
'set gxout contour'
'set clab on'
'd maskout(coll2gr(1,-u),-coll2gr(2,-u))'
'printim F:\pou\dixing.png white'
;