From 8e37b5c4af28fdf594c4fb40abff59995b6fa4c0 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 12 Jul 2024 17:49:47 +0200
Subject: [PATCH] Pierre T. 12/07/2024: bugfix use of JP_DMS even if not
present in BASIC for DMS fluxes
---
src/MNH/ground_paramn.f90 | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/MNH/ground_paramn.f90 b/src/MNH/ground_paramn.f90
index 79aeb4ff4..bfaa10f45 100644
--- a/src/MNH/ground_paramn.f90
+++ b/src/MNH/ground_paramn.f90
@@ -1072,14 +1072,14 @@ END IF
IF (LUSECHEM) THEN
DO JSV=NSV_CHEMBEG,NSV_CHEMEND
PSFSV(:,:,JSV) = ZSFTS(:,:,JSV) * XMD / ( XAVOGADRO * XRHODREF(:,:,IKB))
- IF ((LCHEMDIAG).AND.(CPROGRAM == 'DIAG ')) XCHFLX(:,:,JSV-NSV_CHEMBEG+1) = PSFSV(:,:,JSV)
- END DO
- IF (CPROGRAM == 'MESONH') THEN
+ IF ((LCHEMDIAG).AND.(CPROGRAM == 'DIAG ')) XCHFLX(:,:,JSV-NSV_CHEMBEG+1) = PSFSV(:,:,JSV)
+ IF ((CPROGRAM == 'MESONH').AND.(TRIM(CNAMES(JSV-NSV_CHEMBEG+1))=='DMS')) THEN
! net flux of so2 in kg/m2/s
- XFLX_DMS(:,:) = ZSFTS(:,:,NSV_CHEMBEG-1+JP_DMS) *62.13E-3 / XAVOGADRO
+ XFLX_DMS(:,:) = ZSFTS(:,:,JSV) *62.13E-3 / XAVOGADRO
! total net flux of so2 in kg/m2
- XFLXT_DMS(:,:) = XFLXT_DMS(:,:)+ ZSFTS(:,:,NSV_CHEMBEG-1+JP_DMS)*XTSTEP
+ XFLXT_DMS(:,:) = XFLXT_DMS(:,:)+ ZSFTS(:,:,JSV)*XTSTEP
END IF
+ END DO
ELSE
PSFSV(:,:,NSV_CHEMBEG:NSV_CHEMEND) = 0.
--
GitLab