分享一个由雨滴谱计算降水微物理量的程序

灰色节能君

由雨滴谱计算降水微物理量的程序，程序里输入量d_bin是每个档位中心的直径，w_bin是档位宽度，DSD_bin是m*n的雨滴谱数据矩阵，m个时次，n个档位。输出量包括总数浓度、液态水含量、降水强度、平均直径、质量/体积权重平均直径、归一化截距参数、雷达反射率因子、体积中值直径。
程序计算有错误的地方欢迎指出

mushroom111

这是使用雨滴谱仪观测计算的吗？请问楼主了不了解用双偏振雷达反演雨滴谱呀？

日不不不落

您好 您这个里面计算降水率的时候用到了平均下落速度，程序中我看您使用v=3.78*d_bin.^0.67公式计算的，前期再求雨滴谱DSD_bin其实也用到了这个平均下落速度，因为速度也分了32档，我用的是每一档速度的中值速度，不知道这俩是否一样？您这个求v的公式有参考资料吗

灰色节能君

日不不不落 发表于 2021-3-25 10:28
您好 您这个里面计算降水率的时候用到了平均下落速度，程序中我看您使用v=3.78*d_bin.^0.67公式计算的，前 ...

这个公式在最早出现在Atlas, D., & Ulbrich, C. W. (1977). Path- and Area-Integrated Rainfall Measurement by Microwave Attenuation in the 1–3 cm Band, Journal of Applied Meteorology and Climatology, 16(12), 1322-1331
还有一个常用公式是Brandes E A, Zhang G, Vivekanandan J. Experiments in Rainfall Estimation with a Polarimetric Radar in a Subtropical Environment[J]. Journal of Applied Meteorology, 2002, 41(6): 674–685.
感觉两个算出来差不多吧

DSD_bin是雨滴谱仪输出的，不是我计算的。不太理解计算DSD_bin也用到
平均下落速度的意思。

灰色节能君

mushroom111 发表于 2021-1-26 10:38
这是使用雨滴谱仪观测计算的吗？请问楼主了不了解用双偏振雷达反演雨滴谱呀？

程序DSD_bin就是二维雨滴谱仪输出的。反演雨滴谱了解一点，就是用偏振雷达参量拟合标准化伽马分布里的Nw和D0。能力有限，数据质量控制不好做出来结果也就那样

日不不不落

灰色节能君 发表于 2021-3-25 11:40
这个公式在最早出现在Atlas, D., & Ulbrich, C. W. (1977). Path- and Area-Integrated Rainfall Measure ...

我用的OTT parsivel，输出的变量只有直径，速度和个数，雨滴谱密度还得自己通过公示计算

日不不不落

灰色节能君 发表于 2021-3-25 11:40
这个公式在最早出现在Atlas, D., & Ulbrich, C. W. (1977). Path- and Area-Integrated Rainfall Measure ...

我用您这个代码计算出来的降水率值不对 不知道哪个变量用错了

luzhiang

感谢楼主，先下载看看

MyBeginner

谢谢楼主，下载看看

qq469015280

大佬，请问能分享一个读取雨滴谱文件的程序吗？