Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Finite-amplitude Rossby Wave Diagnostics #566

Open
wants to merge 192 commits into
base: main
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from 181 commits
Commits
Show all changes
192 commits
Select commit Hold shift + click to select a range
561306f
Well-formatted settings.jsonc without comments yet.
csyhuang Jul 31, 2023
75b0f05
start all files
csyhuang Aug 19, 2023
e9722bd
merged main
csyhuang Aug 19, 2023
6e1d37a
try running daily_avg
csyhuang Aug 19, 2023
3b25aa7
daily_avg
csyhuang Aug 19, 2023
2fce674
daily_avg
csyhuang Aug 19, 2023
acf7486
daily_avg
csyhuang Aug 19, 2023
096726b
daily_avg
csyhuang Aug 19, 2023
12d41e8
daily_avg
csyhuang Aug 19, 2023
d45dc46
daily_avg
csyhuang Aug 19, 2023
3d5c410
daily_avg
csyhuang Aug 19, 2023
21ce25e
daily_avg
csyhuang Aug 19, 2023
ccfd952
daily_avg
csyhuang Aug 19, 2023
4897e80
daily_avg
csyhuang Aug 19, 2023
a370a58
add src/conda/env_finite_amplitude_wave_diag.yml
csyhuang Aug 19, 2023
8741055
add back ==
csyhuang Aug 19, 2023
eaa2960
try other dependencies
csyhuang Aug 19, 2023
036b80b
try other dependencies
csyhuang Aug 19, 2023
109cd65
add dependencies
csyhuang Aug 20, 2023
8eb0780
add var_names
csyhuang Aug 20, 2023
61a619f
add var_names
csyhuang Aug 20, 2023
c23077c
add var_names
csyhuang Aug 20, 2023
0502b7e
extract 3 slices to local
csyhuang Aug 20, 2023
f21f126
add
csyhuang Aug 20, 2023
9b7a675
add missing line
csyhuang Aug 20, 2023
4ecc91b
add dask
csyhuang Aug 20, 2023
c709ad1
add missing line
csyhuang Aug 21, 2023
b9c7d74
add missing
csyhuang Aug 21, 2023
04260c4
add missing
csyhuang Aug 21, 2023
957cf00
tried on model data, found unphysical issues
csyhuang Sep 5, 2023
d153a68
Merge branch 'main' into finite_amplitude_wave_diag
csyhuang Sep 5, 2023
b65c391
add all
csyhuang Sep 6, 2023
b04b4bc
season then daily mean
csyhuang Sep 6, 2023
af9310a
open_dataset
csyhuang Sep 6, 2023
373b83d
open_dataset
csyhuang Sep 6, 2023
7e2fa6a
open_dataset
csyhuang Sep 6, 2023
48dfe21
open_dataset
csyhuang Sep 6, 2023
cf3dcce
open_dataset
csyhuang Sep 6, 2023
21e2571
open_dataset
csyhuang Sep 6, 2023
3efa87b
open_dataset
csyhuang Sep 6, 2023
fd4a051
split files
csyhuang Sep 6, 2023
57e6063
retry
csyhuang Sep 6, 2023
1b1c09e
retry
csyhuang Sep 6, 2023
b77f17c
retry
csyhuang Sep 6, 2023
5c7186f
retry
csyhuang Sep 6, 2023
9dc3a02
retry
csyhuang Sep 6, 2023
6fbda41
retry
csyhuang Sep 6, 2023
730716d
retry
csyhuang Sep 6, 2023
d402754
retry
csyhuang Sep 6, 2023
d92ea9d
retry
csyhuang Sep 6, 2023
eae22e1
retry
csyhuang Sep 6, 2023
63818e4
retry
csyhuang Sep 6, 2023
49cfbb6
retry
csyhuang Sep 6, 2023
463d439
test for one month
csyhuang Sep 6, 2023
fe03ea5
print out time
csyhuang Sep 6, 2023
ac56ec0
use datetime
csyhuang Sep 6, 2023
6c24d58
retry
csyhuang Sep 6, 2023
865b115
method="nearest"
csyhuang Sep 6, 2023
69aef56
run in background to time
csyhuang Sep 6, 2023
5cf400a
add new files
csyhuang Sep 6, 2023
cc48fcf
add new files
csyhuang Sep 6, 2023
bd76fd3
update all
csyhuang Sep 7, 2023
91740af
remove unnecessary
csyhuang Sep 7, 2023
afc0699
remove unnecessary
csyhuang Sep 7, 2023
9b55798
fix time str
csyhuang Sep 7, 2023
e7fe7a5
remove unnecessary
csyhuang Sep 7, 2023
c3a357f
remove unnecessary
csyhuang Sep 7, 2023
fd280c9
remove unnecessary
csyhuang Sep 7, 2023
dbdb713
remove unnecessary
csyhuang Sep 7, 2023
3c8898e
remove unnecessary
csyhuang Sep 7, 2023
72184ed
rerun era5_diag_original
csyhuang Sep 7, 2023
23b82c0
move print process time to
csyhuang Sep 7, 2023
0b837dd
otc_model_data_1tstep_extract.py
csyhuang Sep 7, 2023
a60d842
update
csyhuang Sep 7, 2023
fa30a82
update
csyhuang Sep 7, 2023
01b8b2a
update
csyhuang Sep 7, 2023
f864c33
update
csyhuang Sep 7, 2023
f17b08f
send_to_sandro.py
csyhuang Sep 7, 2023
09450c4
update
csyhuang Sep 7, 2023
d488191
update
csyhuang Sep 7, 2023
9fc4a35
update
csyhuang Sep 7, 2023
534a9a1
update
csyhuang Sep 7, 2023
5618389
update
csyhuang Sep 7, 2023
c4b8744
upload file
csyhuang Sep 8, 2023
65f4fe7
Save progress. Main driver: finite_amplitude_wave_diag_zonal_mean.py
csyhuang Sep 8, 2023
01f0b8b
extract 10 timeslices
csyhuang Sep 8, 2023
45a2dfe
extract again
csyhuang Sep 9, 2023
2ffc9f1
extract 10 timeslices
csyhuang Sep 9, 2023
51cc107
try running on OTC
csyhuang Sep 9, 2023
d24cea8
gridfill move to main program
csyhuang Sep 9, 2023
f033ccb
fix environ
csyhuang Sep 9, 2023
2ebfe9a
fix names
csyhuang Sep 9, 2023
84d8fc3
fasten gridfill
csyhuang Sep 9, 2023
6ce07ca
close files when not needed
csyhuang Sep 9, 2023
5511d92
change to new interface to test on otc
csyhuang Nov 4, 2023
6eed37e
save progress
csyhuang Nov 4, 2023
fd5d32f
Merge branch 'main' into finite_amplitude_wave_diag
csyhuang Nov 4, 2023
3f07190
change latitude analysis resolution to 1.5 deg
csyhuang Nov 4, 2023
58acbba
update src/default_tests.jsonc
csyhuang Nov 4, 2023
d376343
Merge branch 'finite_amplitude_wave_diag' of https://github.com/csyhu…
csyhuang Nov 4, 2023
ea4a32a
use env variable in default_tests.jsonc
csyhuang Nov 4, 2023
0441d6e
Print out analysis resolution
csyhuang Nov 4, 2023
e492a7e
sample just first timestamp of a day
csyhuang Nov 4, 2023
8f9736a
'first' with skipna=True is not yet implemented on dask arrays
csyhuang Nov 4, 2023
46287a2
use groupby instead of resample
csyhuang Nov 4, 2023
c16b2a1
do four seasons
csyhuang Nov 4, 2023
1475ecd
remove unused files
csyhuang Nov 4, 2023
a31d402
refactored
csyhuang Nov 5, 2023
a2ad8b3
refactored
csyhuang Nov 5, 2023
15954ec
fix
csyhuang Nov 5, 2023
d6d7247
Try again
csyhuang Nov 5, 2023
9fd76a0
close file
csyhuang Nov 5, 2023
608f3ff
close file
csyhuang Nov 5, 2023
244156a
close file
csyhuang Nov 5, 2023
6310c92
test on cluster
csyhuang Nov 5, 2023
42389db
try running on cluster
csyhuang Nov 5, 2023
9635fc4
masked from second lowest level
csyhuang Nov 5, 2023
ed5c685
add this to test on cluster
csyhuang Nov 5, 2023
d61c63a
Merge branch 'main' into finite_amplitude_wave_diag
csyhuang Dec 5, 2023
3b0d418
update package name
csyhuang Dec 5, 2023
e7095ad
Merge branch 'finite_amplitude_wave_diag' of https://github.com/csyhu…
csyhuang Dec 5, 2023
2803c7b
Merge branch 'main' into finite_amplitude_wave_diag
csyhuang Feb 11, 2024
f53d410
refactored plot scripts
csyhuang Feb 12, 2024
c92d87b
now it runs from end-to-end. Commits before constructing HTML
csyhuang Feb 19, 2024
2b5986d
save progress
csyhuang Mar 2, 2024
405e01b
save progress
csyhuang Mar 3, 2024
14522cb
save progress
csyhuang Mar 13, 2024
67ea924
try to run on OTC
csyhuang Mar 17, 2024
244d07a
add default run file
csyhuang Mar 17, 2024
be6651d
remove typo
csyhuang Mar 17, 2024
59f148d
fix typo
csyhuang Mar 17, 2024
6c26882
fix typo
csyhuang Mar 17, 2024
89563ec
fix typo
csyhuang Mar 17, 2024
7490bbe
freq to 1hr
csyhuang Mar 17, 2024
099e777
freq to 1hr
csyhuang Mar 17, 2024
8725d59
try again after moving files
csyhuang Mar 17, 2024
b9c26bb
try set dry run to true
csyhuang Mar 17, 2024
14c3956
retry
csyhuang Mar 17, 2024
99c977e
try another
csyhuang Mar 17, 2024
faa1ab4
change freq to 'day'
csyhuang Mar 17, 2024
84f1430
try day
csyhuang Mar 17, 2024
0b4558f
switch back to 1hr
csyhuang Mar 17, 2024
04fa13d
add html
csyhuang Mar 17, 2024
1228a11
remove __main__:
csyhuang Mar 17, 2024
3834587
save png
csyhuang Mar 17, 2024
3c5423f
save png in another path
csyhuang Mar 17, 2024
ca36733
temporarily remove clean up
csyhuang Mar 17, 2024
4abd433
turn off test_mode
csyhuang Mar 17, 2024
64533e6
save png first
csyhuang Mar 17, 2024
c5f69e2
Add env var on otc
csyhuang Mar 17, 2024
c4f4ed4
debug
csyhuang Mar 17, 2024
4a52a6b
add print
csyhuang Mar 17, 2024
240ce9e
fix data path
csyhuang Mar 17, 2024
10ffd2b
fix data path
csyhuang Mar 17, 2024
aa7418d
rerun
csyhuang Mar 17, 2024
8af7f7d
add print statements
csyhuang Mar 17, 2024
96b4cb6
add print statements
csyhuang Mar 17, 2024
990f0a6
use first instead of mean
csyhuang Mar 17, 2024
09db6a0
use first instead of mean
csyhuang Mar 17, 2024
acde5b8
skipna=False
csyhuang Mar 17, 2024
273622a
use GFDL-CM3_historical_r1i1p1
csyhuang Mar 17, 2024
285e699
change freq to day
csyhuang Mar 17, 2024
35cea07
change year to 2005
csyhuang Mar 17, 2024
460fe7e
change coord
csyhuang Mar 17, 2024
c6179ef
update env var
csyhuang Mar 17, 2024
a4203dc
change freq to day
csyhuang Mar 17, 2024
090c3c0
plev
csyhuang Mar 17, 2024
e000537
4 seasons
csyhuang Mar 17, 2024
62b7a0c
4 seasons
csyhuang Mar 17, 2024
9a1c893
fix title string
csyhuang Mar 17, 2024
094e71e
use pip in .yml file
csyhuang Mar 19, 2024
f187fa6
Merge remote-tracking branch 'upstream/main' into finite_amplitude_wa…
csyhuang Mar 28, 2024
24dc764
save md
csyhuang Mar 28, 2024
6dd1ee6
fix data path
csyhuang Mar 28, 2024
fcca20a
remove production of png
csyhuang Mar 28, 2024
936da11
clean up branch
csyhuang May 20, 2024
61c2ab3
Add back references
csyhuang May 20, 2024
3bc4373
Added author names
csyhuang May 20, 2024
96138eb
clean up paths
csyhuang May 20, 2024
1d019fe
remove redundant file
csyhuang May 20, 2024
483b9de
Merge branch 'main' into finite_amplitude_wave_diag_fcca20a
wrongkindofdoctor May 20, 2024
6985ff7
merge
csyhuang May 29, 2024
1bf6002
Revert "merge"
csyhuang May 29, 2024
625e606
Merge remote-tracking branch 'upstream/main'
csyhuang May 30, 2024
28cc5b4
Merge remote-tracking branch 'upstream/main'
csyhuang Aug 16, 2024
5029c80
Merge branch 'main' into finite_amplitude_wave_diag_fcca20a
csyhuang Aug 16, 2024
6bb4a88
update src/conda/env_finite_amplitude_wave_diag.yml
csyhuang Aug 16, 2024
9c7b016
update env_finite_amplitude_wave_diag.yml
csyhuang Aug 16, 2024
d8e34f6
Update env file to use PyPI
csyhuang Aug 18, 2024
6d1a94b
move environment file and runtime config file from /src to my POD dir…
csyhuang Aug 18, 2024
31095c0
Merge remote-tracking branch 'upstream/main'
csyhuang Nov 4, 2024
9628cd8
merge
csyhuang Nov 4, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
65 changes: 34 additions & 31 deletions README.md

