Philippe 29/08/2024: bugfix: allocate XIBM_LS to zero size but with last...

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

@@ -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
 !-------------------------------------------------------------------------------
 !