算假相当位温时出现很多缺测值

又又又加零

在进行假相当位温的运算时，计算露点温度下的饱和水汽压es时出现了警告：

C:\Users\user\AppData\Local\Temp\ipykernel_6572\2629662609.py:11: RuntimeWarning: invalid value encountered in double_scalars
  cttd[i,j] = e0*((t0/ttd[i,j])**(cl/Rw))*(np.e**((L0+cl*t0)*(ttd[i,j]-t0)/(Rw*ttd[i,j]*t0)))
E:\Anaconda3\lib\site-packages\xarray\core\computation.py:771: RuntimeWarning: invalid value encountered in log
  result_data = func(*input_data)

我画的是剖面图，发现高层几乎全是缺测值，我想知道怎么才能避免呢？感觉算法没有什么问题，求指导

看了下警告应该是在计算 cttd[i,j] 的过程中可能 ttd[i,j] 的值为零或者负数；第二个是np.log()取对数时出现参数为负，导致了缺测。我看低层基本上都有计算值出来，高层出现很多缺测，不知道怎么解决。

river

印象里应该是求出来的水汽压值有负值，可以加上一个很小的数来解决，我看好多都是直接给水汽压加上1e-10 就可以了。你的原始资料如果有相对湿度，用相对湿度计算露点温度就不会出现这个问题。

又又又加零

解决了一个大乌龙，读错了一个数据

又又又加零

river 发表于 2023-4-15 18:37
印象里应该是求出来的水汽压值有负值，可以加上一个很小的数来解决，我看好多都是直接给水汽压加上1e-10 就 ...

非常谢谢您的帮助，学到了。