Large diffs are not rendered by default.

Original file line number Diff line number Diff line change
@@ -0,0 +1,127 @@
.. This is a comment in RestructuredText format (two periods and a space).

.. Note that all "statements" and "paragraphs" need to be separated by a blank
line. This means the source code can be hard-wrapped to 80 columns for ease
of reading. Multi-line comments or commands like this need to be indented by
exactly three spaces.

.. Underline with '='s to set top-level heading:
https://docutils.sourceforge.io/docs/user/rst/quickref.html#section-structure

Finite Amplitude Rossby Wave Diagnostics Documentation
======================================================

.. rst-class:: center

Clare S. Y. Huang\ |^1|, Christopher Polster |^2| and Noboru Nakamura\ |^1|

.. rst-class:: center

|^1|\ The University of Chicago, Chicago, Illinois

|^2|\ Johannes Gutenberg-Universität Mainz, Germany

.. rst-class:: center

Last update: 03/12/2024

Description
-----------
(to be filled in )

.. Underline with '-'s to make a second-level heading.

Version & Contact info
----------------------

Here you should describe who contributed to the diagnostic, and who should be
contacted for further information:

- Version/revision information: version 1 (03/12/2024)
- PI (name, affiliation, email): Clare S. Y. Huang (The University of Chicago, [email protected])
- Developer/point of contact: Clare S. Y. Huang (The University of Chicago, [email protected])
- Other contributors: Christopher Polster, Noboru Nakamura

