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

PAR.case = '3a';

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.025*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.0008*ones(Ny,Nx);  

% mortality rate of infected [1/day]
PAR.Imort = 8.d-2*ones(Ny,Nx);  

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

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

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

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

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

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

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