diff --git a/srcF/sfx_sheba_noniterative.f90 b/srcF/sfx_sheba_noniterative.f90
index 8d87893e615e1a12c284e0b8ae3ed618c1f0f26e..36fce06ca2c9acae0dc256a6b58b27759f040fe8 100644
--- a/srcF/sfx_sheba_noniterative.f90
+++ b/srcF/sfx_sheba_noniterative.f90
@@ -329,7 +329,7 @@ contains
 
             if (surface_type==3.or.surface_type==6.) then
                 if (Udyn>u_thsnow) then
-            
+            !write(*,*) 'sfx_snow1', Udyn, zeta, surface_type
             call get_sigma(sigma_r, sigma_w, rho_air, rho_s)
             call get_w_snow(w_snow, sigma_w, g, d_s, nu_air)
             call get_h_salt(h_salt, Udyn)
@@ -338,7 +338,7 @@ contains
                       
             Linv=Linv*((1-S_mean)/(1+sigma_w*S_mean))+(g*w_snow*sigma_w*S_mean/(Udyn**3.0))/(1+sigma_w*S_mean)
             zeta = h * Linv 
-            !write(*,*) 'sfx_snow1', Udyn, Linv, zeta, S_mean
+            
                   
                  
              
@@ -351,7 +351,7 @@ contains
             Udyn = kappa * U / (log(h / z0_m) - (psi_m - psi0_m))
             Tdyn = kappa * dT * Pr_t_0_inv / (log(h / z0_t) - (psi_h - psi0_h))
 
-            write(*,*) 'sfx_snow2', Udyn, zeta, S_mean, Linv
+            !write(*,*) 'sfx_snow2', Udyn, zeta, S_mean, Linv
             endif 
             endif