1. Process hydrophone data
Hydrophone data is processed using the script 'hydrophoneArray.m'.
This script reads, calibrates, detrends, and time aligns acoustic data from the two 4-channel recorders. 
A table is provided giving the lat/lon of each source position and the corresponding acoustic data (+/- 2.5 seconds) which was recorded during the time that the boat was at that source position. 
For each source position, acoustic data is averaged across the six hydrophones and TOL analysis performed. The output data is provided as 'ArrayData_11Aug2021_GoatIsland_Koura_+-2_5secs.mat'.

2. Plot hydrophone data per source position (Supplementary Figure S2)
using the script 'SupplementaryFigureS2.m'.

3. Plot audiogram and ambient acoustic data (Supplementary Figure S3)
using the script 'SupplementaryFigureS3.m'.
This script plots the hearing threshold data for each species and each sound modality.
This script also interpolates the hearing thresholds to full octave levels and saves these in the structure 'audiogram.mat', provided.
This structure is necessary for the LSR analysis in subsequent scripts.

3. Model the source level of the boat
This is done using 'modelSL.m' and the output of the hydrophone data processing, 'ArrayData_11Aug2021_GoatIsland_Koura_+-2_5secs.mat', which is organised into transects.
The calculated source levels are used in the remaining scripts to calculate the received levels at a range of different positions.
The source levels are stored in the provided structure 'outMedianSL.mat'.

4. Check cut-off frequency at source/receiver positions, and path between each source position and hydrophone. This requires ArcGIS, and the steps required are outlined in the following scripts:

-The script 'points2depths.m' is used to calculate cut-off frequency at each source position. It is then determined, at which source positions, the cut-off frequency is too high. These positions are then stored in the variable 'pos_2remove' and must be deleted once loading in OUT_F in line 208 of 'LSR_boatRange.m' (see below). 

-The script 'LSR_boatRange_directPath.m' is used to inspect the cut-off frequency on the path between each source position and the central hydrophone array.

-The script 'points2depths.m' is also used to calculate the cut-off frequency at each of the receiver positions surrounding the six randomly selected source positions (presented in Figure 4). The relevant code from this script is duplicated and applied in 'LSR_map.m' (see below). 

4. LSR analysis for all source positions and single listener position (Figures 2 & 3).
This is performed using the script LSR_boatRange.m.
First, PL data for each centre frequency/source position is extracted. This data can be provided, upon request, as 'OUT_F_3000mx72.mat'. This file could not be shared on Figshare due to its large size.
The cut-off frequency at each source position, and on the path between each source position and the hydrophone, is inspected. 
Received levels are then calculated, converted to acceleration, and then to full octave bands.
LSR for a single listener position is then calculated (although there are many listener positions to choose from).
This LSR for a single listener position is then plotted alongside the boat range at each source position.
LSR values for a single listener position are provided as 'LSR_Goatp50_3000mx72.mat'.

5. LSR analysis for six source positions and an array of listener positions (Figure 4).
This is performed using the script LSR_map.m.
As above, PL data for each centre frequency/source position is extracted. This data is provided as 'OUT_F_snapshot_3000mx72.mat'.
Received levels are then calculated and converted from third octave to full octave bands. 
These values are then converted to particle acceleration.
Pressure and particle motion LSR is then calculated for all source positions and the array of receiver positions for each species, and each ambient condition.
Any receiver position where the cut-off frequency is > 80 Hz is removed.
The output is stored in the structure 'LSR_Goatp50_3000mx72_sixPos.mat'.