From e8c79844061992a40e7062bd64d4b94d1013dcc3 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, 13 Feb 2025 21:26:41 +0300 Subject: [PATCH] Snow Ri iter --- srcF/sfx_sheba_noniterative.f90 | 2 +- srcF/sfx_z0m_all_surface.f90 | 4 ++-- srcF/sfx_z0t_all_surface.f90 | 4 ++-- 3 files changed, 5 insertions(+), 5 deletions(-) diff --git a/srcF/sfx_sheba_noniterative.f90 b/srcF/sfx_sheba_noniterative.f90 index f76028f..bde96f6 100644 --- a/srcF/sfx_sheba_noniterative.f90 +++ b/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 diff --git a/srcF/sfx_z0m_all_surface.f90 b/srcF/sfx_z0m_all_surface.f90 index 183403c..249ce69 100644 --- a/srcF/sfx_z0m_all_surface.f90 +++ b/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 ! -------------------------------------------------------------------------------- diff --git a/srcF/sfx_z0t_all_surface.f90 b/srcF/sfx_z0t_all_surface.f90 index 39a570e..c6581f3 100644 --- a/srcF/sfx_z0t_all_surface.f90 +++ b/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 -- GitLab