R/nr_fix_crop.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_fix_crop
#'
#' @param fulldata data to perform actions on
#'
#' @author Jan-Philipp Dietrich
# biological fixation by crops
nr_fix_crop<-function(fulldata) {
    fixrate<-readGDX(fulldata,name="f51_nr_fix_ndfa", format="first_found")
    kcr<-expand.set(fulldata,"kcr")
    out <- fixrate[,,kcr] * (
        as.magpie(nr_harvest(fulldata, nutrient="nr")[,,kcr] 
                  + nr_ag_res(fulldata, nutrient="nr")[,,kcr] 
                  + nr_bg_res(fulldata, nutrient="nr")[,,kcr])
    )
    return(out)
}