Snow Ri iter

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
 

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
      ! --------------------------------------------------------------------------------
 

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