foreSIGHT (version 0.9.81)

modSimulator: modSimulator

Description

Simulates using weather generator models specified using modelTag.

Usage

modSimulator(datStart=NULL,
             datFinish=NULL,
             modelTag=NULL,      
             parS=NULL,
             seed=NULL,
             file=NULL,
             IOmode="suppress"
             )

Arguments

datStart

A date string in an accepted date format e.g. "01-10-1990".

datFinish

A date string in an accepted date format e.g. "01-10-1990". Must occur after datStart.

modelTag

A character vector of which stochastic models to use to create each climate variable. Supported tags are shown in details below.

parS

A list (names must match supplied modelTags) containing numeric vectors of model parameters.

seed

Numeric. Seed value supplied to weather generator.

file

Character. Specifies filename for simulation output.

IOmode

A string that specifies the input-output mode for the time series = "verbose", "dev" or "suppress".

Details

modelTag provides the main function with requested models. modelTag is vector of any of the following supported models:

• "Simple-ann"a simple annual scaling

• "P-ann-wgen"a four parameter annual rainfall model

• "P-seas-wgen"a 16 parameter seasonal rainfall model

• "P-har6-wgen"a harmonic rainfall model with 6 periods

• "P-har12-wgen"a harmonic rainfall model

• "P-har12-wgen-FS"a harmonic rainfall model where seasonality is fixed (phase angles must be specified via modelInfoMod=list("P-har12-wgen-FS"=fixedPars=c(x,x,x,x))

• "P-har26-wgen"a harmonic rainfall model

• "P-2har26-wgen"a double harmonic rainfall model

• "Temp-har26-wgen"a harmonic temperature model not conditional on rainfall

• "Temp-har26-wgen-wd"a harmonic temperature model dependent on wet or dry day

• "Temp-har26-wgen-wdsd"a harmonic temperature model where standard deviation parameters are dependent on wet or dry day

• "PET-har12-wgen"a harmonic potential evapotranspiration model

• "PET-har26-wgen"a harmonic potential evapotranspiration model

• "PET-har26-wgen-wd"a harmonic potential evapotranspiration model dependent on wet or dry day

• "Radn-har26-wgen"a harmonic solar radiation model (MJ/m2)

Examples

# NOT RUN {
data(tankDat); obs=tank_obs                     #Get observed data
modelTag=c("P-har12-wgen","Temp-har26-wgen")    #Select models
pars=modCalibrator(obs=obs,modelTag=modelTag)   #Calibrate models
sim=modSimulator(datStart="1970-01-01",         #Simulate!
                 datFinish="1999-12-31",
                 modelTag=modelTag,
                 parS=pars,
                 seed=123,
                 file=paste0("tester.csv"),
                 IOmode="verbose")
plot(sim$P[1:365])                             #Plot first year of rainfall
# }