From 61033ecf81a4cca83fc17142f8988f080bfcaa06 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@cnrs.fr>
Date: Thu, 29 Aug 2024 09:34:48 +0200
Subject: [PATCH] Philippe 29/08/2024: bugfix: allocate XIBM_LS to zero size
 but with last dimension to 1 if LIBM=F (to prevent crash in PHYS_PARAM_n when
 it is passed to TURB)

---
 src/MNH/ini_modeln.f90 | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/MNH/ini_modeln.f90 b/src/MNH/ini_modeln.f90
index f88d337e5..0ad3212f5 100644
--- a/src/MNH/ini_modeln.f90
+++ b/src/MNH/ini_modeln.f90
@@ -2572,6 +2572,8 @@ IF (LIBM) THEN
   ALLOCATE(XIBM_LS(IIU,IJU,IKU,4)) ; XIBM_LS  = -XIBM_IEPS
   XIBM_LS(:,:,:,1)=ZIBM_LS(:,:,:)
   DEALLOCATE(ZIBM_LS)
+ELSE
+  ALLOCATE( XIBM_LS(0,0,0,1) ) !Last dimension should be at least 1 because it is passed in phys_param_n
 ENDIF
 !-------------------------------------------------------------------------------
 !
-- 
GitLab