diff --git a/srcF/sfx_z0m_all_surface.f90 b/srcF/sfx_z0m_all_surface.f90
index e39a9133ee4ce44702f5813208cfe92269b8d87d..29185386a7c3a7f204f4094aab59ee9ec59cd9f7 100644
--- a/srcF/sfx_z0m_all_surface.f90
+++ b/srcF/sfx_z0m_all_surface.f90
@@ -66,7 +66,7 @@ module sfx_z0m_all_surface
         a = 0.0
         b = 25.0
         c_min = log(h_charnock) / kappa
-        write(*,*) 'sfx_ch'
+        
         do i = 1, maxiters
             f = c1_charnock - 2.0 * log(Uc)
             do j = 1, maxiters
@@ -107,7 +107,7 @@ module sfx_z0m_all_surface
 
         integer :: i, j
         ! ----------------------------------------------------------------------------
-        write(*,*) 'sfx_ow'
+        
         Uc=U
         C_z0=0.007
         betta_u=0.111
@@ -150,7 +150,7 @@ module sfx_z0m_all_surface
 
         integer :: i, j
         ! ----------------------------------------------------------------------------
-        write(*,*) 'sfx_fe'
+        
         Uc = U
         a = 0.0
         b = 25.0
@@ -200,7 +200,7 @@ subroutine get_dynamic_roughness_map(z0_m, u_dyn0, U, h, z0m_map)
         real, intent(in) :: U               !< abs(wind speed) [m/s]
         ! ----------------------------------------------------------------------------
         real :: h0_m
-        write(*,*) 'sfx_map'
+       
         z0_m=z0m_map  
         h0_m = h / z0_m
         u_dyn0 = U * kappa / log(h0_m)
diff --git a/srcF/sfx_z0t_all_surface.f90 b/srcF/sfx_z0t_all_surface.f90
index 751a0c7b668c9405a393e92b73cc55918e661e12..8b9eb71d96b67fce85b188d1a2610362adf3ca80 100644
--- a/srcF/sfx_z0t_all_surface.f90
+++ b/srcF/sfx_z0t_all_surface.f90
@@ -63,7 +63,7 @@ module sfx_z0t_all_surface
         real, intent(in) :: z0_m                !< aerodynamic roughness [m]
         real, intent(in) :: Re                  !< roughness Reynolds number [n/d]
         ! ----------------------------------------------------------------------------
-        write(*,*) 'sfx_kl_water'
+        
          !--- define B = log(z0_m / z0_t)
         if (Re <= Re_rough_min) then
             B = B1_rough * alog(B3_rough * Re) + B2_rough
@@ -91,7 +91,7 @@ module sfx_z0t_all_surface
         real, intent(in) :: Re                  !< roughness Reynolds number [n/d]
        
         
-       write(*,*) 'sfx_cz'
+      
         B=(kappa*10.0**(-0.4*z0_m/0.07))*(Re**0.45)        !Chen and Zhang
         
          ! --- define roughness [thermal]
@@ -111,7 +111,7 @@ module sfx_z0t_all_surface
         real, intent(in) :: Re                  !< roughness Reynolds number [n/d]
        
        
-        write(*,*) 'sfx_zi'
+        
        
         B=0.1*kappa*(Re**0.5)         !6-Zilitinkevich