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

Timestamp:

May 1, 2021, 5:36:56 PM (4 years ago)

Author:

jnaar

Message:

Added a water-cycle related flag : "cap_albedo" (default is false). When true, the albedo on watercaptag is unchanged by seasonal water frost deposition.
JN

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys.h (modified) (2 diffs)
• conf_phys.F (modified) (2 diffs)
• physiq_mod.F (modified) (2 diffs)

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

@@ -17 +17 @@
      &   ,CLFvarying,satindexco2,rdstorm,slpwind,calllott_nonoro        &
      &   ,latentheat_surfwater,gwd_convective_source,startphy_file      &
-     &   ,hdo,hdofrac
+     &   ,hdo,hdofrac,cap_albedo
      
       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 sedimentation
       logical activice,tifeedback,supersat,caps

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

@@ -781 +781 @@
 ! albedo_h2o_frost. Retrocompatible with old
 ! callphys.def with albedo_h2o_ice
-         write(*,*) "water ice albedo ?"
+         write(*,*) "water ice albedo ? Old settings use ",
+     &              "albedo_h2o_ice, new settings use ", 
+     &              "albedo_h2o_cap and albedo_h2o_frost " 
          albedo_h2o_cap=0.35
          albedo_h2o_frost=0.35
@@ -790 +792 @@
          call getin_p("albedo_h2o_frost",albedo_h2o_frost)
          write(*,*) " albedo_h2o_frost = ",albedo_h2o_frost
+
+! Northern polar cap albedo (JN 2021)
+         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
+
 ! inert_h2o_ice
          write(*,*) "water ice thermal inertia ?"

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

@@ -49 +49 @@
      &                     tsurf, co2ice, emis,
      &                     capcal, fluxgrd, qsurf,
-     &                     hmons,summit,base,watercap
+     &                     hmons,summit,base,watercap,watercaptag
       use comsaison_h, only: dist_sol, declin, mu0, fract, local_time
       use slope_mod, only: theta_sl, psi_sl
@@ -2072 +2072 @@
            if ((co2ice(ig).eq.0).and.
      &        (qsurf(ig,igcm_h2o_ice).gt.frost_albedo_threshold)) then
-              albedo(ig,1) = albedo_h2o_frost
-              albedo(ig,2) = albedo_h2o_frost
+             if ((watercaptag(ig)).and.(cap_albedo)) then
+               albedo(ig,1) = albedo_h2o_cap
+               albedo(ig,2) = albedo_h2o_cap
+             else
+               albedo(ig,1) = albedo_h2o_frost
+               albedo(ig,2) = albedo_h2o_frost
+             endif !((watercaptag(ig)).and.(cap_albedo)) then
 c              write(*,*) "frost thickness", qsurf(ig,igcm_h2o_ice)
 c              write(*,*) "physiq.F frost :"