function PAR=PARabalone2
%
%   define model parameters and values
%       define a case identifier
%      2a: local recruitment
%PAR.case = '2a';
%    2b: recruits by exchange with other population
PAR.case = '2b';

switch PAR.case
  case {'2a'}
%       exchange rate from second population  [1/day]
     PAR.Imm = .0;     % case 2a
%       diffusion rate from second population [1/day]
     PAR.Diff = .0;     % case 2a
%       Maximum reproduction rate for susceptable [number]
     PAR.ReproS = .02;   % case 2a
%       Maximum reproduction rate for infected [number]
     PAR.ReproI = .01;   % case 2a
  case {'2b'}
%       exchange rate from second population  [1/day]
     PAR.Imm = .001;  % case 2b
%       diffusion rate from second population [1/day]
     PAR.Diff = .001;  % case 2b
%       Maximum reproduction rate for susceptable [number]
     PAR.ReproS = .0;   % case 2b
%       Maximum reproduction rate for infected [number]
     PAR.ReproI = .0;   % case 2b
end

%  Carrying capacity in first population [number]
PAR.Carry = 100;
%  Carrying capacity in second population [number]
PAR.Carry2 = 10;
%  Maximum reproduction rate in second population [1/day]
PAR.Repro2 = .01;

PAR.depth = 1;   %   vertical dimension of the volume (m)

% infection rate by contact with infections particles
%  [infected produced/infectious particle/day]
PAR.IPinfect = 0.003 ;  
%PAR.IPinfect = 0 ;  

% infection rate by contact with infectious individuals 
%  [infected produced/susceptible animal/day]
%PAR.Iinfect = 0 ;  
PAR.Iinfect = 0.001 ;  

% infection rate by contact with dead infectious individuals 
%  [infected produced/susceptible animal/day]
%PAR.Dinfect = 0 ;  
PAR.Dinfect = 0.0008 ;  

% mortality rate of infected [1/day]
PAR.Imort = 8.d-2;  

% background mortality for susceptable [1/day]
PAR.Bmort = 0;  

% removal rate of dead [1/day]
PAR.DeadDecay = 1.5 ;  

% infectious particles released by infected [particles/animal/day]
PAR.Irelease  = .015;  

% infectious particles released by dead [particles/animal/day]
PAR.Drelease = 1.0;  

% removal rate of IP from the environment [1/day]
PAR.IPremove = 0.001;