Changeset 2378 for trunk/LMDZ.MARS/libf/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

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• hdo_surfex_mod.F (modified) (1 diff)
• phyetat0_mod.F90 (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)
• simpleclouds.F (modified) (1 diff)

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

@@ -105 +105 @@
 
                if (hdofrac) then !do we use fractionation?
-                alpha_c(ig) = exp(16288./zt(ig,1)**2.-9.34e-2)
-c               alpha_c = exp(13525./zt(ig,l)**2.-5.59e-2) !Lamb
+c               alpha_c(ig) = exp(16288./zt(ig,1)**2.-9.34e-2)
+                alpha_c = exp(13525./zt(ig,1)**2.-5.59e-2) !Lamb
                else
                 alpha_c(ig) = 1.

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

@@ -609 +609 @@
          end do
        endif
+
+       if (txt.eq."hdo_ice") then
+         do ig=1,ngrid
+          if (qsurf(ig,iq).lt.0.0) then
+             qsurf(ig,iq) = 0.0
+          end if
+         end do
+       endif
+
       enddo
      endif ! of if (nq.ge.1)

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

@@ -3032 +3032 @@
 
                 do ig=1,ngrid
-                if (qsurf(ig,igcm_h2o_ice).gt.1e-30) then
+                if (qsurf(ig,igcm_h2o_ice).gt.qperemin) then
                     DoH_surf(ig) = 0.5*( qsurf(ig,igcm_hdo_ice)/
      &                  qsurf(ig,igcm_h2o_ice) )/155.76e-6

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

@@ -172 +172 @@
 
                 if (hdofrac) then ! do we use fractionation?
-                alpha_c(ig,l) = exp(16288./zt(ig,l)**2.-9.34d-2)
-c               alpha_c(ig,l) = exp(13525./zt(ig,l)**2.-5.59d-2)  !Lamb
+c               alpha_c(ig,l) = exp(16288./zt(ig,l)**2.-9.34d-2)
+                alpha_c(ig,l) = exp(13525./zt(ig,l)**2.-5.59d-2)  !Lamb
                 else
                 alpha_c(ig,l) = 1.d0