.. Underline with '^'s to make a third-level heading.

Open source copyright agreement
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

The MDTF framework is distributed under the LGPLv3 license (see LICENSE.txt).
Unless you've distributed your script elsewhere, you don't need to change this.

Functionality
-------------

(to be filled in)

Required programming language and libraries
-------------------------------------------

(to be filled in)


Required model output variables
-------------------------------

(to be filled in)

References
----------

References
----------

.. _ref-Nakamura-annual-review:

10241. Nakamura, N. (2024). Large-Scale Eddy-Mean Flow Interaction in the Earth's Extratropical Atmosphere. *Annual Review of Fluid Mechanics*, **56**, 349-377,
`doi:10.1146/annurev-fluid-121021-035602 <https://doi.org/10.1146/annurev-fluid-121021-035602>`__.

.. _ref-Neal-et-al-GRL:

10242. Neal, E., Huang, C. S., & Nakamura, N. (2022). The 2021 Pacific Northwest heat wave and associated blocking: meteorology and the role of an upstream cyclone as a diabatic source of wave activity. *Geophysical Research Letters*, **49(8)**, e2021GL097699. `doi:10.1029/2021GL097699 <https://doi.org/10.1029/2021GL097699>`__.

.. _ref-Nakamura-Science:

10243. Nakamura, N., & Huang, C. S. (2018). Atmospheric blocking as a traffic jam in the jet stream. *Science*, **361(6397)**, 42-47, `doi:10.1126/science.aat0721 <https://doi.org/10.1126/science.aat0721>`__.

