From e9ce520df0ed8a23a5f2046c2a294d1751a7013f Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 3 May 2024 11:39:18 +0200
Subject: [PATCH] Philippe 03/05/2024: Z-split files: manage correctly fields
with 3rd dimension NMNHDIM_NLWB or NMNHDIM_NSWB
---
src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90 | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90 b/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
index 3dc4cc7bd..9cc90f2f1 100644
--- a/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
+++ b/src/LIB/SURCOUCHE/src/mode_io_write_nc4.f90
@@ -2858,7 +2858,7 @@ END SUBROUTINE IO_History_append_nc4
subroutine IO_Select_split_file( tpfile, tpfield, tpfileout, tpfieldout, kvertlevel, kzfile )
-use modd_field, only: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, NMNHDIM_NOTLISTED, NMNHDIM_UNKNOWN
+use modd_field, only: NMNHDIM_LEVEL, NMNHDIM_LEVEL_W, NMNHDIM_NLWB, NMNHDIM_NSWB, NMNHDIM_NOTLISTED, NMNHDIM_UNKNOWN
type(tfiledata), target, intent(in) :: tpfile
class(tfieldmetadata), target, intent(in) :: tpfield
@@ -2890,6 +2890,8 @@ if ( Present( kvertlevel ) ) then
do jdim = 1, size(tpfieldout%ndimlist)
if ( tpfieldout%ndimlist(jdim) == NMNHDIM_LEVEL &
.or. tpfieldout%ndimlist(jdim) == NMNHDIM_LEVEL_W &
+ .or. tpfieldout%ndimlist(jdim) == NMNHDIM_NLWB &
+ .or. tpfieldout%ndimlist(jdim) == NMNHDIM_NSWB &
.or. tpfieldout%ndimlist(jdim) == NMNHDIM_NOTLISTED ) then
if ( jdim < size(tpfieldout%ndimlist) ) &
tpfieldout%ndimlist(jdim:size(tpfieldout%ndimlist)-1) = tpfieldout%ndimlist(jdim+1:)
--
GitLab