Nothing
#-------------------------------------------------------------------------------
# - Functions to ease the work with FLFLeetExt -
#
# - landStock.f: Extract total landings at age of a stock from a FLFleet obj.
# - discStock.f: Extract total discards at age of a stock from a FLFleet obj.
# - catchStock.f: Extract total catch at age of a stock from a FLFleet obj.
# - landWStock.f, discWStock.f, catchWStock.f: The same but in weight.
#
#
# Created: 30/11/2010 13:23:35
# Author: Dorleta Garcia
# Changed: 12/01/2011 15:48:43
#-------------------------------------------------------------------------------
#-------------------------------------------------------------------------------
# landStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
landStock.f <- function(obj, stock){
aux <- 0
res <- FLQuant()
for(m in obj@metiers){
if(!(stock %in% catchNames(m)))
next
if(aux == 0){
aux <- 1
res <- m@catches[[stock]]@landings.n
res[is.na(res)] <- 0
next
}
res <- res + m@catches[[stock]]@landings.n
res[is.na(res)] <- 0
}
return(res)
}
#-------------------------------------------------------------------------------
# discStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
discStock.f <- function(obj, stock){
aux <- 0
res <- FLQuant()
for(m in obj@metiers){
if(!(stock %in% catchNames(m))) next
if(aux == 0){
aux <- 1
res <- m@catches[[stock]]@discards.n
res[is.na(res)] <- 0
next
}
res <- res + m@catches[[stock]]@discards.n
res[is.na(res)] <- 0
}
return(res)
}
#-------------------------------------------------------------------------------
# catchStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
catchStock.f <- function(obj, stock){
return(discStock.f(obj, stock) + landStock.f(obj,stock))
}
#-------------------------------------------------------------------------------
# landWStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
landWStock.f <- function(obj, stock){
aux <- 0
res <- FLQuant()
for(m in obj@metiers){
if(!(stock %in% catchNames(m))) next
if(aux == 0){
aux <- 1
res <- m@catches[[stock]]@landings.n*m@catches[[stock]]@landings.wt
res[is.na(res)] <- 0
next
}
res <- res + m@catches[[stock]]@landings.n*m@catches[[stock]]@landings.wt
res[is.na(res)] <- 0
}
return(res)
}
#-------------------------------------------------------------------------------
# discWStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
discWStock.f <- function(obj, stock){
aux <- 0
res <- FLQuant()
for(m in obj@metiers){
if(!(stock %in% catchNames(m))) next
if(aux == 0){
aux <- 1
res <- m@catches[[stock]]@discards.n*m@catches[[stock]]@discards.wt
res[is.na(res)] <- 0
next
}
res <- res + m@catches[[stock]]@discards.n*m@catches[[stock]]@discards.wt
res[is.na(res)] <- 0
}
return(res)
}
#-------------------------------------------------------------------------------
# catchWStock.f(obj, stock)
# obj : FLFleet object.
# Output: An FLQuant(age, year, unit, n_season,1,iter) with the total catch of
# the 'stock' in the FLFleets objetc.
#-------------------------------------------------------------------------------
catchWStock.f <- function(obj, stock){
return(discWStock.f(obj, stock) + landWStock.f(obj,stock))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.