new ice phy props for SNICAR-ADv4

*********1*********2*********3*********4*********5*********6*********7** functional v of snicar-adv4 in elm, new ice physical props that match measurements better and have better "control" fill albedo physical properties this version was used for simulations in the manuscript
chloewhicker · Sep 6, 2023 · f34ca9e · f34ca9e
1 parent 4247ac3
commit f34ca9e
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 10 deletions.
diff --git a/components/elm/bld/namelist_files/namelist_defaults.xml b/components/elm/bld/namelist_files/namelist_defaults.xml
@@ -587,7 +587,7 @@ this mask will have smb calculated over the entire global land surface
 <!-- SNICAR (SNow, ICe, and Aerosol Radiative model) datasets -->
 <!--  ***********  Resolution independent:   ***********  -->
 <fsnowoptics >lnd/clm2/snicardata/snicar_optics_5bnd_mam_cCAW20230224_nc3.nc</fsnowoptics>
-<ficephyprop >lnd/clm2/snicardata/jja_icephyprop_2000_2021_ValFill_c20230825_nc3.nc</ficephyprop>
+<ficephyprop >lnd/clm2/snicardata/jja_icephyprop_2000_2021_ValFill_c20230826_nc3.nc</ficephyprop>
 <fsnowaging  >lnd/clm2/snicardata/snicar_drdt_bst_fit_60_c070416.nc</fsnowaging>
 <use_snicar_ad>.true.</use_snicar_ad>
 <snow_shape>sphere</snow_shape>

diff --git a/components/elm/src/biogeophys/SnowSnicarMod.F90 b/components/elm/src/biogeophys/SnowSnicarMod.F90
@@ -2639,18 +2639,18 @@ subroutine SNICAR_AD_RT (flg_snw_ice, bounds, num_nourbanc, filter_nourbanc,  &
              elseif(numrad_snw==5) then
                 ! Direct:
                 if (flg_slr_in == 1) then
-                        flx_wgt(1) = 1._r8
-                        flx_wgt(2)  = 0.502936511226529_r8 + 0.05_r8
-                            flx_wgt(3)  = 0.180787216133052_r8 + 0.05_r8
-                            flx_wgt(4)  = 0.116700734417544_r8 - 0.05_r8
-                            flx_wgt(5)  = 0.199547150388836_r8 - 0.05_r8                                            !+CAW
+                            flx_wgt(1) = 1._r8
+                            flx_wgt(2)  = 0.502936511226529_r8 + 0.065_r8
+                            flx_wgt(3)  = 0.180787216133052_r8 + 0.065_r8
+                            flx_wgt(4)  = 0.116700734417544_r8 - 0.065_r8
+                            flx_wgt(5)  = 0.199547150388836_r8 - 0.065_r8                                            !+CAW
                    ! Diffuse:
                 elseif (flg_slr_in == 2) then
                             flx_wgt(1) = 1._r8
-                            flx_wgt(2)  = 0.633177019092330_r8 + 0.05_r8
-                            flx_wgt(3)  = 0.186664819863994_r8 + 0.05_r8
-                            flx_wgt(4)  = 0.093952188565569_r8 - 0.05_r8
-                            flx_wgt(5)  = 0.086205637586152_r8 - 0.05_r8
+                            flx_wgt(2)  = 0.633177019092330_r8 + 0.065_r8
+                            flx_wgt(3)  = 0.186664819863994_r8 + 0.065_r8
+                            flx_wgt(4)  = 0.093952188565569_r8 - 0.065_r8
+                            flx_wgt(5)  = 0.086205637586152_r8 - 0.065_r8
                 endif
              endif ! end if numrad_snw