The correct formulation would look something like this:
` ! rainrt is the rain rate in mm/hr (rainrt= rflx*3600.)
! reflectivity (dum) is in DBz
rthrsh=0.0001_r8/3600._r8
if (rflx(i,k+1).ge.rthrsh) then
dum=log10((rflx(i,k+1)*3600._r8)**6._r8)+16._r8
`
I will work on figuring out how and when to push this back into micro_m3_0.F90. It needs to go in the PUMAS update, plus whatever master version is in the repo
The rain rate (rainrt) used in the analytic reflectivity calculation is never set except to zero.
https://github.com/ESCOMP/PUMAS/blob/master/micro_mg3_0.F90#L3668
The correct formulation would look something like this:
` ! rainrt is the rain rate in mm/hr (rainrt= rflx*3600.) ! reflectivity (dum) is in DBz
`
I will work on figuring out how and when to push this back into micro_m3_0.F90. It needs to go in the PUMAS update, plus whatever master version is in the repo