diff --git a/cicecore/cicedynB/analysis/ice_history_pond.F90 b/cicecore/cicedynB/analysis/ice_history_pond.F90
index f6e4b8737..8818ff94e 100644
--- a/cicecore/cicedynB/analysis/ice_history_pond.F90
+++ b/cicecore/cicedynB/analysis/ice_history_pond.F90
@@ -268,13 +268,8 @@ subroutine accum_hist_pond (iblk)
 
       integer (kind=int_kind) :: &
          nt_apnd, nt_hpnd, nt_alvl, nt_ipnd
-#ifdef UNDEPRECATE_CESMPONDS
-      logical (kind=log_kind) :: &
-         tr_pond_cesm, tr_pond_lvl, tr_pond_topo
-#else
       logical (kind=log_kind) :: &
          tr_pond_lvl, tr_pond_topo
-#endif
 
       real (kind=dbl_kind) :: &
          puny
@@ -289,13 +284,8 @@ subroutine accum_hist_pond (iblk)
       !---------------------------------------------------------------
 
          call icepack_query_parameters(puny_out=puny)
-#ifdef UNDEPRECATE_CESMPONDS
-         call icepack_query_tracer_flags(tr_pond_cesm_out=tr_pond_cesm, &
-              tr_pond_lvl_out=tr_pond_lvl, tr_pond_topo_out=tr_pond_topo)
-#else
          call icepack_query_tracer_flags(tr_pond_lvl_out=tr_pond_lvl, &
               tr_pond_topo_out=tr_pond_topo)
-#endif
          call icepack_query_tracer_indices(nt_apnd_out=nt_apnd, nt_hpnd_out=nt_hpnd, &
               nt_alvl_out=nt_alvl, nt_ipnd_out=nt_ipnd)
          call icepack_warnings_flush(nu_diag)
@@ -303,27 +293,9 @@ subroutine accum_hist_pond (iblk)
             file=__FILE__, line=__LINE__)
 
          if (allocated(a2D)) then
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) then
-         if (f_apond(1:1)/= 'x') &
-             call accum_hist_field(n_apond, iblk, &
-                                   trcr(:,:,nt_apnd,iblk), a2D)
-         if (f_apond_ai(1:1)/= 'x') &
-             call accum_hist_field(n_apond_ai, iblk, &
-                                   aice(:,:,iblk) * trcr(:,:,nt_apnd,iblk), a2D)
-         if (f_hpond(1:1)/= 'x') &
-             call accum_hist_field(n_hpond, iblk, &
-                                   trcr(:,:,nt_apnd,iblk) &
-                                 * trcr(:,:,nt_hpnd,iblk), a2D)
-         if (f_hpond_ai(1:1)/= 'x') &
-             call accum_hist_field(n_hpond_ai, iblk, &
-                                   aice(:,:,iblk) * trcr(:,:,nt_apnd,iblk) &
-                                                  * trcr(:,:,nt_hpnd,iblk), a2D)
 
-         elseif (tr_pond_lvl) then
-#else
          if (tr_pond_lvl) then
-#endif
+
          if (f_apond(1:1)/= 'x') &
              call accum_hist_field(n_apond, iblk, &
                             trcr(:,:,nt_alvl,iblk) * trcr(:,:,nt_apnd,iblk), a2D)
