diff --git a/srcF/sfx_most_snow.f90 b/srcF/sfx_most_snow.f90 index 008fd65742551386c28fe5134941f030d661b55d..06e0026c5a54a07261cca35fd3bc2132f52713dc 100644 --- a/srcF/sfx_most_snow.f90 +++ b/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) - ! --- define relative height + ! 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)