Chunks Long Simulations For Comparison With Hindcasts
This function reorganizes a long run (historical typically) with only one
start date into chunks corresponding to a set of start dates. The expected
input structure is the one output from Load() with 4 to 7 dimensions.
Histo2Hindcast(varin, sdatesin, sdatesout, nleadtimesout)
varin |
Array of model or observational data with dimensions: |
sdatesin |
Start date of the input matrix 'YYYYMMDD'. |
sdatesout |
List of start dates of the output matrix c('YYYYMMDD', 'YYYYMMDD', ...). |
nleadtimesout |
Number of leadtimes in the output matrix. |
An array with the same number of dimensions as varin, the same dimensions 1 and 2 and potentially the same dimensions 5 to 7. Dimensions 3 and 4 are set by the arguments sdatesout and nleadtimesout.
History:
0.1 - 2012-11 (V. Guemas, vguemas@ic3.cat) - Original code
1.0 - 2013-09 (N. Manubens, nicolau.manubens@ic3.cat) - Formatting to CRAN
# See examples on Load() to understand the first lines in this example
## Not run:
data_path <- system.file('sample_data', package = 's2dverification')
exp <- list(
name = 'experiment',
path = file.path(data_path, 'model/$EXP_NAME$/monthly_mean',
'$VAR_NAME$_3hourly/$VAR_NAME$_$START_DATES$.nc')
)
obs <- list(
name = 'observation',
path = file.path(data_path, 'observation/$OBS_NAME$/monthly_mean',
'$VAR_NAME$/$VAR_NAME$_$YEAR$$MONTH$.nc')
)
# Now we are ready to use Load().
startDates <- c('19851101', '19901101', '19951101', '20001101', '20051101')
sampleData <- Load('tos', list(exp), list(obs), startDates,
leadtimemin = 1, leadtimemax = 4, output = 'lonlat',
latmin = 27, latmax = 48, lonmin = -12, lonmax = 40)
## End(Not run)
start_dates_out <- c('19901101', '19911101', '19921101', '19931101', '19941101')
leadtimes_per_startdate <- 12
experimental_data <- Histo2Hindcast(sampleData$mod, startDates[1],
start_dates_out, leadtimes_per_startdate)
observational_data <- Histo2Hindcast(sampleData$obs, startDates[1],
start_dates_out, leadtimes_per_startdate)
PlotAno(experimental_data, observational_data, start_dates_out,
toptitle = paste('anomalies reorganized into shorter chunks'),
ytitle = 'K', fileout='tos_histo2hindcast.eps')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.