diff --git a/columnphysics/icepack_therm_mushy.F90 b/columnphysics/icepack_therm_mushy.F90
index f3c3cc42..79146c34 100644
--- a/columnphysics/icepack_therm_mushy.F90
+++ b/columnphysics/icepack_therm_mushy.F90
@@ -72,7 +72,7 @@ subroutine temperature_changes_salinity(dt,                 &
          lhcoef      , & ! transfer coefficient for latent heat
          Tbot        , & ! ice bottom surfce temperature (deg C)
          sss             ! sea surface salinity (PSU)
-         
+
     real (kind=dbl_kind), intent(inout) :: &
          fswsfc      , & ! SW absorbed at ice/snow surface (W m-2)
          fswint          ! SW absorbed in ice interior below surface (W m-2)
@@ -3139,7 +3139,7 @@ subroutine flushing_velocity(zTin,   phi,   &
 
        ! calculate ocean surface height above bottom of ice
        hocn = (ice_mass + hpond * apond * rhofresh + hsn * rhos) / rhow
-       
+
        ! calculate brine height above bottom of ice
        hbrine = hin + hpond
 
@@ -3151,7 +3151,7 @@ subroutine flushing_velocity(zTin,   phi,   &
 
        ! maximum down flow to drain pond
        w_down_max = (hpond * apond) / dt
-       
+
        ! limit flow
        w = min(w,w_down_max)