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 fix to cldj_interface_mod.F90 to avoid floating-point exception in GCHP integration tests #2111

Merged
merged 2 commits into from
Jan 11, 2024

Conversation

msulprizio
Copy link
Contributor

@msulprizio msulprizio commented Jan 9, 2024

Name and Institution (Required)

Name: Melissa Sulprizio
Institution: Harvard/GCST

Describe the update

The following line in GeosCore/cldj_interface_mod.F90 was causing a floating- point exception error in GCHP integration tests:

    ! Set top of atmosphere relative humidity to 10% of layer below
    RRR(State_Grid%NZ+1) = RRR(State_Grid%NZ) * 1.d-1

This is resolved by moving the lines for setting RRR to before the loop over the number of layers with clouds.

Please provide a clear and concise overview of the update.

Expected changes

@lizziel wrote:

This fix may cause small differences in the results due to a change in relative humidity above level 34 used for retrieving aerosol optical properties. RRR is the relative humidity, which is used only for aerosols in Cloud-J. The way it was before, relative humidity was only assigned up to level 34 which is defined as where there are clouds. In this fix the relative humidity is set for the full column.

Related Github Issue(s)

…n GCHP

The following line in GeosCore/cldj_interface_mod.F90 was causing a floating-
point exception error in GCHP integration tests:

    ! Set top of atmosphere relative humidity to 10% of layer below
    RRR(State_Grid%NZ+1) = RRR(State_Grid%NZ) * 1.d-1

This is resolved by moving the lines for setting RRR to before the loop
over the number of layers with clouds.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio added the category: Bug Fix Fixes a previously-reported bug label Jan 9, 2024
@msulprizio msulprizio added this to the 14.3.0 milestone Jan 9, 2024
@msulprizio msulprizio requested a review from lizziel January 9, 2024 20:47
@msulprizio msulprizio self-assigned this Jan 9, 2024
@msulprizio msulprizio linked an issue Jan 9, 2024 that may be closed by this pull request
@@ -483,6 +483,11 @@ SUBROUTINE Run_CloudJ( Input_Opt, State_Chm, State_Diag, &
CLDF(1:State_Grid%NZ) = State_Met%CLDF(I,J,1:State_Grid%NZ)
CLDF(State_Grid%NZ+1) = CLDF(State_Grid%NZ)

! Set humidity from input meteorology field
! Convert relative humidity [unitless] from percent to fraction
Copy link
Contributor

Choose a reason for hiding this comment

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

Add one more comment: "! Set top of atmosphere relative humidity to 10% of layer below"

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Thanks for catching that. It's now added in 326c63c.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio requested a review from lizziel January 9, 2024 21:15
@msulprizio msulprizio merged commit 648456b into dev/14.3.0 Jan 11, 2024
@msulprizio msulprizio deleted the bugfix/cloudj-rh-fix branch January 11, 2024 20:02
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.

GCHP integration tests failing for 14.3.0-alpha.6
2 participants