Changeset 3291 for trunk/LMDZ.GENERIC
- Timestamp:
- Apr 3, 2024, 8:52:18 PM (9 months ago)
- Location:
- trunk/LMDZ.GENERIC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/changelog.txt
r3270 r3291 1891 1891 1892 1892 - physiq_mod.F90 cleaning + add moist and dry adjustment heating tendencies 1893 1894 == 03/04/2024 == MT 1895 1896 - improvement of parameterization of sea ice albedo vs sea ice thickness in hydrol+photolysis_online 1897 (when ok_slab_ocean activated) -
trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F
r3133 r3291 1173 1173 use chimiedata_h, only: albedo_snow_chim, albedo_co2_ice_chim 1174 1174 ! use slab_ice_h, only: h_alb_ice, alb_ice_min, alb_ice_max 1175 use ocean_slab_mod, only: h_alb_ice 1176 use ocean_slab_mod, only: alb_ice_min 1177 use ocean_slab_mod, only: alb_ice_max 1175 use ocean_slab_mod, only: h_alb_ice,alb_ice_min, snow_min 1178 1176 use tracer_h, only: igcm_h2o_ice, igcm_co2_ice 1179 1177 use callkeys_mod, only: ok_slab_ocean, co2cond, alb_ocean, … … 1222 1220 if(ok_slab_ocean) then 1223 1221 1224 zfra = MAX(0.0,MIN(1.0,qsurf(ig,igcm_h2o_ice)/ 45.0)) ! Snow Fraction (Critical height 45kg/m2~15cm)1225 alb_ice=alb_ice_max-(alb_ice_max-alb_ice_min) ! Ice Albedo 1226 & *exp(-sea_ice(ig)/h_alb_ice) 1222 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. 1223 ! 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). 1224 1227 1225 ! Albedo final calculation : 1228 1226 do iw=1,nw - 1 1227 alb_ice=albedo_snow_chim(iw) 1228 & -(albedo_snow_chim(iw)-alb_ice_min) 1229 & *exp(-sea_ice(ig)/h_alb_ice) ! this replaces the formulation from BC2014 1230 ! More details on the parameterization of sea ice albedo vs thickness is provided in the wiki : 1231 ! https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Slab_ocean_model 1232 ! 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 1229 1233 albedo_chim(iw) = pctsrf_sic(ig)* 1230 1234 & (albedo_snow_chim(iw)*zfra -
trunk/LMDZ.GENERIC/libf/phystd/hydrol.F90
r3270 r3291 15 15 USE tracer_h 16 16 ! use slab_ice_h 17 USE ocean_slab_mod, ONLY: alb_ice_m ax,alb_ice_min,h_alb_ice,snow_min17 USE ocean_slab_mod, ONLY: alb_ice_min,h_alb_ice,snow_min 18 18 use callkeys_mod, only: albedosnow,alb_ocean,albedoco2ice,ok_slab_ocean,Tsaldiff,maxicethick,co2cond 19 19 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. 222 222 ! 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). 223 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)223 224 224 ! Albedo final calculation : 225 225 do nw=1,L_NSPECTV 226 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 227 ! More details on the parameterization of sea ice albedo vs thickness is provided in the wiki : 228 ! https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Slab_ocean_model 229 ! 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 226 230 albedo(ig,nw) = pctsrf_sic(ig)* & 227 231 (albedo_snow_SPECTV(nw)*frac_snow + alb_ice*(1.0-frac_snow)) & -
trunk/LMDZ.GENERIC/libf/phystd/ocean_slab_mod.F90
r3268 r3291 131 131 REAL, PARAMETER, PUBLIC :: capcalsno=2.3867e+06*0.15 132 132 133 REAL, PARAMETER, PUBLIC :: h_alb_ice=0. 5*ice_den ! height for full ice albedo134 REAL, PARAMETER, PUBLIC :: h_sno_alb=0.02*sno_den ! height for control of snow fraction135 136 REAL, PARAMETER, PUBLIC :: alb_ice_min=0.2 137 REAL, PARAMETER, PUBLIC :: alb_ice_m ax=0.65133 REAL, PARAMETER, PUBLIC :: h_alb_ice=0.3*ice_den ! height (in kg/m2) used in the calculation of sea ice albedo vs thickness 134 ! (changed from 50cm to 30cm based on comparisons with Brandt et al. 2005) ; more info in the slab ocean wiki page 135 REAL, PARAMETER, PUBLIC :: h_sno_alb=0.02*sno_den ! height (in kg/m2) for control of snow fraction 136 137 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) 138 138 139 139 ! Horizontal transport parameters
Note: See TracChangeset
for help on using the changeset viewer.