From 02c313d23b37f60a82b87f1ab4284002a94f4d5d Mon Sep 17 00:00:00 2001
From: Victor Stepanenko <stepanen@srcc.msu.ru>
Date: Wed, 3 Feb 2021 16:02:41 +0300
Subject: [PATCH] Bug fixed in bottom oxygen flux computation using logarithmic
 law

---
 source/model/oxygen_mod.f90 | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/source/model/oxygen_mod.f90 b/source/model/oxygen_mod.f90
index 999fc4d..f56db45 100644
--- a/source/model/oxygen_mod.f90
+++ b/source/model/oxygen_mod.f90
@@ -39,6 +39,8 @@ use PHYS_CONSTANTS, only : &
 use T_SOLVER_MOD, only : &
 & DIFF_COEF
 
+use ATMOS, only : velfrict_bot
+
 implicit none
 
 !> Grid size group
@@ -102,7 +104,7 @@ allocate (a(1:vector_length),b(1:vector_length),c(1:vector_length), &
 a(:) = 0.; b(:) = 0.; c(:) = 0.; f(:) = 0.; y(:) = 0.
 
 Foxyg1 = sodbot
-!Foxyg1 = kappa*oxyg(gs%M+1,1)/log(0.25*gsp%ddz(gs%M)*ls%h1/z0_bot) !note: explicit scheme for flux
+!Foxyg1 = kappa*oxyg(gs%M+1,1)*velfrict_bot/log(0.25*gsp%ddz(gs%M)*ls%h1/z0_bot) !note: explicit scheme for flux
 
 if (ls%h1 > 0) then
 ! 1-st step of splitting-up scheme - diffusion, explicit scheme for flux
-- 
GitLab