Changeset 2378 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Jun 18, 2020, 5:12:08 PM (5 years ago)

Author:

lrossi

Message:

MARS GCM

HDO
Correction of an error in newstart for inihdo.
Other minor corrections for HDO cycle.
Transition from fractionation coefficients from Merlivat et al. 1967 to Lamb et al. 2017

LR

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -1010 +1010 @@
         !D/H (SMOW) = 155.76e-6 so HDO/H2O is twice that
 
-          DO l=1,llm
-            DO j=1,jjp1
-              DO i=1,iip1-1
-              q(i,j,l,igcm_hdo_vap)= q(i,j,l,igcm_h2o_vap) * DoverH
-              q(i,j,l,igcm_hdo_ice)= q(i,j,l,igcm_h2o_ice) * DoverH
-              ENDDO
-              ! We want to have the very same value at lon -180 and lon 180
-              q(iip1,j,l,igcm_hdo_vap) = q(1,j,l,igcm_hdo_vap)
-              q(iip1,j,l,igcm_hdo_ice) = q(1,j,l,igcm_hdo_ice)
-            ENDDO
-          ENDDO
-
            do ig=1,ngridmx
            qsurf(ig,igcm_h2o_ice)=max(0.,qsurf(ig,igcm_h2o_ice))
-           qsurf(ig,igcm_hdo_ice)= DoverH * qsurf(ig,igcm_h2o_ice)
            end do
+
+        ! Update the hdo tracers
+         q(1:iip1,1:jjp1,1:llm,igcm_hdo_vap)
+     &      =q(1:iip1,1:jjp1,1:llm,igcm_h2o_vap)* DoverH
+         q(1:iip1,1:jjp1,1:llm,igcm_hdo_ice)
+     &      =q(1:iip1,1:jjp1,1:llm,igcm_h2o_ice)* DoverH
+
+         qsurf(1:ngridmx,igcm_hdo_ice)
+     &      =qsurf(1:ngridmx,igcm_h2o_ice)*DoverH