From 3b211c6bcf16dfbc376c4d98d3a27e33424799f9 Mon Sep 17 00:00:00 2001 From: apcraig Date: Tue, 19 Jan 2021 11:36:26 -0700 Subject: [PATCH] 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 | 4 ++-- configuration/scripts/ice_in | 4 ++-- configuration/scripts/options/set_nml.zsal | 1 + doc/source/cice_index.rst | 2 +- doc/source/user_guide/ug_case_settings.rst | 2 +- icepack | 2 +- 6 files changed, 8 insertions(+), 7 deletions(-) diff --git a/cicecore/cicedynB/general/ice_init.F90 b/cicecore/cicedynB/general/ice_init.F90 index a16cd7649..3670e174d 100644 --- a/cicecore/cicedynB/general/ice_init.F90 +++ b/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' diff --git a/configuration/scripts/ice_in b/configuration/scripts/ice_in index 64060d245..1721541d4 100644 --- a/configuration/scripts/ice_in +++ b/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. diff --git a/configuration/scripts/options/set_nml.zsal b/configuration/scripts/options/set_nml.zsal index 7e8ed47ff..c099fed8c 100644 --- a/configuration/scripts/options/set_nml.zsal +++ b/configuration/scripts/options/set_nml.zsal @@ -1,3 +1,4 @@ +nblyr = 7 ktherm = 1 sw_redist = .true. tfrz_option = 'linear_salt' diff --git a/doc/source/cice_index.rst b/doc/source/cice_index.rst index 8ea16261d..caaf87d5a 100644 --- a/doc/source/cice_index.rst +++ b/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`" diff --git a/doc/source/user_guide/ug_case_settings.rst b/doc/source/user_guide/ug_case_settings.rst index 3cf271386..a14e571dc 100644 --- a/doc/source/user_guide/ug_case_settings.rst +++ b/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``" diff --git a/icepack b/icepack index af7e38f81..349b3f49a 160000 --- a/icepack +++ b/icepack @@ -1 +1 @@ -Subproject commit af7e38f8107b2ec8680cf5a49d7274b43fe34248 +Subproject commit 349b3f49af407be06f9f37014883ecae0daed5d1