Debug build

.zenodo.json

@@ -1,22 +1,22 @@
 {
     "license": "BSD-3-Clause",
-    "copyright": "Copyright (c) 1998, 2017 Triad National Security, LLC",
+    "copyright": "Copyright 1998-2026 Triad National Security, LLC",
     "description": "View detailed release notes at https://github.com/CICE-Consortium/CICE/releases", 
     "language": "eng", 
-    "title": "CICE-Consortium/CICE: CICE Version 6.6.1", 
+    "title": "CICE-Consortium/CICE: CICE Version 6.6.3", 
     "keywords": [
         "sea ice model",
         "CICE",
         "Icepack"
     ], 
-    "version": "6.6.1", 
+    "version": "6.6.3", 
     "upload_type": "software", 
     "communities": [
         {
             "identifier": "cice-consortium"
         }
     ], 
-    "publication_date": "2025-07-21",
+    "publication_date": "2026-01-22", 
     "creators": [
         {
             "affiliation": "Los Alamos National Laboratory", 

cicecore/cicedyn/analysis/ice_history_bgc.F90

@@ -3398,7 +3398,7 @@ subroutine init_history_bgc
 
       use ice_arrays_column, only: PP_net, grow_net, hbri, &
           ice_bio_net, snow_bio_net, fbio_snoice, fbio_atmice, &
-          zfswin
+          upNO, upNH, bTiz, bphi, zfswin, iDi, iki
       use ice_flux_bgc, only: flux_bio, flux_bio_ai, fnit, fsil, &
           famm, fdmsp, fdms, fhum, fdust, falgalN, fdoc, fdic, &
           fdon, ffep, ffed
@@ -3407,14 +3407,20 @@ subroutine init_history_bgc
 
       PP_net        (:,:,:) = c0
       grow_net      (:,:,:) = c0
+      upNO          (:,:,:) = c0
+      upNH          (:,:,:) = c0
       hbri          (:,:,:) = c0
       flux_bio    (:,:,:,:) = c0
       flux_bio_ai (:,:,:,:) = c0
       ice_bio_net (:,:,:,:) = c0
       snow_bio_net(:,:,:,:) = c0
       fbio_snoice (:,:,:,:) = c0
       fbio_atmice (:,:,:,:) = c0
+      bTiz      (:,:,:,:,:) = c0
+      bphi      (:,:,:,:,:) = c0
       zfswin    (:,:,:,:,:) = c0
+      iDi       (:,:,:,:,:) = c0
+      iki       (:,:,:,:,:) = c0
       fnit          (:,:,:) = c0
       fsil          (:,:,:) = c0
       famm          (:,:,:) = c0

cicecore/cicedyn/analysis/ice_history_mechred.F90

@@ -37,6 +37,11 @@ module ice_history_mechred
            f_aredistn  = 'x', f_vredistn   = 'x', &
            f_araftn    = 'x', f_vraftn     = 'x'
 
+      ! CMIP ridging variables.
+
+      character (len=max_nstrm), public :: &
+           f_sirdgconc = 'm', f_sirdgthick = 'm'
+
       !---------------------------------------------------------------
       ! namelist variables
       !---------------------------------------------------------------
@@ -51,13 +56,14 @@ module ice_history_mechred
            f_dvirdgndt, &
            f_aparticn,  f_krdgn    , &
            f_aredistn,  f_vredistn , &
-           f_araftn,    f_vraftn
+           f_araftn,    f_vraftn   , &
+           f_sirdgconc, f_sirdgthick
 
       !---------------------------------------------------------------
       ! field indices
       !---------------------------------------------------------------
 
-      integer (kind=int_kind), dimension(max_nstrm) :: &
+      integer (kind=int_kind), dimension(max_nstrm), public :: &
            n_ardg       , n_vrdg       , &
            n_alvl       , n_vlvl       , &
            n_dardg1dt   , n_dardg2dt   , &
@@ -69,6 +75,9 @@ module ice_history_mechred
            n_aredistn   , n_vredistn   , &
            n_araftn     , n_vraftn
 
+      integer (kind=int_kind), dimension(max_nstrm) :: &
+           n_sirdgconc, n_sirdgthick
+
 !=======================================================================
 
       contains
@@ -152,6 +161,8 @@ subroutine init_hist_mechred_2D
          f_vrdgn = 'x'
          f_araftn = 'x'
          f_vraftn = 'x'
+         f_sirdgconc  = 'x'
+         f_sirdgthick = 'x'
       endif
       if (f_araftn /= 'x' .or. f_vraftn /= 'x') f_ardgn = f_araftn
 
@@ -174,6 +185,8 @@ subroutine init_hist_mechred_2D
       call broadcast_scalar (f_vredistn, master_task)
       call broadcast_scalar (f_araftn, master_task)
       call broadcast_scalar (f_vraftn, master_task)
+      call broadcast_scalar (f_sirdgconc, master_task)
+      call broadcast_scalar (f_sirdgthick, master_task)
 
       ! 2D variables
 
@@ -182,49 +195,63 @@ subroutine init_hist_mechred_2D
 
       if (f_alvl(1:1) /= 'x') &
          call define_hist_field(n_alvl,"alvl","1",tstr2D, tcstr, &
-             "level ice area fraction",                            &
-             "none", c1, c0,                                       &
+             "level ice area fraction",                          &
+             "none", c1, c0,                                     &
              ns, f_alvl)
       if (f_vlvl(1:1) /= 'x') &
          call define_hist_field(n_vlvl,"vlvl","m",tstr2D, tcstr, &
-             "level ice volume",                           &
-             "grid cell mean level ice thickness", c1, c0, &
+             "level ice volume",                                 &
+             "grid cell mean level ice thickness", c1, c0,       &
              ns, f_vlvl)
       if (f_ardg(1:1) /= 'x') &
          call define_hist_field(n_ardg,"ardg","1",tstr2D, tcstr, &
-             "ridged ice area fraction",                           &
-             "none", c1, c0,                                       &
+             "ridged ice area fraction",                         &
+             "none", c1, c0,                                     &
              ns, f_ardg)
       if (f_vrdg(1:1) /= 'x') &
          call define_hist_field(n_vrdg,"vrdg","m",tstr2D, tcstr, &
-             "ridged ice volume",                          &
-             "grid cell mean level ridged thickness", c1, c0, &
+             "ridged ice volume",                                &
+             "grid cell mean level ridged thickness", c1, c0,    &
              ns, f_vrdg)
 
       if (f_dardg1dt(1:1) /= 'x') &
          call define_hist_field(n_dardg1dt,"dardg1dt","%/day",tstr2D, tcstr, &
-             "ice area ridging rate",                                      &
-             "none", secday*c100, c0,                                      &
+             "ice area ridging rate",                                        &
+             "none", secday*c100, c0,                                        &
              ns, f_dardg1dt)
 
       if (f_dardg2dt(1:1) /= 'x') &
          call define_hist_field(n_dardg2dt,"dardg2dt","%/day",tstr2D, tcstr, &
-             "ridge area formation rate",                                  &
-             "none", secday*c100, c0,                                      &
+             "ridge area formation rate",                                    &
+             "none", secday*c100, c0,                                        &
              ns, f_dardg2dt)
 
       if (f_dvirdgdt(1:1) /= 'x') &
          call define_hist_field(n_dvirdgdt,"dvirdgdt","cm/day",tstr2D, tcstr, &
-             "ice volume ridging rate",                                     &
-             "none", mps_to_cmpdy, c0,                                      &
+             "ice volume ridging rate",                                       &
+             "none", mps_to_cmpdy, c0,                                        &
              ns, f_dvirdgdt)
 
       if (f_opening(1:1) /= 'x') &
          call define_hist_field(n_opening,"opening","%/day",tstr2D, tcstr, &
-             "lead area opening rate",                                   &
-             "none", secday*c100, c0,                                    &
+             "lead area opening rate",                                     &
+             "none", secday*c100, c0,                                      &
              ns, f_opening)
 
+      if (f_sirdgconc(1:1) /= 'x') &
+         call define_hist_field(n_sirdgconc,"sirdgconc","%",tstr2D, tcstr, &
+             "percentage of ridged sea ice",                               &
+             "area percentage of sea ice surface that is ridged sea ice",  &
+             c100, c0,                                                     &
+             ns, f_sirdgconc, avg_ice_present='none', mask_ice_free_points=.false.)
+
+      if (f_sirdgthick(1:1) /= 'x') &
+         call define_hist_field(n_sirdgthick,"sirdgthick","m",tstr2D, tcstr, &
+             "ridged ice thickness",                                         &
+             "total volume of ridged sea ice divided by area of ridges",     &
+             c1, c0,                                                         &
+             ns, f_sirdgthick, avg_ice_present='ridge', mask_ice_free_points=.true.)
+
       endif ! histfreq(ns) /= 'x'
       enddo ! nstreams
 
@@ -331,6 +358,7 @@ end subroutine init_hist_mechred_3Dc
 
       subroutine accum_hist_mechred (iblk)
 
+      use ice_blocks, only: nx_block, ny_block
       use ice_history_shared, only: n2D, a2D, a3Dc, ncat_hist, &
           accum_hist_field
       use ice_state, only: aice, vice, trcr, aicen, vicen, trcrn
@@ -380,6 +408,15 @@ subroutine accum_hist_mechred (iblk)
          if (f_opening(1:1) /= 'x') &
              call accum_hist_field(n_opening, iblk, opening(:,:,iblk), a2D)
 
+         if (f_sirdgconc(1:1)/= 'x') &
+             call accum_hist_field(n_sirdgconc,   iblk, &
+                             aice(:,:,iblk) * (c1 - trcr(:,:,nt_alvl,iblk)), a2D)
+
+         if (f_sirdgthick(1:1)/= 'x') then
+             call accum_hist_field(n_sirdgthick,  iblk, &
+                             vice(:,:,iblk) * (c1 - trcr(:,:,nt_vlvl,iblk)), a2D)
+         endif
+
          endif ! allocated(a2D)
 
          ! 3D category fields