plot_full_length_timeseries {StrathE2E2} | R Documentation |
Multi-panel time series plots of state variable data for the full duration of a model run generated by the StrathE2E() function
plot_full_length_timeseries(model, results)
model |
R-list object defining the model configuration compiled by the read_model() function |
results |
R-list object containing model results generated by the StrathE2E() function |
The function plots a multi-panel page of time series plots of daily state variable values aggregated over the inshore and offshore zones of the domain and over sediment classes, for the full duration of a model run. Currently the masses of macrophytes, corpses and discards are not included in the plots due to space constraints.
Be warned that if the run is more than about 10 years then the plot becomes extremely compressed and messy. It is not intended to be of publication quality. The intention is to provide a quick-look diagnostic of trends in the state variables. This is useful to assess whether the model is close to stationary state or not.
Units of the plotted varaiables mMN ie. mass in the model domain without any scaling to zone-area or layer thickness
Graphical display in a new graphics window
read_model
, StrathE2E
, plot_final_year
, plot_final_year_HTL_inshore_vs_offshore
, plot_final_year_LTL_inshore_vs_offshore
, plot_sediment_habitats_final_year
# Load the 2003-2013 version of the North Sea model supplied with the package: model <- read_model("North_Sea", "2003-2013") #Run the model and generate the results object results <- StrathE2E(model,nyears=10) # Plot the time series of outputs plot_full_length_timeseries(model, results) # Time series plot of state varaiables in a new window leaving any existing windows open dev.new() plot_full_length_timeseries(model, results) #Direct the graphics output to a file (Windows OS)... pdf("C:/Users/username/Documents/Foldername/plot.pdf",width=8,height=6) # or jpeg("plot.jpg"), png("plot.png") plot_full_length_timeseries(model, results) dev.off() # Read the North Sea/1970-1999 model and set the identifier for output files to "Baseline", run the model and plot the full length output model <- read_model("North_Sea", "1970-1999",model.ident="Baseline") results <- StrathE2E(model,nyears=10) plot_full_length_timeseries(model, results) # Create a new scenario version of the North Sea/1970-1999 model by increasing the activity rate of gear 1 (pelagic trawls and seines) by a factor of 2 model$data$fleet.model$gear_mult[1] <- 2 model$setup$model.ident <- "gear1x2" # Set a new identifier for the outputs results <- StrathE2E(model,nyears=10) dev.new() # Open a new graphics window plot_full_length_timeseries(model, results)