From 3aa0e789f000d8435f0a15f64537a280887f6ada Mon Sep 17 00:00:00 2001 From: Alistair Adcroft Date: Tue, 29 Nov 2016 08:25:54 -0500 Subject: [PATCH] Fix "edges" attribute for dimension variable in z-diag output - The old z-diagnostics had "edge positions" in the output file in addition to the level centers. This was missed in the generalized diagnostic coordinate code. --- src/framework/MOM_diag_remap.F90 | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/src/framework/MOM_diag_remap.F90 b/src/framework/MOM_diag_remap.F90 index da814487bb..98b899750f 100644 --- a/src/framework/MOM_diag_remap.F90 +++ b/src/framework/MOM_diag_remap.F90 @@ -169,12 +169,13 @@ subroutine diag_remap_configure_axes(remap_cs, GV, param_file) interfaces(:) = getCoordinateInterfaces(remap_cs%regrid_cs) layers(:) = 0.5 * ( interfaces(1:remap_cs%nz) + interfaces(2:remap_cs%nz+1) ) - remap_cs%layer_axes_id = diag_axis_init(lowercase(trim(remap_cs%diag_coord_name))//'_l', & - layers, trim(units), 'z', & - trim(longname)//' at cell center', direction=-1) remap_cs%interface_axes_id = diag_axis_init(lowercase(trim(remap_cs%diag_coord_name))//'_i', & interfaces, trim(units), 'z', & trim(longname)//' at interface', direction=-1) + remap_cs%layer_axes_id = diag_axis_init(lowercase(trim(remap_cs%diag_coord_name))//'_l', & + layers, trim(units), 'z', & + trim(longname)//' at cell center', direction=-1, & + edges=remap_cs%interface_axes_id) ! Axes have now been configured. remap_cs%configured = .true.