diff --git a/cicecore/cicedynB/dynamics/ice_transport_driver.F90 b/cicecore/cicedynB/dynamics/ice_transport_driver.F90
index 43fe5af13..30fe546e0 100644
--- a/cicecore/cicedynB/dynamics/ice_transport_driver.F90
+++ b/cicecore/cicedynB/dynamics/ice_transport_driver.F90
@@ -1534,13 +1534,8 @@ subroutine state_to_work (nx_block, ny_block,        &
       integer (kind=int_kind) :: &
          nt_alvl, nt_apnd, nt_fbri
 
-#ifdef UNDEPRECATE_CESMPONDS
-      logical (kind=log_kind) :: &
-         tr_pond_cesm, tr_pond_lvl, tr_pond_topo
-#else
       logical (kind=log_kind) :: &
          tr_pond_lvl, tr_pond_topo
-#endif
 
       integer (kind=int_kind) ::      &
          i, j, n, it, & ! counting indices
@@ -1548,13 +1543,8 @@ subroutine state_to_work (nx_block, ny_block,        &
 
       character(len=*), parameter :: subname = '(state_to_work)'
 
-#ifdef UNDEPRECATE_CESMPONDS
-      call icepack_query_tracer_flags(tr_pond_cesm_out=tr_pond_cesm, &
-           tr_pond_lvl_out=tr_pond_lvl, tr_pond_topo_out=tr_pond_topo)
-#else
       call icepack_query_tracer_flags(tr_pond_lvl_out=tr_pond_lvl, &
            tr_pond_topo_out=tr_pond_topo)
-#endif
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_apnd_out=nt_apnd, &
            nt_fbri_out=nt_fbri)
       call icepack_warnings_flush(nu_diag)
@@ -1613,13 +1603,8 @@ subroutine state_to_work (nx_block, ny_block,        &
                                         * trcrn(i,j,it     ,n)
                enddo
                enddo
-#ifdef UNDEPRECATE_CESMPONDS
-            elseif (trcr_depend(it) == 2+nt_apnd .and. &
-                    tr_pond_cesm .or. tr_pond_topo) then
-#else
             elseif (trcr_depend(it) == 2+nt_apnd .and. &
                     tr_pond_topo) then
-#endif
                do j = 1, ny_block
                do i = 1, nx_block
                   works(i,j,narrays+it) = aicen(i,j        ,n) &
diff --git a/cicecore/cicedynB/general/ice_init.F90 b/cicecore/cicedynB/general/ice_init.F90
index a6c4d7941..45ae58d8b 100644
--- a/cicecore/cicedynB/general/ice_init.F90
+++ b/cicecore/cicedynB/general/ice_init.F90
@@ -73,15 +73,9 @@ subroutine input_data
                               npt, dt, ndtd, days_per_year, use_leap_years, &
                               write_ic, dump_last, npt_unit
       use ice_arrays_column, only: oceanmixed_ice
-#ifdef UNDEPRECATE_CESMPONDS
-      use ice_restart_column, only: restart_age, restart_FY, restart_lvl, &
-          restart_pond_cesm, restart_pond_lvl, restart_pond_topo, restart_aero, &
-          restart_fsd, restart_iso, restart_snow
-#else
       use ice_restart_column, only: restart_age, restart_FY, restart_lvl, &
           restart_pond_lvl, restart_pond_topo, restart_aero, &
           restart_fsd, restart_iso, restart_snow
-#endif
       use ice_restart_shared, only: &
           restart, restart_ext, restart_coszen, restart_dir, restart_file, pointer_file, &
           runid, runtype, use_restart_time, restart_format, lcdf64
@@ -165,23 +159,12 @@ subroutine input_data
 
       logical (kind=log_kind) :: tr_iage, tr_FY, tr_lvl, tr_pond
       logical (kind=log_kind) :: tr_iso, tr_aero, tr_fsd, tr_snow
-#ifdef UNDEPRECATE_CESMPONDS
-      logical (kind=log_kind) :: tr_pond_cesm, tr_pond_lvl, tr_pond_topo
-#else
       logical (kind=log_kind) :: tr_pond_lvl, tr_pond_topo
-#endif
       integer (kind=int_kind) :: numin, numax  ! unit number limits
 
-#ifdef UNDEPRECATE_CESMPONDS
-      integer (kind=int_kind) :: rpcesm, rplvl, rptopo
-#else
       integer (kind=int_kind) :: rplvl, rptopo
-#endif
       real (kind=dbl_kind) :: Cf, ksno, puny
       character (len=char_len) :: abort_list
-#ifdef CESMCOUPLED
-      character (len=64) :: tmpstr
-#endif
       character (len=128) :: tmpstr2
 
       character(len=*), parameter :: subname='(input_data)'
@@ -221,9 +204,6 @@ subroutine input_data
         tr_iage, restart_age,                                           &
         tr_FY, restart_FY,                                              &
         tr_lvl, restart_lvl,                                            &
-#ifdef UNDEPRECATE_CESMPONDS
-        tr_pond_cesm, restart_pond_cesm,                                &
-#endif
         tr_pond_lvl, restart_pond_lvl,                                  &
         tr_pond_topo, restart_pond_topo,                                &
         tr_snow, restart_snow,                                          &
@@ -445,11 +425,7 @@ subroutine input_data
       conserv_check = .false.  ! tracer conservation check
       shortwave = 'ccsm3'      ! 'ccsm3' or 'dEdd' (delta-Eddington)
       albedo_type = 'ccsm3'    ! 'ccsm3' or 'constant'
-#ifdef UNDEPRECATE_0LAYER
-      ktherm = 1               ! -1 = OFF, 0 = 0-layer, 1 = BL99, 2 = mushy thermo
-#else
       ktherm = 1               ! -1 = OFF, 1 = BL99, 2 = mushy thermo
-#endif
       conduct = 'bubbly'       ! 'MU71' or 'bubbly' (Pringle et al 2007)
       coriolis = 'latitude'    ! latitude dependent, or 'constant'
       ssh_stress = 'geostrophic'  ! 'geostrophic' or 'coupled'
@@ -559,10 +535,6 @@ subroutine input_data
       restart_FY   = .false. ! ice age restart
       tr_lvl       = .false. ! level ice
       restart_lvl  = .false. ! level ice restart
-#ifdef UNDEPRECATE_CESMPONDS
-      tr_pond_cesm = .false. ! CESM melt ponds
-      restart_pond_cesm = .false. ! melt ponds restart
-#endif
       tr_pond_lvl  = .false. ! level-ice melt ponds
       restart_pond_lvl  = .false. ! melt ponds restart
       tr_pond_topo = .false. ! explicit melt ponds (topographic)
@@ -1034,10 +1006,6 @@ subroutine input_data
       call broadcast_scalar(restart_FY,           master_task)
       call broadcast_scalar(tr_lvl,               master_task)
       call broadcast_scalar(restart_lvl,          master_task)
-#ifdef UNDEPRECATE_CESMPONDS
-      call broadcast_scalar(tr_pond_cesm,         master_task)
-      call broadcast_scalar(restart_pond_cesm,    master_task)
-#endif
       call broadcast_scalar(tr_pond_lvl,          master_task)
       call broadcast_scalar(restart_pond_lvl,     master_task)
       call broadcast_scalar(tr_pond_topo,         master_task)
@@ -1130,9 +1098,6 @@ subroutine input_data
             restart_age =  .false.
             restart_fy =  .false.
             restart_lvl =  .false.
-#ifdef UNDEPRECATE_CESMPONDS
-            restart_pond_cesm =  .false.
-#endif
             restart_pond_lvl =  .false.
             restart_pond_topo =  .false.
             restart_snow = .false.
@@ -1249,29 +1214,15 @@ subroutine input_data
          endif
       endif
 
-#ifdef UNDEPRECATE_CESMPONDS
-      rpcesm = 0
-#endif
       rplvl  = 0
       rptopo = 0
-#ifdef UNDEPRECATE_CESMPONDS
-      if (tr_pond_cesm) rpcesm = 1
-#endif
       if (tr_pond_lvl ) rplvl  = 1
       if (tr_pond_topo) rptopo = 1
 
       tr_pond = .false. ! explicit melt ponds
-#ifdef UNDEPRECATE_CESMPONDS
-      if (rpcesm + rplvl + rptopo > 0) tr_pond = .true.
-#else
       if (rplvl + rptopo > 0) tr_pond = .true.
-#endif
 
-#ifdef UNDEPRECATE_CESMPONDS
-      if (rpcesm + rplvl + rptopo > 1) then
-#else
       if (rplvl + rptopo > 1) then
-#endif
          if (my_task == master_task) then
             write(nu_diag,*) subname//' ERROR: Must use only one melt pond scheme'
          endif
@@ -1495,13 +1446,6 @@ subroutine input_data
             abort_list = trim(abort_list)//":16"
          endif
 
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) then
-            if (my_task == master_task) write(nu_diag,*) subname//' ERROR: formdrag=T and frzpnd=cesm'
-            abort_list = trim(abort_list)//":17"
-         endif
-#endif
-
          if (.not. tr_lvl) then
             if (my_task == master_task) write(nu_diag,*) subname//' ERROR: formdrag=T and tr_lvl=F'
             abort_list = trim(abort_list)//":18"
@@ -1925,10 +1869,6 @@ subroutine input_data
             tmpstr2 = ' : Bitz and Lipscomb 1999 thermo'
          elseif (ktherm == 2) then
             tmpstr2 = ' : mushy-layer thermo'
-#ifdef UNDEPRECATE_0LAYER
-         elseif (ktherm == 0) then
-            tmpstr2 = ' : zero-layer thermo'
-#endif
          elseif (ktherm < 0) then
             tmpstr2 = ' : Thermodynamics disabled'
          else
@@ -2109,14 +2049,7 @@ subroutine input_data
          write(nu_diag,*) ' '
          write(nu_diag,*) ' Melt ponds'
          write(nu_diag,*) '--------------------------------'
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) then
-            write(nu_diag,1010) ' tr_pond_cesm     = ', tr_pond_cesm,' : CESM pond formulation'
-            write(nu_diag,1002) ' pndaspect        = ', pndaspect,' : ratio of pond depth to area fraction'
-         elseif (tr_pond_lvl) then
-#else
          if (tr_pond_lvl) then
-#endif
             write(nu_diag,1010) ' tr_pond_lvl      = ', tr_pond_lvl,' : level-ice pond formulation'
             write(nu_diag,1002) ' pndaspect        = ', pndaspect,' : ratio of pond depth to area fraction'
             write(nu_diag,1000) ' dpscale          = ', dpscale,' : time scale for flushing in permeable ice'
@@ -2229,9 +2162,6 @@ subroutine input_data
          if (tr_lvl)       write(nu_diag,1010) ' tr_lvl           = ', tr_lvl,' : ridging related tracers'
          if (tr_pond_lvl)  write(nu_diag,1010) ' tr_pond_lvl      = ', tr_pond_lvl,' : level-ice pond formulation'
          if (tr_pond_topo) write(nu_diag,1010) ' tr_pond_topo     = ', tr_pond_topo,' : topo pond formulation'
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) write(nu_diag,1010) ' tr_pond_cesm     = ', tr_pond_cesm,' : CESM pond formulation'
-#endif
          if (tr_snow)      write(nu_diag,1010) ' tr_snow          = ', tr_snow,' : advanced snow physics'
          if (tr_iage)      write(nu_diag,1010) ' tr_iage          = ', tr_iage,' : chronological ice age'
          if (tr_FY)        write(nu_diag,1010) ' tr_FY            = ', tr_FY,' : first-year ice area'
@@ -2356,9 +2286,6 @@ subroutine input_data
          write(nu_diag,1011) ' restart_age      = ', restart_age
          write(nu_diag,1011) ' restart_FY       = ', restart_FY
          write(nu_diag,1011) ' restart_lvl      = ', restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-         write(nu_diag,1011) ' restart_pond_cesm= ', restart_pond_cesm
