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

Timestamp:

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

Author:

mturbet

Message:

improved parameterization of sea ice albedo vs thickness

File:

• trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F

@@ -1173 +1173 @@
       use chimiedata_h,  only: albedo_snow_chim, albedo_co2_ice_chim
 !      use slab_ice_h,    only: h_alb_ice, alb_ice_min, alb_ice_max
-      use ocean_slab_mod, only: h_alb_ice
-      use ocean_slab_mod, only: alb_ice_min
-      use ocean_slab_mod, only: alb_ice_max
+      use ocean_slab_mod, only: h_alb_ice,alb_ice_min, snow_min
       use tracer_h,      only: igcm_h2o_ice, igcm_co2_ice
       use callkeys_mod,  only: ok_slab_ocean, co2cond, alb_ocean,
@@ -1222 +1220 @@
           if(ok_slab_ocean) then
         
-            zfra = MAX(0.0,MIN(1.0,qsurf(ig,igcm_h2o_ice)/45.0))     ! Snow Fraction (Critical height 45kg/m2~15cm)
-            alb_ice=alb_ice_max-(alb_ice_max-alb_ice_min)            ! Ice Albedo
-     &                          *exp(-sea_ice(ig)/h_alb_ice)
+            zfra = MAX(0.0,MIN(1.0,qsurf(ig,igcm_h2o_ice)/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).
+
             ! Albedo final calculation :
             do iw=1,nw - 1
+               alb_ice=albedo_snow_chim(iw)
+     &         -(albedo_snow_chim(iw)-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_chim(iw) = pctsrf_sic(ig)*
      &                          (albedo_snow_chim(iw)*zfra