Set useKPP = .false. (default) and comment call to diag_grid_storage_end

src/parameterizations/vertical/MOM_diabatic_driver.F90

@@ -151,7 +151,7 @@ module MOM_diabatic_driver
   real    :: evap_CFL_limit = 0.8    !< The largest fraction of a layer that can be
                                      !! evaporated in one time-step (non-dim).
 
-  logical :: useKPP                  !< use CVmix/KPP diffusivities and non-local transport
+  logical :: useKPP = .false.        !< use CVmix/KPP diffusivities and non-local transport
   logical :: salt_reject_below_ML    !< If true, add salt below mixed layer (layer mode only)
   logical :: KPPisPassive            !< If true, KPP is in passive mode, not changing answers.
   logical :: debug                   !< If true, write verbose checksums for debugging purposes.
@@ -2450,9 +2450,15 @@ subroutine diabatic_driver_end(CS)
     deallocate(CS%optics)
   endif
 
-  call diag_grid_storage_end(CS%diag_grids_prev)
+  ! GMM, the following is commented out because arrays in
+  ! CS%diag_grids_prev are neither pointers or allocatables
+  ! and, therefore, cannot be deallocated.
+
+  !call diag_grid_storage_end(CS%diag_grids_prev)
+
   if (associated(CS)) deallocate(CS)
 
+
 end subroutine diabatic_driver_end