Skip to content
Snippets Groups Projects
Commit 41b5e811 authored by Виктория Суязова's avatar Виктория Суязова Committed by Anna Shestakova
Browse files

no_write

parent 10a558a8
Branches
No related tags found
No related merge requests found
......@@ -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)
......
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment