Changeset 3291 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Apr 3, 2024, 8:52:18 PM (11 months ago)

Author:

mturbet

Message:

improved parameterization of sea ice albedo vs thickness

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• hydrol.F90 (modified) (2 diffs)
• ocean_slab_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/hydrol.F90

@@ -15 +15 @@
   USE tracer_h
 !  use slab_ice_h
-  USE ocean_slab_mod, ONLY: alb_ice_max,alb_ice_min,h_alb_ice,snow_min
+  USE ocean_slab_mod, ONLY: alb_ice_min,h_alb_ice,snow_min
   use callkeys_mod, only: albedosnow,alb_ocean,albedoco2ice,ok_slab_ocean,Tsaldiff,maxicethick,co2cond
   use radinc_h, only : L_NSPECTV
@@ -221 +221 @@
                frac_snow = MAX(0.0,MIN(1.0,zqsurf(ig,iice)/snow_min)) ! Critical snow height (in kg/m2) from ocean_slab_ice routine. 
                ! Standard value should be 15kg/m2 (i.e. about 5 cm). Note that in the previous ocean param. (from BC2014), this value was 45kg/m2 (i.e. about 15cm).
-               alb_ice=alb_ice_max-(alb_ice_max-alb_ice_min)*exp(-sea_ice(ig)/h_alb_ice) ! this is the old formulation from BC2014 (earth-centric so will be replaced)
+               
                ! Albedo final calculation :
                do nw=1,L_NSPECTV
+                  alb_ice=albedo_snow_SPECTV(nw)-(albedo_snow_SPECTV(nw)-alb_ice_min)*exp(-sea_ice(ig)/h_alb_ice) ! this replaces the formulation from BC2014
+                  ! More details on the parameterization of sea ice albedo vs thickness is provided in the wiki : 
+                  ! https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Slab_ocean_model
+                  ! sea_ice is the ice thickness (calculated in ocean_slab routine) in kg/m2 ; h_alb_ice is fixed to 275.1kg/m2 i.e. 30cm based on comparisons with Brandt et al. 2005
                   albedo(ig,nw) = pctsrf_sic(ig)*                                        &
                                  (albedo_snow_SPECTV(nw)*frac_snow + alb_ice*(1.0-frac_snow))      &

trunk/LMDZ.GENERIC/libf/phystd/ocean_slab_mod.F90

@@ -131 +131 @@
     REAL, PARAMETER, PUBLIC :: capcalsno=2.3867e+06*0.15
 
-    REAL, PARAMETER, PUBLIC :: h_alb_ice=0.5*ice_den ! height for full ice albedo
-    REAL, PARAMETER, PUBLIC :: h_sno_alb=0.02*sno_den ! height for control of snow fraction
-
-    REAL, PARAMETER, PUBLIC :: alb_ice_min=0.2
-    REAL, PARAMETER, PUBLIC :: alb_ice_max=0.65
+    REAL, PARAMETER, PUBLIC :: h_alb_ice=0.3*ice_den ! height (in kg/m2) used in the calculation of sea ice albedo vs thickness
+    ! (changed from 50cm to 30cm based on comparisons with Brandt et al. 2005) ; more info in the slab ocean wiki page
+    REAL, PARAMETER, PUBLIC :: h_sno_alb=0.02*sno_den ! height (in kg/m2) for control of snow fraction
+
+    REAL, PARAMETER, PUBLIC :: alb_ice_min=0.08 ! minimum sea ice albedo used for calculation of albedo as a function of sea ice thickness (https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Slab_ocean_model)
 
 ! Horizontal transport parameters