Vorticity Calculation in vorticity.F

[cmoffat]

Hello,

I was very glad to see vorticity added to average file output. However, I believe there are a few issues with the way it is being calculated:

• The last term in pvor (line 378) should read -dVdz_pr*dRdx_pr (instead of -dVdz_pr*dRde_pr)

• It looks like rvor is being treated as a scalar, so that all three components are added together. Shouldn't this either be the magnitude of the relative vorticity, or just its vertical component (the first term)?

• I believe there's an extra d(pden)/d(z) in the preamble of the code, where the expressions for pvor and rvor are.

I might be missing something, of course. I've attached a diff file in case it is useful.

Best Regards, Carlos

[cmoffat]

Diff File