Function to compute maximun eddy growth rate

Open avatar101 opened this issue 3 years ago • 1 comments

What should we add?

The maximum Eady growth rate (EGR) is a measure of baroclinic instability. I couldn't find a function in MetPy which implements that. NCL has a function implemented to do that here. The formula is not complicated to implement.

However, does MetPy also considers making the NCL functionality available or should people try using python-NCL instead? I see that it has computations for a wide range of mesoscale to synoptic scale variables. So in some ways, MetPy does implement some of the existing NCL functions. I don't know if python-NCL now offers xarray compatibility and thus, I have been reluctant to use it so far.

If the developers decide to implement EGR function then the formula and the appropriate academic references are available on the NCL webpage linked above.

Thank you

Reference

References and the formula are given on the NCL webpage

Mar 17 '22 08:03 avatar101

I think it makes sense to include this in MetPy. Not sure when we'll have cycles internally to implement this, but would be happy to review a PR contributing this with a set of tests (probably using values from NCL or the publication) to validate the correctness of the implementation.

Mar 17 '22 23:03 dopplershift