function PAR=PARabalone3c
%
%   define model parameters and values
%    some parameters are independent of space, others are arrays

PAR.case = '3c';

PAR.Nx = 3;   %   number of cross shore grid boxes (columns)
PAR.Ny = 4;   %   number of along shore grid boxes (rows)

Nx=PAR.Nx;Ny=PAR.Ny;

PAR.depth = ones(Ny,Nx);   %   vertical dimension of the volume (m)

% infection rate by contact with infections particles
%  [infected produced/infectious particle/day]
PAR.IPinfect = 0.003*ones(Ny,Nx);  
%PAR.IPinfect = zeros(Ny,Nx);  

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

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

% rate of symptomatic infection from asymptomatic infection
PAR.Einfect = 0.077*ones(Ny,Nx);
%PAR.Einfect = zeros(Ny,Nx);

% mortality rate of infected [1/day]
%PAR.Imort = 0.0095*ones(Ny,Nx);  
%PAR.Imort(4,3) = 0.1;
PAR.Imort = [0.0095 .0095 .0095; .0025 .0025 .0025; .0005 .0005 .0005; .0005 .0005 .0005];

% background mortality for susceptable [1/day]
%PAR.Bmort = 0.0009*ones(Ny,Nx);  
PAR.Bmort = zeros(Ny,Nx);  

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

% infectious particles released by infected [particles/animal/day]
%PAR.Irelease  = 1.d3*ones(Ny,Nx); 
PAR.Irelease  = 1.d1*ones(Ny,Nx); 

% infectious particles released by dead [particles/animal/day]
PAR.Drelease = zeros(Ny,Nx);  

% removal rate of IP from the environment [1/day]
%PAR.IPremove = 1.5d2*ones(Ny,Nx);  
PAR.IPremove = 2.3d-1*ones(Ny,Nx);  

% rate of transfer of IP in the E/W direction [1/day]
%PAR.Uex = zeros(Ny,Nx-1);  
PAR.Uex = [0.01 0.01 ; 0.01 0.01; -0.01 -0.01; -0.01 -0.01 ];  
%PAR.Uex = [0.05 0.05 ; -0.05 -0.05; -0.05 -0.05; -0.05 -0.05 ];  
%PAR.Uex = [0.005 0.005 ; 0 0; 0 0; -0.005 -0.005 ];  

% rate of transfer of IP in the N/S direction [1/day]
%PAR.Vex = zeros(Ny-1,Nx);  
%PAR.Vex = [-0.1 0  0.1; -0.1 0  0.1; -0.1 0  0.1];
%PAR.Vex = [-0.01 0  0.01; -0.01 0  0.01; -0.01 0  0.01];
PAR.Vex = [-0.1 0  0.1; -0.1 0  0.1; -0.1 0  0.1];