R/nr_feed.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_feed
#'
#' @param fulldata data to perform actions on
#'
#' @author Jan-Philipp Dietrich
nr_feed<-function(fulldata){
    kcr<-expand.set(fulldata,"kcr")
    kcr<-kcr[-which(kcr=="foddr")]
    kli<-expand.set(fulldata,"kli")
    concentrate_feed<-setNames(dimSums(nr_use(fulldata=fulldata)[,,"feed"][,,kcr],dim=c(3.1,3.2)),"concentrate")
    fodder_feed<-setNames(dimSums(nr_use(fulldata=fulldata)[,,"feed"][,,"foddr"],dim=c(3.1,3.2)),"fodder")
    pasture_feed<-setNames(dimSums(collapseNames(readGDX(fulldata, "ov_dem_feed", format="first_found")[,,"pasture"][,,"level"]),dim=3.1)*collapseNames(readGDX(fulldata, "im_attributes_harvest", format="first_found")[,,"nr"][,,"pasture"]),"pasture")
    residues_feed<-setNames(dimSums(nr_ag_res_use(fulldata=fulldata)[,,"feed"],dim=c(3.1,3.2)),"residues")
    convby_feed<-setNames(dimSums(nr_convby_use(fulldata=fulldata)[,,"feed"],dim=3.1),"convby")
    animal_feed<-setNames(dimSums(nr_use(fulldata=fulldata)[,,"feed"][,,kli],dim=c(3.1,3.2)),"animal_feed")
    scavenging_feed<-setNames(dimSums(collapseNames(readGDX(fulldata,"ov_dem_feed", format="first_found")[,,"scavenging"][,,"level"]),dim=3.1)*collapseNames(readGDX(fulldata, "im_attributes_harvest", format="first_found")[,,"nr"][,,"pasture"]),"scavenging")
    out<-mbind(concentrate_feed,fodder_feed,pasture_feed,residues_feed,convby_feed,animal_feed,scavenging_feed)
    return(out)                           
}