NCAR
SRH calculation values overestimated
Using the SRH calculation in wrf-python and also NCL (I think they are equivalent), I have found that the SRH values that I am getting are too large when compared to values computed in MetPy and those on observed soundings and in operational model output. After parsing through the source code (wrf_relhl.f90), I think that this overestimation is due to the estimation of the storm motion, which is assumed to be 0.75 * the mean wind over the 3-10 km layer and 30 degrees to the right of this vector. Since most other methods of computing SRH rely on the Bunkers storm motion, which is related to the mean wind over the ~0-6 km layer, I think that the values output in this routine would be quite a bit larger than what is standard as they accumulate through the integration. Would it be possible to update this code using the Bunkers calculation?
