R/a13.eq.ac.R

In EcoTroph: An Implementation of the EcoTroph Ecosystem Modelling Approach

#'  function used within the multiplier effort analysis
#' @param compteur counter
#' @param ET_Main output of the create.ETmain function
#' @param biom.mf Effort multiplier applied on biomass
#' @param fish.m.ac Effort multiplier applied on natural mortality
#' @param TopD Parameters of the formula
#' @param FormD Parameters of the formula
#' @param range.TLpred Range of predators trophic level
#' @export
a13.eq.ac=function(compteur,ET_Main,biom.mf,fish.m.ac,TopD,FormD,range.TLpred){
  kin.ac=(ET_Main[compteur, "Kin_acc"] - ET_Main[compteur, "Fish_mort_acc"]) * (1 + TopD[compteur] * (sum(biom.mf[(compteur + range.TLpred[1]):(compteur + range.TLpred[2])])^(FormD[compteur]) - sum(ET_Main[(compteur+ range.TLpred[1]):(compteur + range.TLpred[2]), "B"])^(FormD[compteur]))/(sum(ET_Main[(compteur + range.TLpred[1]):(compteur + range.TLpred[2]), "B"])^(FormD[compteur]))) +  fish.m.ac[compteur]
  return(kin.ac)
}