.. _ref-Nakamura-Solomon-JAS-2010:

10244. Nakamura, N., & Solomon, A. (2010). Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part I: Quasigeostrophic theory and analysis. *Journal of the atmospheric sciences*, **67(12)**, 3967-3983, `doi:10.1175/2010JAS3503.1 <https://doi.org/10.1175/2010JAS3503.1>`__.

.. _ref-Nakamura-Solomon-JAS-2011:

10245. Nakamura, N., & Solomon, A. (2011). Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part II: Analysis in the isentropic coordinate. Journal of the atmospheric sciences, 68(11), 2783-2799, `doi:10.1175/2011JAS3685.1 <https://doi.org/10.1175/2011JAS3685.1>`__.

.. _ref-Huang-Nakamura-JAS-2016:

10246. Huang, C. S., & Nakamura, N. (2016). Local finite-amplitude wave activity as a diagnostic of anomalous weather events. Journal of the Atmospheric Sciences, 73(1), 211-229, `doi:10.1175/JAS-D-15-0194.1 <https://doi.org/10.1175/JAS-D-15-0194.1>`__.

.. _ref-Huang-Nakamura-GRL-2017:

10247. Huang, C. S., & Nakamura, N. (2017). Local wave activity budgets of the wintertime Northern Hemisphere: Implication for the Pacific and Atlantic storm tracks. Geophysical Research Letters, 44(11), 5673-5682, `doi:10.1002/2017GL073760 <https://doi.org/10.1002/2017GL073760>`__.