-#endif
          write(nu_diag,1011) ' restart_pond_lvl = ', restart_pond_lvl
          write(nu_diag,1011) ' restart_pond_topo= ', restart_pond_topo
          write(nu_diag,1011) ' restart_snow     = ', restart_snow
@@ -2457,11 +2384,7 @@ subroutine input_data
       call icepack_init_tracer_flags(tr_iage_in=tr_iage, tr_FY_in=tr_FY, &
          tr_lvl_in=tr_lvl, tr_iso_in=tr_iso, tr_aero_in=tr_aero, &
          tr_fsd_in=tr_fsd, tr_snow_in=tr_snow, tr_pond_in=tr_pond, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_pond_cesm_in=tr_pond_cesm, tr_pond_lvl_in=tr_pond_lvl, tr_pond_topo_in=tr_pond_topo)
-#else
          tr_pond_lvl_in=tr_pond_lvl, tr_pond_topo_in=tr_pond_topo)
-#endif
       call icepack_init_tracer_sizes(ncat_in=ncat, nilyr_in=nilyr, nslyr_in=nslyr, nblyr_in=nblyr, &
          nfsd_in=nfsd, n_algae_in=n_algae, n_iso_in=n_iso, n_aero_in=n_aero, &
          n_DOC_in=n_DOC, n_DON_in=n_DON, &
@@ -2518,18 +2441,10 @@ subroutine init_state
          it          , & ! tracer index
          iblk            ! block index
 
-#ifdef UNDEPRECATE_0LAYER
-      logical (kind=log_kind) :: &
-         heat_capacity   ! from icepack
-#endif
 
       integer (kind=int_kind) :: ntrcr
       logical (kind=log_kind) :: tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero
-#ifdef UNDEPRECATE_CESMPONDS
-      logical (kind=log_kind) :: tr_pond_cesm, tr_pond_lvl, tr_pond_topo
-#else
       logical (kind=log_kind) :: tr_pond_lvl, tr_pond_topo
-#endif
       logical (kind=log_kind) :: tr_snow, tr_fsd
       integer (kind=int_kind) :: nt_Tsfc, nt_sice, nt_qice, nt_qsno, nt_iage, nt_FY
       integer (kind=int_kind) :: nt_alvl, nt_vlvl, nt_apnd, nt_hpnd, nt_ipnd
@@ -2543,17 +2458,10 @@ subroutine init_state
 
       !-----------------------------------------------------------------
 
-#ifdef UNDEPRECATE_0LAYER
-      call icepack_query_parameters(heat_capacity_out=heat_capacity)
-#endif
       call icepack_query_tracer_sizes(ntrcr_out=ntrcr)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
         tr_lvl_out=tr_lvl, tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#ifdef UNDEPRECATE_CESMPONDS
-        tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, tr_pond_topo_out=tr_pond_topo, &
-#else
         tr_pond_lvl_out=tr_pond_lvl, tr_pond_topo_out=tr_pond_topo, &
-#endif
         tr_snow_out=tr_snow, tr_fsd_out=tr_fsd)
       call icepack_query_tracer_indices(nt_Tsfc_out=nt_Tsfc, nt_sice_out=nt_sice, &
         nt_qice_out=nt_qice, nt_qsno_out=nt_qsno, nt_iage_out=nt_iage, nt_fy_out=nt_fy, &
@@ -2587,25 +2495,6 @@ subroutine init_state
                file=__FILE__, line=__LINE__)
          endif
 
-#ifdef UNDEPRECATE_0LAYER
-         if (.not.heat_capacity) then
-
-            if (nilyr > 1) then
-               write(nu_diag,*) subname//' ERROR: Must have nilyr = 1 if heat_capacity=F'
-               write(nu_diag,*) subname//' ERROR:   nilyr =', nilyr
-               call abort_ice(error_message=subname//' Too many ice layers', &
-                  file=__FILE__, line=__LINE__)
-            endif
-
-            if (nslyr > 1) then
-               write(nu_diag,*) subname//' ERROR: Must have nslyr = 1 if heat_capacity=F'
-               write(nu_diag,*) subname//' ERROR:  nslyr =', nslyr
-               call abort_ice(error_message=subname//' Too many snow layers', &
-                  file=__FILE__, line=__LINE__)
-            endif
-
-         endif   ! heat_capacity = F
-#endif
       endif      ! my_task
 
       !-----------------------------------------------------------------
@@ -2624,12 +2513,6 @@ subroutine init_state
       if (tr_FY)   trcr_depend(nt_FY)    = 0   ! area-weighted first-year ice area
       if (tr_lvl)  trcr_depend(nt_alvl)  = 0   ! level ice area
       if (tr_lvl)  trcr_depend(nt_vlvl)  = 1   ! level ice volume
-#ifdef UNDEPRECATE_CESMPONDS
-      if (tr_pond_cesm) then
-                   trcr_depend(nt_apnd)  = 0           ! melt pond area
-                   trcr_depend(nt_hpnd)  = 2+nt_apnd   ! melt pond depth
-      endif
-#endif
       if (tr_pond_lvl) then
                    trcr_depend(nt_apnd)  = 2+nt_alvl   ! melt pond area
                    trcr_depend(nt_hpnd)  = 2+nt_apnd   ! melt pond depth
@@ -2691,12 +2574,6 @@ subroutine init_state
          nt_strata   (it,2) = 0
       enddo
 
-#ifdef UNDEPRECATE_CESMPONDS
-      if (tr_pond_cesm) then
-         n_trcr_strata(nt_hpnd)   = 1       ! melt pond depth
-         nt_strata    (nt_hpnd,1) = nt_apnd ! on melt pond area
-      endif
-#endif
       if (tr_pond_lvl) then
          n_trcr_strata(nt_apnd)   = 1       ! melt pond area
          nt_strata    (nt_apnd,1) = nt_alvl ! on level ice area
