露点温度无法画出垂直剖面

summerthere

借鉴了论坛露点温度的计算方法 但在画垂直方向的剖面时提示错误

如果改为区域范围内某一层就可以正常出图了 求大神解答

gs内容如下：
'reinit'
'sdfopen C:\opengrads\Classic\data\air.2019.nc'
'sdfopen C:\opengrads\Classic\data\rhum.2019.nc'
'set lon 60 120'
'set lat 40'
'set t 1360'
'set lev 1000 200'
'define T=air.1-273.16'
'define RH=rhum.2'
'define td=T-((14.55+0.114*T)*(1-0.01*RH) + pow((2.5+0.007*T)*(1-0.01*RH),3) + (15.9+0.117*T)*pow((1-0.01*RH),14))'
'set gxout shaded'
'd td'
'run C:\opengrads\Contents\Resources\Scripts\cbarn.gs'
'q dims'
'printim C:\opengrads\Classic\data\td.png white'
'c'

伽蓝鸟

你用的资料，湿度是不是只到300 hPa？200 hPa上没有湿度，当然没法定义，算出来也就有问题了

伽蓝鸟

哦对了，一般大家用的是opengrads，这里头有内置函数来算露点温度啊。额论坛上算露点温度的帖子太多，结果又不太一样，晕。用dewpt(T,RH)这个函数就好了。
This function computes the dew-point temperature given the Temperature T in Kelvin and the relative humidity RH in percent.

summerthere

伽蓝鸟 发表于 2020-3-5 08:10
哦对了，一般大家用的是opengrads，这里头有内置函数来算露点温度啊。额论坛上算露点温度的帖子太多，结果 ...

确实是因为上层没有湿度的原因，受教了，万分感谢。