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

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/dynphy_lonlat/phymars/newstart.F (modified) (4 diffs)

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

@@ -24 +24 @@
      &                      igcm_ccn_number, igcm_ccn_mass,
      &                      igcm_h2o_vap, igcm_h2o_ice, igcm_co2,
+     &                      igcm_hdo_vap, igcm_hdo_ice,
      &                      igcm_n2, igcm_ar, igcm_o2, igcm_co,
      &                      igcm_o, igcm_h2
@@ -157 +158 @@
 
       INTEGER :: itau
+      real DoverH !D/H ratio
       
       INTEGER :: numvanle
@@ -514 +516 @@
       write(*,*) 'composition  : change atm main composition: CO2,N2,Ar,
      $ O2,CO'
+      write(*,*) 'inihdo       : initialize HDO'
       write(*,*) 'ini_h2osurf  : reinitialize surface water ice '
       write(*,*) 'noglacier    : Remove tropical H2O ice if |lat|<45'
@@ -986 +989 @@
           write(*,*) 'inichim_newstart: chemical species initialised
      $ (except water vapour)'
+
+c      inihdo : initialize HDO with user D/H value
+c      --------------------------------------------------------
+       else if (trim(modif) .eq. 'inihdo') then
+        ! check that there is indeed a water vapor tracer
+        if (igcm_h2o_vap.eq.0) then
+          write(*,*) "No water vapour tracer! Can't use this option"
+          stop
+        endif
+
+         write(*,*)'Input D/H ratio (in SMOW)'
+         write(*,*)'If value is <0 then HDO=H2O'
+ 303     read(*,*, iostat=ierr) DoverH
+         if(ierr.ne.0) goto 303
+
+        DoverH = DoverH * 2 * 155.76e-6 
+        if (DoverH.lt.0.0) then
+           DoverH = 1.
+        endif
+        !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
+
+
 
 c      composition : change main composition: CO2,N2,Ar,O2,CO (FF 03/2014)