Data used in 'Slumping regime in lock-release turbidity currents'

This repository contains the data used in the paper:

Gadal, C., Mercier, M., Rastello, M., & Lacaze, L. (2023). Slumping regime in lock-release turbidity currents. Journal of Fluid Mechanics, 974, A4. doi:10.1017/jfm.2023.762
where the slumping regime of turbidity currents is studied with respect to the initial volume fraction, the bottom slope and the particle settling velocity. The folder 'runs' contains 169 netcdf4 files corresponding to each experimental run used in the paper. For each run, the structure of the NetCDF file is the following:

attributes:

particle_type: particle type used (silica sand, glass beads or saline water)
run_number: NetCDF file name
expe_type: always lock-release here
surface_type: can be 'open surface' or 'rigid lid'
set_up: can be 'set-up 1' or 'set-up 2'
run_oldID: run name corresponding to the experimental notebook
groups:

initial_parameters:

dimensions(sizes):
variables(dimensions):

Bottom slope(): bottom slope
Current density(): initial average (fluid + particle) lock density
Grain density(): particle density (not measured, estimated)
Grain diameter(): particle diameter
Initial Reynolds number(): initial Reynolds number, [rho_0 * u_0 * h_0 / mu]
Initial Rouse number(): initial Rouse number, [v_s / u_0]
Initial volume fraction(): initial lock particle volume fraction
Reduced gravity(): reduced gravity, [g*(rho_0 - rho_f)/rho_f]
Settling velocity(): particle settling velocity
Temperature(): water temperature (not measured)
V0 (lock volume)(): lock suspension volume
Water density(): water density
Water dynamic viscosity(): water dynamic viscosity (not measured)
h0 (lock height)(): suspension height inside lock
u0 (velocity scale)():  velocity scale, [sqrt(g'*h_0)]
w0 (tank width)(): lock crosstream width
x0 (lock length)(): lock streamwise length
scalar_variables:

dimensions(sizes): tuples(2), x(1181)
variables(dimensions):

Av. shape('x',): current average shape
Av. shape head volume(): Volume per unit of width of the head part of current average shape
Av. shape tail volume(): Volume per unit of width of the tail part of current average shape
Av. shape volume(): Volume per unit of width of current average shape
Bulk entrainment coefficient(): Bulk entrainment coefficient during slumping
Current Froude number(): Current Froude number, [u_c/sqrt(g' * h_b)]
Current Reynolds number(): Current Reynolds number, [rho_0 * u_c * h_b / mu]
Current Rouse number(): Current Rouse number, [v_s / u_c]
Current head height (log fit)(): current height h_h coming from log fit
Current height (benjamin fit)(): current height h_b coming from fit of Benjamin's shape
Current nose height (benjamin fit)(): current nose height h_n coming from fit of Benjamin's shape
Current nose height (log fit)(): current nose h_n coming from log fit
Geometrical Froude number(): Current Geometrical Froude number, [u_c/sqrt(g' * h0)]
Geometrical Reynolds number(): Current Geometrical Reynolds number [rho_0 * u_c * h0 / mu]
Times lock opening (tstart, tend)('tuples',): Start and end times of lock opening
Times slumping regime (tstart, tend)('tuples',): Start and end times of constant velocity regime
Velocity (slumping regime)(): Current velocity during slumping
time_series:
dimensions(sizes): time(3071)
variables(dimensions):
Volume('time',): current volume per unit of width
contour time series (x)('time',): x coordinate time series of the current contours
contour time series (y)('time',): y coordinate time series of the current contours
position('time',): front position
time('time',): time vector
velocity('time',): front velocity
Most variables possess the following attributes:

unit: corresponding unit
std: error(s) on the given quantity, calculated by error propagation from measurement uncertainties using the `uncertainties` module (https://pythonhosted.org/uncertainties/) in Python.
comments: comments on the given quantity (definition, formulas, etc ..)
Note that all variables related to the current shape are not available for experimental runs carried out in set-up 2.The script ReadPlotData.py shows how to display the structure of a NetCDF file, and gives examples of how to load some variables and plot them by reproducing some of the paper's figures.The CSV file 'dataset_summary.csv' offers a summary of all runs and corresponding experimental parameters, allowing for easier access for testing purposes. *Note that errors are not given in this file.*
OpenData License: licence-ouverte-v2.0
If you use this open data in your work (research or other), please cite in your bibliography the following reference doi:https://doi.org/10.1017/jfm.2023.762