This repository contains scripts and data for reproducing the figures from the manuscript:
Fajardo-Urbina, J.M., Arts, G., Gräwe, U., Clercx, H.J.H., Gerkema, T., & Duran-Matute, M. (2023). Atmospherically Driven Seasonal and Interannual Variability in the Lagrangian Transport Time Scales of a Multiple-Inlet Coastal System. Journal of Geophysical Research: Oceans, 128, e2022JC019522. https://doi.org/10.1029/2022JC019522
There are in total 9 NetCDF files which can be used to reproduce Figure2 to Figure10 from the main manuscript. They are provided inside the folder data.
The environment employed for the analysis is based on Python v3.8 and can be found in the file environment.yml. The wavelet analysis is based on the Python package Pycwt v0.3.0a22, but we added a script to perform the bias correction (Liu et al., 2007) and a wavelet filter (Torrence & Compo, 1998). This script (called get_results.py) and the core functions from Pycwt are in /functions/pycwt_wavelet_package.
All the scripts (notebooks) necessary for reproducing the figures are in the folder notebooks. They can be run using any of the following instructions:
- Mybinder.org: click on the binder icon to open a jupyter-lab session.
- Google Colab: follow the instructions from the notebook clone_repo_using_google_colab.ipynb
- Clone or download the repository in your PC: install the packages of the file environment.yml.
Jup-lab:
The 9 netCDF files used to generate the figures were obtained from the following raw data:
- Eulerian data from the GETM/GOTM model, and its set-up is described in:
- Duran-Matute et al. (2014): https://doi.org/10.5194/os-10-611-2014
- Grawe et al. (2016): https://doi.org/10.1002/2016JC011655
- Lagrangian data from the model Parcels v2.1.1, which can be installed from:
- Monthly mean sea level pressure from NCEP-NCAR Reanalysis 1 (used for the large-scale analytical model):