diff --git a/physics/PBL/SATMEDMF/satmedmfvdifq.F b/physics/PBL/SATMEDMF/satmedmfvdifq.F
index 17db798da..9738c009a 100644
--- a/physics/PBL/SATMEDMF/satmedmfvdifq.F
+++ b/physics/PBL/SATMEDMF/satmedmfvdifq.F
@@ -95,9 +95,9 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
 !     &     canopylaixy, canopyfchxy, canopyfrtxy, canopycluxy,        & 
      &     ntqv,dtend,dtidx,index_of_temperature,index_of_x_wind,       &
      &     index_of_y_wind,index_of_process_pbl,gen_tend,ldiag3d,       &
-     &     errmsg,errflg,                                               &
+     &     errmsg,errflg)                                               
 !IVAI: aux arrays
-     &     naux2d,naux3d,aux2d,aux3d)
+!     &     naux2d,naux3d,aux2d,aux3d)
 
 !
       use machine  , only : kind_phys
@@ -321,9 +321,9 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
 !----------------------------------------------
 
 !IVAI
-      integer, intent(in) :: naux2d,naux3d
-      real(kind_phys), intent(inout) :: aux2d(:,:)
-      real(kind_phys), intent(inout) :: aux3d(:,:,:)
+!      integer, intent(in) :: naux2d,naux3d
+!      real(kind_phys), intent(inout) :: aux2d(:,:)
+!      real(kind_phys), intent(inout) :: aux3d(:,:,:)
 !IVAI
 
 !!
@@ -1383,14 +1383,14 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
 !        print*, 'SATMEDMFVDIFQ_RUN: CCLU = ' , cclu (:)
 !        print*, 'SATMEDMFVDIFQ_RUN: CPOPU= ' , cpopu(:)
 ! 2D aux arrays: canopy data in diffusion
-        aux2d(:,1) = cfch (:)
-        aux2d(:,2) = claie(:)
-        aux2d(:,3) = cfrt(:)
+!        aux2d(:,1) = cfch (:)
+!        aux2d(:,2) = claie(:)
+!        aux2d(:,3) = cfrt(:)
 
 ! 3D aux arrays: before canopy correction
-        aux3d(:,:,1) = dkq(:,:)
-        aux3d(:,:,2) = dkt(:,:)
-        aux3d(:,:,3) = dku(:,:)
+!        aux3d(:,:,1) = dkq(:,:)
+!        aux3d(:,:,2) = dkt(:,:)
+!        aux3d(:,:,3) = dku(:,:)
 !IVAI
         do k = 1, km1-1
           do i = 1, im
@@ -1464,7 +1464,7 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
             ELSE ! There is a contiguous forest canopy, apply correction over canopy layers
 
 ! Output contiguous canopy mask
-              if (kount .EQ. 0 ) aux2d(i,5) = aux2d(i,5) + 1
+!              if (kount .EQ. 0 ) aux2d(i,5) = aux2d(i,5) + 1
 
 !Raupauch M. R. A Practical Lagrangian method for relating scalar
 !concentrations to
@@ -1552,9 +1552,9 @@ subroutine satmedmfvdifq_run(im,km,ntrac,ntcw,ntrw,               &
               dku(i,k)= (dku(i,k)/EDDYVEST1) * EDDYVEST_INT !Scale dku to resolved eddy diffusivity
 
 !IVAI: Output contiguos canopy correction bottom layer and 3D
-              if ( kount .EQ.  0)
-     &           aux2d(i,4) = 1./EDDYVEST1  * EDDYVEST_INT
-              aux3d(i,k,4) = 1./EDDYVEST1  * EDDYVEST_INT
+!              if ( kount .EQ.  0)
+!     &           aux2d(i,4) = 1./EDDYVEST1  * EDDYVEST_INT
+!              aux3d(i,k,4) = 1./EDDYVEST1  * EDDYVEST_INT
 !IVAI
 
             END IF   ! contigous canopy conditions
diff --git a/physics/PBL/SATMEDMF/satmedmfvdifq.meta b/physics/PBL/SATMEDMF/satmedmfvdifq.meta
index 64a529bf0..6cc86ad05 100644
--- a/physics/PBL/SATMEDMF/satmedmfvdifq.meta
+++ b/physics/PBL/SATMEDMF/satmedmfvdifq.meta
@@ -755,35 +755,33 @@
   dimensions = ()
   type = integer
   intent = out
-#IVAI
-[naux2d]
-  standard_name = number_of_xy_dimensioned_auxiliary_arrays
-  long_name = number of 2d auxiliary arrays to output (for debugging)
-  units = count
-  dimensions = ()
-  type = integer
-  intent = out
-[naux3d]
-  standard_name  = number_of_xyz_dimensioned_auxiliary_arrays
-  long_name = number of 3d auxiliary arrays to output (for debugging)
-  units = count
-  dimensions = ()
-  type = integer
-  intent = out
-[aux2d]
-  standard_name = auxiliary_2d_arrays
-  long_name = auxiliary 2d arrays to output (for debugging)
-  units = none
-  dimensions = (horizontal_loop_extent,number_of_3d_auxiliary_arrays)
-  type = real
-  kind = kind_phys
-  intent = out
-[aux3d]
-  standard_name = auxiliary_3d_arrays
-  long_name = auxiliary 3d arrays to output (for debugging)
-  units = none
-  dimensions = (horizontal_loop_extent,vertical_layer_dimension,number_of_3d_auxiliary_arrays)
-  type = real
-  kind = kind_phys
-  intent = out
-#IVAI
+#[naux2d]
+#  standard_name = number_of_xy_dimensioned_auxiliary_arrays
+#  long_name = number of 2d auxiliary arrays to output (for debugging)
+#  units = count
+#  dimensions = ()
+#  type = integer
+#  intent = out
+#[naux3d]
+#  standard_name  = number_of_xyz_dimensioned_auxiliary_arrays
+#  long_name = number of 3d auxiliary arrays to output (for debugging)
+#  units = count
+#  dimensions = ()
+#  type = integer
+#  intent = out
+#[aux2d]
+#  standard_name = auxiliary_2d_arrays
+#  long_name = auxiliary 2d arrays to output (for debugging)
+#  units = none
+#  dimensions = (horizontal_loop_extent,number_of_3d_auxiliary_arrays)
+#  type = real
+#  kind = kind_phys
+#  intent = out
+#[aux3d]
+#  standard_name = auxiliary_3d_arrays
+#  long_name = auxiliary 3d arrays to output (for debugging)
+#  units = none
+#  dimensions = (horizontal_loop_extent,vertical_layer_dimension,number_of_3d_auxiliary_arrays)
+#  type = real
+#  kind = kind_phys
+#  intent = out