inst/tests2/test_readNcdf.R

In REddyProc: Data processing and plotting utilities of (half-)hourly eddy-covariance measurements

#require(REddyProc)

context('fLoadFluxNCIntoDataframe')

FileName.s <- 'Example_DE-Tha.1996.1998.hourly.nc' 
Dir.s <- 'inst/examples'
VarList.V.s <- c('NEE', 'Rg', 'rH', 'Tair', 'NEE_f')
nRec <- 1000L


expectDefaults <- function(ds){
	expect_that( nrow(ds), equals(nRec) )
	expect_true( all( VarList.V.s %in% colnames(ds)) )
	expect_true( all( c('DateTime','year','month','day','hour') %in% colnames(ds)) )
	expect_true( is.character(sapply(ds[,-1], function(var){attributes(var)$units})) )
	expect_that( as.vector(sapply(ds[,-1], function(var){attributes(var)$varnames})), equals(colnames(ds)[-1]) )
}

test_that('ncdf4',{
	ds <- fLoadFluxNCIntoDataframe(VarList.V.s, FileName.s, Dir.s, count=c(1L,1L,nRec) 
		,NcPackage.s='ncdf4'
	)
	expectDefaults(ds)
})

test_that('RNetCDF',{
			ds <- fLoadFluxNCIntoDataframe(VarList.V.s, FileName.s, Dir.s, count=c(1L,1L,nRec)
				,NcPackage.s='RNetCDF'
			)
			expectDefaults(ds)
		})

test_that('RNetCDF missing col',{
			expect_warning(
				ds <- fLoadFluxNCIntoDataframe(c(VarList.V.s,"bla"), FileName.s, Dir.s, count=c(1L,1L,nRec)
						,NcPackage.s='RNetCDF'
				)
			)
			expectDefaults(ds)
		})

test_that('ncdf4 missing col',{
			expect_warning(
					ds <- fLoadFluxNCIntoDataframe(c(VarList.V.s,"bla"), FileName.s, Dir.s, count=c(1L,1L,nRec)
							,NcPackage.s='ncdf4'
					)
			)
			expectDefaults(ds)
		})

.tmp.f <- function(){
	# file still missing units and currently not in inst/examples
	test_that('RNetCDF Berkeley',{
				path <- "./tmp"
				fname <- "ES-Ln2.HH.2009.2009.nc"
				ds <-  fLoadFluxNCIntoDataframe(c('SW_IN_F', 'SW_IN_POT'), fname, path, 'RNetCDF', count=nRec)			
				expect_true( all( c('SW_IN_F', 'SW_IN_POT') %in% colnames(ds)) )
				expect_true( all( c('DateTime','year','month','day','hour') %in% colnames(ds)) )
				#expect_true( is.character(sapply(ds[,-1], function(var){attributes(var)$units})) )
				expect_that( as.vector(sapply(ds[,-1], function(var){attributes(var)$varnames})), equals(colnames(ds)[-1]) )
			})
}