Changeset 2312 for trunk/LMDZ.MARS/libf/phymars/phyredem.F90

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/phyredem.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/phyredem.F90

@@ -189 +189 @@
   integer :: i_h2o_vap=0
   integer :: i_h2o_ice=0
+  integer :: i_hdo_vap=0
+  integer :: i_hdo_ice=0
 
   
@@ -245 +247 @@
       i_h2o_ice=iq
     endif
+
+  ! look for HDO vapour & HDO ice tracers (if any)
+    if (noms(iq).eq."hdo_vap") then
+      i_hdo_vap=iq
+    endif
+    if (noms(iq).eq."hdo_ice") then
+      i_hdo_ice=iq
+    endif
   enddo
+
   
   if (nq.gt.0) then
@@ -264 +275 @@
         endif ! of if (igcm_h2o_ice.eq.igcm_h2o_vap)
       endif ! of if (txt.eq."h2o_vap")
+
+      if (txt.eq."hdo_vap") then
+        write(*,*)"physdem1: skipping HDO vapour tracer"
+        if (i_hdo_ice.eq.i_hdo_vap) then
+          ! then there is no "water ice" tracer; but still
+          ! there is some water ice on the surface
+          write(*,*)"          writing HDO ice instead"
+          txt="hdo_ice"
+        else
+          ! there is a "water ice" tracer which has been / will be
+          ! delt with in due time
+          cycle
+        endif ! of if (igcm_hdo_ice.eq.igcm_hdo_vap)
+      endif ! of if (txt.eq."hdo_vap")
+
       call put_field(trim(txt),"tracer on surface",qsurf(:,iq),time)
     enddo