R/N2O_tier1.R

In emissions: What the package does (short line)

.N2O_tier2 <- function(x){
	
	index <- 1
	 for (EF in x$EF){
		 if ( (EF <  x$EF_L[index]) || (EF > x$EF_H[index])){
			 print.noquote("Warning: emission factor outside uncertainty bounds.")
		 }
		index <- index + 1
	 }
		
	
	sum(x$amount_kg * x$EF, na.rm=TRUE) + sum(x$area_ha * x$EF, na.rm=TRUE)
	
}

	

N2O_tier1 <- function(v, p){
#v is a vector, matrix, or data.frame
#p is a data.frame or list
	if (missing(p)) {
		thisenvir <- new.env()
		path <- system.file("data", package="emissions")
		d <- paste(path, '/EF_N2Otier1.RData', sep='')
		p <- get(load(d, thisenvir), thisenvir)
	}

	if (is.vector(v)) {
		v <- t(as.matrix(v))
	}
	
	stopifnot(all(colnames(v) %in% names(p)))
	varnames <- colnames(p)
	pp <- v[varnames]
	apply(p, 1, function(x, ...)sum(x*pp), na.rm=TRUE)
}