polyDat: Simulation data extracted for 51 south Florida indicator...

polyDatR Documentation

Simulation data extracted for 51 south Florida indicator regions used in development of the Combined Operational Plan

Description

A list containing the output from running extractSimData on an EverForeCast simulation dataset beginning 20201201. The variables are as follows:

  • listOfSims a list of simulation data for each indicator region. The list is composed of a dataframe for each simulation (n = 100); each dataframe has a column for each indicator region (n = 51) and a row for each day (n = 182). Values are mean water depth within the indicator region.

  • listOfLocations a list of simulation data, re-organized. The list is composed of a dataframe for each indicator region (n = 51) with a row for each simulation (n = 100) and a column for each day (n = 181; first day is removed because it is identical in each simulation, creating a problem in cluster analysis)

  • traceDataWide. a somewhat-wide dataframe, with a column for each indicator region (n = 51) and a row for each simulation (n = 100) and date (n = 182). Additional rows specify the simulation number and date. In total there are 18200 rows and 53 columns.

  • traceDataLong a long dataframe, with one row for each unique combination of indicator region (n = 51), day (n = 182), and simulation (n = 100). There are 51x182x100 = 928,200 rows and five columns.

Format

A list with four elements

Examples

summary(polyDat)

## Not run: 
### code used to generate object 
### step 1: identify region(s) of interest
locs   <- IRMap[[2]]
loc    <- locs[locs$INDICATOR %in% 118:119, ]
### step 2: homogenize projections (or let the function do this internally)
# loc  <- spTransform(loc, crs(edenDEM))

### step 3: load EverForeCast data from netCDFs
### available here: 
### https://s3.amazonaws.com/jem.models.headless/Ever4Cast_2020_12/simulations_interpolated.zip
simDat  <- loadSims()

# step 4: extract simulation data at each ROIs
polyDat <- extractSimData(simulationData      = simDat, 
                          targetLocations     = loc, 
                          targetLocationNames = loc$NAME)
# save("polyDat", file =  paste0(here::here(), "/data/polyDat.RData"))


## End(Not run)

troyhill/EvergladesEBM documentation built on Nov. 19, 2023, 8:39 a.m.