From ea432eb5425858ce07b563e1f6b877203d6fb2dc Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?=D0=92=D0=B8=D0=BA=D1=82=D0=BE=D1=80=D0=B8=D1=8F=20=D0=A1?= =?UTF-8?q?=D1=83=D1=8F=D0=B7=D0=BE=D0=B2=D0=B0?= <viktoriasuazova@MacBook-Pro-Viktoria.local> Date: Thu, 31 Oct 2024 01:38:01 +0300 Subject: [PATCH] h0_m u_dyn0 after z0m call --- srcF/sfx_most_snow.f90 | 25 +++++++++++++------------ 1 file changed, 13 insertions(+), 12 deletions(-) diff --git a/srcF/sfx_most_snow.f90 b/srcF/sfx_most_snow.f90 index 008fd65..06e0026 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) -- GitLab