Merge PR #341 (Fix missing headers in aerosol table)

This merge brings PR #341 (Fix missing headers in aerosol table, by @yuanjianz) into the GCPy 1.6.0 development stream. This PR makes sure that the line containing the header for the benchmark aerosol burden table is written to file. Signed-off-by: Bob Yantosca <[email protected]>
geoschem · Sep 13, 2024 · 5958edf · 5958edf
2 parents cfe6b55 + 05a394d
commit 5958edf
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -99,8 +99,9 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 - Import error in `gcpy/examples/diagnostics/compare_diags.py`
 - Added missing `n_cores` to `gcpy/examples/diagnostics/compare_diags.yml`
 - Added missing `plot_drydep` option to `gcpy/gcpy/benchmark/config/1yr_ch4_benchmark.yml`
-- Add `docs/requirements.txt` symbolic link to `docs/environment_files/read_the_docs_requirements.txt` for RTD builds
+- Added `docs/requirements.txt` symbolic link to `docs/environment_files/read_the_docs_requirements.txt` for RTD builds
 - `gcpy/file_regrid.py` now tests if `lon_bnds`, `lat_bnds` are in the dataset before trying to drop them
+- Ensured that the header line in the aerosol burden table is written to file
 
 ### Removed
 - Example script `gcpy/examples/plotting/mda8_o3_timeseries.py`

diff --git a/gcpy/benchmark/modules/benchmark_funcs.py b/gcpy/benchmark/modules/benchmark_funcs.py
@@ -4175,7 +4175,7 @@ def print_aerosol_metrics(data, species_list, filename, title, label):
 
             # Print top header
             print("%" * 79, file=f)
-            print(f" {title} for {year} in {devstr}")
+            print(f" {title} for {year} in {devstr}", file=f)
             print(" (weighted by the number of days per month)", file=f)
             print("%" * 79, file=f)
             line = "\n" + " " * 40 + "Strat         Trop         Strat+Trop\n"