From b7f18ae97be149a4cbf00e99c5692582223c99bf Mon Sep 17 00:00:00 2001
From: Juan Escobar <juan.escobar@aero.obs-mip.fr>
Date: Thu, 6 Nov 2014 14:56:13 +0000
Subject: [PATCH] Maud & Juan : 6/11/2014 Bug in pH computing
---
src/MNH/ch_monitorn.f90 | 4 +---
src/MNH/read_field.f90 | 6 +++---
2 files changed, 4 insertions(+), 6 deletions(-)
diff --git a/src/MNH/ch_monitorn.f90 b/src/MNH/ch_monitorn.f90
index 34084ab9f..65972ab68 100644
--- a/src/MNH/ch_monitorn.f90
+++ b/src/MNH/ch_monitorn.f90
@@ -106,6 +106,7 @@ END MODULE MODI_CH_MONITOR_n
!! + supress line transfer H2SO4 from AP to gas phase
!! imply transfer H2SO4 AP in aqueous phase if aq.chem.
!! 04/2014 (C.Lac) Remove GCENTER with FIT temporal scheme
+!! 06/11/14 (M Leriche) Bug in pH computing
!!
!! EXTERNAL
!! --------
@@ -1085,7 +1086,6 @@ DO JL=1,ISVECNMASK
IF (LUSECHAQ.AND.LCH_PH) THEN
SELECT CASE(NRRL)
CASE(1)
- XPHC(:,:,:) = 0.0
DO JM=0,ISVECNPT-1
JI=JM-IDTI*(JM/IDTI)+ISVECMASK(1,JL)
JJ=JM/IDTI-IDTJ*(JM/(IDTI*IDTJ))+ISVECMASK(3,JL)
@@ -1093,8 +1093,6 @@ DO JL=1,ISVECNMASK
XPHC(JI,JJ,JK) = ZPH(JM+1,1)
END DO
CASE(2)
- XPHC(:,:,:) = 0.0
- XPHR(:,:,:) = 0.0
DO JM=0,ISVECNPT-1
JI=JM-IDTI*(JM/IDTI)+ISVECMASK(1,JL)
JJ=JM/IDTI-IDTJ*(JM/(IDTI*IDTJ))+ISVECMASK(3,JL)
diff --git a/src/MNH/read_field.f90 b/src/MNH/read_field.f90
index 8c70c5322..4e4ee8654 100644
--- a/src/MNH/read_field.f90
+++ b/src/MNH/read_field.f90
@@ -224,6 +224,7 @@ END MODULE MODI_READ_FIELD
!! P. Peyrille 06/12 2D west african monsoon : add reading of ADV forcing and addy fluxes
!! C.Lac 03/13 add prognostic supersaturation for C2R2/KHKO
!! Bosseur & Filippi 07/13 Adds Forefire
+!! M. Leriche 11/14 correct bug in pH initialization
!!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -243,7 +244,6 @@ USE MODD_PASPOL
USE MODD_FOREFIRE
#endif
USE MODD_CH_AEROSOL
-USE MODD_CH_MNHC_n, ONLY : XCH_PHINIT
!
USE MODD_RAIN_C2R2_DESCR, ONLY: C2R2NAMES
USE MODD_ICE_C1R3_DESCR, ONLY: C1R3NAMES
@@ -987,7 +987,7 @@ SELECT CASE(HGETPHC) ! pH in cloud water
CALL FMREAD(HINIFILE,YRECFM,HLUOUT,YDIR,PPHC,IGRID,ILENCH, &
YCOMMENT,IRESP)
CASE('INIT')
- PPHC(:,:,:)=XCH_PHINIT
+ PPHC(:,:,:)=0.
END SELECT
!
SELECT CASE(HGETPHR) ! pH in rainwater
@@ -997,7 +997,7 @@ SELECT CASE(HGETPHR) ! pH in rainwater
CALL FMREAD(HINIFILE,YRECFM,HLUOUT,YDIR,PPHR,IGRID,ILENCH, &
YCOMMENT,IRESP)
CASE('INIT')
- PPHR(:,:,:)=XCH_PHINIT
+ PPHR(:,:,:)=0.
END SELECT
!
IRESP=0
--
GitLab