Commit 70a3afeb authored by Dany Dumont's avatar Dany Dumont

Nouveau modele bio_npzd4

parent 8465f2a2
......@@ -16,8 +16,9 @@ bio_sed.F90 \
bio_fasham.F90 \
bio_ismer.F90 \
bio_gsj.F90 \
bio_nocera.F90 \
bio_save.F90 \
bio_nocera.F90 \
bio_npzd4.F90 \
nitrate.F90 \
ammonium.F90 \
hydrocarbon.F90
......@@ -33,6 +34,7 @@ ${LIB}(bio_ismer.o) \
${LIB}(bio_gsj.o) \
${LIB}(bio_fasham.o) \
${LIB}(bio_nocera.o) \
${LIB}(bio_npzd4.o) \
${LIB}(bio_sed.o) \
${LIB}(bio_save.o) \
${LIB}(nitrate.o) \
......
......@@ -40,6 +40,9 @@
use bio_nocera, only : init_bio_nocera,init_var_nocera,var_info_nocera
use bio_nocera, only : light_nocera,do_bio_nocera
use bio_npzd4, only : init_bio_npzd4,init_var_npzd4,var_info_npzd4
use bio_npzd4, only : light_npzd4,do_bio_npzd4
use bio_sed, only : init_bio_sed,init_var_sed,var_info_sed
use bio_mab, only : init_bio_mab,init_var_mab,var_info_mab
......@@ -340,6 +343,16 @@
call var_info_nocera()
case (9) ! The model with four detrital compartments
call init_bio_npzd4(namlst,'bio_npzd4.nml',unit)
call allocate_memory(nlev)
call init_var_npzd4(nlev)
call var_info_npzd4()
case default
stop "bio: no valid biomodel specified in bio.nml !"
end select
......@@ -726,6 +739,9 @@
case (8)
call light_nocera(nlev,bioshade_feedback)
call ode_solver(ode_method,numc,nlev,dt_eff,cc,do_bio_nocera)
case (9)
call light_npzd4(nlev,bioshade_feedback)
call ode_solver(ode_method,numc,nlev,dt_eff,cc,do_bio_npzd4)
end select
end do
......
......@@ -609,6 +609,7 @@
!Diurnal migration of zooplancton as a function of light and phytoplancton
!(food) concentration
!ws(z,ci) = -1.0*w_zmax*tanh(bertha*(par(ci)-parcrit))
if (cc(p,ci) .lt. pmin) then
ws(z,ci) = -1.0*w_zmax*tanh(bertha*(par(ci)-parcrit))
else
......
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