open_efts: Creates a EftsDataSet for access to a netCDF EFTS data set

In jmp75/efts: High-Level Functions to Read and Write Ensemble Forecast Time Series in netCDF

Creates a EftsDataSet for access to a netCDF EFTS data set

Description

Creates a EftsDataSet for access to a netCDF EFTS data set

Usage

open_efts(ncfile, writein = FALSE)

Arguments

ncfile	name of the netCDF file, or an object of class 'ncdf4'
writein	if TRUE the data set is opened in write mode

Value

A EftsDataSet object

Examples

library(efts)
ext_data <- system.file('extdata', package='efts')
ens_fcast_file <- file.path(ext_data, 'Upper_Murray_sample_ensemble_rain_fcast.nc')
stopifnot(file.exists(ens_fcast_file))
snc <- open_efts(ens_fcast_file)
(variable_names <- snc$get_variable_names())
(stations_ids <- snc$get_values('station_id'))
nEns <- snc$get_ensemble_size()
nLead <- snc$get_lead_time_count()
td <- snc$get_time_dim()
stopifnot('rain_fcast_ens' %in% variable_names)

ens_fcast_rainfall <- snc$get_ensemble_forecasts('rain_fcast_ens',
  stations_ids[1], start_time=td[2])
names(ens_fcast_rainfall) <- as.character(1:ncol(ens_fcast_rainfall))
plot(ens_fcast_rainfall, legend.loc='right')

snc$close()

jmp75/efts documentation built on Feb. 3, 2023, 2:44 p.m.