Changeset 2378 for trunk/LMDZ.MARS/libf

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

Files:

• dynphy_lonlat/phymars/newstart.F (modified) (1 diff)
• phymars/hdo_surfex_mod.F (modified) (1 diff)
• phymars/phyetat0_mod.F90 (modified) (1 diff)
• phymars/physiq_mod.F (modified) (1 diff)
• phymars/simpleclouds.F (modified) (1 diff)

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

@@ -1010 +1010 @@
         !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
+
+        ! Update the hdo tracers
+         q(1:iip1,1:jjp1,1:llm,igcm_hdo_vap)
+     &      =q(1:iip1,1:jjp1,1:llm,igcm_h2o_vap)* DoverH
+         q(1:iip1,1:jjp1,1:llm,igcm_hdo_ice)
+     &      =q(1:iip1,1:jjp1,1:llm,igcm_h2o_ice)* DoverH
+
+         qsurf(1:ngridmx,igcm_hdo_ice)
+     &      =qsurf(1:ngridmx,igcm_h2o_ice)*DoverH
 
 

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