plot_final_year_LTL_inshore_vs_offshore {StrathE2E2} | R Documentation |
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
plot_final_year_LTL_inshore_vs_offshore(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 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
Graphical display in a new graphics window
read_model
, StrathE2E
, plot_final_year
, plot_final_year_HTL_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 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)