%  PlotParam.m
%
%  plot the parameters as a function of body mass.
%
mass=10.^[-3:3];
Nvar=4;  %  number of variables in the model
Nspecies=2;  %  number of species in the model
PAR=allometricPAR(Nspecies,Nvar);  %   define model parameters


%  define parameters
   fecund=PAR.fecundA*mass.^PAR.fecundB;
   mort1=PAR.mort1A*mass.^PAR.mort1B;
   mort2=PAR.mort2A*mass.^PAR.mort2B;
   imort=PAR.ImortA*mass.^PAR.ImortB;
   Ddecay=PAR.DeadDecayA*mass.^PAR.DeadDecayB;
   release=PAR.releaseA*mass.^PAR.releaseB;
   beta=PAR.betaA*mass.^PAR.betaB;

figure
loglog(mass,fecund,'k-',mass,mort1,'k--',mass,mort2,'k-.')
title('Fecundity (-), mort1(--), mort2(-.)')
xlabel('mass (arbitrary units)')
ylabel('Rates (per day)')

figure
loglog(mass,beta,'k-',mass,release,'k--',...
mass,Ddecay,'k-.',mass,imort,'k:')
title('beta (-) release(--) DeadDecay(-.) Imort(..) ')
xlabel('mass (arbitrary units)')
ylabel('Rates (per day)')