diff --git a/physics/SFC_Layer/UFS/sfc_diff.f b/physics/SFC_Layer/UFS/sfc_diff.f
index f4a102c91..eb5bd7b5c 100644
--- a/physics/SFC_Layer/UFS/sfc_diff.f
+++ b/physics/SFC_Layer/UFS/sfc_diff.f
@@ -168,7 +168,6 @@ subroutine sfc_diff_run (im,rvrdm1,eps,epsm1,grav,                &  !intent(in)
       errmsg = ''
       errflg = 0
 
-
 !  initialize variables. all units are supposedly m.k.s. unless specified
 !  ps is in pascals, wind is wind speed,
 !  surface roughness length is converted to m from cm
@@ -282,11 +281,11 @@ subroutine sfc_diff_run (im,rvrdm1,eps,epsm1,grav,                &  !intent(in)
 !
             call stability
 !  ---  inputs:
-     &         (z1(i), zvfun(i), gdx, tv1, thv1, wind(i),
-     &         z0max, ztmax_lnd(i), tvs, grav, thsfc_loc,
+     &       (z1(i), zvfun(i), gdx, tv1, thv1, wind(i),
+     &        z0max, ztmax_lnd(i), tvs, grav, thsfc_loc,
 !  ---  outputs:
-     &         rb_lnd(i), fm_lnd(i), fh_lnd(i), fm10_lnd(i), fh2_lnd(i),
-     &         cm_lnd(i), ch_lnd(i), stress_lnd(i), ustar_lnd(i))
+     &        rb_lnd(i), fm_lnd(i), fh_lnd(i), fm10_lnd(i), fh2_lnd(i),
+     &        cm_lnd(i), ch_lnd(i), stress_lnd(i), ustar_lnd(i))
           endif ! Dry points
 
           if (icy(i)) then ! Some ice
@@ -336,11 +335,11 @@ subroutine sfc_diff_run (im,rvrdm1,eps,epsm1,grav,                &  !intent(in)
 !
             call stability
 !  ---  inputs:
-     &       (z1(i), zvfun(i), gdx, tv1, thv1, wind(i),
-     &        z0max, ztmax_ice(i), tvs, grav, thsfc_loc,
+     &     (z1(i), zvfun(i), gdx, tv1, thv1, wind(i),
+     &      z0max, ztmax_ice(i), tvs, grav, thsfc_loc,
 !  ---  outputs:
-     &        rb_ice(i), fm_ice(i), fh_ice(i), fm10_ice(i), fh2_ice(i),
-     &        cm_ice(i), ch_ice(i), stress_ice(i), ustar_ice(i))
+     &      rb_ice(i), fm_ice(i), fh_ice(i), fm10_ice(i), fh2_ice(i),
+     &      cm_ice(i), ch_ice(i), stress_ice(i), ustar_ice(i))
       endif ! Icy points
 
 ! BWG: Everything from here to end of subroutine was after
diff --git a/physics/SFC_Layer/UFS/sfc_nst.f90 b/physics/SFC_Layer/UFS/sfc_nst.f90
index 1dd9d6117..9c3804211 100644
--- a/physics/SFC_Layer/UFS/sfc_nst.f90
+++ b/physics/SFC_Layer/UFS/sfc_nst.f90
@@ -27,7 +27,7 @@ module sfc_nst
   subroutine sfc_nst_run                                          &
        ( im, hvap, cp, hfus, jcal, eps, epsm1, rvrdm1, rd, rhw0,  &  ! --- inputs:
        pi, tgice, sbc, ps, u1, v1, usfco, vsfco, icplocn2atm, t1, &
-       q1, tref, cm, ch, lseaspray, fm, fm10,                                       &
+       q1, tref, cm, ch, lseaspray, fm, fm10,                     &
        prsl1, prslki, prsik1, prslk1, wet, use_lake_model, xlon,  &
        sinlat, stress,                                            &
        sfcemis, dlwflx, sfcnsw, rain, timestep, kdt, solhr,xcosz, &
diff --git a/physics/SFC_Models/Ocean/UFS/sfc_ocean.F b/physics/SFC_Models/Ocean/UFS/sfc_ocean.F
index 505476510..88d23a7aa 100644
--- a/physics/SFC_Models/Ocean/UFS/sfc_ocean.F
+++ b/physics/SFC_Models/Ocean/UFS/sfc_ocean.F
@@ -179,7 +179,6 @@ subroutine sfc_ocean_run                                          &
 !           rho is density, qss is sat. hum. at surface
 
         if ( flag(i) ) then
-          
           if (use_med_flux) then
             q0(i)    = max( q1(i), qmin )
             rho(i)   = prsl1(i) / (rd*t1(i)*(one + rvrdm1*q0(i)))