Analysis and plotting scripts for first manuscript submission of "Temporal and spatial scale dependency of air-sea interactions via the vertical mixing mechanism"

Analysis and plotting scripts for figures in first manuscript submission of "Temporal and spatial scale dependency of air-sea interactions via the vertical mixing mechanism" by Dian A. Putrasahan and Jin-Song von Storch.

Summary: (25 analysis scripts + 6 function scripts)

0)  Collection of python functions pertaining to this analysis

• • edgevertexcell.py

  • extractICONdata.py

  • getSpectraPlot_r2b9.py

  • plotspectra_r2b9.py

  • r2b9ocnrgn.py

  • spatialspec_r2b9.py

1)  Extract daily and monthly fields (SST and windstress) on ocean grid

• • extractR2B9ocn_monthly.job

  • extractR2B9ocn_daily.job

2)  Compute the SST gradients and the wind stress divergence

• • create_daily_SSTgrad_TAUdivcurl_ngc2009.py

  • create_monthly_SSTgrad_TAUdivcurl_ngc2009.py

  • create_hourly_SSTgrad_TAUdivcurl_ocn_ngc2009.py

3a)  Spatially filter full fields (2deg) = high-pass filtered + smoothed fields 

• • sm_monthly_SSTgradTAUdivcurl_ngc2009.job 

  • submit_smmonthly_SSTgradTAUdivcurl_ngc2009.sh

  • sm_daily_SSTgradTAUdivcurl_ngc2009.job

  • submit_smdaily_SSTgradTAUdivcurl_ngc2009.sh

  • sm_hourly_SSTgradTAUdivcurl_ngc2009.job

  • submit_smhourly_SSTgradTAUdivcurl_ngc2009.sh

3b)  Compute global distribution of coupling coefficient and pattern correlation [Figure 1]

• • Fig1.ipynb

4)  Use full daily fields and interpolate onto regular grid 3600x1800

• • create_daily_0.05deg_downSST.py

  • create_daily_0.05deg_TAUdiv.py

5) For each ocean basin, compute FFT for each latitude.

• • FFT_monthly_0.05deg_downSST_TAUdiv.py

  • FFTthres_daily_0.05deg.py

  • execFFT_thres_005deg.job

  • submit_FFTrgn_thres_daily_005.sh

6)  Compute and plot spatial spectra for each ocean basin  [Figure 2]

• • Fig2.py

7) Compute and plot frequency-wavenumber spectra for N. Atlantic, tropical Pacific and Southern Ocean [Figure 3]

• • Fig3a.py

  • Fig3b.py

  • Fig3c.py

8)  Global distribution of coupling coefficient and pattern correlation based on 10deg smoothing [Figure S1]

• • FigS1.ipynb

9)  Spatial power spectra of SST and SST gradients [Figure S2]

• • FigS2.py

10)  Spatial coherence of downwind SST gradient and windstress divergence as a function of latitude based on daily data  [Figure S3]

• • FigS3.ipynb