From 41b5e81184461a10922e70632837f031ca9b8cb2 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: Tue, 4 Feb 2025 06:21:19 +0300
Subject: [PATCH] no_write
---
srcF/sfx_z0m_all_surface.f90 | 8 ++++----
srcF/sfx_z0t_all_surface.f90 | 6 +++---
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/srcF/sfx_z0m_all_surface.f90 b/srcF/sfx_z0m_all_surface.f90
index e39a913..2918538 100644
--- a/srcF/sfx_z0m_all_surface.f90
+++ b/srcF/sfx_z0m_all_surface.f90
@@ -66,7 +66,7 @@ module sfx_z0m_all_surface
a = 0.0
b = 25.0
c_min = log(h_charnock) / kappa
- write(*,*) 'sfx_ch'
+
do i = 1, maxiters
f = c1_charnock - 2.0 * log(Uc)
do j = 1, maxiters
@@ -107,7 +107,7 @@ module sfx_z0m_all_surface
integer :: i, j
! ----------------------------------------------------------------------------
- write(*,*) 'sfx_ow'
+
Uc=U
C_z0=0.007
betta_u=0.111
@@ -150,7 +150,7 @@ module sfx_z0m_all_surface
integer :: i, j
! ----------------------------------------------------------------------------
- write(*,*) 'sfx_fe'
+
Uc = U
a = 0.0
b = 25.0
@@ -200,7 +200,7 @@ subroutine get_dynamic_roughness_map(z0_m, u_dyn0, U, h, z0m_map)
real, intent(in) :: U !< abs(wind speed) [m/s]
! ----------------------------------------------------------------------------
real :: h0_m
- write(*,*) 'sfx_map'
+
z0_m=z0m_map
h0_m = h / z0_m
u_dyn0 = U * kappa / log(h0_m)
diff --git a/srcF/sfx_z0t_all_surface.f90 b/srcF/sfx_z0t_all_surface.f90
index 751a0c7..8b9eb71 100644
--- a/srcF/sfx_z0t_all_surface.f90
+++ b/srcF/sfx_z0t_all_surface.f90
@@ -63,7 +63,7 @@ module sfx_z0t_all_surface
real, intent(in) :: z0_m !< aerodynamic roughness [m]
real, intent(in) :: Re !< roughness Reynolds number [n/d]
! ----------------------------------------------------------------------------
- write(*,*) 'sfx_kl_water'
+
!--- define B = log(z0_m / z0_t)
if (Re <= Re_rough_min) then
B = B1_rough * alog(B3_rough * Re) + B2_rough
@@ -91,7 +91,7 @@ module sfx_z0t_all_surface
real, intent(in) :: Re !< roughness Reynolds number [n/d]
- write(*,*) 'sfx_cz'
+
B=(kappa*10.0**(-0.4*z0_m/0.07))*(Re**0.45) !Chen and Zhang
! --- define roughness [thermal]
@@ -111,7 +111,7 @@ module sfx_z0t_all_surface
real, intent(in) :: Re !< roughness Reynolds number [n/d]
- write(*,*) 'sfx_zi'
+
B=0.1*kappa*(Re**0.5) !6-Zilitinkevich
--
GitLab