diff --git a/columnphysics/icepack_therm_mushy.F90 b/columnphysics/icepack_therm_mushy.F90
index d593bae8..fc99668b 100644
--- a/columnphysics/icepack_therm_mushy.F90
+++ b/columnphysics/icepack_therm_mushy.F90
@@ -3160,7 +3160,7 @@ subroutine flushing_velocity(zTin,   phi,   &
        if (tr_pond_lvl) then
           w_down_max = (hpond * apnd * alvl) / dt
        else
-         w_down_max = (hpond * apnd) / dt
+          w_down_max = (hpond * apnd) / dt
        endif
 
        ! limit flow
diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90
index f5eac5f2..9737e0f6 100644
--- a/columnphysics/icepack_therm_vertical.F90
+++ b/columnphysics/icepack_therm_vertical.F90
@@ -84,6 +84,7 @@ subroutine thermo_vertical (dt,          aicen,     &
                                   Tsf,         zSin,      &
                                   zqin,        zqsn,      &
                                   apnd,        hpond,     &
+                                  alvl,                   &
                                   flw,         potT,      &
                                   Qa,          rhoa,      &
                                   fsnow,       fpond,     &
@@ -107,7 +108,7 @@ subroutine thermo_vertical (dt,          aicen,     &
                                   congel,      snoice,    &
                                   mlt_onset,   frz_onset, &
                                   yday,        dsnow,     &
-                                  prescribed_ice, alvl)
+                                  prescribed_ice)
 
       real (kind=dbl_kind), intent(in) :: &
          dt      , & ! time step
@@ -2695,6 +2696,7 @@ subroutine icepack_step_therm1(dt,                        &
                                  Tsf=Tsfc        (n), zSin=zSin         (:,n), &
                                  zqin=zqin     (:,n), zqsn=zqsn         (:,n), &
                                  apnd=apnd       (n), hpond=hpnd          (n), &
+                                 alvl=alvl       (n),                          &
                                  flw=flw,             potT=potT,               &
                                  Qa=Qa,               rhoa=rhoa,               &
                                  fsnow=fsnow,         fpond=fpond,             &
@@ -2718,7 +2720,7 @@ subroutine icepack_step_therm1(dt,                        &
                                  congel=congeln  (n), snoice=snoicen      (n), &
                                  mlt_onset=mlt_onset, frz_onset=frz_onset,     &
                                  yday=yday,           dsnow=dsnown        (n), &
-                                 prescribed_ice=prescribed_ice, alvl=alvl (n))
+                                 prescribed_ice=prescribed_ice)
 
             if (icepack_warnings_aborted(subname)) then
                write(warnstr,*) subname, ' ice: Vertical thermo error, cat ', n