From d1237bddf8a70bb8763ce07a0abfde4fb35e11c5 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Wed, 8 Feb 2023 13:13:49 +0100
Subject: [PATCH] Philippe 08/02/2023: bugfix: IBM: forgot to write XIBM_XMUT

---
 src/MNH/write_lfin.f90 | 15 ++++++++++++++-
 1 file changed, 14 insertions(+), 1 deletion(-)

diff --git a/src/MNH/write_lfin.f90 b/src/MNH/write_lfin.f90
index f20b46cef..154aeff86 100644
--- a/src/MNH/write_lfin.f90
+++ b/src/MNH/write_lfin.f90
@@ -229,7 +229,7 @@ USE MODD_GRID_n
 USE MODD_HURR_CONF, ONLY: LFILTERING,CFILTERING,NDIAG_FILT
 USE MODD_HURR_FIELD_n
 USE MODD_IBM_LSF,         ONLY: LIBM_LSF
-USE MODD_IBM_PARAM_n,     ONLY: LIBM, XIBM_LS
+USE MODD_IBM_PARAM_n,     ONLY: LIBM, XIBM_LS, XIBM_XMUT
 USE MODD_IO, ONLY: TFILEDATA
 USE MODD_LATZ_EDFLX
 USE MODD_LIMA_PRECIP_SCAVENGING_n
@@ -498,6 +498,19 @@ IF (LIBM .OR. LIBM_LSF) THEN
   !
   CALL IO_Field_write(TPFILE,TZFIELD,XIBM_LS(:,:,:,1))
   !
+  TZFIELD = TFIELDMETADATA( &
+    CMNHNAME   = 'XMUT',    &
+    CLONGNAME  = 'XMUT',    &
+    CSTDNAME   = '',        &
+    CUNITS     = 'm2 s-1',  &
+    CDIR       = 'XY',      &
+    NGRID      = 1,         &
+    NTYPE      = TYPEREAL,  &
+    NDIMS      = 3,         &
+    LTIMEDEP   = .TRUE.     )
+   !
+   CALL IO_Field_write(TPFILE,TZFIELD,XIBM_XMUT)
+   !
 ENDIF
 !
 IF (LRECYCL) THEN
-- 
GitLab