画相当位温垂直剖面图怎样一层显示一层的风向杆

┏_(゜ω゜)=

'reinit'
'open e:/fnl/fnl_20200519_12_00.ctl'
'enable print e:/fnl/se.gmf'
'set lon 112.475'
'set lat 22.991'  

'set lev 1000 100'
'set t 1 18'
'define prs=lev';*获得某层高度的气压
'set ccolor 1'
'set grid off'
'define tc=TMPprs-273.16' ;*某高度层的摄氏温度C
'define tk=TMPprs' ;*某高度层的开氏温度K
'define rh=RHprs' ;*某高度层的相对湿度Relative humidity%
*求饱和水汽压Tetens经验公式
*水面es,tk开氏温度,tc摄氏温度
if(tk>273.16)
'define es=6.1078*exp(17.2693882*tc/(tk-35.86))'
endif
*冰面es,tk开氏温度,tc摄氏温度
if(tk<=273.16)
'define es=6.1078*exp(21.8745584*tc/(tk-7.66))'
endif
*饱和比湿
'define qs=0.622*es/(prs-0.378*es)'
*用相对湿度等求比湿
'define q=rh*qs/100'
*水汽压
'define e=prs*q/(0.622+q)'
*凝结高度的绝对温度,tk起始面上绝对温度K，
'define tlcl=55.0+2840.0/(3.5*log(tk)-log(e)-4.805)'
*求假相当位温Bolton公式,se为开氏温度K
'define theta=tk*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define se=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set parea 1 9.4 1 8'
'set xlopts 1 4 0.14'
'set ylopts 1 4 0.14'
'set ylevs 1000 925 850 700 500 200 100'
'set grads off'
'set zlog on'
'set gxout shaded'
'set cint 2'

'set rgb 20 0 89 162'
'set rgb 21 10 105 199'
'set rgb 22 15 114 255'
'set rgb 23 8 109 241'
'set rgb 24 21 161 255'
'set rgb 25 143 213 255'
'set rgb 26 215 224 255'
'set rgb 27 255 229 230'
'set rgb 28 252 198 172'
'set rgb 29 241 178 145'
'set rgb 30 236 146 122'
'set rgb 31 232 122 99'
'set rgb 32 255 88 21'
'set rgb 33 254 155 25'
'set rgb 34 255 209 39'
'set rgb 35 244 243 117'
'set rgb 36 255 252 101'
'set rgb 37 254 254 158'

'set rbcols 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 '
'd se'
*'cbar.gs'

'run cbar.gs'
'set gxout barb'


'd UGRDprs;VGRDprs'

'draw ylab Pa'
'print'
'disable print'

┏_(゜ω゜)=

最后的效果
'reinit'
'open e:/fnl/fnl_20200519_12_00.ctl'
'enable print e:/fnl/se.gmf'
'set lon 112.475'
'set lat 22.991'  

'set lev 1000 100'
'set t 1 18'
'define prs=lev';*获得某层高度的气压
'set ccolor 1'
'set grid off'
'define tc=TMPprs-273.16' ;*某高度层的摄氏温度C
'define tk=TMPprs' ;*某高度层的开氏温度K
'define rh=RHprs' ;*某高度层的相对湿度Relative humidity%
*求饱和水汽压Tetens经验公式
*水面es,tk开氏温度,tc摄氏温度
if(tk>273.16)
'define es=6.1078*exp(17.2693882*tc/(tk-35.86))'
endif
*冰面es,tk开氏温度,tc摄氏温度
if(tk<=273.16)
'define es=6.1078*exp(21.8745584*tc/(tk-7.66))'
endif
*饱和比湿
'define qs=0.622*es/(prs-0.378*es)'
*用相对湿度等求比湿
'define q=rh*qs/100'
*水汽压
'define e=prs*q/(0.622+q)'
*凝结高度的绝对温度,tk起始面上绝对温度K，
'define tlcl=55.0+2840.0/(3.5*log(tk)-log(e)-4.805)'
*求假相当位温Bolton公式,se为开氏温度K

'define theta=tk*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define se=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set parea 1 9.4 1 8'
*设置画布尺寸
'set xlopts 1 4 0.14'
*设置X轴标注颜色 粗细 大小
'set ylopts 1 4 0.14'
'set ylevs 1000 925 850 700 500 200 100'
'set grads off'
'set zlog on'
'set gxout shaded'
'set cint 2'

'set rgb 20 0 89 162'
'set rgb 21 10 105 199'
'set rgb 22 15 114 255'
'set rgb 23 8 109 241'
'set rgb 24 21 161 255'
'set rgb 25 143 213 255'
'set rgb 26 215 224 255'
'set rgb 27 255 229 230'
'set rgb 28 252 198 172'
'set rgb 29 241 178 145'
'set rgb 30 236 146 122'
'set rgb 31 232 122 99'
'set rgb 32 255 88 21'
'set rgb 33 254 155 25'
'set rgb 34 255 209 39'
'set rgb 35 244 243 117'
'set rgb 36 255 252 101'
'set rgb 37 254 254 158'
*色标采用气象家园调色盘克隆功能。
'set rbcols 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 '
'set clevs 300 305 310 315 320 325 330 335 340 345 350 355 360 365 370 375 380'
'd se'

*'cbar.gs'

'run cbarn.gs'
'set gxout barb'
'set bargap 50'
'set digsiz 0.05'
*设置箭头大小
'set csmooth on'
'set arrowhead 0.05'
*设置箭头尺寸
*'set lev 1000 925 850 700 500 200 100'
'd skip(UGRDprs*2.5,2,1);skip(VGRDprs*2.5,2,1)'
*风向杆的换算*2.5   http://bbs.06climate.com/forum.p ... ght=barb&page=1
'draw ylab Pa'
'print'
'disable print'

'reinit'

┏_(゜ω゜)=

经过修改，最后的效果，还是差点

ELLAZH

要画业务网的图的话，你改成网格经度2试试看，我是画经度二的

┏_(゜ω゜)=

ELLAZH 发表于 2020-10-20 13:00
要画业务网的图的话，你改成网格经度2试试看，我是画经度二的

改哪个网格经度

ELLAZH

┏_(゜ω゜)= 发表于 2020-10-20 14:44
改哪个网格经度

d skip(UGRDprs,1);skip(VGRDprs,1)用这个语句

Itsljw

谢谢楼主分享！