Changeset 2965 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

May 17, 2023, 12:31:18 PM (20 months ago)

Author:

abierjon

Message:

Mars GCM:
Fixed r2963 which was preventing to compile the model without XIOS
+ changed the computation of variables rhowater_* so that they are real densities (factor 1/rvap missing ; this doesn't affect the previous PEM results as these densities were only compared between each other)
+ added comments and units for ice table variables in physiq_mod.F
+ made Clapeyron coefficient names in physiq_mod.F coherent with how they are defined in the PEM
+ fixed a reference in constants_marspem_mod.F90
+ fixed unit attribute of surface/soil water densities in field_def_physics_mars.xml

AB

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (3 diffs)

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

@@ -546 +546 @@
 
 ! Variable for ice table
-      REAL :: rhowater_surf(ngrid,nslope)
-      REAL :: rhowater_surf_sat(ngrid,nslope)
-      REAL :: rhowater_soil(ngrid,nsoilmx,nslope)
-      REAL,PARAMETER  :: alph_clap = -6143.7
-      REAL,PARAMETER :: beta_clap = 28.9074
-      REAL :: pvap_surf(ngrid)
-      REAL,PARAMETER :: m_co2 = 44.01E-3  ! CO2 molecular mass (kg/mol)   
-      REAL,PARAMETER :: m_noco2 = 33.37E-3  ! Non condensible mol mass (kg/mol)   
-      REAL :: ztmp1,ztmp2
+      REAL :: rhowater_surf(ngrid,nslope)         ! Water density at the surface [kg/m^3]
+      REAL :: rhowater_surf_sat(ngrid,nslope)     ! Water density at the surface at saturation [kg/m^3]
+      REAL :: rhowater_soil(ngrid,nsoilmx,nslope) ! Water density in soil layers [kg/m^3]
+      REAL,PARAMETER :: alpha_clap_h2o = 28.9074  ! Coeff for Clapeyron law [/]
+      REAL,PARAMETER :: beta_clap_h2o = -6143.7   ! Coeff for Clapeyron law [K]
+      REAL :: pvap_surf(ngrid)                    ! Water vapor partial pressure in first layer [Pa]
+      REAL,PARAMETER :: m_co2 = 44.01E-3          ! CO2 molecular mass [kg/mol]
+      REAL,PARAMETER :: m_noco2 = 33.37E-3        ! Non condensible mol mass [kg/mol]
+      REAL :: ztmp1,ztmp2                         ! intermediate variables to compute the mean molar mass of the layer
 
 c=======================================================================
@@ -3847 +3847 @@
 
             DO islope = 1,nslope
+            ! Clapeyron law for psat (psat = exp(beta/Th2o+alpha)),following Murphy and Koop 2005
              rhowater_surf_sat(ig,islope)  =
-     &         exp(alph_clap/tsurf(ig,islope)+beta_clap)
+     &         exp(beta_clap_h2o/tsurf(ig,islope)+alpha_clap_h2o)
      &         / tsurf(ig,islope)
+     &         * mmol(igcm_h2o_vap)/(mugaz*r)
  
              if(qsurf(ig,igcm_h2o_ice,islope).gt.(1.e-4)) then
-               rhowater_surf(ig,islope) =
-     &         exp(alph_clap/tsurf(ig,islope)+beta_clap)
-     &         / tsurf(ig,islope)
+               ! we consider to be at saturation above 1.e-4 kg.m-2
+               rhowater_surf(ig,islope) = rhowater_surf_sat(ig,islope)
              else
+              ! otherwise, use vapor partial pressure
               rhowater_surf(ig,islope) = pvap_surf(ig)
      &         / tsurf(ig,islope)
+     &         * mmol(igcm_h2o_vap)/(mugaz*r)
              endif
              DO isoil = 1,nsoilmx
              rhowater_soil(ig,isoil,islope) =
-     &         exp(alph_clap/tsoil(ig,isoil,islope)+beta_clap)
+     &         exp(beta_clap_h2o/tsoil(ig,isoil,islope)+alpha_clap_h2o)
      &         / tsoil(ig,isoil,islope)
+     &         * mmol(igcm_h2o_vap)/(mugaz*r)
              ENDDO
           ENDDO
@@ -3869 +3873 @@
       DO islope = 1,nslope
         write(str2(1:2),'(i2.2)') islope
-        CALL send_xios_field("Waterdensity_soil_slope"//str2,
+        CALL write_output("Waterdensity_soil_slope"//str2,
+     &     "rhowater_soil_slope"//str2,'kg.m-3',
      &     rhowater_soil(:,:,islope))
-        CALL send_xios_field("Waterdensity_surface"//str2,
+        CALL write_output("Waterdensity_surface"//str2,
+     &     "rhowater_surface"//str2,'kg.m-3',
      &     rhowater_surf(:,islope))
       ENDDO
 
-      CALL send_xios_field("h2o_layer1",zq(:,1,igcm_h2o_vap))
-      CALL send_xios_field("co2_layer1",zq(:,1,igcm_co2))
+      CALL write_output("h2o_layer1","h2o mass mr in the first layer",
+     &     'kg/kg',zq(:,1,igcm_h2o_vap))
+      CALL write_output("co2_layer1","co2 mass mr in the first layer",
+     &     'kg/kg',zq(:,1,igcm_co2))
 
 !PREVIOUSLY IN 1D ONLY