diff --git a/cicecore/cicedynB/general/ice_step_mod.F90 b/cicecore/cicedynB/general/ice_step_mod.F90
index 5742bb2b9..8dd6fe49a 100644
--- a/cicecore/cicedynB/general/ice_step_mod.F90
+++ b/cicecore/cicedynB/general/ice_step_mod.F90
@@ -274,11 +274,7 @@ subroutine step_therm1 (dt, iblk)
          nt_isosno, nt_isoice, nt_rsnw, nt_smice, nt_smliq
 
       logical (kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_iso, tr_aero, tr_pond, tr_pond_cesm, &
-#else
          tr_iage, tr_FY, tr_iso, tr_aero, tr_pond, &
-#endif
          tr_pond_lvl, tr_pond_topo, calc_Tsfc, highfreq, tr_snow
 
       real (kind=dbl_kind) :: &
@@ -304,11 +300,7 @@ subroutine step_therm1 (dt, iblk)
       call icepack_query_tracer_sizes(ntrcr_out=ntrcr)
       call icepack_query_tracer_flags( &
          tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_iso_out=tr_iso, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_aero_out=tr_aero, tr_pond_out=tr_pond, tr_pond_cesm_out=tr_pond_cesm, &
-#else
          tr_aero_out=tr_aero, tr_pond_out=tr_pond, &
-#endif
          tr_pond_lvl_out=tr_pond_lvl, tr_pond_topo_out=tr_pond_topo, &
          tr_snow_out=tr_snow)
       call icepack_query_tracer_indices( &
diff --git a/cicecore/cicedynB/infrastructure/io/io_binary/ice_restart.F90 b/cicecore/cicedynB/infrastructure/io/io_binary/ice_restart.F90
index b2b438ebe..3e7abe3a3 100644
--- a/cicecore/cicedynB/infrastructure/io/io_binary/ice_restart.F90
+++ b/cicecore/cicedynB/infrastructure/io/io_binary/ice_restart.F90
@@ -58,11 +58,7 @@ subroutine init_restart_read(ice_ic)
 
       logical (kind=log_kind) :: &
          solve_zsal, tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, tr_pond_cesm, &
-#else
          tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, &
-#endif
          tr_pond_topo, tr_pond_lvl, tr_brine, tr_snow
 
       character(len=char_len_long) :: &
@@ -87,11 +83,7 @@ subroutine init_restart_read(ice_ic)
          nbtrcr_out=nbtrcr)
       call icepack_query_tracer_flags( &
          tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_lvl_out=tr_lvl, tr_fsd_out=tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iso_out=tr_iso, tr_aero_out=tr_aero, tr_pond_cesm_out=tr_pond_cesm, &
-#else
          tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_pond_lvl_out=tr_pond_lvl, &
          tr_snow_out=tr_snow, tr_brine_out=tr_brine)
       call icepack_warnings_flush(nu_diag)
@@ -236,28 +228,6 @@ subroutine init_restart_read(ice_ic)
          endif
       endif
 
-#ifdef UNDEPRECATE_CESMPONDS
-      if (tr_pond_cesm) then
-         if (my_task == master_task) then
-            n = index(filename0,trim(restart_file))
-            if (n == 0) call abort_ice(subname//'ERROR: pond_cesm restart: filename discrepancy')
-            string1 = trim(filename0(1:n-1))
-            string2 = trim(filename0(n+lenstr(restart_file):lenstr(filename0)))
-            write(filename,'(a,a,a,a)') &
-               string1(1:lenstr(string1)), &
-               restart_file(1:lenstr(restart_file)),'.pond_cesm', &
-               string2(1:lenstr(string2))
-            if (restart_ext) then
-               call ice_open_ext(nu_restart_pond,filename,0)
-            else
-               call ice_open(nu_restart_pond,filename,0)
-            endif
-            read (nu_restart_pond) iignore,rignore,rignore
-            write(nu_diag,*) 'Reading ',filename(1:lenstr(filename))
-         endif
-      endif
-#endif
-
       if (tr_pond_lvl) then
          if (my_task == master_task) then
             n = index(filename0,trim(restart_file))
@@ -424,11 +394,7 @@ subroutine init_restart_write(filename_spec)
 
       logical (kind=log_kind) :: &
          solve_zsal, tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, tr_pond_cesm, &
-#else
          tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, &
-#endif
          tr_pond_topo, tr_pond_lvl, tr_brine, tr_snow
 
       integer (kind=int_kind) :: &
@@ -444,11 +410,7 @@ subroutine init_restart_write(filename_spec)
          nbtrcr_out=nbtrcr)
       call icepack_query_tracer_flags( &
          tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_lvl_out=tr_lvl, tr_fsd_out=tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iso_out=tr_iso, tr_aero_out=tr_aero, tr_pond_cesm_out=tr_pond_cesm, &
-#else
          tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_pond_lvl_out=tr_pond_lvl, &
          tr_snow_out=tr_snow, tr_brine_out=tr_brine)
       call icepack_warnings_flush(nu_diag)
@@ -581,28 +543,6 @@ subroutine init_restart_write(filename_spec)
 
       endif
 
-#ifdef UNDEPRECATE_CESMPONDS
-      if (tr_pond_cesm) then
-
-         write(filename,'(a,a,a,i4.4,a,i2.2,a,i2.2,a,i5.5)') &
-              restart_dir(1:lenstr(restart_dir)), &
-              restart_file(1:lenstr(restart_file)),'.pond_cesm.', &
-              myear,'-',mmonth,'-',mday,'-',msec
-
-         if (restart_ext) then
-            call ice_open_ext(nu_dump_pond,filename,0)
-         else
-            call ice_open(nu_dump_pond,filename,0)
-         endif
-
-         if (my_task == master_task) then
-           write(nu_dump_pond) istep1,timesecs,time_forc
-           write(nu_diag,*) 'Writing ',filename(1:lenstr(filename))
-         endif
-
-      endif
-#endif
-
       if (tr_pond_lvl) then
 
          write(filename,'(a,a,a,i4.4,a,i2.2,a,i2.2,a,i5.5)') &
@@ -871,11 +811,7 @@ subroutine final_restart()
 
       logical (kind=log_kind) :: &
          solve_zsal, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, tr_pond_cesm, &
-#else
          tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, &
-#endif
          tr_pond_topo, tr_pond_lvl, tr_brine, tr_snow
 
       integer (kind=int_kind) :: &
@@ -889,11 +825,7 @@ subroutine final_restart()
          nbtrcr_out=nbtrcr)
       call icepack_query_tracer_flags( &
          tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_lvl_out=tr_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iso_out=tr_iso, tr_aero_out=tr_aero, tr_pond_cesm_out=tr_pond_cesm, &
-#else
          tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_pond_lvl_out=tr_pond_lvl, &
          tr_snow_out=tr_snow, tr_brine_out=tr_brine)
       call icepack_warnings_flush(nu_diag)
@@ -908,9 +840,6 @@ subroutine final_restart()
          if (tr_iage)      close(nu_dump_age)
          if (tr_FY)        close(nu_dump_FY)
          if (tr_lvl)       close(nu_dump_lvl)
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) close(nu_dump_pond)
-#endif
          if (tr_pond_lvl)  close(nu_dump_pond)
          if (tr_pond_topo) close(nu_dump_pond)
          if (tr_snow)      close(nu_dump_snow)
diff --git a/cicecore/cicedynB/infrastructure/io/io_netcdf/ice_restart.F90 b/cicecore/cicedynB/infrastructure/io/io_netcdf/ice_restart.F90
index f647bd96b..ed49a48f5 100644
--- a/cicecore/cicedynB/infrastructure/io/io_netcdf/ice_restart.F90
+++ b/cicecore/cicedynB/infrastructure/io/io_netcdf/ice_restart.F90
@@ -145,11 +145,7 @@ subroutine init_restart_write(filename_spec)
 
       logical (kind=log_kind) :: &
          solve_zsal, skl_bgc, z_tracers, tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, tr_pond_cesm, &
-#else
          tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, &
-#endif
          tr_pond_topo, tr_pond_lvl, tr_brine, tr_snow, &
          tr_bgc_N, tr_bgc_C, tr_bgc_Nit, &
          tr_bgc_Sil, tr_bgc_DMS, &
@@ -185,11 +181,7 @@ subroutine init_restart_write(filename_spec)
          nbtrcr_out=nbtrcr)
       call icepack_query_tracer_flags( &
          tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_lvl_out=tr_lvl, tr_fsd_out=tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iso_out=tr_iso, tr_aero_out=tr_aero, tr_pond_cesm_out=tr_pond_cesm, &
-#else
          tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_pond_lvl_out=tr_pond_lvl, &
          tr_snow_out=tr_snow, tr_brine_out=tr_brine, &
          tr_bgc_N_out=tr_bgc_N, tr_bgc_C_out=tr_bgc_C, tr_bgc_Nit_out=tr_bgc_Nit, &
@@ -416,13 +408,6 @@ subroutine init_restart_write(filename_spec)
             call define_rest_field(ncid,'vlvl',dims)
          end if
 
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) then
-            call define_rest_field(ncid,'apnd',dims)
-            call define_rest_field(ncid,'hpnd',dims)
-         end if
-#endif
-
          if (tr_pond_topo) then
             call define_rest_field(ncid,'apnd',dims)
             call define_rest_field(ncid,'hpnd',dims)
