Maud & Juan : 6/11/2014 Bug in pH computing

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)

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