diff --git a/source/model/oxygen_mod.f90 b/source/model/oxygen_mod.f90
index 999fc4db35011f6fe945df500c87154deb9fcd48..f56db4569d2e1dcc39bd0829d4a3386b398be75c 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