function PAR=PARabalone3b(Nspecies,Nvar)
%
%   define model parameters and values
%    some parameters are independent of space, others are arrays
PAR.case='3b';

PAR.Nvar=Nvar;
PAR.Nspecies=Nspecies;
PAR.iS=1;PAR.iI=2;PAR.iD=3;PAR.iP=4;

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

% infection rate by contact with infections particles
%  [infected produced/infectious particle/day]
%  the convention is IPinfect(source species, infected species)
PAR.IPinfect = zeros(Nspecies,Nspecies);  
%  particles from species 1 infect species 2
PAR.IPinfect(1,2) = 0.005;
PAR.IPinfect(2,1) = 0.0002;
%PAR.IPinfect(1,1) = 0.025;
%PAR.IPinfect(2,2) = 0.025;
%PAR.IPinfect = 0.025*ones(Nspecies,Nspecies);  

% infection rate by contact with infectious individuals 
%  [infected produced/susceptible animal/day]
PAR.Iinfect = zeros(Nspecies,Nspecies);  
%PAR.Iinfect = 0.02*ones(Nspecies,Nspecies);  

% infection rate by contact with dead infected individuals 
%  [infected produced/susceptible animal/day]
PAR.Dinfect = zeros(Nspecies,Nspecies);  
%PAR.Dinfect = 0.0008*ones(Nspecies,Nspecies);  

%  reproduction rate for susceptibles
PAR.Srepro = zeros(Nspecies,1);  %
PAR.Srepro(1) = .09;
PAR.Srepro(2) = .05;
%  reproduction rate for infected
PAR.Irepro = zeros(Nspecies,1);  
PAR.Irepro(1) = .005;
PAR.Irepro(2) = .005;
%   carrying capacity for the population
PAR.Carry = zeros(Nspecies,1);  
PAR.Carry(1)=150;
PAR.Carry(2)=150;

% mortality rate of infected [1/day]
PAR.Imort = 0.04*ones(Nspecies,1);  

% background mortality rate of infected [1/day]
%PAR.Bmort = .08*ones(Nspecies,1);  
PAR.Bmort = zeros(Nspecies,1);  

% removal rate of dead [1/day]
PAR.DeadDecay = 0.5*ones(Nspecies,1);  

% infectious particles released by infected [particles/animal/day]
%PAR.Irelease  = .015*ones(Nspecies,1); 
PAR.Irelease  = zeros(Nspecies,1); 
PAR.Irelease(1)=.001;
PAR.Irelease(2)=.0001;
% infectious particles released by dead [particles/animal/day]
PAR.Drelease = ones(Nspecies,1);  

% removal rate of IP from the environment [1/day]
%PAR.IPremove = 0.001*ones(Nspecies,1);  
PAR.IPremove = zeros(Nspecies,1);  
PAR.IPremove(1) = 0.1;
PAR.IPremove(2) = 0.01;