Clean up some comments

CICE-Consortium · Jul 30, 2024 · 13f3be6 · 13f3be6
1 parent 7e37c6d
commit 13f3be6
Show file tree

Hide file tree

Showing 3 changed files with 2 additions and 22 deletions.
diff --git a/columnphysics/icepack_fsd.F90 b/columnphysics/icepack_fsd.F90
@@ -564,9 +564,6 @@ subroutine fsd_lateral_growth (ncat,      nfsd,         &
          vi0new_lat = vi0new * lead_area / (c1 + aice/latsurf_area)
       end if
 
-      ! for history/diagnostics
-!     frazil = vi0new - vi0new_lat
-
       ! lateral growth increment
       if (vi0new_lat > puny) then
          G_radial = vi0new_lat/dt
@@ -588,7 +585,6 @@ subroutine fsd_lateral_growth (ncat,      nfsd,         &
 
       ! Use remaining ice volume as in standard model,
       ! but ice cannot grow into the area that has grown laterally
-!     vi0new = vi0new - vi0new_lat
       tot_latg = SUM(d_an_latg(:))
 
       end subroutine fsd_lateral_growth

diff --git a/columnphysics/icepack_therm_itd.F90 b/columnphysics/icepack_therm_itd.F90
@@ -1632,11 +1632,6 @@ subroutine add_new_ice (ncat,      nilyr,      &
             vi0new = c0
          endif               ! aice0 > puny
 
-         ! volume added to each category from lateral growth
-!        do n = 1, ncat
-!           if (aicen(n) > c0) vin0new(n) = d_an_latg(n) * vicen(n)/aicen(n)
-!        end do
-
          ! combine areal change from new ice growth and lateral growth
          d_an_newi(1)     = ai0new
          d_an_tot(2:ncat) = d_an_latg(2:ncat)

diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
@@ -650,9 +650,6 @@ subroutine frzmlt_bottom_lateral (dt,       ncat,     &
                call icepack_warnings_add(warnstr)
             endif
 
-            ! CMB there was a lot of wasted calcs here, particularly the multiplcataion and then division by aicen(n)
-            ! that has been removed. The results aren't quite bfb without the extra calc but results are
-            ! negligeably different
             bin1_arealoss = -afsdn(1,n)  / floe_binwidth(1) ! when scaled by *G_radialn(n)*dt*aicen(n)
 
             delta_an(n) = c0
@@ -705,11 +702,6 @@ subroutine frzmlt_bottom_lateral (dt,       ncat,     &
       !-----------------------------------------------------------------
       ! Limit bottom and lateral heat fluxes if necessary.
       ! Limit rside so we don't melt laterally more ice than frzmlt permits
-      ! FYI: fside is not yet correct for fsd, furthermore rside is
-      ! what matters in lateral melt. fside is simply used here to limit rside
-      ! and fsd code disregards rside anyway. So new idea is to make fside an upper limit
-      ! on what the lateral heat flux could be given fzmlt and fbot, which gave bfb same
-      ! answers until additional mods were made in icepack_therm_itd to utlize it.
       !-----------------------------------------------------------------
 
          xtmp = frzmlt/(fbot + fside - puny)
@@ -2318,7 +2310,7 @@ subroutine icepack_step_therm1(dt, ncat, nilyr, nslyr,    &
          strocnxT    , & ! ice-ocean stress, x-direction
          strocnyT    , & ! ice-ocean stress, y-direction
          fbot        , & ! ice-ocean heat flux at bottom surface (W/m^2)
-         frzmlt      , & ! freezing/melting potential         (W/m^2)       !CMB notes this is not changing so should be in only I think
+         frzmlt      , & ! freezing/melting potential         (W/m^2)
          sst         , & ! sea surface temperature                (C)
          Tf          , & ! freezing temperature                   (C)
          Tbot        , & ! ice bottom surface temperature     (deg C)
@@ -2600,12 +2592,9 @@ subroutine icepack_step_therm1(dt, ncat, nilyr, nslyr,    &
       endif
 
       !-----------------------------------------------------------------
-      ! Adjust frzmlt to account for ice-ocean heat fluxes since last
+      ! Use frzmlt to account for ice-ocean heat fluxes since last
       !  call to coupler.
       ! Compute lateral and bottom heat fluxes.
-      ! CMB notes this routine does not adust frzmlt! instead fhocn is variable
-      ! CMB sent back to ocn with "actual" heat flux used. It is computed as energy residual
-      ! CMB of Etot change minus top surface fluxes so does not even depend on fbot of fside
       !-----------------------------------------------------------------
 
       call frzmlt_bottom_lateral (dt,        ncat,      &