diff --git a/cicecore/cicedynB/infrastructure/io/io_pio2/ice_restart.F90 b/cicecore/cicedynB/infrastructure/io/io_pio2/ice_restart.F90
index 1124cc048..679a2b6e6 100644
--- a/cicecore/cicedynB/infrastructure/io/io_pio2/ice_restart.F90
+++ b/cicecore/cicedynB/infrastructure/io/io_pio2/ice_restart.F90
@@ -151,11 +151,7 @@ subroutine init_restart_write(filename_spec)
           solve_zsal, skl_bgc, z_tracers
 
       logical (kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, tr_pond_cesm, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_iso, tr_aero, &
-#endif
           tr_pond_topo, tr_pond_lvl, tr_brine, tr_snow, &
           tr_bgc_N, tr_bgc_C, tr_bgc_Nit, &
           tr_bgc_Sil, tr_bgc_DMS, &
@@ -191,11 +187,7 @@ subroutine init_restart_write(filename_spec)
       call icepack_query_tracer_sizes(nbtrcr_out=nbtrcr)
       call icepack_query_tracer_flags( &
           tr_iage_out=tr_iage, tr_FY_out=tr_FY, tr_lvl_out=tr_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iso_out=tr_iso, tr_aero_out=tr_aero, tr_pond_cesm_out=tr_pond_cesm, &
-#else
           tr_iso_out=tr_iso, tr_aero_out=tr_aero, &
-#endif
           tr_pond_topo_out=tr_pond_topo, tr_pond_lvl_out=tr_pond_lvl, &
           tr_snow_out=tr_snow, tr_brine_out=tr_brine, &
           tr_bgc_N_out=tr_bgc_N, tr_bgc_C_out=tr_bgc_C, tr_bgc_Nit_out=tr_bgc_Nit, &
@@ -420,13 +412,6 @@ subroutine init_restart_write(filename_spec)
             call define_rest_field(File,'vlvl',dims)
          end if
 
-#ifdef UNDEPRECATE_CESMPONDS
-         if (tr_pond_cesm) then
-            call define_rest_field(File,'apnd',dims)
-            call define_rest_field(File,'hpnd',dims)
-         end if
-#endif
-
          if (tr_pond_topo) then
             call define_rest_field(File,'apnd',dims)
             call define_rest_field(File,'hpnd',dims)
diff --git a/cicecore/drivers/direct/hadgem3/CICE_InitMod.F90 b/cicecore/drivers/direct/hadgem3/CICE_InitMod.F90
index 87dc8d9a1..0b8ed689e 100644
--- a/cicecore/drivers/direct/hadgem3/CICE_InitMod.F90
+++ b/cicecore/drivers/direct/hadgem3/CICE_InitMod.F90
@@ -235,17 +235,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_fsd, read_restart_fsd, &
@@ -260,11 +253,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_fsd, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -283,11 +272,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_fsd_out=tr_fsd)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -346,21 +331,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/direct/hadgem3/CICE_RunMod.F90 b/cicecore/drivers/direct/hadgem3/CICE_RunMod.F90
index 91f7985bd..b67e1a223 100644
--- a/cicecore/drivers/direct/hadgem3/CICE_RunMod.F90
+++ b/cicecore/drivers/direct/hadgem3/CICE_RunMod.F90
@@ -145,11 +145,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_bgc, write_restart_hbrine
       use ice_restart_driver, only: dumpfile
@@ -172,11 +168,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -185,11 +177,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_fsd_out=tr_fsd)
       call icepack_warnings_flush(nu_diag)
@@ -333,9 +321,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_fsd)       call write_restart_fsd
diff --git a/cicecore/drivers/direct/nemo_concepts/CICE_InitMod.F90 b/cicecore/drivers/direct/nemo_concepts/CICE_InitMod.F90
index 87dc8d9a1..0b8ed689e 100644
--- a/cicecore/drivers/direct/nemo_concepts/CICE_InitMod.F90
+++ b/cicecore/drivers/direct/nemo_concepts/CICE_InitMod.F90
@@ -235,17 +235,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_fsd, read_restart_fsd, &
@@ -260,11 +253,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_fsd, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -283,11 +272,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_fsd_out=tr_fsd)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -346,21 +331,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/direct/nemo_concepts/CICE_RunMod.F90 b/cicecore/drivers/direct/nemo_concepts/CICE_RunMod.F90
index ea6a65165..c9875d769 100644
--- a/cicecore/drivers/direct/nemo_concepts/CICE_RunMod.F90
+++ b/cicecore/drivers/direct/nemo_concepts/CICE_RunMod.F90
@@ -145,11 +145,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_bgc, write_restart_hbrine
       use ice_restart_driver, only: dumpfile
@@ -172,11 +168,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -185,11 +177,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_fsd_out=tr_fsd)
       call icepack_warnings_flush(nu_diag)
@@ -333,9 +321,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_fsd)       call write_restart_fsd
diff --git a/cicecore/drivers/mct/cesm1/CICE_InitMod.F90 b/cicecore/drivers/mct/cesm1/CICE_InitMod.F90
index cfc5bece9..a8bf96ad2 100644
--- a/cicecore/drivers/mct/cesm1/CICE_InitMod.F90
+++ b/cicecore/drivers/mct/cesm1/CICE_InitMod.F90
@@ -264,17 +264,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, init_snowtracers, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_snow, read_restart_snow, &
@@ -291,11 +284,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_snow, tr_fsd, tr_iso, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -315,11 +304,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_snow_out=tr_snow, tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -381,21 +366,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/mct/cesm1/CICE_RunMod.F90 b/cicecore/drivers/mct/cesm1/CICE_RunMod.F90
index b96086c6d..5836479b4 100644
--- a/cicecore/drivers/mct/cesm1/CICE_RunMod.F90
+++ b/cicecore/drivers/mct/cesm1/CICE_RunMod.F90
@@ -137,11 +137,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_iso, write_restart_bgc, write_restart_hbrine, &
           write_restart_snow
@@ -166,11 +162,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, tr_snow, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -188,11 +180,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_iso_out=tr_iso, tr_fsd_out=tr_fsd, tr_snow_out=tr_snow)
       call icepack_warnings_flush(nu_diag)
@@ -403,9 +391,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_snow)      call write_restart_snow
diff --git a/cicecore/drivers/nuopc/cmeps/CICE_InitMod.F90 b/cicecore/drivers/nuopc/cmeps/CICE_InitMod.F90
index 3d5e5cc2a..5fbde9cce 100644
--- a/cicecore/drivers/nuopc/cmeps/CICE_InitMod.F90
+++ b/cicecore/drivers/nuopc/cmeps/CICE_InitMod.F90
@@ -222,17 +222,10 @@ subroutine init_restart()
     use ice_grid, only: tmask
     use ice_init, only: ice_ic
     use ice_init_column, only: init_age, init_FY, init_lvl, init_snowtracers, &
-#ifdef UNDEPRECATE_CESMPONDS
-         init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
          init_meltponds_lvl, init_meltponds_topo, &
-#endif
          init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
     use ice_restart_column, only: restart_age, read_restart_age, &
          restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-         restart_pond_cesm, read_restart_pond_cesm, &
-#endif
          restart_pond_lvl, read_restart_pond_lvl, &
          restart_pond_topo, read_restart_pond_topo, &
          restart_snow, read_restart_snow, &
@@ -249,11 +242,7 @@ subroutine init_restart()
          i, j        , & ! horizontal indices
          iblk            ! block index
     logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
          tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
          tr_pond_topo, tr_fsd, tr_iso, tr_aero, tr_brine, tr_snow, &
          skl_bgc, z_tracers, solve_zsal
     integer(kind=int_kind) :: &
