%
%  Parameters values for the HTVEP model
%

% Phosphorus to another nutrient ratio.

gammaH = 1;
gammaT = 1;
gammaV = 1;
gammaE = 1;

gamma = [gammaH, gammaT, gammaV, gammaE];

% Maximum phosphorus comsumption rates. 

cH = 3; % Kuang04
cT = 6^(2/3); % surface area argument.
cV = 0.05; % Kuang04
cE = 3; % no reference
cVT = 6; % Kuang04
aET = .25; % Kirshner97, range [0,.5]

c = [cH, cT, cV,cE,cVT,aET];

% Semi-saturation constants.

sH = 4; % capacityH = 10 kg % Kuang04
sT = .002; % capacityT = .006 kg guessing % capacityVT = 3 kg % Kuang04
sV = .12; % capacityV = .3 kg guess
sE = .0012; % 100th of the vien
sVT = 1; % capacityT = 3 kg % Kuang04
sET =.0012; %guessed
%sH=6;
s = [sH, sT, sV,sE,sVT,sET];

% Intrinsic death rates;

dH = 1;  % Kuang04
dT = 1;  % Kuang04
dV = 0.2;   % Kuang04
dE = 0.03; % Kirshner97

dH = .1;  % our guess
dT = .1;  % ours guess
dV = 0.01;   % ours guess
dE = 0.03; % Kirshner97
%dE=.003;

d = [dH, dT, dV,dE];

% Intra-, inter-specific competition rates.

mH = .2; 
mT = 5/3; 
mV = .12; 
mE = 12;

m = [mH, mT, mV,mE];

% Dimensionless-Phusphorus flow rate

D =  20.16; %Vlandimir keeps the record.


D = 1008; % .7*60*24
% Incoming (usable) phusphorus mass.
% D=57600;

P0 = .000035; % 1\% phosphorus Kuang04
%P0=.0095;

% Chemothreshold

tauC =1;

dosageC=1;
frequencyC=[7;10;2];
dosageI=0;
frequencyI=[0;0;0];

dosage=[dosageC,dosageI];
frequency=[frequencyC frequencyI];

qH = 0; 
qT = 0; 
qV = 0; 
qE = 0;
qVT =0;

q=[qH,qT,qV,qE,qVT];

QH = 1; 
QT = 1; 
QV = 1;
QE = 1;
QVT = 1;

Q=[QH,QT,QV,QE,QVT];

% Percentage of recycling.

r=1;


