Skip to content
Snippets Groups Projects
Commit dc061ac3 authored by Evgeny Mortikov's avatar Evgeny Mortikov
Browse files

fixing model init

parent 083f3a52
Branches
No related tags found
No related merge requests found
......@@ -370,6 +370,8 @@ contains
#ifdef USE_CONFIG_PARSER
block
!> just a temporaty to set smth
call set_initial_conditions_kato(grid)
end block
#endif
......
......@@ -103,7 +103,7 @@ module obl_k_epsilon
real, intent(out) :: TKE(nz) !< TKE, [J/kg]
integer :: k !< counter
do k = 2, nz-1
do k = 1, nz
TKE(k) = 1.0 / (param%C_mu)**(1.0/2.0) * 0.001 * 0.001
end do
end subroutine
......@@ -116,7 +116,7 @@ module obl_k_epsilon
real, intent(out) :: eps(nz) !< TKE dissipation rate, [W/kg]
integer :: k !< counter
do k = 2, nz-1
do k = 1, nz
eps(k) = 0.001 * 0.001 * 0.001 / (param%kappa * z) !not full z, but current z - must be fixed
end do
end subroutine
......
......@@ -243,7 +243,11 @@ program obl_main
! ----------------------------------------------------------------------------
if (closure_mode.eq.3 .or. closure_mode.eq.4) then
call TKE_init(TKE, param_k_epsilon, grid%cz)
call eps_init(EPS, param_k_epsilon, grid%cz, grid%height)
call EPS_init(EPS, param_k_epsilon, grid%cz, grid%height)
!< have to define Km, Kh at init
Km = 0.0
Kh = 0.0
endif
! ----------------------------------------------------------------------------
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment