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

test_zt/zm

parent ef5678a8
Branches
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
srcF/sfx_surface.f90
+ 8
− 8
View file @ fe14a2fa
...@@ -98,20 +98,20 @@ module sfx_surface ...@@ -98,20 +98,20 @@ module sfx_surface
integer, public, parameter :: ocean_z0m_id = z0m_ch !< ocean surface integer, public, parameter :: ocean_z0m_id = z0m_ch !< ocean surface
integer, public, parameter :: land_z0m_id = z0m_map_id !< land surface integer, public, parameter :: land_z0m_id = z0m_map_id !< land surface
integer, public, parameter :: lake_z0m_id = z0m_map_id !< lake surface integer, public, parameter :: lake_z0m_id = z0m_fe !< lake surface
integer, public, parameter :: snow_z0m_id = z0m_map_id !< snow covered surface integer, public, parameter :: snow_z0m_id = z0m_ow !< snow covered surface
integer, public, parameter :: forest_z0m_id = z0m_map_id !< forest csurface integer, public, parameter :: forest_z0m_id = z0m_map_id !< forest csurface
integer, public, parameter :: usersf_z0m_id = z0m_ch !< user surface integer, public, parameter :: usersf_z0m_id = z0m_ch !< user surface
integer, public, parameter :: ice_z0m_id = z0m_ch !< ice surface integer, public, parameter :: ice_z0m_id = z0m_map_id !< ice surface
integer, public, parameter :: ocean_z0t_id = z0t_kl_water !< ocean surface integer, public, parameter :: ocean_z0t_id = z0t_kl_water !< ocean surface
integer, public, parameter :: land_z0t_id = z0t_kl_land !< land surface integer, public, parameter :: land_z0t_id = z0t_zi !< land surface
integer, public, parameter :: lake_z0t_id = z0t_kl_land !< lake surface integer, public, parameter :: lake_z0t_id = z0t_cz !< lake surface
integer, public, parameter :: snow_z0t_id = z0t_kl_land !< snow covered surface integer, public, parameter :: snow_z0t_id = z0t_cz !< snow covered surface
integer, public, parameter :: forest_z0t_id = z0t_kl_land !< forest csurface integer, public, parameter :: forest_z0t_id = z0t_zi !< forest csurface
integer, public, parameter :: usersf_z0t_id = z0t_kl_water !< user surface integer, public, parameter :: usersf_z0t_id = z0t_kl_water !< user surface
integer, public, parameter :: ice_z0t_id = z0t_kl_water !< user surface integer, public, parameter :: ice_z0t_id = z0t_ca !< user surface
! -------------------------------------------------------------------------------- ! --------------------------------------------------------------------------------
real, parameter, private :: kappa = 0.40 !< von Karman constant [n/d] real, parameter, private :: kappa = 0.40 !< von Karman constant [n/d]
......
srcF/sfx_z0m_all_surface.f90
+ 7
− 6
View file @ fe14a2fa
...@@ -66,7 +66,7 @@ module sfx_z0m_all_surface ...@@ -66,7 +66,7 @@ module sfx_z0m_all_surface
a = 0.0 a = 0.0
b = 25.0 b = 25.0
c_min = log(h_charnock) / kappa c_min = log(h_charnock) / kappa
write(*,*) 'sfx_ch'
do i = 1, maxiters do i = 1, maxiters
f = c1_charnock - 2.0 * log(Uc) f = c1_charnock - 2.0 * log(Uc)
do j = 1, maxiters do j = 1, maxiters
...@@ -107,6 +107,7 @@ module sfx_z0m_all_surface ...@@ -107,6 +107,7 @@ module sfx_z0m_all_surface
integer :: i, j integer :: i, j
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
write(*,*) 'sfx_ow'
Uc=U Uc=U
C_z0=0.007 C_z0=0.007
betta_u=0.111 betta_u=0.111
...@@ -149,7 +150,7 @@ module sfx_z0m_all_surface ...@@ -149,7 +150,7 @@ module sfx_z0m_all_surface
integer :: i, j integer :: i, j
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
write(*,*) 'sfx_fe'
Uc = U Uc = U
a = 0.0 a = 0.0
b = 25.0 b = 25.0
...@@ -199,7 +200,7 @@ subroutine get_dynamic_roughness_map(z0_m, u_dyn0, U, h, z0m_map) ...@@ -199,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, intent(in) :: U !< abs(wind speed) [m/s]
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
real :: h0_m real :: h0_m
write(*,*) 'sfx_map'
z0_m=z0m_map z0_m=z0m_map
h0_m = h / z0_m h0_m = h / z0_m
u_dyn0 = U * kappa / log(h0_m) u_dyn0 = U * kappa / log(h0_m)
......
srcF/sfx_z0t_all_surface.f90
+ 3
− 3
View file @ fe14a2fa
...@@ -63,7 +63,7 @@ module sfx_z0t_all_surface ...@@ -63,7 +63,7 @@ module sfx_z0t_all_surface
real, intent(in) :: z0_m !< aerodynamic roughness [m] real, intent(in) :: z0_m !< aerodynamic roughness [m]
real, intent(in) :: Re !< roughness Reynolds number [n/d] real, intent(in) :: Re !< roughness Reynolds number [n/d]
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
write(*,*) 'sfx_kl_water'
!--- define B = log(z0_m / z0_t) !--- define B = log(z0_m / z0_t)
if (Re <= Re_rough_min) then if (Re <= Re_rough_min) then
B = B1_rough * alog(B3_rough * Re) + B2_rough B = B1_rough * alog(B3_rough * Re) + B2_rough
...@@ -91,7 +91,7 @@ module sfx_z0t_all_surface ...@@ -91,7 +91,7 @@ module sfx_z0t_all_surface
real, intent(in) :: Re !< roughness Reynolds number [n/d] 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 B=(kappa*10.0**(-0.4*z0_m/0.07))*(Re**0.45) !Chen and Zhang
! --- define roughness [thermal] ! --- define roughness [thermal]
...@@ -111,7 +111,7 @@ module sfx_z0t_all_surface ...@@ -111,7 +111,7 @@ module sfx_z0t_all_surface
real, intent(in) :: Re !< roughness Reynolds number [n/d] real, intent(in) :: Re !< roughness Reynolds number [n/d]
write(*,*) 'sfx_zi'
B=0.1*kappa*(Re**0.5) !6-Zilitinkevich B=0.1*kappa*(Re**0.5) !6-Zilitinkevich
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment