kuchaale / grl_2020

Code for the paper submitted to GRL

Diverse dynamical response to orographic gravity wave drag hotspots - a zonal mean perspective

P. Sacha, A. Kuchar, R. Eichinger, P. Pisoft, Ch. Jacobi and H. Rieder

Published in Geophysical Research Letters.

Code used to process and visualise the model and other data outputs in order to reproduce figures in the manuscript. Model data are available here. All datasets already preprocessed can be found here.

Notebooks for each individual figure as well as for two data tables are in the code/ directory, while the figures themselves are in the plots/ directory.

Figures

#	Figure	Notebook	Dependencies
1	Decomposition into leading zonal wavenumbers of Eliassen-Palm flux composite average at lag=0 for the HI, EA and WA hotspot, respectively	GRL_reproduce_Fig1.ipynb	TEM_calculation-w-GetWaves.py, tem_calculation.py
2	Zonal wind composite at lag=0 representing the HI, EA and WA hotspot, respectively	GRL_reproduce_Fig2.ipynb	ssw_composite_cmam_optimized2-wss_hotspots.py, montecarlo_composites_script.py, resampling.py, climatology_CMAM_woSSW.py, calc_zmzw_tendency.py
3	Zonal mean FAWA composite averages at lag=0 representing the HI, EA and WA hotspot, respectively	GRL_reproduce_Fig3.ipynb	lwa_calc.py, lwa_tendency_calc.py

Supplementary figures

#	Figure	Notebook	Dependencies
S1	Residual anomalies during boreal winter representing the HI, EA and WA hotspot, respectively	GRL_reproduce_FigS1.ipynb
S2	{Surface downward eastward wind stress anomalies during boreal winter representing the HI, EA and WA hotspot, respectively	GRL_reproduce_FigS2.ipynb
S3	EPFD composite average at lag=0 for the HI, EA and WA hotspot, respectively	GRL_reproduce_FigS3.ipynb
S4	Zonal mean of zonal wind anomalies (shading; units: m/s) at lag=0 representing the Himalaya hotspot composites differentiated according to QBO phases	GRL_reproduce_FigS4.ipynb	ssw_composite_cmam_optimized2-wss_hotspots_QBO.py, GRL_QBO_timeseries4composites_CMAM.ipynb

Supplementary tables

#	Table	Notebook	Dependencies
S1	Number of detected peak events	GRL_QBO_timeseries4composites_CMAM.ipynb

Review figures

Notebook
GRL_OGWD_absolute_composite.ipynb
OGWD_correlation.ipynb

Required package installation

pip install -r requirements.txt

lwa_calc.py requires hn2016_falwa to be installed. TEM_calculation-w-GetWaves.py was adapted using aostools.

References

• Kuchar, A., Sacha, P., Eichinger, R., Jacobi, C., Pisoft, P., and Rieder, H. E.: On the intermittency of orographic gravity wave hotspots and its importance for middle atmosphere dynamics, Weather Clim. Dynam., 1, 481-495, https://doi.org/10.5194/wcd-2020-21, 2020.
• Sacha, P., Kuchar, A., Eichinger, R., Pisoft, P., Jacobi, C., and Rieder, H. E.: Diverse dynamical response to orographic gravity wave drag hotspots — a zonal mean perspective, Geophysical Research Letters, 48, e2021GL093305, https://doi.org/10.1029/2021GL093305, 2021.