nvogtvincent/WIO_Marine_Debris

Repository for scripts used to reproduce figures in "Sources of marine debris for Seychelles and other remote islands in the western Indian Ocean" # WIO_Marine_Debris Repository for scripts used to reproduce figures in "Sources of marine debris for Seychelles and other remote islands in the western Indian Ocean" If you use our data and/or scripts for a publication or other application, it would be greatly appreciated if you could let us know! ## Directory overview ``` ANALYSIS/ Scripts to reproduce figures FIGURES/ Figures are saved here FISHERIES/ Fisheries input scripts and postprocessed data GRID_DATA/ Contains analysis grid (with source/sink sites) MATRICES/ Postprocessed terrestrial and marine matrices PLASTIC_DATA/ List of input and sink countries used in analysis REFERENCE/ Reference datasets for plotting/analyses SHIPPING/ Reference shipping density file used for Figure 5 SIM/ Scripts to rerun particle-tracking analyses (relies on ocean/wind/wave data as described in manuscript) TRAJ/ Postprocessing scripts to generate matrices from model output ``` ## User guide Firstly, set up a conda environment with the `environment.yml` file in the root directory. Alternatively, the packages needed for these analyses are numpy, matplotlib, xarray, xskillscore, pandas, cmasher, scipy, scikit-image, gdal, cartopy, resterio, geopandas, and tqdm. ### I want to generate results using a new combination of beaching and sinking rate, and windage 1. Retrieve the trajectory files for the physical scenario of interest from the BODC repository, and put them in in the relevant subdirectory in `TRAJ`. For instance, for the terrestrial CS3 scenario, you should put all 1056 terrestrial CS3 netcdf files in `TRAJ/TER/CS3/`; for the marine C0 scenario, you should put all 960 marine C0 netcdf files in `TRAJ/MAR/C0`, etc. 2. Run the relevant script to generate a processed matrix (note that this may take several hours): 1. **Terrestrial**: Run `land_process.py [us_d] [ub_d] [scenario]` where `[us_d]` is the sinking timescale in days, `[us_b]` is the …