From dc061ac3a5d87cc88beec0d0345cc7cf8b54794d Mon Sep 17 00:00:00 2001
From: Evgeny Mortikov <evgeny.mortikov@gmail.com>
Date: Wed, 18 Dec 2024 13:43:30 +0300
Subject: [PATCH] fixing model init

---
 obl_config.f90    | 2 ++
 obl_k_epsilon.f90 | 4 ++--
 obl_main.f90      | 6 +++++-
 obl_run_kato.f90  | 2 +-
 4 files changed, 10 insertions(+), 4 deletions(-)

diff --git a/obl_config.f90 b/obl_config.f90
index 12c9a46..b5e2808 100644
--- a/obl_config.f90
+++ b/obl_config.f90
@@ -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
diff --git a/obl_k_epsilon.f90 b/obl_k_epsilon.f90
index 26f54b6..78b2874 100644
--- a/obl_k_epsilon.f90
+++ b/obl_k_epsilon.f90
@@ -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
diff --git a/obl_main.f90 b/obl_main.f90
index b88073d..0930760 100644
--- a/obl_main.f90
+++ b/obl_main.f90
@@ -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
     ! ---------------------------------------------------------------------------- 
 
diff --git a/obl_run_kato.f90 b/obl_run_kato.f90
index 8b441f2..650aaa2 100644
--- a/obl_run_kato.f90
+++ b/obl_run_kato.f90
@@ -25,7 +25,7 @@ module obl_run_kato
     ! --------------------------------------------------------------------------------
     ! --------------------------------------------------------------------------------
     
-    
+
     contains
 
     ! --------------------------------------------------------------------------------
-- 
GitLab