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 https://github.com/CESM-Development/cime/issues/486, or mute the thread https://github.com/notifications/unsubscribe-auth/AKO3mcztPNPGwr2bi4_VTZ-eh49mwQ-hks5ryRTDgaJpZM4NEw5e .

-- Alice Bertini NCAR - CSEG Software Engineer