Commit 4fa08a35 authored by Gwenaelle Gremion's avatar Gwenaelle Gremion
Browse files

Accumulation size_msn au fond

parent dae7253b
......@@ -676,7 +676,7 @@
else if (j==14 .and. bio_model ==10) then
call adv_center(nlev,dt,h,h,ws(j,:),flux, &
flux,_ZERO_,_ZERO_,w_adv_discr,adv_mode_3,cc(j,:))
zeroDivergence3,_ZERO_,_ZERO_,w_adv_discr,adv_mode_3,cc(j,:))
! Specific advection scheme for Bio_model =10 , for nitrate (5) and ammonium (6)
else if (j ==5 .and. bio_model==10) then
do i=1,nlev
......
......@@ -418,6 +418,8 @@ endif
end if
case (zeroDivergence)
cu(0) = cu(1)
case (zeroDivergence3)
cu(0) = ww(1)*Y(1)
case default
FATAL 'unkown lower boundary condition type in adv_center()'
stop
......
......@@ -53,7 +53,7 @@
integer,parameter :: value = 2
integer,parameter :: oneSided = 3
integer,parameter :: zeroDivergence = 4
integer,parameter :: zeroDivergence3 = 5
!
! !REVISION HISTORY:
! Original author(s): Lars Umlauf
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment