Commit cbbc906a authored by Karsten Bolding's avatar Karsten Bolding
Browse files

prey-predator: updated yaml file - likely problems with results of jacob_monod...

prey-predator: updated yaml file - likely problems with results of jacob_monod model - now 2 LV models
parent 5ca44c5e
......@@ -27,8 +27,6 @@
type (type_state_variable_id) :: id_prey,id_predator
!
! Model parameters
! b,p,d,r are the parameters in Lotka-Volterra model: b is the natural growth rate of prey in the absence of predation;p measures the impact of predation on prey;
! d is the death (or emigration) rate of predator in the absence of interaction with prey; r is the efficiency of turning predated prey into pred.
real(rk) :: b ! natural growth rate of prey
real(rk) :: p ! impact of predation
real(rk) :: r ! efficiency
......@@ -69,10 +67,10 @@
!-----------------------------------------------------------------------
!BOC
! Register model parameters
call self%get_parameter(self%b,'b','d-1','??',default=1.00_rk,scale_factor=d_per_s)
call self%get_parameter(self%p,'p','d-1','??',default=0.05_rk,scale_factor=d_per_s)
call self%get_parameter(self%r,'r','d-1','??',default=0.02_rk,scale_factor=d_per_s)
call self%get_parameter(self%d,'d','d-1','??',default=0.50_rk,scale_factor=d_per_s)
call self%get_parameter(self%b,'b','d-1','growth rate of prey', default=1.00_rk,scale_factor=d_per_s)
call self%get_parameter(self%p,'p','d-1','impact of predation', default=0.05_rk,scale_factor=d_per_s)
call self%get_parameter(self%r,'r','d-1','growth efficiency rate',default=0.02_rk,scale_factor=d_per_s)
call self%get_parameter(self%d,'d','d-1','death rate', default=0.50_rk,scale_factor=d_per_s)
! Register state variables
call self%register_state_variable(self%id_prey,'prey','mmol/m**3','nutrient', &
......
instances:
lotka_volterra:
LV1:
long_name: lotka_volterra
model: au/lotka_volterra
use: true
parameters:
b: 1.0
p: 0.05
r: 0.02
d: 0.5
b: 1.0
p: 0.05
r: 0.02
d: 0.5
initialization:
prey: 10.
predator: 10.
prey: 10.
predator: 10.
LV2:
long_name: lotka_volterra
model: au/lotka_volterra
use: true
parameters:
b: 1.0
p: 0.05
r: 0.02
d: 0.01
initialization:
prey: 10.
predator: 10.
jacob_monod:
long_name: jacob_monod
model: au/jacob_monod
use: true
use: false
parameters:
V: 1.0
K: 1.0
Y: 0.5
V: 1.0
K: 1.0
Y: 0.5
initialization:
prey: 10.
predator: 10.
prey: 10.
predator: 10.
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