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

Timestamp:

Sep 13, 2021, 12:19:38 PM (3 years ago)

Author:

jnaar

Message:

Water cycle flag update & addition :

• "cap_albedo" renamed "cst_cap_albedo" (default false) : if true, water ice cap albedo remains constant even when frost with higher albedo condensates on it
• "refill_watercap" added (default false) : turns h2o_ice_s into watercap when above a given threshold
• "frost_metam_threshold" added (default 0.05 km.m-2) : threshold used by "refill_watercap"

JN

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys.h (modified) (2 diffs)
• conf_phys.F (modified) (2 diffs)
• nuclea.F (modified) (1 diff)
• physiq_mod.F (modified) (4 diffs)
• surfdat_h.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -17 +17 @@
      &   ,CLFvarying,satindexco2,rdstorm,slpwind,calllott_nonoro        &
      &   ,latentheat_surfwater,gwd_convective_source,startphy_file      &
-     &   ,hdo,hdofrac,cap_albedo,temp_dependant_m
+     &   ,hdo,hdofrac,cst_cap_albedo,temp_dependant_m,refill_watercap
      
       COMMON/callkeys_i/iradia,iaervar,ilwd,ilwb,ilwn,ncouche           &
@@ -69 +69 @@
       logical slpwind ! entrainment by slope wind parametrization
       logical latentheat_surfwater ! latent heat release from ground water ice sublimation/condensation
-      logical cap_albedo ! polar cap albedo remains unchanged by water frost deposition
+      logical cst_cap_albedo ! polar cap albedo remains unchanged by water frost deposition
       logical temp_dependant_m ! temperature-dependant water contact parameter
+      logical refill_watercap ! h2o_ice_s is converted to watercap when above threshold
       logical sedimentation
       logical activice,tifeedback,supersat,caps

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

@@ -37 +37 @@
      &                       nuice_ref,nuiceco2_ref
       use surfdat_h, only: albedo_h2o_cap,albedo_h2o_frost,
-     &                     frost_albedo_threshold, inert_h2o_ice
+     &                     frost_albedo_threshold, inert_h2o_ice,
+     &                     frost_metam_threshold
       use time_phylmdz_mod, only: ecritphy,day_step,iphysiq,ecritstart,
      &                            daysec,dtphys
@@ -805 +806 @@
          write(*,*)"Watercaptag albedo is unchanged by water frost",
      &              " deposition (default is false)"
-         cap_albedo=.false. ! default value
-         call getin_p("cap_albedo",cap_albedo)
-         write(*,*)"cap_albedo = ",cap_albedo
+         cst_cap_albedo=.false. ! default value
+         call getin_p("cst_cap_albedo",cst_cap_albedo)
+         write(*,*)"cst_cap_albedo = ",cst_cap_albedo
+
+! Watercap evolution & refill (with discriminated albedos) (JN 2021)
+         write(*,*)"Watercap is replenished by water frost",
+     &              " accumulation (default is false)"
+         refill_watercap=.false. ! default value
+         call getin_p("refill_watercap",refill_watercap)
+         write(*,*)"refill_watercap= ",refill_watercap
+! frost thickness threshold for refill_watercap (ice metamorphism)
+         write(*,*) "frost thickness threshold for metamorphism ?",
+     &              "ie converted into watercap",
+     &              "only if refill_watercap is .true."
+         frost_metam_threshold=0.05 !  (i.e 0.05 kg.m-2)
+         call getin_p("frost_metam_threshold",
+     &    frost_metam_threshold)
+         write(*,*) " frost_metam_threshold = ",
+     &            frost_metam_threshold
+
 
 ! inert_h2o_ice

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

@@ -57 +57 @@
 
       IF (temp_dependant_m) THEN
+c     Simple linear parametrisation from Maattaanen 2014
+c     Smectite sample
+c     Maxed out at 0.97 for physical realism
          mtetalocal = min(0.0044*temp + 0.1831,0.97)
       ENDIF ! (temp_dependant_m) THEN

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

@@ -46 +46 @@
       use surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam,
      &                     zthe, z0, albedo_h2o_cap,albedo_h2o_frost,
-     &                     frost_albedo_threshold,
+     &                     frost_albedo_threshold,frost_metam_threshold,
      &                     tsurf, co2ice, emis,
      &                     capcal, fluxgrd, qsurf,
@@ -2078 +2078 @@
            if ((co2ice(ig).eq.0).and.
      &        (qsurf(ig,igcm_h2o_ice).gt.frost_albedo_threshold)) then
-             if ((watercaptag(ig)).and.(cap_albedo)) then
+             if ((watercaptag(ig)).and.(cst_cap_albedo)) then
                albedo(ig,1) = albedo_h2o_cap
                albedo(ig,2) = albedo_h2o_cap
@@ -2084 +2084 @@
                albedo(ig,1) = albedo_h2o_frost
                albedo(ig,2) = albedo_h2o_frost
-             endif !((watercaptag(ig)).and.(cap_albedo)) then
+             endif !((watercaptag(ig)).and.(cst_cap_albedo)) then
 c              write(*,*) "frost thickness", qsurf(ig,igcm_h2o_ice)
 c              write(*,*) "physiq.F frost :"
@@ -2151 +2151 @@
          watercap(ig)=watercap(ig)+ptimestep*dwatercap(ig)
       ENDDO
+
+      IF (refill_watercap) THEN
+
+        DO ig=1,ngrid
+           if (watercaptag(ig).and.
+     &        (qsurf(ig,igcm_h2o_ice).gt.frost_metam_threshold)) then
+
+                watercap(ig)=watercap(ig)+qsurf(ig,igcm_h2o_ice)
+     &                  - frost_metam_threshold
+                qsurf(ig,igcm_h2o_ice) = frost_metam_threshold
+           endif ! (watercaptag(ig).and.
+        ENDDO
+
+      ENDIF ! (refill_watercap) THEN
+
 
 c-----------------------------------------------------------------------

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

@@ -17 +17 @@
   real,save :: inert_h2o_ice ! water ice thermal inertia
   real,save :: frost_albedo_threshold ! water frost thickness on the ground (kg.m^-2, ie mm)
+  real,save :: frost_metam_threshold ! water frost threshold before conversion to ice (kg.m^-2, ie mm)
   real,save :: TESice_Ncoef ! coefficient for TES ice albedo in Northern hemisphere
   real,save :: TESice_Scoef ! coefficient for TES ice albedo in Southern hemisphere