diff --git a/srcF/sfx_sheba_noniterative.f90 b/srcF/sfx_sheba_noniterative.f90
index f76028ff54ff52b29b0c7fe93cfa48e0aa6e71b2..bde96f6cf831e1088c599f8979636159cec1d112 100644
--- a/srcF/sfx_sheba_noniterative.f90
+++ b/srcF/sfx_sheba_noniterative.f90
@@ -337,7 +337,7 @@ contains
             Ri_sn = (g * sigma_r * deltaS * h) / U**2
             
             end do
-            !write(*,*) 'sfx_snow1', Ri_sn, Rib, Udyn, S_salt, S_mean
+            write(*,*) 'sfx_snow', Ri_sn, Rib, Udyn, surface_type
         else if (Rib < Rib_conv_lim) then
             ! --- strong instability block
 
diff --git a/srcF/sfx_z0m_all_surface.f90 b/srcF/sfx_z0m_all_surface.f90
index 183403cfaff2ff31ff8133af9c7c90f0a6a6df86..249ce699569617e050b827cfec42e57894756620 100644
--- a/srcF/sfx_z0m_all_surface.f90
+++ b/srcF/sfx_z0m_all_surface.f90
@@ -84,7 +84,7 @@ module sfx_z0m_all_surface
 
         ! --- define dynamic velocity in neutral conditions
         u_dyn0 = Uc / c
-        !write(*,*) 'ch', z0_m, u_dyn0
+        write(*,*) 'ch', z0_m, u_dyn0
     end subroutine
     ! --------------------------------------------------------------------------------
     
@@ -204,7 +204,7 @@ subroutine get_dynamic_roughness_map(z0_m, u_dyn0, U, h, z0m_map)
         z0_m=z0m_map  
         h0_m = h / z0_m
         u_dyn0 = U * kappa / log(h0_m)
-        !write(*,*) 'map', z0_m, u_dyn0
+        write(*,*) 'map', z0_m, u_dyn0
         end subroutine
      ! --------------------------------------------------------------------------------
 
diff --git a/srcF/sfx_z0t_all_surface.f90 b/srcF/sfx_z0t_all_surface.f90
index 39a570e3567fd20b8af52870bdc229b697257f07..c6581f3b82f83bed460de75bfed918af5a7be352 100644
--- a/srcF/sfx_z0t_all_surface.f90
+++ b/srcF/sfx_z0t_all_surface.f90
@@ -49,7 +49,7 @@ module sfx_z0t_all_surface
             B = min(B, B_max_land)
        
        z0_t = z0_m / exp(B)
-       !write(*,*) 'kl_land', z0_t, B
+       write(*,*) 'kl_land', z0_t, B
      end subroutine
     ! --------------------------------------------------------------------------------  
         
@@ -75,7 +75,7 @@ module sfx_z0t_all_surface
             B = min(B, B_max_ocean)
        
        z0_t = z0_m / exp(B)
-       !write(*,*) 'kl_water', z0_t, B
+       write(*,*) 'kl_water', z0_t, B
      end subroutine