@@ -274,11 +263,7 @@ subroutine init_restart()
     call icepack_query_parameters(skl_bgc_out=skl_bgc, &
          z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
     call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
          tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
          tr_snow_out=tr_snow, tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
     call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -340,21 +325,6 @@ subroutine init_restart()
           enddo ! iblk
        endif
     endif
-#ifdef UNDEPRECATE_CESMPONDS
-    ! CESM melt ponds
-    if (tr_pond_cesm) then
-       if (trim(runtype) == 'continue') &
-            restart_pond_cesm = .true.
-       if (restart_pond_cesm) then
-          call read_restart_pond_cesm
-       else
-          do iblk = 1, nblocks
-             call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                  trcrn(:,:,nt_hpnd,:,iblk))
-          enddo ! iblk
-       endif
-    endif
-#endif
     ! level-ice melt ponds
     if (tr_pond_lvl) then
        if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/nuopc/cmeps/CICE_RunMod.F90 b/cicecore/drivers/nuopc/cmeps/CICE_RunMod.F90
index c68583648..e908f509f 100644
--- a/cicecore/drivers/nuopc/cmeps/CICE_RunMod.F90
+++ b/cicecore/drivers/nuopc/cmeps/CICE_RunMod.F90
@@ -121,11 +121,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_iso, write_restart_bgc, write_restart_hbrine, &
           write_restart_snow
@@ -150,11 +146,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, tr_snow, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -172,11 +164,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_iso_out=tr_iso, tr_fsd_out=tr_fsd, tr_snow_out=tr_snow)
       call icepack_warnings_flush(nu_diag)
@@ -382,9 +370,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_snow)      call write_restart_snow
diff --git a/cicecore/drivers/nuopc/dmi/CICE_InitMod.F90 b/cicecore/drivers/nuopc/dmi/CICE_InitMod.F90
index 78d462d4c..22596429d 100644
--- a/cicecore/drivers/nuopc/dmi/CICE_InitMod.F90
+++ b/cicecore/drivers/nuopc/dmi/CICE_InitMod.F90
@@ -271,17 +271,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, init_snowtracers, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_snow, read_restart_snow, &
@@ -298,11 +291,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_snow, tr_fsd, tr_iso, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -322,11 +311,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_snow_out=tr_snow, tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -388,21 +373,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/nuopc/dmi/CICE_RunMod.F90 b/cicecore/drivers/nuopc/dmi/CICE_RunMod.F90
index 6e799723e..77bb7738e 100644
--- a/cicecore/drivers/nuopc/dmi/CICE_RunMod.F90
+++ b/cicecore/drivers/nuopc/dmi/CICE_RunMod.F90
@@ -157,11 +157,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_iso, write_restart_bgc, write_restart_hbrine, &
           write_restart_snow
@@ -184,11 +180,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, tr_snow, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -206,11 +198,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_iso_out=tr_iso, tr_fsd_out=tr_fsd, tr_snow_out=tr_snow)
       call icepack_warnings_flush(nu_diag)
@@ -412,9 +400,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_snow)      call write_restart_snow
diff --git a/cicecore/drivers/standalone/cice/CICE_InitMod.F90 b/cicecore/drivers/standalone/cice/CICE_InitMod.F90
index 07a151a01..9c30b15a3 100644
--- a/cicecore/drivers/standalone/cice/CICE_InitMod.F90
+++ b/cicecore/drivers/standalone/cice/CICE_InitMod.F90
@@ -266,17 +266,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, init_snowtracers, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_snow, read_restart_snow, &
@@ -293,11 +286,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_snow, tr_fsd, tr_iso, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -317,11 +306,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_snow_out=tr_snow, tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -383,21 +368,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/standalone/cice/CICE_RunMod.F90 b/cicecore/drivers/standalone/cice/CICE_RunMod.F90
index 00c7921d1..ae7f7ab1f 100644
--- a/cicecore/drivers/standalone/cice/CICE_RunMod.F90
+++ b/cicecore/drivers/standalone/cice/CICE_RunMod.F90
@@ -151,11 +151,7 @@ subroutine ice_step
       use ice_history_bgc, only: init_history_bgc
       use ice_restart, only: final_restart
       use ice_restart_column, only: write_restart_age, write_restart_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-          write_restart_lvl, write_restart_pond_cesm, write_restart_pond_lvl, &
-#else
           write_restart_lvl, write_restart_pond_lvl, &
-#endif
           write_restart_pond_topo, write_restart_aero, write_restart_fsd, &
           write_restart_iso, write_restart_bgc, write_restart_hbrine, &
           write_restart_snow
@@ -178,11 +174,7 @@ subroutine ice_step
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, tr_fsd, tr_snow, &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_pond_cesm, tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#else
           tr_pond_lvl, tr_pond_topo, tr_brine, tr_iso, tr_aero, &
-#endif
           calc_Tsfc, skl_bgc, solve_zsal, z_tracers, wave_spec
 
       character(len=*), parameter :: subname = '(ice_step)'
@@ -200,11 +192,7 @@ subroutine ice_step
            solve_zsal_out=solve_zsal, z_tracers_out=z_tracers, ktherm_out=ktherm, &
            wave_spec_out=wave_spec)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_aero_out=tr_aero, &
            tr_iso_out=tr_iso, tr_fsd_out=tr_fsd, tr_snow_out=tr_snow)
       call icepack_warnings_flush(nu_diag)
@@ -408,9 +396,6 @@ subroutine ice_step
             if (tr_iage)      call write_restart_age
             if (tr_FY)        call write_restart_FY
             if (tr_lvl)       call write_restart_lvl
-#ifdef UNDEPRECATE_CESMPONDS
-            if (tr_pond_cesm) call write_restart_pond_cesm
-#endif
             if (tr_pond_lvl)  call write_restart_pond_lvl
             if (tr_pond_topo) call write_restart_pond_topo
             if (tr_snow)      call write_restart_snow
diff --git a/cicecore/drivers/unittest/gridavgchk/CICE_InitMod.F90 b/cicecore/drivers/unittest/gridavgchk/CICE_InitMod.F90
index ad355d783..84d1a3a60 100644
--- a/cicecore/drivers/unittest/gridavgchk/CICE_InitMod.F90
+++ b/cicecore/drivers/unittest/gridavgchk/CICE_InitMod.F90
@@ -241,17 +241,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_fsd, read_restart_fsd, &
@@ -267,11 +260,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_fsd, tr_iso, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -290,11 +279,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -354,21 +339,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/drivers/unittest/sumchk/CICE_InitMod.F90 b/cicecore/drivers/unittest/sumchk/CICE_InitMod.F90
index ad355d783..84d1a3a60 100644
--- a/cicecore/drivers/unittest/sumchk/CICE_InitMod.F90
+++ b/cicecore/drivers/unittest/sumchk/CICE_InitMod.F90
@@ -241,17 +241,10 @@ subroutine init_restart
       use ice_grid, only: tmask
       use ice_init, only: ice_ic
       use ice_init_column, only: init_age, init_FY, init_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          init_meltponds_cesm,  init_meltponds_lvl, init_meltponds_topo, &
-#else
           init_meltponds_lvl, init_meltponds_topo, &
-#endif
           init_isotope, init_aerosol, init_hbrine, init_bgc, init_fsd
       use ice_restart_column, only: restart_age, read_restart_age, &
           restart_FY, read_restart_FY, restart_lvl, read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-          restart_pond_cesm, read_restart_pond_cesm, &
-#endif
           restart_pond_lvl, read_restart_pond_lvl, &
           restart_pond_topo, read_restart_pond_topo, &
           restart_fsd, read_restart_fsd, &
@@ -267,11 +260,7 @@ subroutine init_restart
          i, j        , & ! horizontal indices
          iblk            ! block index
       logical(kind=log_kind) :: &
-#ifdef UNDEPRECATE_CESMPONDS
-          tr_iage, tr_FY, tr_lvl, tr_pond_cesm, tr_pond_lvl, &
-#else
           tr_iage, tr_FY, tr_lvl, tr_pond_lvl, &
