diff --git a/CHANGELOG.md b/CHANGELOG.md
index 59dd3e7e..fca6a5dd 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -155,7 +155,8 @@ sequential patch number (starting from `01`).
     were written by default.
   - FMS2_io does not do this. Users can specify real long_names and units by calling register_variable_attribute.
 ### Removed
-- FMS_io was almost completely removed from FMScoupler and replaced with fms2_io.
+- FMS_io was almost completely removed from FMScoupler and replaced with fms2_io
+	- The only remaining usage is `fms_io_exit` calls. These calls are made at the end of each coupler driver and only called if the CPP macro `use_deprecated_io` is set.
 ### Tag Commit Hashes
 - 2021.02-alpha1 (c1c8044a6c3efb8ddbbd01a3769bbf2610b34937)
 - 2021.02-alpha2 (c1c8044a6c3efb8ddbbd01a3769bbf2610b34937)
diff --git a/simple/flux_exchange.F90 b/simple/flux_exchange.F90
index 97083b64..307a45b5 100644
--- a/simple/flux_exchange.F90
+++ b/simple/flux_exchange.F90
@@ -617,14 +617,6 @@ subroutine flux_up_to_atmos (Time, Land, Ice, Boundary )
      Boundary%dt_tr(:,:,isphum) = f_q_delt_n  + dt_t_surf*e_q_n
   endwhere
 
-#ifdef DEBUG_COUPLER_FLUX_TO_ATMOS
-print *, 'PE,dt_t(L)(mn,mx)=',fms_mpp_pe(),minval(Boundary%dt_t,mask=Land%mask(:,:,1)), &
-         maxval(Boundary%dt_t,mask=Land%mask(:,:,1))
-print *, 'PE,dt_q(L)(mn,mx)=',fms_mpp_pe(),minval(Boundary%dt_q,mask=Land%mask(:,:,1)), &
-         maxval(Boundary%dt_q,mask=Land%mask(:,:,1))
-print *, 'PE,dt_t(I)(mn,mx)=',fms_mpp_pe(),minval(Boundary%dt_t,mask=Ice%mask),maxval(Boundary%dt_t,mask=Ice%mask)
-print *, 'PE,dt_q(I)(mn,mx)=',fms_mpp_pe(),minval(Boundary%dt_q,mask=Ice%mask),maxval(Boundary%dt_q,mask=Ice%mask)
-#endif
 
 !=======================================================================
 !-------------------- diagnostics section ------------------------------