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

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/dyn1d/testphys1d.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -117 +117 @@
       integer :: nq=1 ! number of tracers
       real :: latitude(1), longitude(1), cell_area(1)
+      integer iqh2ovap
+      integer iqh2oice
 
       character*2 str2
@@ -294 +296 @@
           ! WATER VAPOUR
           if (txt.eq."h2o_vap") then
+             iqh2ovap=iq !remember index for water vap
             !look for a "profile_h2o_vap" input file
             open(91,file='profile_h2o_vap',status='old',
@@ -309 +312 @@
           ! WATER ICE
           if (txt.eq."h2o_ice") then
+                iqh2oice = iq !remember index for water ice
             !look for a "profile_h2o_vap" input file
             open(91,file='profile_h2o_ice',status='old',
@@ -322 +326 @@
             close(91)
           endif ! of if (txt.eq."h2o_ice")
+
+          ! HDO VAPOUR
+          if (txt.eq."hdo_vap") then
+            !look for a "profile_hdo_vap" input file
+            open(91,file='profile_hdo_vap',status='old',
+     &       form='formatted',iostat=ierr)
+            if (ierr.eq.0) then
+              read(91,*) qsurf(iq)
+              do ilayer=1,nlayer
+                read(91,*) q(ilayer,iq)
+              enddo
+            else
+              write(*,*) "No profile_hdo_vap file!"
+              do ilayer=1,nlayer
+                q(ilayer,iq) =  q(ilayer,iqh2ovap)*2*155.76e-6*5 
+              enddo
+            endif
+            close(91)
+          endif ! of if (txt.eq."hdo_ice")
+          ! HDO ICE
+          if (txt.eq."hdo_ice") then
+            !look for a "profile_hdo_vap" input file
+            open(91,file='profile_hdo_ice',status='old',
+     &       form='formatted',iostat=ierr)
+            if (ierr.eq.0) then
+              read(91,*) qsurf(iq)
+              do ilayer=1,nlayer
+                read(91,*) q(ilayer,iq)
+              enddo
+            else
+              write(*,*) "No profile_hdo_ice file!"
+                qsurf(iq) = qsurf(iqh2oice) * 2*155.76e-6*5
+               do ilayer=1,nlayer
+                q(ilayer,iq) = q(ilayer,iqh2oice) * 2*155.76e-6*5 
+              enddo
+            endif
+            close(91)
+          endif ! of if (txt.eq."hdo_ice")
+
+
           ! DUST
           !if (txt(1:4).eq."dust") then