-#endif
           tr_pond_topo, tr_fsd, tr_iso, tr_aero, tr_brine, &
           skl_bgc, z_tracers, solve_zsal
       integer(kind=int_kind) :: &
@@ -290,11 +279,7 @@ subroutine init_restart
       call icepack_query_parameters(skl_bgc_out=skl_bgc, &
            z_tracers_out=z_tracers, solve_zsal_out=solve_zsal)
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
-#ifdef UNDEPRECATE_CESMPONDS
-           tr_lvl_out=tr_lvl, tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
            tr_lvl_out=tr_lvl, tr_pond_lvl_out=tr_pond_lvl, &
-#endif
            tr_pond_topo_out=tr_pond_topo, tr_aero_out=tr_aero, tr_brine_out=tr_brine, &
            tr_fsd_out=tr_fsd, tr_iso_out=tr_iso)
       call icepack_query_tracer_indices(nt_alvl_out=nt_alvl, nt_vlvl_out=nt_vlvl, &
@@ -354,21 +339,6 @@ subroutine init_restart
             enddo ! iblk
          endif
       endif
-#ifdef UNDEPRECATE_CESMPONDS
-      ! CESM melt ponds
-      if (tr_pond_cesm) then
-         if (trim(runtype) == 'continue') &
-              restart_pond_cesm = .true.
-         if (restart_pond_cesm) then
-            call read_restart_pond_cesm
-         else
-            do iblk = 1, nblocks
-               call init_meltponds_cesm(trcrn(:,:,nt_apnd,:,iblk), &
-                                        trcrn(:,:,nt_hpnd,:,iblk))
-            enddo ! iblk
-         endif
-      endif
-#endif
       ! level-ice melt ponds
       if (tr_pond_lvl) then
          if (trim(runtype) == 'continue') &
diff --git a/cicecore/shared/ice_init_column.F90 b/cicecore/shared/ice_init_column.F90
index 5339aa6ec..06ab79cdb 100644
--- a/cicecore/shared/ice_init_column.F90
+++ b/cicecore/shared/ice_init_column.F90
@@ -44,11 +44,7 @@ module ice_init_column
       private
       public :: init_thermo_vertical, init_shortwave, &
                 init_age, init_FY, init_lvl, init_fsd, &
-#ifdef UNDEPRECATE_CESMPONDS
-                init_meltponds_cesm, init_meltponds_lvl, init_meltponds_topo, &
-#else
                 init_meltponds_lvl, init_meltponds_topo, &
-#endif
                 init_aerosol, init_bgc, init_hbrine, init_zbgc, input_zbgc, &
                 count_tracers, init_isotope, init_snowtracers
 
@@ -547,23 +543,6 @@ subroutine init_lvl(iblk, alvl, vlvl)
 
       end subroutine init_lvl
 
-#ifdef UNDEPRECATE_CESMPONDS
-!=======================================================================
-
-!  Initialize melt ponds.
-
-      subroutine init_meltponds_cesm(apnd, hpnd)
-
-      real(kind=dbl_kind), dimension(:,:,:), intent(out) :: &
-         apnd , & ! melt pond area fraction
-         hpnd     ! melt pond depth
-      character(len=*),parameter :: subname='(init_meltponds_cesm)'
-
-      apnd(:,:,:) = c0
-      hpnd(:,:,:) = c0
-
-      end subroutine init_meltponds_cesm
-#endif
 !=======================================================================
 
 !  Initialize melt ponds.
@@ -1818,11 +1797,7 @@ subroutine count_tracers
       integer (kind=int_kind) :: ntrcr
       logical (kind=log_kind) :: tr_iage, tr_FY, tr_lvl, tr_pond, tr_aero, tr_fsd
       logical (kind=log_kind) :: tr_snow
-#ifdef UNDEPRECATE_CESMPONDS
-      logical (kind=log_kind) :: tr_iso, tr_pond_cesm, tr_pond_lvl, tr_pond_topo
-#else
       logical (kind=log_kind) :: tr_iso, tr_pond_lvl, tr_pond_topo
-#endif
       integer (kind=int_kind) :: nt_Tsfc, nt_sice, nt_qice, nt_qsno, nt_iage, nt_FY
       integer (kind=int_kind) :: nt_alvl, nt_vlvl, nt_apnd, nt_hpnd, nt_ipnd, nt_aero
       integer (kind=int_kind) :: nt_fsd, nt_isosno, nt_isoice
@@ -1907,11 +1882,7 @@ subroutine count_tracers
 
       call icepack_query_tracer_flags(tr_iage_out=tr_iage, tr_FY_out=tr_FY, &
          tr_lvl_out=tr_lvl, tr_aero_out=tr_aero, tr_pond_out=tr_pond, &
-#ifdef UNDEPRECATE_CESMPONDS
-         tr_pond_cesm_out=tr_pond_cesm, tr_pond_lvl_out=tr_pond_lvl, &
-#else
          tr_pond_lvl_out=tr_pond_lvl, &
-#endif
          tr_pond_topo_out=tr_pond_topo, tr_brine_out=tr_brine, tr_fsd_out=tr_fsd, &
          tr_snow_out=tr_snow, tr_iso_out=tr_iso, &
          tr_bgc_Nit_out=tr_bgc_Nit, tr_bgc_Am_out =tr_bgc_Am,  tr_bgc_Sil_out=tr_bgc_Sil,   &
diff --git a/cicecore/shared/ice_restart_column.F90 b/cicecore/shared/ice_restart_column.F90
index a15f9d2c1..86ff170c7 100644
--- a/cicecore/shared/ice_restart_column.F90
+++ b/cicecore/shared/ice_restart_column.F90
@@ -29,9 +29,6 @@ module ice_restart_column
       public ::  write_restart_age,       read_restart_age, &
                  write_restart_FY,        read_restart_FY, &
                  write_restart_lvl,       read_restart_lvl, &
-#ifdef UNDEPRECATE_CESMPONDS
-                 write_restart_pond_cesm, read_restart_pond_cesm, &
-#endif
                  write_restart_pond_lvl,  read_restart_pond_lvl, &
                  write_restart_pond_topo, read_restart_pond_topo, &
                  write_restart_snow,      read_restart_snow, &
@@ -45,9 +42,6 @@ module ice_restart_column
          restart_age      , & ! if .true., read age tracer restart file
          restart_FY       , & ! if .true., read FY tracer restart file
          restart_lvl      , & ! if .true., read lvl tracer restart file
-#ifdef UNDEPRECATE_CESMPONDS
-         restart_pond_cesm, & ! if .true., read meltponds restart file
-#endif
          restart_pond_lvl , & ! if .true., read meltponds restart file
          restart_pond_topo, & ! if .true., read meltponds restart file
          restart_snow     , & ! if .true., read snow tracer restart file
@@ -260,74 +254,6 @@ subroutine read_restart_lvl()
 
       end subroutine read_restart_lvl
 
-#ifdef UNDEPRECATE_CESMPONDS
-!=======================================================================
-!
-! Dumps all values needed for restarting
-!
-! authors Elizabeth C. Hunke, LANL
-!         David A. Bailey, NCAR
-
-      subroutine write_restart_pond_cesm()
-
-      use ice_fileunits, only: nu_dump_pond
-      use ice_state, only: trcrn
-
-      ! local variables
-
-      logical (kind=log_kind) :: diag
-      integer (kind=int_kind) :: nt_apnd, nt_hpnd
-      character(len=*),parameter :: subname='(write_restart_pond_cesm)'
-
-      call icepack_query_tracer_indices(nt_apnd_out=nt_apnd, nt_hpnd_out=nt_hpnd)
-      call icepack_warnings_flush(nu_diag)
-      if (icepack_warnings_aborted()) call abort_ice(error_message=subname, &
-         file=__FILE__, line=__LINE__)
-
-      diag = .true.
-
-      call write_restart_field(nu_dump_pond,0,trcrn(:,:,nt_apnd,:,:),'ruf8', &
-                               'apnd',ncat,diag)
-      call write_restart_field(nu_dump_pond,0,trcrn(:,:,nt_hpnd,:,:),'ruf8', &
-                               'hpnd',ncat,diag)
-
-      end subroutine write_restart_pond_cesm
-
-!=======================================================================
-
-! Reads all values needed for a meltpond volume restart
-!
-! authors Elizabeth C. Hunke, LANL
-!         David A. Bailey, NCAR
-
-      subroutine read_restart_pond_cesm()
-
-      use ice_fileunits, only: nu_restart_pond
-      use ice_state, only: trcrn
-
-      ! local variables
-
-      logical (kind=log_kind) :: &
-         diag
-      integer (kind=int_kind) :: nt_apnd, nt_hpnd
-      character(len=*),parameter :: subname='(read_restart_pond_cesm)'
-
-      call icepack_query_tracer_indices(nt_apnd_out=nt_apnd, nt_hpnd_out=nt_hpnd)
-      call icepack_warnings_flush(nu_diag)
-      if (icepack_warnings_aborted()) call abort_ice(error_message=subname, &
-         file=__FILE__, line=__LINE__)
-
-      diag = .true.
-
-      if (my_task == master_task) write(nu_diag,*) subname,'min/max cesm ponds'
-
-      call read_restart_field(nu_restart_pond,0,trcrn(:,:,nt_apnd,:,:),'ruf8', &
-                              'apnd',ncat,diag,field_loc_center,field_type_scalar)
-      call read_restart_field(nu_restart_pond,0,trcrn(:,:,nt_hpnd,:,:),'ruf8', &
-                              'hpnd',ncat,diag,field_loc_center,field_type_scalar)
-
-      end subroutine read_restart_pond_cesm
-#endif
 !=======================================================================
 !
 ! Dumps all values needed for restarting
diff --git a/configuration/tools/cice4_restart_conversion/convert_restarts.f90 b/configuration/tools/cice4_restart_conversion/convert_restarts.f90
index 30c952510..3717f7d86 100644
--- a/configuration/tools/cice4_restart_conversion/convert_restarts.f90
+++ b/configuration/tools/cice4_restart_conversion/convert_restarts.f90
@@ -39,12 +39,8 @@ program convert_restarts
       logical (kind=log_kind), parameter :: &
          oceanmixed_ice = .true., & ! if true, read/write ocean mixed layer fields
          heat_capacity  = .true., & ! if true, ice has nonzero heat capacity
-#ifdef UNDEPRECATE_0LAYER
-                                    ! if false, use zero-layer thermodynamics
-#else
                                     ! heat_capacity = .false. (zero-layer thermodynamics) 
                                     ! has been deprecated in CICE and Icepack
-#endif
          diag = .true.              ! write min/max diagnostics for fields
 
       ! file names
diff --git a/doc/source/cice_index.rst b/doc/source/cice_index.rst
index 763475992..a8b9d08f1 100644
--- a/doc/source/cice_index.rst
+++ b/doc/source/cice_index.rst
@@ -9,7 +9,7 @@ This index defines many of the symbols used frequently in the CICE model
 code.  Namelist variables are partly included here but also documented
 elsewhere, see Section :ref:`tabnamelist`. All
 quantities in the code are expressed in MKS units (temperatures may take
-either Celsius or Kelvin units).  
+either Celsius or Kelvin units).  Deprecated parameters are listed at the end.
 
 .. csv-table:: *Alphabetical Index*
    :header: " ", " ", " "
@@ -330,7 +330,7 @@ either Celsius or Kelvin units).
    "hpmin", "minimum melt pond depth (shortwave)", "0.005 m"
    "hpondn", "melt pond depth", "m"
    "hs_min", "minimum thickness for which :math:`T_s` is computed", "1.\ :math:`\times`\ 10\ :math:`^{-4}` m"
