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

This repository contains the data used in the paper:<br> "Gadal C., Mercier M. and Lacaze L. Slumping regime in lock-release turbidity currents. Submitted to Journal of FLuid Mechanics.",<br> where the slumping regime of turbidity currents is studied with respect to the initial volume fraction,<br> the bottom slope and the particle settling velocity. The folder 'runs' contains 169 netcdf4 files corresponding to each experimental run used in the paper.<br> For each run, the structure of the netcdf file is the following:<br> - attributes:<br> - particle_type: particle type used (silica sand, glass beads or saline water)<br> - run_number: netcdf file name<br> - expe_type: always lock-release here<br> - surface_type: can be 'open surface' or 'rigid lid'<br> - set_up: can be 'set-up 1' or 'set-up 2'<br> - run_oldID: run name corresponding to the experimental notebook - groups:<br> - initial_parameters:<br> - dimensions(sizes):<br> - variables(dimensions):<br> - Bottom slope(): bottom slope<br> - Current density(): initial average (fluid + particle) lock density<br> - Grain density(): particle density (not measured, estimated)<br> - Grain diameter(): particle diameter<br> - Initial Reynolds number(): initial Reynolds number, [rho_0 * u_0 * h_0 / mu]<br> - Initial Rouse number(): initial Rouse number, [v_s / u_0]<br> - Initial volume fraction(): initial lock particle volume fraction<br> - Reduced gravity(): reduced gravity, [g*(rho_0 - rho_f)/rho_f]<br> - Settling velocity(): particle settling velocity<br> - Temperature(): water temperature (not measured)<br> - V0 (lock volume)(): lock suspension volume<br> - Water density(): water density<br> - Water dynamic viscosity(): water dynamic viscosity (not measured)<br> - h0 (lock height)(): suspension height inside lock<br> - u0 (velocity scale)(): velocity scale, [sqrt(g'*h_0)]<br> - w0 (tank width)(): lock crossstream width<br> - x0 (lock length)(): lock streamwise length - scalar_variables:<br> - dimensions(sizes): tuples(2), x(1181)<br> - variables(dimensions):<br> - Av. shape('x',): current average shape<br> - Av. shape head volume(): Volume per unit of width of the head part of current average shape<br> - Av. shape tail volume(): Volume per unit of width of the tail part of current average shape<br> - Av. shape volume(): Volume per unit of width of current average shape<br> - Current Froude number(): Current Froude number, [u_c/sqrt(g' * h_b)]<br> - Current Reynolds number(): Current Reynolds number, [rho_0 * u_c * h_b / mu]<br> - Current Rouse number(): Current Rouse number, [v_s / u_c]<br> - Current head height (log fit)(): current height h_h coming from log fit<br> - Current height (benjamin fit)(): current height h_b coming from fit of Benjamin's shape<br> - Current nose height (benjamin fit)(): current nose height h_n coming from fit of Benjamin's shape<br> - Current nose height (log fit)(): current nose h_n coming from log fit<br> - Geometrical Froude number(): Current Geometrical Froude number, [u_c/sqrt(g' * h0)]<br> - Geometrical Reynolds number(): Current Geometrical Reynolds number [rho_0 * u_c * h0 / mu]<br> - Times lock opening (tstart, tend)('tuples',): Start and end times of lock opening<br> - Times slumping regime (tstart, tend)('tuples',): Start and end times of constant velocity regime<br> - Velocity (slumping regime)(): Current velocity during slumping - time_series:<br> - dimensions(sizes): time(3071)<br> - variables(dimensions):<br> - Volume('time',): current volume per unit of width<br> - contour time series (x)('time',): x coordinate time series of the current contours<br> - contour time series (y)('time',): y coordinate time series of the current contours<br> - position('time',): front position<br> - time('time',): time vector<br> - velocity('time',): front velocity <br> Most variables possess the following attributes:<br> - unit: corresponding unit<br> - std: error(s) on the given quantity, calculated by error propagation from measurement uncertainties using the `uncertainties` module (https://pythonhosted.org/uncertainties/) in python.<br> - comments: comments on the given quantity (definition, formulaes, etc ..) Note that all variables related to the current shape are no available for experimental runs carried in set-up 2. The script ReadPlotData.py shows how to display the structure of a netcdf file,<br> and gives exemples on 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,<br> allowing for an easier access for testing purposes. *Note that errors are not given in this file.*