diff --git a/config_src/nuopc_driver/mom_cap.F90 b/config_src/nuopc_driver/mom_cap.F90
index 12d70fe4e7..cbc4e2d82c 100644
--- a/config_src/nuopc_driver/mom_cap.F90
+++ b/config_src/nuopc_driver/mom_cap.F90
@@ -1597,7 +1597,6 @@ subroutine ModelAdvance(gcomp, rc)
                 "ocn_stoch.res.", year, month, day, hour, minute, seconds,".nc"
         endif
         call ESMF_LogWrite("MOM_cap: Writing restart :  "//trim(restartname), ESMF_LOGMSG_INFO)
-        if (is_root_pe()) print*,'calling write_stoch_restart_ocn ',trim(restartname)
         call write_stoch_restart_ocn('RESTART/'//trim(restartname))
 #endif
      endif
diff --git a/config_src/solo_driver/MOM_driver.F90 b/config_src/solo_driver/MOM_driver.F90
index 81fb79207f..ba52d9c02a 100644
--- a/config_src/solo_driver/MOM_driver.F90
+++ b/config_src/solo_driver/MOM_driver.F90
@@ -74,9 +74,6 @@ program MOM_main
 
   use MOM_wave_interface, only: wave_parameters_CS, MOM_wave_interface_init
   use MOM_wave_interface, only: MOM_wave_interface_init_lite, Update_Surface_Waves
-#ifdef UFS
-  use get_stochy_pattern_mod, only: write_stoch_restart_ocn
-#endif
 
   implicit none
 
diff --git a/src/core/MOM.F90 b/src/core/MOM.F90
index 5143ffbf77..2e749f5240 100644
--- a/src/core/MOM.F90
+++ b/src/core/MOM.F90
@@ -2370,8 +2370,6 @@ subroutine initialize_MOM(Time, Time_init, param_file, dirs, CS, restart_CSp, &
   me=PE_here()
   master=root_PE()
 
-  !print*,'callling init_stochastic_physics_ocn',maxval(G%geoLatT)
-  if (master) print*,'about to call init_stochastic_physics'
   call init_stochastic_physics_ocn(CS%dt_therm,G%geoLonT,G%geoLatT,G%ied-G%isd+1,G%jed-G%jsd+1,nz,do_epbl,do_sppt,master,mom_comm,iret)
 #endif
 
diff --git a/src/parameterizations/vertical/MOM_diabatic_driver.F90 b/src/parameterizations/vertical/MOM_diabatic_driver.F90
index 78068ba44d..f1dac73e7a 100644
--- a/src/parameterizations/vertical/MOM_diabatic_driver.F90
+++ b/src/parameterizations/vertical/MOM_diabatic_driver.F90
@@ -311,11 +311,7 @@ subroutine diabatic(u, v, h, tv, Hml, fluxes, visc, ADp, CDp, dt, Time_end, &
      t_in=tv%T
      s_in=tv%S
    endif
-   print*,'calling run_stochastic_physics'
    call run_stochastic_physics_ocn(t_rp,sppt_wts)
-   !print*,'in diabatic',CS%do_sppt,size(t_in,1),size(t_in,2),size(t_in,3),size(sppt_wts,1),size(sppt_wts,2)
-   print*,'in diabatic',CS%do_sppt,minval(sppt_wts),maxval(sppt_wts)
-   print*,'in diabatic',CS%do_sppt,minval(t_rp),maxval(t_rp)
   if (CS%id_t_rp1 > 0) then
      call post_data(CS%id_t_rp1, t_rp(:,:,1), CS%diag)
   endif