plot_final_year {StrathE2E2} | R Documentation |
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 final year of a model run. Currently the masses of macrophytes, corpses and discards are not included in the plots due to space constraints.
plot_final_year(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 |
Units of the plotted varaiables are mass concentrations, scaled to seas surface area or layer volume, ie. mMN/m2 or mMN/m3
Graphical display in a new graphics window
read_model
, StrathE2E
, plot_full_length_timeseries
, 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 final year of outputs plot_final_year(model, results) # Time series plot of state varaiables in a new window leaving any existing windows open dev.new() plot_final_year(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_final_year(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 final year output model <- read_model("North_Sea", "1970-1999",model.ident="Baseline") results <- StrathE2E(model,nyears=10) plot_final_year(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_final_year(model, results)