Bug in AMWG diags: Modification of the AMWG scripts related to Taylor diagrams

[bertinia]

From @cecilehannay on April 21, 2017 20:42

A script for Taylor diagrams in the AMWG diagnostics package uses the pressure data for
vertical interpolation for 300-mb zonal wind, and vertical weighted average of relative humidity
and temperature. In the current script, “lev”, lev = po*(hyam+hybm), is used for pressure, but lev
is not the actual pressure at each grid point.

Thus, pressure should be computed using the NCL function, pres_hybrid_ccm, which
calculates pressure at hybrid levels [pres(i,j)=po*hyam+ps(i,j)*hybm]. Also, a variable at a
specific pressure level should be calculated by a vertical interpolation function (vinth2p) and
should not be a value at a closest pressure level.

20160826-AMWG_Taylor_debug.pdf

Copied from original issue: CESM-Development/cime#486

[bertinia]

Hi Cecile,

For these issues associated with bugs in the AMWG diags did you want to
post
them to the CESM_postprocessing github issues (
https://github.com/NCAR/CESM_postprocessing/issues)
instead of the CESM-Development/CIME issues?

Thanks,
Alice

On Fri, Apr 21, 2017 at 2:42 PM, Cecile Hannay [email protected]
wrote:

A script for Taylor diagrams in the AMWG diagnostics package uses the
pressure data for
vertical interpolation for 300-mb zonal wind, and vertical weighted
average of relative humidity
and temperature. In the current script, “lev”, lev = po*(hyam+hybm), is
used for pressure, but lev
is not the actual pressure at each grid point.

Thus, pressure should be computed using the NCL function, pres_hybrid_ccm,
which
calculates pressure at hybrid levels [pres(i,j)=po*hyam+ps(i,j)*hybm].
Also, a variable at a
specific pressure level should be calculated by a vertical interpolation
function (vinth2p) and
should not be a value at a closest pressure level.

20160826-AMWG_Taylor_debug.pdf
https://github.com/CESM-Development/cime/files/948280/20160826-AMWG_Taylor_debug.pdf

—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
CESM-Development/cime#486, or mute the thread
https://github.com/notifications/unsubscribe-auth/AKO3mcztPNPGwr2bi4_VTZ-eh49mwQ-hks5ryRTDgaJpZM4NEw5e
.

--
Alice Bertini
NCAR - CSEG Software Engineer