fnl资料中水汽通量（散度）、垂直速度等作图

蔡彬彬

一些作图仅供参考！！！！！！！

zlh

你好 可以请问下你的ctl文件是怎么生成的？

蔡彬彬

zlh 发表于 2021-4-16 17:54
你好 可以请问下你的ctl文件是怎么生成的？

主要是套用她们画t-lnp的ctl命令

maxia

你好   用了你的gs文件绘制垂直速度剖面图，但是总是提示第一行错误，我搞不清楚自己哪儿出了错误？

maxia

这是文件'reinit'
'open G:\t\fnl_2021.ctl'
'set lon 100 110'
'set lat 35 40'
'set t 1 5'
'set lev 1000 200'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define t=TMPprs+0.608*q*TMPprs+1e-10'
'define w=-(pre*9.8*VVELprs)/(287.05*t)'
*--------GRADS设置-----------
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
*---------垂直速度------------
'set lev 1000 200'
'set grads off'
'set lon 106.2'
'set lat 35.3'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200'
'set xlabs 00|06|12|18|00'
'set gxout shaded'
'set cterp off'
'set csmooth off'
'set clap on'
'set strsiz 0.3'
'set ccolor rainbow'
'd w'
'set gxout contour'
'set cint 0.3'
'set cthick 5'
'set clskip 1'
'draw ylab pressure(hPa)'
'draw xlab time(UTC)'
'd w'
'print'

'printim g:\t\out\w.jpg'
;

lhhyjn

非常感谢分享谢谢

我陪你一世流年

请问楼主，最后一幅图 windsher.gs这样报错怎么解决呢？Warning: Contour levels are not strictly increasing
         This may lead to errors or undesired results

sicuo151

maxia 发表于 2021-12-4 20:28
这是文件'reinit'
'open G:\t\fnl_2021.ctl'
'set lon 100 110'

解决了嘛
大神

云南小黑哥

66666