Update emissivity to 0.985 and set nblyr default to 1

Change emissivity default to 0.985 (from 0.95)
Change nblyr default to 1 (from 7)
Add nblyr=7 to set_nml.zsal
Update Icepack to include emissivity default changes
Update documentation as needed

cicecore/cicedynB/general/ice_init.F90

@@ -368,7 +368,7 @@ subroutine input_data
       calc_Tsfc = .true.     ! calculate surface temperature
       update_ocn_f = .false. ! include fresh water and salt fluxes for frazil
       ustar_min = 0.005      ! minimum friction velocity for ocean heat flux (m/s)
-      emissivity = 0.95      ! emissivity of snow and ice
+      emissivity = 0.985     ! emissivity of snow and ice
       l_mpond_fresh = .false.     ! logical switch for including meltpond freshwater
                                   ! flux feedback to ocean model
       fbot_xfer_type = 'constant' ! transfer coefficient type for ocn heat flux
@@ -1436,7 +1436,7 @@ subroutine input_data
                   write(nu_diag,1002) ' ahmax            = ', ahmax,' : albedo is constant above this thickness'
                endif
             endif
-            write(nu_diag,1002) ' emissivity       = ', emissivity,' : emissivity of snow and ice'
+            write(nu_diag,1000) ' emissivity       = ', emissivity,' : emissivity of snow and ice'
             write(nu_diag,1010) ' sw_redist        = ', sw_redist,' : redistribute internal shortwave to surface'
             if (sw_redist) then
                write(nu_diag,1002) ' sw_frac          = ', sw_frac,' : fraction redistributed'

configuration/scripts/ice_in

@@ -63,7 +63,7 @@
     nfsd         = 1
     nilyr        = 7
     nslyr        = 1
-    nblyr        = 7
+    nblyr        = 1
     orca_halogrid = .false.
 /
 
@@ -195,7 +195,7 @@
     natmiter        = 5
     atmiter_conv    = 0.0d0
     ustar_min       = 0.0005
-    emissivity      = 0.95
+    emissivity      = 0.985
     fbot_xfer_type  = 'constant'
     update_ocn_f    = .false.
     l_mpond_fresh   = .false.

configuration/scripts/options/set_nml.zsal

@@ -1,3 +1,4 @@
+nblyr           = 7
 ktherm          = 1
 sw_redist       = .true.
 tfrz_option     = 'linear_salt'

doc/source/cice_index.rst

@@ -188,7 +188,7 @@ either Celsius or Kelvin units).
    "e11, e12, e22", "strain rate tensor components", ""
    "ecci", "yield curve minor/major axis ratio, squared", "1/4"
    "eice(n)", "energy of melting of ice per unit area (in category n)", "J/m\ :math:`^2`"
-   "emissivity", "emissivity of snow and ice", "0.95"
+   "emissivity", "emissivity of snow and ice", "0.985"
    "eps13", "a small number", "10\ :math:`^{-13}`"
    "eps16", "a small number", "10\ :math:`^{-16}`"
    "esno(n)", "energy of melting of snow per unit area (in category n)", "J/m\ :math:`^2`"

doc/source/user_guide/ug_case_settings.rst

@@ -492,7 +492,7 @@ forcing_nml
    "``calc_Tsfc``", "logical", "calculate surface temperature", "``.true.``"
    "``default_season``", "``summer``", "forcing initial summer values", "``winter``"
    "", "``winter``", "forcing initial winter values", ""
-   "``emissivity``", "real", "emissivity of snow and ice", "0.95"
+   "``emissivity``", "real", "emissivity of snow and ice", "0.985"
    "``fbot_xfer_type``", "``Cdn_ocn``", "variabler ocean heat transfer coefficient scheme", "``constant``"
    "", "``constant``", "constant ocean heat transfer coefficient", ""
    "``fe_data_type``", "``clim``", "ocean climatology forcing value for iron", "``default``"