## This file is written and formatted in notepad in the normal text language. 
## Please use this program to view this file.

## Author - 	Simon sharp, Rose Johnson, Edwin van der sar
## Project - 	How debris-flow composition affects bed erosion quantity and mechanisms - an experimental assessment
## Date - 	2 December 2021

------------------------------------------------------------------------------------------------------------ 

General
=======
This readme.txt file briefly describes the contents of the datapackage required to process, 
plot and analyse data of the experiments in the debris flow flume. The experiments simulate the erosion of an unconsolidated 
loosely packed bed by debris flows of different compositions and volumes. 

A few experiments are not used in the paper, either due to slumping of the bed or failing of the data collection, these 
experiments are noted down in the supplementary table added to the paper. 

Funded by the Netherlands Organisation for Scientific Research (NWO).


Folder Structure
================
Two folders are provided: one with the raw and processed DEMs before DF erosion (T0) and after erosion (T1)  and one with 
the raw sensor data.

- DataSupplement
	-- DEMs
	   -- exp0..
	      -- T0
		-- Processed
	        -- Raw	
	      -- T1
		-- Processed
	        -- Raw	
	-- SensorData
	-- readmeYoda
	-- Supplement_Table1_ExperimentalInformation




Folder contents
===============
Raw DEMs are given in .pct format, processed DEMs are in .mat format.

Raw sensor data is given in .csv format. In the .csv files the order of the stored data is as follows: 
Sample Number, Date/Time, OADM_01 (laser) , OADM_02 (laser), PORE_01, PORE_02, WEIGHT (load cell), SHEAR, TILT (of the flume), 
GEOPHONE_V, GEOPHONE_H, FADK_01 (laser), FADK_02 (laser), FADK_03 (laser).



Software version required
=========================
Matlab 2020 or later.



Abbreviations
=============
DEM: Digital Elevation Model



Data processing and analyses 
=========
For calculations related to velocity (frontal flow velocity, momentum and dimensionless numbers), the data from the laser 
OADM_02 (at  290cm, right before the erodible bed). 
In the paper the data from the other laser distance sensors is not used (OADM_01, FADK_01-03), as well as the data from the 
load cell and the shear sensor as they proved to be unreliable. 

The following calibrations are used for the laser distance sensors
FADK_01 = 0.03457 * Raw [V] + 0.21183 [m]		[at 138.5 cm]
FADK_02 = 0.03564 * Raw [V] + 0.25961 [m]	    	[at 281.4 cm]
OADM_01 = -0.05037 * Raw [V] + 0.19897 [m]    		[at 298.1 cm]
OADM_02 = -0.05005 * Raw [V] + 0.19686 [m]    		[at 290 cm]
FADK_03 = 0.03483 * Raw [V] + 0.25715 [m]	    	[at 535.7 cm]

For the geophone data the following procedure is used to calibrate and process the raw data
1. Apply band-stop filter - MATLAB code: Gbs = filter1('bs',G,'fs',9500,'fc',[2.5 50],'order',1; - [G=raw (vertical/horizontal) geophone signal]
2. Convert data to amplitude (make data absolute) -> GVa or GHa - [GVa=vertical geophone amplitude, GHa=horizontal geophone amplitude]
3. Convert data to (m/s) - GVa=GVa/20
4. Calculate total seismic energy (m/s)^2 by numerical integration - MATALB code: Energy(i,2) = trapz(GVa(1:end-1,1), GVa(1:end-1,2).^2);


