Repository under construction!

Codes used for the paper: Gomez-Navarro et al. (2026)

Reference: Gómez-Navarro, L.; van Sebille, E.; Morales-Márquez, V.; Hernández-Carrasco, I.; Albert, A.; Ubelmann, C.; Molines, J.M.; Le Sommer, J. and Brodeau, L. Impact of Tidal Forcing on Surface Particle Transport Properties: Insights From Twin Ocean Simulations. Journal of Advances in Modeling Earth Systems 2026, 10, 599. (edit) www.doi.org/10.1029/2024MS004805

venv_??

Use? https://github.com/Parcels-code/Azores_TidalForcing/blob/main/Code/Plotting/plots_source_code.md

Figures:

• fig_01.ipynb

  • Output:
    Figure 1: Spatial domain of the eNATL60 simulation (except the Gulf of Mexico, Black Sea and eastern Mediterranean Sea domains). Black box shows the region of this study where virtual surface particles are released. Red box shows the subregion used for some of the analyses.

• fig_02.ipynb

  • In:
    

    • xx.nc Obtained from: xx.py

  • Out:
    Figure 02: Box plots of mean cumulative distance [km] (top) and mean absolute distance [km] (bottom) travelled by the virtual particles per month.

  • Pendings:
    

    • thicker lines? to see better bottom part????

• fig_03_04.ipynb

  • In:

    • .nc

  • Out:
    

    • Figure 03: Gaussian Kernel Density Estimation (GKDE) comparison between non-tidal (top) and tidal (bottom) simulations. Maximum GKDE value (top) and percentage of particles with a high GKDE value (greater than 0.008) (bottom) are shown in the top left textbox.
    • Figure 04: Comparison of the percentage of particles with a high GKDE value (≥0.008) per month. Non-tidal results are shown in blue and tidal in red.

  • Pendings:
    

    • add dot in titles?? Apr. 2010
    • equation needed somewhere? (fig 4 caption??)
    • fig. 4 : Date [dd/mm/YYYY] , mm or MM, dd or DD?

• fig_05

  • In: List of input files

  • Out: Figure 05: Probability density functions of the surface particle accumulation (from the 2D, 0.1◦ binned histograms) after 28 days of advection. Results of the non-tidal simulation are shown in blue and from the tidal in red. Vertical, dashed lines indicate the maximum value.

• fig_06.ipynb

  • In:

    • *.nc

  • Out:
    Figure 06: Skewness temporal evolution in time for each month for non-tidal (blue) and tidal (red) simulations. Values in the text box are the variance of the skewness from the beginning of the month until the vertical line. The vertical line indicates the moment in time when any of the particles released at the boundary of the region first enters the subregion (see Section 2.3.2.2). Bottom plot shows the skewness value at the vertical line.

  • Pendings:
    

    • add dot in titles?? Apr. 2010
    • colour vertical line

• fig_07_08_09.ipynb

  • In: *

  • Out:

    • Figure 07: Comparison of attracting LCS structures (bFTLEs) on day 1 of each month from July to December 2009, for the non-tidal (left) and the tidal simulation (right).
    • Figure 08: Comparison of attracting LCS structures (bFTLEs) on day 1 of each month from January to June 2010, for the non-tidal (left) and the tidal simulation (right).
    • Figure 09: Top: Percentage of virtual particles with backward FTLE >0.5 days−1. Bottom: Skewness values of the backward FTLE fields.

  • Pendings:

    • Date [dd/MM/YYYY] , MM and dd or DD?

• fig_10.ipynb

  • In: List of input files

  • Out: Figure 10: Percentage difference with tidal forcing per month for each diagnostic calculated. From top to bottom: cumulative distance (CD), absolute distance (AD), percentage of particles with high Gaussian Kernel Density Estimation (GKDE ≥ 0.008) and percentage of particles with high backward Finite Time Lyapunov Exponents (bFTLE ≥ 0.5 days−1).

    • Pendings:
      • Date [dd/MM/YYYY] , MM and dd or DD?

Appendix figures:

• fig_A1.ipynb

  • In:

    • wnoises_loop_try2_cut.mat

  • Out:
    Figure A1: Probability density function of the backward FTLE fields on day 1 of each month for the no tidal forcing (blue) and tidal forcing (red) simulations. The vertical black line shows the threshold used at 0.5 days−1.

Supplementary figures:

• fig_S1_S2.ipynb

  • In:

    • wnoises_loop_try2_cut.mat
    • wnoises_loop_try2.mat

  • Out:
    

    • Figure S1: Gaussian Kernel Density Estimation (GKDE) comparison between non-tidal (top) and tidal (bottom) simulations. Maximum GKDE value (top) and percentage of particles with a high GKDE value (greater than 0.008) (bottom) are shown in the top left textbox.
    • Figure S2: Comparison of the percentage of particles with a high GKDE value (≥0.008) per month. Non-tidal results are shown in blue and tidal in red.

  • pENDING:

    • UPDATE AS 7 and 8??
    • cbar label
    • . in title Apr.

Calculations: