Multi-panel time series plots of state variable data for the final year of a model run generated by the StrathE2E() function

Description

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.

Usage

plot_final_year(model, results)

Arguments

Details

Units of the plotted varaiables are mass concentrations, scaled to seas surface area or layer volume, ie. mMN/m2 or mMN/m3

Value

Graphical display in a new graphics window

See Also

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

Examples

# 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)