From 08f213457d69da709d32d1bb35cd1b757541523d Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 15 Jan 2021 13:42:59 +0100
Subject: [PATCH] Philippe 15/01/2021: bugfix: SV_FLUX for flyers were not
computed correctly and overwritten in files
---
src/MNH/aircraft_balloon_evol.f90 | 1 +
src/MNH/write_aircraft_balloon.f90 | 2 +-
2 files changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/MNH/aircraft_balloon_evol.f90 b/src/MNH/aircraft_balloon_evol.f90
index dbb4d6609..9642bcc4e 100644
--- a/src/MNH/aircraft_balloon_evol.f90
+++ b/src/MNH/aircraft_balloon_evol.f90
@@ -632,6 +632,7 @@ IF ( TPFLYER%FLY) THEN
END DO
ZTHW_FLUX(:,:,:) = ZRHO(:,:,:)*XCPD *XTHW_FLUX(II:II+1,IJ:IJ+1,:)
ZRCW_FLUX(:,:,:) = ZRHO(:,:,:)*XLVTT*XRCW_FLUX(II:II+1,IJ:IJ+1,:)
+ ZSVW_FLUX(:,:,:,:) = XSVW_FLUX(II:II+1,IJ:IJ+1,:,:)
END IF
!
diff --git a/src/MNH/write_aircraft_balloon.f90 b/src/MNH/write_aircraft_balloon.f90
index 822a7b2a7..c4036ee46 100644
--- a/src/MNH/write_aircraft_balloon.f90
+++ b/src/MNH/write_aircraft_balloon.f90
@@ -398,7 +398,7 @@ ZWORK6 (1,1,1,:,1,JPROC) = TPFLYER%RCW_FLUX(:)
!
DO JSV=1,SIZE(TPFLYER%SVW_FLUX,2)
JPROC = JPROC + 1
- YTITLE (JPROC) = 'SV_FLUX'
+ WRITE ( YTITLE(JPROC), FMT = '( A7, I3.3 )' ) 'SV_FLUX', JSV
YUNIT (JPROC) = 'SVUNIT m s-1'
YCOMMENT (JPROC) = 'scalar flux'
ZWORK6 (1,1,1,:,1,JPROC) = TPFLYER%SVW_FLUX(:,JSV)
--
GitLab