More about this diagnostic
--------------------------

(to be filled in)

Links to external sites
^^^^^^^^^^^^^^^^^^^^^^^

(to be filled in)

More references and citations
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

(to be filled in)

Figures
^^^^^^^

(to be filled in)

Equations
^^^^^^^^^

(to be filled in)
Original file line number Diff line number Diff line change
@@ -0,0 +1,120 @@
<!-- This file is part of the EOF_500hPa module of the MDTF code package (see LICENSE.txt) -->

<TITLE> Finite-amplitude Rossby wave diagnostics</TITLE>
<img src="../mdtf_diag_banner.png">
<h3>Finite-amplitude Rossby wave diagnostics</h3>

<td class="d-block comment-body markdown-body js-comment-body">
<p dir="auto">The goal to incorporate Finite-amplitude wave activity (FAWA) formalism into MDTF repo is to provide diagnostic tools for model comparison. Each diagnostic shall associate with an interpretation related to the physical process.</p>
<p dir="auto">Below are some proposed candidates and associated physical interpretations:</p>
<table role="table">
<thead>
<tr>
<th>Diagnostic</th>
<th>Interpretation</th>
</tr>
</thead>
<tbody>
<tr>
<td>Seasonal climatology of zonal mean FAWA</td>
<td>It quantifies the strength of vertically propagating planetary waves, and also relative strength of synoptic eddies in troposphere</td>
</tr>
<tr>
<td>Seasonal climatology of Uref</td>
<td>To quantify strength of mean-flow interaction, i.e., adiabatic adjustment of flow from an eddy-free reference state is given by Uref-Ubar(zonal mean wind)</td>
</tr>
<tr>
<td>Seasonal climatology of &lt;LWA&gt; (&lt;...&gt; = vertically averaged)</td>
<td>To quantify geographical distribution (and amplitude) of eddies</td>
</tr>
<tr>
<td>Seasonal climatology of temporal Covariance of &lt;LWA&gt; and &lt;U&gt; (as in NH18 Fig. 2A)</td>
<td>A measure of strength of nonlinear wave-mean flow interaction (via nonlinear zonal wave activity flux)</td>
</tr>
</tbody>
</table>
</td>

