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

Add fixes for transport tracers simulation #1890

Merged
merged 5 commits into from Jul 28, 2023

Conversation

msulprizio
Copy link
Contributor

Name and Institution (Required)

Name: Melissa Sulprizio
Institution: Harvard / GCST

Confirm you have reviewed the following documentation

Describe the update

This pull request fixes several remaining issues in the transport tracers simulation following #1816. These fixes primarily involve proper calculation and application of masks for applying sources and sinks to latitude zones and vertical regions. An additional fix that was also addressed in the dev/14.2.1 branch is to define MW_g: 1.0 for species with no current molecular weight defined to avoid negative or zero concentrations for those species.

Expected changes

These updates only impact the TransportTracers simulation and will change concentrations of the aoa*, st80_25, and e90(* species.

Transport tracer species aoa*, nh*, St80_25, and StOX did not have molecular
weights defined because they're not traditional chemical species. This causes
issues with unit conversions in the code and results in negative and/or
zero SpeciesConc output. To avoid this, we define MW_g: 1.0 for those
species.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
1. Move initialization of mask arrays into loop over species.
2. Make sure to do calculations for maintaining the mixing ratio only in the defined region (where Mask is 1)
3. Fix typo in the logic for determining masks of latitude zones.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio added this to the 14.2.0 milestone Jul 27, 2023
@msulprizio msulprizio self-assigned this Jul 27, 2023
The stOX tracer is defined in GMAO's TR_GridComp but is not commonly utilized
and isn't needed in GEOS-Chem. It is removed here to avoid confusion. The
fields needed for stOX (GLOBAL_O3 and O3_LOSS) are also removed from
tracer_mod.F90 and TransportTracer configuration files.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
Instead of multiplying the source or sink application by the Mask array,
the code now checks if Mask > 0 before applying the source or sink.
Multiplying by the Mask array may have inadvertently set species concentrations
to zero instead of not implying the source or sink.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
The Budget diagnostics collection is now needed for TransportTracer benchmarking.
It is now disabled by default to avoid saving out unnecessary data.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio marked this pull request as ready for review July 27, 2023 20:01
@msulprizio msulprizio requested a review from lizziel July 27, 2023 20:01
Copy link
Contributor

@lizziel lizziel left a comment

Choose a reason for hiding this comment

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

This looks good to me. We can bring back stOX at a later time if we need to.

@msulprizio
Copy link
Contributor Author

All GCClassic and GCHP integration tests for this update have passed.

@msulprizio msulprizio merged commit 112c4b6 into dev/14.2.0 Jul 28, 2023
@msulprizio msulprizio deleted the bugfix/transport-tracers branch July 28, 2023 14:30
@yantosca yantosca added category: Bug Fix Fixes a previously-reported bug and removed category: Bug Something isn't working labels Aug 1, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Bug Fix Fixes a previously-reported bug
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants