Plots of inshore and offshore low trophic level state variables for the final year of a model run

Description

Multi-panel time series plots of inshore and offshore low trophic level state variables for the final year of a model run generated by the StrathE2E() function

Usage

plot_final_year_LTL_inshore_vs_offshore(model, results)

Arguments

Details

The function plots a multi-panel page of time series plots of daily values of nutrients, detritus, macrophyte, phytoplankton, zooplankton and suspension/deposit feeding benthos (including larvae) for the inshore zone and upper and lower levels of the offshore zone, for the final year of a model run.

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_final_year , plot_final_year_HTL_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 low trophic level outputs
plot_final_year_LTL_inshore_vs_offshore(model, results)

# Time series plot of state varaiables in a new window leaving any existing windows open
dev.new()
plot_final_year_LTL_inshore_vs_offshore(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_LTL_inshore_vs_offshore(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 low trophic level output
model <- read_model("North_Sea", "1970-1999",model.ident="Baseline")
results <- StrathE2E(model,nyears=10)
plot_final_year_LTL_inshore_vs_offshore(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_LTL_inshore_vs_offshore(model, results)