From daff9222394570a52937fca6e08a5fc800a9f8d8 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Thu, 11 Jul 2024 11:56:00 +0200
Subject: [PATCH] Philippe 11/07/2024: outputs: fix for 3D variables in boxes
---
src/LIB/SURCOUCHE/src/mode_io_tools.f90 | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/LIB/SURCOUCHE/src/mode_io_tools.f90 b/src/LIB/SURCOUCHE/src/mode_io_tools.f90
index 0bc64e813..baa85c944 100644
--- a/src/LIB/SURCOUCHE/src/mode_io_tools.f90
+++ b/src/LIB/SURCOUCHE/src/mode_io_tools.f90
@@ -226,7 +226,7 @@ SUBROUTINE IO_Dim_box_bounds( TPBOX, TPFIELD, KIINF, KISUP, KJINF, KJSUP, KKINF,
! Subroutine to return the box boundaries in global coordinates
USE MODD_BAKOUT, ONLY: LOUT_PHYSICAL_SIMPLIFIED
USE MODD_DIM_n, ONLY: NIMAX_ll, NJMAX_ll, NKMAX
- USE MODD_FIELD, ONLY: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, TFIELDMETADATA
+ USE MODD_FIELD, ONLY: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, NMNHDIM_BOX_LEVEL, NMNHDIM_BOX_LEVEL_W, TFIELDMETADATA
USE MODD_OUT_n, ONLY: CMAINDOMAINNAME, TOUTBOXMETADATA
USE MODD_PARAMETERS, ONLY: JPHEXT, JPVEXT
@@ -267,7 +267,8 @@ SUBROUTINE IO_Dim_box_bounds( TPBOX, TPFIELD, KIINF, KISUP, KJINF, KJSUP, KKINF,
IF ( .NOT.PRESENT(KKINF) .OR. .NOT.PRESENT(KKSUP) ) &
call Print_msg( NVERB_FATAL, 'IO', 'IO_Dim_box_bounds', trim(tpfield%cmnhname) // ': missing needed dummy arguments' )
- IF ( TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL_W ) THEN
+ IF ( TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_LEVEL_W .OR. &
+ TPFIELD%NDIMLIST(3) == NMNHDIM_BOX_LEVEL .OR. TPFIELD%NDIMLIST(3) == NMNHDIM_BOX_LEVEL_W ) THEN
KKINF = JPVEXT + TPBOX%NKINF
IF ( TPBOX%CNAME == CMAINDOMAINNAME .AND. TPFIELD%NGRID == 4 .AND. .NOT.LOUT_PHYSICAL_SIMPLIFIED ) THEN
! There is one more PHYSICAL vertical layer for w points (but the same number for physical + unphysical borders)
--
GitLab