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

Restore consideration of sea salt alkalinity to heterogenerous acid-catalyzed reactions of halogens on sea salt aerosols (now fixing typos) #1952

Merged
merged 3 commits into from Sep 11, 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 is the corresponding pull request to #1547, now fixing typos in the calls to Ars_L1K. In several calls to that routine the gamma and SrMW arguments were swapped, leading to integration errors in KPP. Fixing these, seems to get rid of the errors and allows us to properly restore consideration of sea salt alkalinity to heterogenerous acid-catalyzed reactions of halogens on sea salt aerosols.

Tagging @beckyalexander @viral211 @yantosca.

Expected changes

See the original description of the issue and fix in #1547.

Reference(s)

If this is a science update, please provide a literature citation.

Related Github Issue(s)

…atalyzed reactions of halogens on sea salt aerosols

These are the fixes implemented originally by Becky Alexander in #1702 as
discussed in issue #1547.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
…cid* functions

In fullchem_RateLawFuncs.F90, the arguments passed to function Ars_L1k were
not in the correct order - area, radius, gamma, srMw. In the following functions,
gamma and srMw were swapped:

  - IuptkbyAlkSALA1stOrd
  - IuptkByAlkSALC1stOrd
  - IbrkdnbyAcidBrSALA
  - IbrkdnbyAcidBrSALC
  - IbrkdnbyAcidSALACl
  - IbrkdnbyAcidSALCCl

These typos were introduced in the original PRs (#1548, #1702) for these
bug fixes, possible because the fixes were implemented off of an older
version of GEOS-Chem (13.4.1).

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio added category: Bug Something isn't working topic: Chemical Mechanisms Related to KPP and/or GEOS-Chem chemistry mechanisms Attn: Chemistry WG Attention: Chemistry Working Group labels Sep 11, 2023
@msulprizio msulprizio added this to the 14.2.1 milestone Sep 11, 2023
@msulprizio msulprizio self-assigned this Sep 11, 2023
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.

Hi @msulprizio. Thanks for this. I don't quite see where the incorrect argument sequences were. Were they in the IuptkBySALA1stOrd and IuptkBySALC1stOrd functions? Thanks.

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.

Never mind @msulprizio. I understand now. Thanks for the fix! Good to merge.

@msulprizio
Copy link
Contributor Author

msulprizio commented Sep 11, 2023

Thanks @yantosca . For clarity for others reading this, the sequence of arguments passed to Ars_L1K are corrected in commit 76004d9.

…#1933

Arguments for id_SAL* needed to be added to fullchem_SetStateHet in
stub_fullchem_HetStateFuncs.F90 for consistency with fullchem_HetStateFuncs.F90.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
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.

Oops, sorry I forgot to mention the stub file issue. Now good to go.

@msulprizio msulprizio merged commit 5218fe8 into dev/14.2.0 Sep 11, 2023
@msulprizio
Copy link
Contributor Author

All GCClassic and GCHP integration tests have passed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Attn: Chemistry WG Attention: Chemistry Working Group category: Bug Something isn't working topic: Chemical Mechanisms Related to KPP and/or GEOS-Chem chemistry mechanisms
Projects
None yet
2 participants