Update documentation including interfaces

columnphysics/icepack_parameters.F90

@@ -117,8 +117,8 @@ module icepack_parameters
          Timelt    = 0.0_dbl_kind     ,&! melting temperature, ice top surface  (C)
          Tsmelt    = 0.0_dbl_kind     ,&! melting temperature, snow top surface (C)
          ice_ref_salinity =4._dbl_kind,&! (ppt)
-                                        ! kice is not used for mushy thermo
          kice      = 2.03_dbl_kind    ,&! thermal conductivity of fresh ice(W/m/deg)
+                                        ! kice is only used with ktherm=1 (BL99) and conduct='MU71'
          ksno      = 0.30_dbl_kind    ,&! thermal conductivity of snow  (W/m/deg)
          hs_min    = 1.e-4_dbl_kind   ,&! min snow thickness for computing zTsn (m)
          snowpatch = 0.02_dbl_kind    ,&! parameter for fractional snow area (m)

columnphysics/icepack_therm_shared.F90

@@ -216,7 +216,7 @@ end subroutine dsurface_heat_flux_dTsf
 
 !=======================================================================
 !autodocument_start icepack_init_salinity
-! Initialize the vertical profile of ice salinity and melting temperature.
+! Initialize the vertical profile of ice salinity.
 ! This subroutine was renamed from icepack_init_thermo in Oct 2024
 !
 ! authors: C. M. Bitz, UW
@@ -245,7 +245,7 @@ subroutine icepack_init_salinity(sprofile)
       if (saltmax > min_salin) l_brine = .true.
 
       !-----------------------------------------------------------------
-      ! Prescibe vertical profile of salinity and melting temperature.
+      ! Prescibe vertical profile of salinity.
       ! Note this profile is only used for BL99 thermodynamics.
       !-----------------------------------------------------------------
 
@@ -298,7 +298,7 @@ end function icepack_salinity_profile
 !autodocument_start icepack_init_enthalpy
 ! This subroutine was renamed from icepack_init_trcr in Oct 2024
 !
-      subroutine icepack_init_enthalpy(Tair,     Tf,       &
+      subroutine icepack_init_enthalpy(Tair, Tf,      &
                                   Sprofile, Tprofile, &
                                   Tsfc,               &
                                   qin,      qsn)

doc/source/icepack_index.rst

@@ -249,7 +249,7 @@ section :ref:`tabnamelist`.
    "kappav", "visible extinction coefficient in ice, wavelength\ :math:`<`\ 700nm", "1.4 m\ :math:`^{-1}`"
    "kcatbound", ":math:`\bullet` category boundary formula", ""
    "kg_to_g", "kg to g conversion factor", "1000."
-   "kice", "thermal conductivity of fresh ice (:cite:`Bitz99`)", "2.03 W/m/deg"
+   "kice", "thermal conductivity of fresh ice with BL99+MU71 (:cite:`Bitz99`)", "2.03 W/m/deg"
    "kitd", ":math:`\bullet` type of itd conversions (0 = delta function, 1 = linear remap)", "1"
    "krdg_partic", ":math:`\bullet` ridging participation function", "1"
    "krdg_redist", ":math:`\bullet` ridging redistribution function", "1"