function PAR=allometricPAR(Nspecies,Nvar)
%
%   define model parameters and values
%    some parameters are independent of space, others are arrays
PAR.case='1';
 
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)

PAR.mass=[1 10];

PAR.betaA = .00002;PAR.betaB = -.25;

PAR.fecundA = (.2/365);PAR.fecundB = -.25;

PAR.mort1A = (.1/365);PAR.mort1B = -.25;

PAR.mort2A = .00001;PAR.mort2B = -.25;

PAR.ImortA = .01;PAR.ImortB = -.25;

PAR.DeadDecayA = .1;PAR.DeadDecayB = .25;

PAR.releaseA = 10 ;PAR.releaseB = 0 ;

PAR.remove = .1;

%PAR.A = ;PAR.B = ;