From 846b8de9e14230f08c4162f08ebd7c8a783294b6 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 5 Apr 2024 16:19:32 +0200
Subject: [PATCH] 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 | 5 ++++-
src/MNH/write_lfin.f90 | 26 ++++++++++++++------------
2 files changed, 18 insertions(+), 13 deletions(-)
diff --git a/src/MNH/read_field.f90 b/src/MNH/read_field.f90
index 877437197..00593e541 100644
--- a/src/MNH/read_field.f90
+++ b/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
!
diff --git a/src/MNH/write_lfin.f90 b/src/MNH/write_lfin.f90
index 43e3f1e22..5d6517151 100644
--- a/src/MNH/write_lfin.f90
+++ b/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
!
--
GitLab