transect exceedance module

The transect exceedance plotting module is also a variant of the transect plotting module, and able to perform comparison of model predicted water quality against environmental guidelines or thresholds, such as comparing the predicted water quality concentrations to the ANZECC guidelines.

Summary of transect exceedance module configuration

Field name	Description	Options	Comments
start_plot_ID	select which variable to start plotting		refer to the setting in master.varname
end_plot_ID	select which variable to finish plotting		refer to the setting in master.varname
plotvalidation	option to add field data	1: add field data; 0: no field data	only model data to be plotted in this module, need to be 0
plotmodel	option to add model results	1: add model results; 0: no model results	only model data to be plotted in this module, need to be 1
isHTML	option to add all plots into one HTML page	1: add to HTML; 0: no HTML	0 by default
polygon_file	define the line GIS file for transect
outputdirectory	define directory to save plots
htmloutput	define directory to save HTML files
istitled	option to add title	1: add; 0: not add	1 by default
isylabel	option to add y label	1: add; 0: not add	1 by default
islegend	option to add legend	1: add; 0: not add	1 by default
isYlim	option to define Y axis limits	1: add; 0: not add	0 by default
isGridon	option to add grid on	1: add; 0: not add	1 by default
isSurf	option to choose surface or bottom layer	1: surface layer; 0: bottom layer	1 by default
isSpherical	option to define the coordinate type	1: spherical; 0: Euler	0 by default
thresh.value	threshold values to be added onto plots		user can define multiple thresholds
thresh.legend	legends for thresholds		string array match the value size
boxlegend	define field data box plot legend location		'southeast' by default
rangelegend	define model result legend location		'northwest' by default
dimensions	define figure dimensions in centimeters	[20 10] by default
pdates.value	time period for transect plots		can define multiple periods
binfielddata	option to include nearby field data	1: include; 0: not include	1 by default
binradius	define the radius of bins for field data search	number with unit of km	0.5 by default
linedist	distance from model polyline to be considered	number with unit of m	500 by default
xlim	define limits in X axis
xticks	define ticks in X axis
xlabel	define label on X axis
cAxis.value	define limits of Y axis	empty [] by default, matlab will automatically adjust the y limit
filetype	define figure file type	'png' or 'eps'	'png' to save figures to PNG format only; 'eps' to save figures in both EPS and 300dpi JPG formats

Configuration Example

transectExc.start_plot_ID = 13;
transectExc.end_plot_ID = 13;

transectExc.polygon_file = 'E:\database\AED-MARVl-v0.2\Examples\Cockburn\GIS\Curtain_polyline_100m_QC.shp';
% Add field data to figure
transectExc.plotvalidation = 0; % 1 or 0
transectExc.pred_lims = [0.05,0.25,0.5,0.75,0.95];

transectExc.isRange = 1;
transectExc.istitled = 1;
transectExc.isylabel = 1;
transectExc.islegend = 0;
transectExc.isYlim = 1;
transectExc.isHTML = 1;
transectExc.isSurf = 1; %plot surface (1) or bottom (0)
transectExc.isSpherical = 0;
%transectExc.use_matfiles = 0;
%transectExc.add_obs_num = 1;
%config.boxon = 1;

transectExc.thresh(13).value = [0.5 1];
transectExc.thresh(13).legend = {'%time > 0.5NTU',...
    '%time > 1NTU'};
% ___
transectExc.outputdirectory = 'plotting/transect_exceedance/RAW/';
transectExc.htmloutput = 'plotting/transect_exceedance/HTML/';

% plotting configuration
transectExc.dimc = [0.9 0.9 0.9]; % dimmest (lightest) color
transectExc.boxlegend = 'southeast';
transectExc.rangelegend = 'northwest';
transectExc.dimensions = [20 10]; % Width & Height in cm

i=1;
transectExc.pdates(1).value = [datenum(2021,06,16) datenum(2021,07,01)];i=i+1;
transectExc.pdates(2).value = [datenum(2021,07,01) datenum(2021,07,15)];i=i+1;
transectExc.pdates(3).value = [datenum(2021,07,16) datenum(2021,08,01)];i=i+1;

transectExc.binfielddata = 1;
% radius distance to include field data. Used to bin data where number of
% sites is higher, but the frequency of sampling is low. The specified
% value will also make where on the line each polygon will be created. So
% if radius == 5, then there will be a search polygon found at r*2, so 0km, 10km, 20km etc. In windy rivers these polygons may overlap.
transectExc.binradius = 0.5;% in km;


%distance from model polyline to be consided.
%Field data further than specified distance won't be included.
%Even if found with search radius. This is to attempt to exclude data
%sampled outside of the domain.
transectExc.linedist = 1500;%  in m

transectExc.xlim = [0 45];% xlim in KM
transectExc.xticks = [0:5:45];
transectExc.xlabel = 'Distance from Southern CS (km)';

% ylim
for vvvv=1:size(MARVLs.master.varname,1)
    transectExc.cAxis(vvvv).value = [0 100];
end

transectExc.ncfile(1).symbol = {'-'};
transectExc.ncfile(1).translate = 1;
transectExc.ncfile(1).colour = [166,86,40]./255;% Surface and Bottom
transectExc.ncfile(1).edge_color = [166,86,40]./255;
transectExc.ncfile(1).col_pal_color =[[176 190 197]./255;[162 190 197]./255;[150 190 197]./255;[150 190 197]./255];

Example outputs

Example transect exceedance plot for suspended solids concentration against guidelines in Gladstone.