h0_m u_dyn0 after z0m call

srcF/sfx_most_snow.f90

@@ -172,19 +172,20 @@ contains
             surface_type = surface_land
         end if
 
-        if (surface_type == surface_ocean) then
-            ! --- define surface roughness [momentum] & dynamic velocity in neutral conditions
-            call get_charnock_roughness(z0_m, u_dyn0, U, h, numerics%maxiters_charnock)
+       ! if (surface_type == surface_ocean) then
+       !     ! --- define surface roughness [momentum] & dynamic velocity in neutral conditions
+        !    call get_charnock_roughness(z0_m, u_dyn0, U, h, numerics%maxiters_charnock)
+        !    ! --- define relative height
+        !    h0_m = h / z0_m
+        !endif
+        !if (surface_type == surface_land) then
+        !    ! --- define relative height
+        !    h0_m = h / z0_m
+        !    ! --- define dynamic velocity in neutral conditions
+        !    u_dyn0 = U * kappa / log(h0_m)
+        !end if
        ! --- define relative height
             h0_m = h / z0_m
-        endif
-        if (surface_type == surface_land) then
-            ! --- define relative height
-            h0_m = h / z0_m
-            ! --- define dynamic velocity in neutral conditions
-            u_dyn0 = U * kappa / log(h0_m)
-        end if
-
         ! --- define thermal roughness & B = log(z0_m / z0_h)
         Re = u_dyn0 * z0_m / nu_air
         call get_thermal_roughness(z0_t, B, z0_m, Re, surface_type)