Quentin + Philippe 05/04/2024: write XMUT only if program=MESONH + set it to...

Quentin + Philippe 05/04/2024: write XMUT only if program=MESONH + set it to 0. by default if not found at reading (and remove warning if not found)

src/MNH/read_field.f90

@@ -660,7 +660,10 @@ IF (LIBM .AND. CPROGRAM=='MESONH') THEN
     NDIMS      = 3,         &
     LTIMEDEP   = .TRUE.     )
    !
-   CALL IO_Field_read(TPINIFILE,TZFIELD,PIBM_XMUT)
+   CALL IO_Field_read( TPINIFILE, TZFIELD, PIBM_XMUT, IRESP )
+   !
+  !If field not found (only available in restarts) => set to 0.
+   IF ( IRESP /= 0 ) PIBM_XMUT = 0.
    !
 ENDIF
 !

src/MNH/write_lfin.f90

@@ -507,18 +507,20 @@ 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)
+  IF ( CPROGRAM == 'MESONH' ) THEN
+    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)
+  END IF
    !
 ENDIF
 !