Changeset 2312 for trunk/LMDZ.MARS/libf/phymars/watercloud_mod.F

Timestamp:

May 6, 2020, 4:46:33 PM (5 years ago)

Author:

lrossi

Message:

MARS GCM:
Implementation of HDO in the physics
New tracers hdo_vap and hdo_ice are added. Cf. traceur.def.isotopes in deftank for exemple of traceur.def.
All HDO related computations are activated by flag hdo=.true. in callphys.def. (see callphys.def.hdo in deftank for an example)
Additional option hdofrac (true by default) allows for turning on/off fractionation (for tests).
For now, only functional with simplified cloud scheme (so microphys=.false. and activice=.false. also recommended).
Initialisation of start files can be done with option 'inihdo' in newstart.
LR

File:

• trunk/LMDZ.MARS/libf/phymars/watercloud_mod.F (modified) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/watercloud_mod.F

@@ -19 +19 @@
       USE watersat_mod, ONLY: watersat
       use tracer_mod, only: nqmx, igcm_h2o_vap, igcm_h2o_ice,
+     &                      igcm_hdo_vap, igcm_hdo_ice,
      &                      igcm_dust_mass, igcm_dust_number,
      &                      igcm_ccn_mass, igcm_ccn_number,
@@ -124 +125 @@
       REAL :: mincloud ! min cloud frac
       LOGICAL, save :: flagcloud=.true. 
+
 c    ** un petit test de coherence
 c       --------------------------
@@ -312 +314 @@
      &          sum_subpdq(ig,l,igcm_h2o_vap)
      &        + pdq(ig,l,igcm_h2o_vap)
+            IF (hdo) THEN
+            sum_subpdq(ig,l,igcm_hdo_ice) = 
+     &          sum_subpdq(ig,l,igcm_hdo_ice)
+     &        + pdq(ig,l,igcm_hdo_ice)
+            sum_subpdq(ig,l,igcm_hdo_vap) = 
+     &          sum_subpdq(ig,l,igcm_hdo_vap)
+     &        + pdq(ig,l,igcm_hdo_vap)
+            ENDIF !hdo
           ENDDO
         ENDDO      
@@ -361 +371 @@
      &          sum_subpdq(ig,l,igcm_h2o_vap)
      &        + subpdqcloud(ig,l,igcm_h2o_vap)
+
+            IF (hdo) THEN
+            sum_subpdq(ig,l,igcm_hdo_ice) =
+     &          sum_subpdq(ig,l,igcm_hdo_ice)
+     &        + subpdqcloud(ig,l,igcm_hdo_ice)
+            sum_subpdq(ig,l,igcm_hdo_vap) =
+     &          sum_subpdq(ig,l,igcm_hdo_vap)
+     &        + subpdqcloud(ig,l,igcm_hdo_vap)
+            ENDIF ! hdo
+
           ENDDO
         ENDDO
@@ -397 +417 @@
      &        sum_subpdq(ig,l,igcm_h2o_vap)/real(imicro) 
      &       - pdq(ig,l,igcm_h2o_vap)
+            IF (hdo) THEN
+            pdqcloud(ig,l,igcm_hdo_ice) = 
+     &        sum_subpdq(ig,l,igcm_hdo_ice)/real(imicro) 
+     &       - pdq(ig,l,igcm_hdo_ice)
+            pdqcloud(ig,l,igcm_hdo_vap) = 
+     &        sum_subpdq(ig,l,igcm_hdo_vap)/real(imicro) 
+     &       - pdq(ig,l,igcm_hdo_vap)
+            ENDIF !hdo
          ENDDO
        ENDDO
@@ -478 +506 @@
         DO l=1,nlay
          DO ig=1,ngrid
+
           IF (pq(ig,l,igcm_h2o_ice) + ptimestep*
      &       (pdq(ig,l,igcm_h2o_ice) + pdqcloud(ig,l,igcm_h2o_ice)) 
@@ -484 +513 @@
      &     - pq(ig,l,igcm_h2o_ice)/ptimestep - pdq(ig,l,igcm_h2o_ice)
            pdqcloud(ig,l,igcm_h2o_vap) = -pdqcloud(ig,l,igcm_h2o_ice)
+             if (hdo) then
+           pdqcloud(ig,l,igcm_hdo_ice) = 
+     &     - pq(ig,l,igcm_hdo_ice)/ptimestep - pdq(ig,l,igcm_hdo_ice)
+           pdqcloud(ig,l,igcm_hdo_vap) = -pdqcloud(ig,l,igcm_hdo_ice)
+             endif
           ENDIF 
           IF (pq(ig,l,igcm_h2o_vap) + ptimestep*
@@ -491 +525 @@
      &     - pq(ig,l,igcm_h2o_vap)/ptimestep - pdq(ig,l,igcm_h2o_vap)
            pdqcloud(ig,l,igcm_h2o_ice) = -pdqcloud(ig,l,igcm_h2o_vap)
+             if (hdo) then
+           pdqcloud(ig,l,igcm_hdo_vap) = 
+     &     - pq(ig,l,igcm_hdo_vap)/ptimestep - pdq(ig,l,igcm_hdo_vap)
+           pdqcloud(ig,l,igcm_hdo_ice) = -pdqcloud(ig,l,igcm_hdo_vap)
+             endif
           ENDIF
-         ENDDO
-        ENDDO
-
+
+         ENDDO
+        ENDDO
 
 c------Update the ice and dust particle size "rice" for output or photochemistry
@@ -627 +666 @@
      &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
      &      igcm_h2o_vap))
+      if (hdo) then
+      call WRITEDIAGFI(ngrid,"pdqiceD","pdqiceD apres microphysique"
+     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,igcm_hdo_ice))
+      call WRITEDIAGFI(ngrid,"pdqvapD","pdqvapD apres microphysique"
+     &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,
+     &      igcm_hdo_vap))
+      endif
       call WRITEDIAGFI(ngrid,"pdqccn2","pdqcloudccn apres microphysique"
      &      ,"kg/kg.s-1",3,pdqcloud(1:ngrid,1:nlay,