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

Optional mass accumulation table in 1-month benchmark #216

Merged
merged 2 commits into from
May 23, 2023

Conversation

lizziel
Copy link
Contributor

@lizziel lizziel commented May 12, 2023

Name: Lizzie Lundgren
Institution: Harvard University

Description:

The PR adds a mass accumulation table (full column only) as a 1-month benchmark option. The mass accumulation table compares the difference in mass change for each species across two different runs, ref and dev. The mass change per species is computed as the difference between start and end based on values in the restart files (end restart mass minus start restart mass). Concentrations are converted to mass in the same way as done for the benchmark mass tables.

Here is an example of what the accumulation table looks like.

#########################################################################################
### Global mass accumulation (Gg) at end of simulation (Trop + Strat)                 ###
###                                                                                   ###
### Computed as change in instantaneous mass across period                            ###
###                                                                                   ###
### Ref = GCC 14.2.0                                                                  ###
### Dev = GCHP 14.2.0                                                                 ###
###                                                                                   ###
### Ref period: 2019-07-01T00:00:00 - 2019-07-01T01:00:00                             ###
### Dev period: 2019-07-01T00:00:00 - 2019-07-01T01:00:00                             ###
###                                                                                   ###
### Dev and Ref show 320 differences                                                  ###
#########################################################################################
                                    Ref                 Dev     Dev - Ref    % diff diffs
A3O2               :           0.010931            0.011537      0.000606     5.545   *
ACET               :          -0.826196           -1.264580     -0.438384    53.061   *
ACTA               :           2.208566            2.235849      0.027283     1.235   *
AERI               :           0.007904            0.008448      0.000544     6.880   *
ALD2               :           2.804017            2.823217      0.019200     0.685   *
ALK4               :           1.063624            0.743184     -0.320439   -30.127   *
AONITA             :          -0.027543           -0.015672      0.011872   -43.103   *
AROMP4             :           0.052310            0.056877      0.004567     8.731   *
AROMP5             :           0.048970            0.053550      0.004579     9.351   *
AROMRO2            :           0.008096            0.009819      0.001723    21.287   *

This table can be used with the operations budget table to look at changes in mass across a run, and to validate the operations budget diagnostics. Its values correspond to the ACCUMULATION entries in the operations budget table and is an alternative (and cheaper) way to get those values.

The table is currently only produced for full column and is limited to the 1-month benchmark. It can be expanded in a future PR if there is demand.

Expected differences:

None. This table is off by default in the 1-month benchmark yaml configuration file.

… option

The mass accumulation table compares the difference in mass change for each
species across two different runs, ref and dev. The mass change per species
is computed as the difference between start and end based on values
in the restart files (end restart mass minus start restart mass).
Concentrations are converted to mass in the same way as done for the benchmark
mass tables.

This table can be used with the operations budget table to look at changes
in mass across a run, and to validate the operations budget diagnostics. Its
values correspond to the ACCUMULATION entries in the operations budget table
and is an alternative (and cheaper) way to get those values.

The table is currently only produced for full column and is limited to the
1-month benchmark.

Signed-off-by: Lizzie Lundgren <elundgren@seas.harvard.edu>

initial commit

Signed-off-by: Lizzie Lundgren <elundgren@seas.harvard.edu>
…hmark

This table is optional and is off by default.

Signed-off-by: Lizzie Lundgren <elundgren@seas.harvard.edu>
@lizziel lizziel force-pushed the feature/add_mass_accumulation_table_to_1mo_benchmark branch from d944f43 to 51b4e8b Compare May 16, 2023 15:30
@lizziel lizziel marked this pull request as ready for review May 16, 2023 15:30
@lizziel lizziel requested a review from yantosca May 16, 2023 15:30
@yantosca yantosca changed the base branch from main to dev May 23, 2023 15:24
Copy link
Contributor

@yantosca yantosca left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks @lizziel! This is good to merge.

@yantosca yantosca merged commit dbd5fe2 into dev May 23, 2023
@yantosca yantosca deleted the feature/add_mass_accumulation_table_to_1mo_benchmark branch May 23, 2023 21:23
@yantosca yantosca added this to the 1.4.0 milestone May 23, 2023
hannahnesser pushed a commit to hannahnesser/gcpy that referenced this pull request Jan 22, 2024
This merge brings PR geoschem#216 (by @lizziel) into the GCPy 1.4.0
development stream.

PR geoschem#216 adds an option for a mass accumulation table (full column
only) as a 1-month benchmark option. The mass accumulation table
compares the difference in mass change for each species across two
different runs, ref and dev. The mass change per species is computed
as the difference between start and end based on values in the restart
files (end restart mass minus start restart mass). Concentrations are
converted to mass in the same way as done for the benchmark mass tables.

Signed-off-by: Bob Yantosca <yantosca@seas.harvard.edu>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Feature Request New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

[FEATURE REQUEST] Include species delta mass from restart files in full column budget tables
2 participants