E3SM-Project · eclare108213 · Nov 26, 2024 · Dec 17, 2024 · Jan 17, 2025 · Feb 13, 2025
diff --git a/.zenodo.json b/.zenodo.json
@@ -1,19 +1,21 @@
 {
-    "license": "other-open", 
-    "description": "No description provided", 
+    "license": "BSD-3-Clause",
+    "copyright": "Copyright 1998-2025 Triad National Security, LLC",
+    "description": "View detailed release notes at https://github.com/CICE-Consortium/Icepack/releases", 
     "language": "eng", 
-    "title": "CICE-Consortium/Icepack: Icepack m.n.p", 
+    "title": "CICE-Consortium/Icepack: Icepack 1.5.2", 
     "keywords": [
-        "sea ice model, Icepack"
+        "sea ice model",
+        "Icepack"
     ], 
-    "version": "m.n.p", 
+    "version": "1.5.2", 
     "upload_type": "software", 
     "communities": [
         {
             "identifier": "cice-consortium"
         }
     ], 
-    "publication_date": "2019-06-20", 
+    "publication_date": "2025-08-15", 
     "creators": [
         {
             "affiliation": "Los Alamos National Laboratory", 
@@ -31,6 +33,10 @@
             "affiliation": "Environment and Climate Change Canada", 
             "name": "Philippe Blain"
         }, 
+        {
+            "affiliation": "National Center for Atmospheric Research", 
+            "name": "David Clemens-Sewall"
+        }, 
         {
             "affiliation": "National Oceanic and Atmospheric Administration", 
             "name": "Anthony Craig"
@@ -75,20 +81,33 @@
             "affiliation": "Danish Meteorological Institute", 
             "name": "Mads Ribergaard"
         }, 
+        {
+            "affiliation": "Alfred-Wegener-Institut Helmholz-Zentrum fur Polar- und Meersforschung", 
+            "name": "Lettie Roach"
+        }, 
         {
             "affiliation": "Los Alamos National Laboratory", 
             "name": "Andrew Roberts"
         }, 
+        {
+            "affiliation": "Australian Earth-System Simulator", 
+            "name": "Anton Steketee"
+        }, 
         {
             "affiliation": "Los Alamos National Laboratory", 
             "name": "Matthew Turner"
         }, 
         {
             "affiliation": "Geophysical Fluid Dynamics Laboratory", 
             "name": "Michael Winton"
+        },
+        {
+            "affiliation": "National Aeronautics and Space Administration", 
+            "name": "Bin Zhao"
         }
     ], 
     "access_right": "open", 
+    "repository_url": "https://github.com/CICE-Consortium/Icepack", 
     "related_identifiers": [
         {
             "identifier": "https://github.com/CICE-Consortium/Icepack/tree/Icepack1.0.0", 

diff --git a/COPYRIGHT.pdf b/COPYRIGHT.pdf
diff --git a/DistributionPolicy.pdf b/DistributionPolicy.pdf
diff --git a/LICENSE.pdf b/LICENSE.pdf
diff --git a/columnphysics/icepack_algae.F90 b/columnphysics/icepack_algae.F90
@@ -394,8 +394,7 @@ subroutine zbio   (dt,                        &
                                bio_index,                &
                                aicen,                    &
                                vicen,        vsnon,      &
-                               iphin,                    &! ntrcr
-                               trcrn,        aice_old,   &!aice_old
+                               iphin,        trcrn,      &
                                flux_bion,    flux_bio,   &
                                upNOn,        upNHn,      &
                                upNO,         upNH,       &

diff --git a/columnphysics/icepack_atmo.F90 b/columnphysics/icepack_atmo.F90
@@ -507,8 +507,7 @@ end subroutine atmo_boundary_const
 !
 ! changes: Andrew Roberts, NPS (RASM/CESM coupling and documentation)
 
-      subroutine neutral_drag_coeffs (apnd,     hpnd,     &
-                                      ipnd,               &
+      subroutine neutral_drag_coeffs (apondn,               &
                                       alvl,     vlvl,     &
                                       aice,     vice,     &
                                       vsno,     aicen,    &
@@ -526,9 +525,7 @@ subroutine neutral_drag_coeffs (apnd,     hpnd,     &
       use icepack_tracers, only: tr_pond
 
       real (kind=dbl_kind), dimension (:), intent(in) :: &
-         apnd     ,& ! melt pond fraction of sea ice
-         hpnd     ,& ! mean melt pond depth over sea ice
-         ipnd     ,& ! mean ice pond depth over sea ice in cat n
+         apondn   ,& ! melt pond fraction of sea ice category
          alvl     ,& ! level ice area fraction (of grid cell ?)
          vlvl        ! level ice mean thickness
 
@@ -666,7 +663,7 @@ subroutine neutral_drag_coeffs (apnd,     hpnd,     &
          if (tr_pond) then
             do n = 1,ncat
                ! area of pond per unit area of grid cell
-               apond = apond+apnd(n)*aicen(n)
+               apond = apond+apondn(n)*aicen(n)
             enddo
          endif
 

diff --git a/columnphysics/icepack_brine.F90 b/columnphysics/icepack_brine.F90
@@ -9,8 +9,8 @@ module icepack_brine
       use icepack_kinds
       use icepack_parameters, only: p01, p001, p5, c0, c1, c2, c1p5, puny, p25
       use icepack_parameters, only: gravit, rhoi, rhow, rhos, depressT
-      use icepack_parameters, only: salt_loss, min_salin, rhosi
-      use icepack_parameters, only: dts_b, l_sk
+      use icepack_parameters, only: min_salin, rhosi
+      use icepack_parameters, only: l_sk
       use icepack_tracers, only: nilyr, nblyr, ntrcr, nt_qice, nt_sice
       use icepack_tracers, only: nt_Tsfc
       use icepack_zbgc_shared, only: k_o, exp_h, Dm, Ra_c, viscos_dynamic, thinS
@@ -28,8 +28,7 @@ module icepack_brine
                 compute_microS_mushy, &
                 update_hbrine, &
                 calculate_drho, &
-                icepack_init_hbrine, &
-                icepack_init_zsalinity   ! deprecated
+                icepack_init_hbrine
 
       real (kind=dbl_kind), parameter :: &
          maxhbr  = 1.25_dbl_kind  , & ! brine overflows if hbr > maxhbr*hin
@@ -51,7 +50,7 @@ module icepack_brine
 
 !=======================================================================
 ! Computes the top and bottom brine boundary changes for flushing
-! works for zsalinity and tr_salinity
+! works for tr_salinity
 !
 ! NOTE: In this subroutine, trcrn(nt_fbri) is the volume fraction of ice with
 ! dynamic salinity or the height ratio = hbr/vicen*aicen, where hbr is the
@@ -680,7 +679,7 @@ subroutine icepack_init_hbrine(bgrid_out, igrid_out, cgrid_out, &
       ! Calculate bio gridn: 0 to 1 corresponds to ice top to bottom
       !-----------------------------------------------------------------
 
-      bgrid(:)       = c0 ! zsalinity grid points
+      bgrid(:)       = c0 ! biology nondimensional vertical grid points
       bgrid(nblyr+2) = c1 ! bottom value
       igrid(:)       = c0 ! bgc interface grid points
       igrid(1)       = c0 ! ice top
@@ -738,36 +737,6 @@ subroutine icepack_init_hbrine(bgrid_out, igrid_out, cgrid_out, &
 
       end subroutine icepack_init_hbrine
 
-!=======================================================================
-!autodocument_start icepack_init_zsalinity
-!  **DEPRECATED**, all code removed
-!  Interface provided for backwards compatibility
-
-      subroutine icepack_init_zsalinity(Rayleigh_criteria, &
-               Rayleigh_real, trcrn_bgc, sss)
-
-      logical (kind=log_kind), intent(inout) :: &
-       Rayleigh_criteria
-
-      real (kind=dbl_kind), intent(inout):: &
-       Rayleigh_real
-
-      real (kind=dbl_kind), intent(in):: &
-       sss
-
-      real (kind=dbl_kind), dimension(:,:), intent(inout):: &
-       trcrn_bgc  ! bgc subset of trcrn
-
-!autodocument_end
-
-      ! local variables
-
-      character(len=*),parameter :: subname='(icepack_init_zsalinity)'
-
-      call icepack_warnings_add(subname//' DEPRECATED, do not use')
-!      call icepack_warnings_setabort(.true.,__FILE__,__LINE__)
-
-      end subroutine icepack_init_zsalinity
 !=======================================================================
 
       end module icepack_brine

diff --git a/columnphysics/icepack_flux.F90 b/columnphysics/icepack_flux.F90
@@ -12,7 +12,7 @@ module icepack_flux
       use icepack_parameters, only: c1, emissivity, snwgrain
       use icepack_warnings, only: warnstr, icepack_warnings_add
       use icepack_warnings, only: icepack_warnings_setabort, icepack_warnings_aborted
-      use icepack_tracers, only: tr_iso
+      use icepack_tracers, only: tr_iso, tr_pond
 
       implicit none
       private
@@ -43,6 +43,8 @@ subroutine merge_fluxes (aicen,                &
                                fhocnn,   fswthrun,   &
                                fswthrun_vdr, fswthrun_vdf,&
                                fswthrun_idr, fswthrun_idf,&
+                               fswthrun_uvrdr, fswthrun_uvrdf,&
+                               fswthrun_pardr, fswthrun_pardf,&
                                strairxT, strairyT,   &
                                Cdn_atm_ratio,        &
                                fsurf,    fcondtop,   &
@@ -56,6 +58,8 @@ subroutine merge_fluxes (aicen,                &
                                fhocn,    fswthru,    &
                                fswthru_vdr, fswthru_vdf,&
                                fswthru_idr, fswthru_idf,&
+                               fswthru_uvrdr, fswthru_uvrdf,&
+                               fswthru_pardr, fswthru_pardf,&
                                melttn, meltsn, meltbn, congeln, snoicen, &
                                meltt,  melts,        &
                                meltb,  dsnow, dsnown,&
@@ -64,12 +68,18 @@ subroutine merge_fluxes (aicen,                &
                                Uref,     Urefn,      &
                                Qref_iso, Qrefn_iso,  &
                                fiso_ocn, fiso_ocnn,  &
-                               fiso_evap, fiso_evapn)
-
-      ! single category fluxes
+                               fiso_evap, fiso_evapn,&
+                               dpnd_flush,   dpnd_flushn,   &
+                               dpnd_expon,   dpnd_exponn,   &
+                               dpnd_freebd,  dpnd_freebdn,  &
+                               dpnd_initial, dpnd_initialn, &
+                               dpnd_dlid,    dpnd_dlidn)
+
+      ! concentration is aicen_init in call to subroutine
       real (kind=dbl_kind), intent(in) :: &
           aicen       ! concentration of ice
 
+      ! single category fluxes
       real (kind=dbl_kind), optional, intent(in) :: &
           flw     , & ! downward longwave flux          (W/m**2)
           strairxn, & ! air/ice zonal  strss,           (N/m**2)
@@ -98,10 +108,19 @@ subroutine merge_fluxes (aicen,                &
           dsnown  , & ! change in snow depth            (m)
           congeln , & ! congelation ice growth          (m)
           snoicen , & ! snow-ice growth                 (m)
+          dpnd_flushn , & ! pond flushing rate due to ice permeability (m/step)
+          dpnd_exponn , & ! exponential pond drainage rate (m/step)
+          dpnd_freebdn, & ! pond drainage rate due to freeboard constraint (m/step)
+          dpnd_initialn,& ! runoff rate due to rfrac (m/step)
+          dpnd_dlidn  , & ! pond loss/gain due to ice lid (m/step)
           fswthrun_vdr, & ! vis dir sw radiation through ice bot    (W/m**2)
           fswthrun_vdf, & ! vis dif sw radiation through ice bot    (W/m**2)
           fswthrun_idr, & ! nir dir sw radiation through ice bot    (W/m**2)
           fswthrun_idf, & ! nir dif sw radiation through ice bot    (W/m**2)
+          fswthrun_uvrdr, & !   > 700nm vis uvr dir sw radiation through ice bot (W/m**2)
+          fswthrun_uvrdf, & !   > 700nm vis uvr dif sw radiation through ice bot (W/m**2)
+          fswthrun_pardr, & ! 400-700nm vis par dir sw radiation through ice bot (W/m**2)
+          fswthrun_pardf, & ! 400-700nm vis par dif sw radiation through ice bot (W/m**2)
           Urefn       ! air speed reference level       (m/s)
 
       ! cumulative fluxes
@@ -131,10 +150,19 @@ subroutine merge_fluxes (aicen,                &
           meltsliq, & ! mass of snow melt               (kg/m^2)
           congel  , & ! congelation ice growth          (m)
           snoice  , & ! snow-ice growth                 (m)
+          dpnd_flush , & ! pond flushing rate due to ice permeability (m/step)
+          dpnd_expon , & ! exponential pond drainage rate (m/step)
+          dpnd_freebd, & ! pond drainage rate due to freeboard constraint (m/step)
+          dpnd_initial,& ! runoff rate due to rfrac (m/step)
+          dpnd_dlid  , & ! pond loss/gain (+/-) to ice lid freezing/melting (m/step)
           fswthru_vdr, & ! vis dir sw radiation through ice bot    (W/m**2)
           fswthru_vdf, & ! vis dif sw radiation through ice bot    (W/m**2)
           fswthru_idr, & ! nir dir sw radiation through ice bot    (W/m**2)
           fswthru_idf, & ! nir dif sw radiation through ice bot    (W/m**2)
+          fswthru_uvrdr, & !   > 700nm vis uvr dir sw radiation through ice bot (W/m**2)
+          fswthru_uvrdf, & !   > 700nm vis uvr dif sw radiation through ice bot (W/m**2)
+          fswthru_pardr, & ! 400-700nm vis par dir sw radiation through ice bot (W/m**2)
+          fswthru_pardf, & ! 400-700nm vis par dif sw radiation through ice bot (W/m**2)
           dsnow,    & ! change in snow depth            (m)
           Uref        ! air speed reference level       (m/s)
 
@@ -228,6 +256,15 @@ subroutine merge_fluxes (aicen,                &
       if (present(fswthrun_idf) .and. present(fswthru_idf)) &
          fswthru_idf = fswthru_idf + fswthrun_idf  * aicen
 
+      if (present(fswthrun_uvrdr) .and. present(fswthru_uvrdr)) &
+         fswthru_uvrdr   = fswthru_uvrdr   + fswthrun_uvrdr  * aicen
+      if (present(fswthrun_uvrdf) .and. present(fswthru_uvrdf)) &
+         fswthru_uvrdf   = fswthru_uvrdf   + fswthrun_uvrdf  * aicen
+      if (present(fswthrun_pardr) .and. present(fswthru_pardr)) &
+         fswthru_pardr   = fswthru_pardr   + fswthrun_pardr  * aicen
+      if (present(fswthrun_pardf) .and. present(fswthru_pardf)) &
+         fswthru_pardf   = fswthru_pardf   + fswthrun_pardf  * aicen
+
       ! ice/snow thickness
 
       if (present(melttn) .and. present(meltt)) &
@@ -247,6 +284,19 @@ subroutine merge_fluxes (aicen,                &
          congel    = congel    + congeln   * aicen
       if (present(snoicen) .and. present(snoice)) &
          snoice    = snoice    + snoicen   * aicen
+      ! Meltwater fluxes
+      if (tr_pond) then
+         if (present(dpnd_flushn)  .and. present(dpnd_flush))   &
+            dpnd_flush   = dpnd_flush   + dpnd_flushn   * aicen
+         if (present(dpnd_exponn)  .and. present(dpnd_expon))   &
+            dpnd_expon   = dpnd_expon   + dpnd_exponn   * aicen
+         if (present(dpnd_freebdn) .and. present(dpnd_freebd))  &
+            dpnd_freebd  = dpnd_freebd  + dpnd_freebdn  * aicen
+         if (present(dpnd_initialn).and. present(dpnd_initial)) &
+            dpnd_initial = dpnd_initial + dpnd_initialn * aicen
+         if (present(dpnd_dlidn)   .and. present(dpnd_dlid))    &
+            dpnd_dlid    = dpnd_dlid    + dpnd_dlidn    * aicen
+      endif
 
       end subroutine merge_fluxes
 

diff --git a/columnphysics/icepack_fsd.F90 b/columnphysics/icepack_fsd.F90
@@ -180,7 +180,7 @@ subroutine icepack_init_fsd_bounds( &
          floe_area_c         (nfsd), & ! fsd area at bin centre (m^2)
          floe_area_binwidth  (nfsd), & ! floe area bin width (m^2)
          floe_binwidth       (nfsd), & ! floe bin width (m)
-         c_fsd_range         (nfsd), & ! 
+         c_fsd_range         (nfsd), & !
          iweld         (nfsd, nfsd), & ! fsd categories that can weld
          stat=ierr)
       if (ierr/=0) then
@@ -620,7 +620,6 @@ subroutine fsd_add_new_ice (n, &
                                   wave_sig_ht,               &
                                   wave_spectrum,             &
                                   wavefreq,                  &
-                                  dwavefreq,                 &
                                   d_afsd_latg,               &
                                   d_afsd_newi,               &
                                   afsdn,      aicen_init,    &
@@ -640,8 +639,7 @@ subroutine fsd_add_new_ice (n, &
                         ! power spectral density of surface elevation, E(f) (units m^2 s)
 
       real(kind=dbl_kind), dimension(:), intent(in) :: &
-         wavefreq   , & ! wave frequencies (s^-1)
-         dwavefreq      ! wave frequency bin widths (s^-1)
+         wavefreq       ! wave frequencies (s^-1)
 
       real (kind=dbl_kind), dimension(:), intent(in) :: &
          d_an_latg  , & ! change in aicen due to lateral growth
@@ -747,8 +745,7 @@ subroutine fsd_add_new_ice (n, &
                if (wave_spec) then
                   if (wave_sig_ht > puny) then
                      call wave_dep_growth (wave_spectrum, wave_sig_ht, &
-                                           wavefreq, dwavefreq, &
-                                           new_size)
+                                           wavefreq, new_size)
                      if (icepack_warnings_aborted(subname)) return
                   end if
 
@@ -775,8 +772,7 @@ subroutine fsd_add_new_ice (n, &
                if (wave_spec) then
                   if (wave_sig_ht > puny) then
                      call wave_dep_growth (wave_spectrum, wave_sig_ht, &
-                                           wavefreq, dwavefreq, &
-                                           new_size)
+                                           wavefreq, new_size)
                      if (icepack_warnings_aborted(subname)) return
                   end if
 
@@ -817,20 +813,18 @@ end subroutine fsd_add_new_ice
 !  authors: Lettie Roach, NIWA/VUW
 !
       subroutine wave_dep_growth (local_wave_spec, wave_height, &
-                                  wavefreq, dwavefreq, &
-                                  new_size)
+                                  wavefreq, new_size)
 
       real (kind=dbl_kind), dimension(:), intent(in) :: &
          local_wave_spec ! ocean surface wave spectrum as a function of frequency
                          ! power spectral density of surface elevation, E(f) (units m^2 s)
                          ! dimension set in ice_forcing
 
       real(kind=dbl_kind), dimension(:), intent(in) :: &
-         wavefreq,     & ! wave frequencies (s^-1)
-         dwavefreq       ! wave frequency bin widths (s^-1)
+         wavefreq        ! wave frequencies (s^-1)
 
-      real(kind=dbl_kind), intent(in), optional :: &
-         wave_height ! significant wave height (m)
+      real(kind=dbl_kind), intent(in) :: &
+         wave_height     ! significant wave height (m)
 
       integer (kind=int_kind), intent(out) :: &
          new_size        ! index of floe size category in which new floes will grow
@@ -847,11 +841,7 @@ subroutine wave_dep_growth (local_wave_spec, wave_height, &
 
       integer (kind=int_kind) :: k
 
-      if (present(wave_height)) then
-         w_amp = p5 * wave_height  ! amplitude is 1/2 sig wave hight
-      else 
-         w_amp = c2* SQRT(SUM(local_wave_spec*dwavefreq))   ! sig wave amplitude
-      endif
+      w_amp = p5 * wave_height  ! amplitude is 1/2 sig wave height
       f_peak = wavefreq(MAXLOC(local_wave_spec, DIM=1))  ! peak frequency
 
       ! tensile failure