<a href="https://mdtf-diagnostics.readthedocs.io/en/latest/sphinx_pods/finite_amplitude_wave_diag_zonal_mean.html">Full Documentation and Contact Information</a>

<h4>Zonal-mean finite-amplitude wave diagnostics for {{CASENAME}} (Climatologies)</h4>

<table class="tg">
<thead>
<tr>
<th>Name of Diagnostic</th>
<th>DJF</th>
<th>MAM</th>
<th>JJA</th>
<th>SON</th>
</tr>
</thead>
<tbody>
<tr>
<td>zonal mean zonal wind</td>
<td><a href=model/DJF_zonal_mean_u.png>plot</a></td>
<td><a href=model/MAM_zonal_mean_u.png>plot</a></td>
<td><a href=model/JJA_zonal_mean_u.png>plot</a></td>
<td><a href=model/SON_zonal_mean_u.png>plot</a></td>
</tr>
<tr>
<td>zonal mean wave activity (FAWA)</td>
<td><a href=model/DJF_zonal_mean_lwa.png>plot</a></td>
<td><a href=model/MAM_zonal_mean_lwa.png>plot</a></td>
<td><a href=model/JJA_zonal_mean_lwa.png>plot</a></td>
<td><a href=model/SON_zonal_mean_lwa.png>plot</a></td>
</tr>
<tr>
<td>zonal mean reference states (Uref)</td>
<td><a href=model/DJF_zonal_mean_uref.png>plot</a></td>
<td><a href=model/MAM_zonal_mean_uref.png>plot</a></td>
<td><a href=model/JJA_zonal_mean_uref.png>plot</a></td>
<td><a href=model/SON_zonal_mean_uref.png>plot</a></td>
</tr>
<tr>
<td>zonal mean wind adjustment (\Delta U)</td>
<td><a href=model/DJF_zonal_mean_delta_u.png>plot</a></td>
<td><a href=model/MAM_zonal_mean_delta_u.png>plot</a></td>
<td><a href=model/JJA_zonal_mean_delta_u.png>plot</a></td>
<td><a href=model/SON_zonal_mean_delta_u.png>plot</a></td>
</tr>
</tbody>
</table>


<h4>Vertically-averaged finite-amplitude wave diagnostics for {{CASENAME}} (Climatologies)</h4>

<table class="tg">
<thead>
<tr>
<th>Name of Diagnostic</th>
<th>DJF</th>
<th>MAM</th>
<th>JJA</th>
<th>SON</th>
</tr>
</thead>
<tbody>
<tr>
<td>Barotropic zonal mean zonal wind</td>
<td><a href=model/DJF_u_baro.png>plot</a></td>
<td><a href=model/MAM_u_baro.png>plot</a></td>
<td><a href=model/JJA_u_baro.png>plot</a></td>
<td><a href=model/SON_u_baro.png>plot</a></td>
</tr>
<tr>
<td>Barotropic local wave activity</td>
<td><a href=model/DJF_lwa_baro.png>plot</a></td>
<td><a href=model/MAM_lwa_baro.png>plot</a></td>
<td><a href=model/JJA_lwa_baro.png>plot</a></td>
<td><a href=model/SON_lwa_baro.png>plot</a></td>
</tr>
<tr>
<td>Covariance between barotropic zonal wind and LWA</td>
<td><a href=model/DJF_u_lwa_covariance.png>plot</a></td>
<td><a href=model/MAM_u_lwa_covariance.png>plot</a></td>
<td><a href=model/JJA_u_lwa_covariance.png>plot</a></td>
<td><a href=model/SON_u_lwa_covariance.png>plot</a></td>
</tr>
</tbody>
</table>
Loading
Loading