search

josephhardinee / PyDSD

Python Library for working with disdrometer data.
GNU Lesser General Public License v2.1
45 stars 33 forks source link

rain rate estimation #5

Closed bdolan44 closed 10 years ago

bdolan44 commented 10 years ago

It would be nice to have the double-moment rain estimators from the disdrometer data as well, R(Kdp,Zdr) and R(Zh, Zdr). Thanks!

josephhardinee commented 10 years ago

Commit c7ff1d5f1ac07fdb6e00c12fae87bf7d2ba5f234 on the 2moment branch starts to address this. It has all three dual moment estimators, as well as all three single moment estimators. It uses the Levenberg-Marquardt algorithm to estimate the parameters based upon nonlinear least squares. It requires a rain rate exist(and I added a calc_RR function that uses estimated terminal velocity to reconstruct rain rate if the disdrometer file did not provide it). I'm leaving this open until some testing has been done on this branch and it gets merged into master.

josephhardinee commented 10 years ago

Commit 9a14b148828043bdd960305640da1578c6e52b18 adds support for the different pairs of 2 moment estimators. This requires that the dsd object has a rain_rate. I've added a calc_RR() function as well to support data files that don't have the rain rate embedded in them. This closes this issue. I'll work on merging the two nasa GV formats for disdrometers next to bring in the flux measurement for rain rate.