Changeset 2508

Timestamp:

Apr 28, 2021, 5:38:26 PM (4 years ago)

Author:

jnaar

Message:

The water ice albedo is now conceptually decoupled between old perennial ice
and seasonal frost. For now, "old ice" is only on watercaptag = .true.
The variable albedo_h2o_ice has been decomposed into albedo_h2o_cap and
albedo_h2o_frost. Default values are both 0.35 and retrocompatible with
old callphys.def and albedo_h2o_ice if needed.
JN

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/albedocaps.F90 (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/physiq_mod.F (modified) (3 diffs)
• libf/phymars/surfdat_h.F90 (modified) (1 diff)
• libf/phymars/surfini.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3344 +3344 @@
 of the simulation and the remaining time is in Time (often=0), if we write 
 multiple restart nothing changes 
+
+== 22/04/2021 == JN
+ - The variable albedo_h2o_ice is now decomposed into albedo_h2o_cap 
+and albedo_h2o_frost, in order to discriminate between perennial ice
+and seasonnal frost. Retrocompatible with old callphys.def and 
+albedo_h2o_ice. Current default values are 0.35 for both albedos.
+

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

@@ -7 +7 @@
 use geometry_mod, only: latitude ! grid point latitudes (rad)
 use surfdat_h, only: TESicealbedo, TESice_Ncoef, TESice_Scoef, &
-                     emisice, albedice, watercaptag, albedo_h2o_ice, &
+                     emisice, albedice, watercaptag, albedo_h2o_cap, &
                      emissiv, albedodat
 implicit none
@@ -72 +72 @@
     ! to do : emissivity
     emisref(ig) = 1
-    psolaralb(ig,1)=albedo_h2o_ice
-    psolaralb(ig,2)=albedo_h2o_ice
+    psolaralb(ig,1)=albedo_h2o_cap
+    psolaralb(ig,2)=albedo_h2o_cap
   else
     ! set emissivity of surface to be bare ground emissivity

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

@@ -36 +36 @@
       use tracer_mod, only : nuice_sed, ccn_factor, nuiceco2_sed,
      &                       nuice_ref,nuiceco2_ref
-      use surfdat_h, only: albedo_h2o_ice, inert_h2o_ice,
-     &                     frost_albedo_threshold
+      use surfdat_h, only: albedo_h2o_cap,albedo_h2o_frost,
+     &                     frost_albedo_threshold, inert_h2o_ice
       use time_phylmdz_mod, only: ecritphy,day_step,iphysiq,ecritstart,
      &                            daysec,dtphys
@@ -778 +778 @@
          write(*,*) " caps = ",caps
 
-! albedo_h2o_ice
+! JN : now separated between albedo_h2o_cap and 
+! albedo_h2o_frost. Retrocompatible with old
+! callphys.def with albedo_h2o_ice
          write(*,*) "water ice albedo ?"
-         albedo_h2o_ice=0.45
-         call getin_p("albedo_h2o_ice",albedo_h2o_ice)
-         write(*,*) " albedo_h2o_ice = ",albedo_h2o_ice
+         albedo_h2o_cap=0.35
+         albedo_h2o_frost=0.35
+         call getin_p("albedo_h2o_ice",albedo_h2o_cap)
+         albedo_h2o_frost=albedo_h2o_cap
+         call getin_p("albedo_h2o_cap",albedo_h2o_cap)
+         write(*,*) " albedo_h2o_cap = ",albedo_h2o_cap
+         call getin_p("albedo_h2o_frost",albedo_h2o_frost)
+         write(*,*) " albedo_h2o_frost = ",albedo_h2o_frost
 ! inert_h2o_ice
          write(*,*) "water ice thermal inertia ?"

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

@@ -45 +45 @@
       use comgeomfi_h, only: sinlon, coslon, sinlat, coslat
       use surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam,
-     &                     zthe, z0, albedo_h2o_ice,
+     &                     zthe, z0, albedo_h2o_cap,albedo_h2o_frost,
      &                     frost_albedo_threshold,
      &                     tsurf, co2ice, emis,
@@ -680 +680 @@
 
         IF (tracer.AND.water.AND.(ngrid.NE.1)) THEN
-          write(*,*)"physiq: water_param Surface water ice albedo:", 
-     .                  albedo_h2o_ice
+          write(*,*)"physiq: water_param Surface water frost albedo:", 
+     .                  albedo_h2o_frost
+          write(*,*)"physiq: water_param Surface watercap albedo:", 
+     .                  albedo_h2o_cap
         ENDIF
 
@@ -2070 +2072 @@
            if ((co2ice(ig).eq.0).and.
      &        (qsurf(ig,igcm_h2o_ice).gt.frost_albedo_threshold)) then
-              albedo(ig,1) = albedo_h2o_ice
-              albedo(ig,2) = albedo_h2o_ice
+              albedo(ig,1) = albedo_h2o_frost
+              albedo(ig,2) = albedo_h2o_frost
 c              write(*,*) "frost thickness", qsurf(ig,igcm_h2o_ice)
 c              write(*,*) "physiq.F frost :"

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

@@ -13 +13 @@
   logical,save :: temptag !temp tag for water caps
       
-  real,save :: albedo_h2o_ice ! water ice albedo
+  real,save :: albedo_h2o_cap ! water cap albedo
+  real,save :: albedo_h2o_frost ! water frost albedo
   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)

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

@@ -7 +7 @@
      &                     cell_area ! for watercaptag diagnosis
       use surfdat_h, only: watercaptag, frost_albedo_threshold,
-     &                     albedo_h2o_ice, inert_h2o_ice, albedodat,
+     &                     albedo_h2o_cap, inert_h2o_ice, albedodat,
      &                     albedice, dryness
 #ifndef MESOSCALE
@@ -84 +84 @@
                  watercaptag(ig)  = .true.
                  dryness(ig)      = 1.
-                 albedodat(ig)    = albedo_h2o_ice  !! pour output 
+                 albedodat(ig)    = albedo_h2o_cap  !! pour output 
          else
                  watercaptag(ig)  = .false.