-   "hs0", "snow depth at which transition to ice occurs (dEdd)", "0.03 m"
+   "hs0", "snow depth at which transition to ice occurs (dEdd)", "m"
    "hs1", "snow depth of transition to pond ice", "0.03 m"
    "hs_ssl", "snow surface scattering layer thickness", "0.04 m"
    "Hstar", "determines mean thickness of ridged ice", "25. m"
@@ -699,7 +699,6 @@ either Celsius or Kelvin units).
    "tr_FY", "if true, use first-year area tracer", ""
    "tr_iage", "if true, use ice age tracer", ""
    "tr_lvl", "if true, use level ice area and volume tracers", ""
-   "tr_pond_cesm", "if true, use CESM melt pond scheme", ""
    "tr_pond_lvl", "if true, use level-ice melt pond scheme", ""
    "tr_pond_topo", "if true, use topo melt pond scheme", ""
    "trcr", "ice tracers", ""
@@ -772,9 +771,11 @@ either Celsius or Kelvin units).
    "zref", "reference height for stability", "10. m"
    "zTrf", "reference height for :math:`T_{ref}`, :math:`Q_{ref}`, :math:`U_{ref}`", "2. m"
    "zvir", "gas constant (water vapor)/gas constant (air) - 1", "0.606"
-
-..
-   ktherm=0 has been deprecated
+   "**Deprecated options and parameters**", "", ""
    "heat_capacity", "if true, use salinity-dependent thermodynamics", "T"
    "kseaice", "thermal conductivity of ice for zero-layer thermodynamics", "2.0 W/m/deg"
    "ktherm", "thermodynamic formulation (0 = zero-layer, 1 = :cite:`Bitz99`, 2 = mushy)", ""        
+   "tr_pond_cesm", "if true, use CESM melt pond scheme", ""
+
+.. 
+   new deprecation comments
diff --git a/doc/source/user_guide/ug_case_settings.rst b/doc/source/user_guide/ug_case_settings.rst
index b2dd54c33..a34c69822 100644
--- a/doc/source/user_guide/ug_case_settings.rst
+++ b/doc/source/user_guide/ug_case_settings.rst
@@ -403,9 +403,6 @@ thermo_nml
    "``Rac_rapid_mode``", "real", "critical Rayleigh number", "10.0"
    "", "", "", ""
 
-..
-   ktherm=0 has been deprecated
-   "", "``0``", "zero-layer thermodynamic model", ""
 
 .. _dynamics_nml:
 
@@ -538,7 +535,7 @@ ponds_nml
    "``frzpnd``", "``cesm``", "CESM pond refreezing forumulation", "``cesm``"
    "", "``hlid``", "Stefan refreezing with pond ice thickness", ""
    "``hp1``", "real", "critical ice lid thickness for topo ponds in m", "0.01"
-   "``hs0``", "real", "snow depth of transition to bare sea ice in m", "0.03"
+   "``hs0``", "real", "snow depth of transition to bare sea ice in m", ""
    "``hs1``", "real", "snow depth of transition to pond ice in m", "0.03"
    "``pndaspect``", "real", "aspect ratio of pond changes (depth:area)", "0.8"
    "``rfracmax``", ":math:`0 \le r_{max} \le 1`", "maximum melt water added to ponds", "0.85"
diff --git a/icepack b/icepack
index 460ddf8de..18fc1c9b7 160000
--- a/icepack
+++ b/icepack
@@ -1 +1 @@
-Subproject commit 460ddf8de6cf9f55572e8bf0e7672f24d6a7ec09
+Subproject commit 18fc1c9b79d81604eafdb1fac4ddd039b78ad390