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

HFIP 2022 Hannah Update 2 #162

Open
wants to merge 6 commits into
base: HFIP2022
Choose a base branch
from
Open

Conversation

hannahcbarnes
Copy link

@hannahcbarnes hannahcbarnes commented Sep 14, 2022

Enables GF to use AOD analysis instead of AOD from MERRA2 climatology
Fixes GF - radiation coupling
Fixes dry / moist mixing ratios in GF
Additional tuning to GF
Reverts previous change to sgscloud_radpre

The regression test logs are in /scratch1/BMC/gsd-fv3-dev/hcbarnes/CommitTest/sorc/ufs_model.fd/tests/

@hannahcbarnes
Copy link
Author

@haiqinli Can you look double check these modifications?

@haiqinli
Copy link
Collaborator

@hannahcbarnes The variable of aod_da should be added into cu_gf_driver.meta. The others look good. Thanks.

@hannahcbarnes
Copy link
Author

@haiqinli thanks. aod_da should now be a part of cu_gf_driver.meta

- Minor fixes to pass regression tests
@@ -90,7 +92,8 @@ subroutine sgscloud_radpre_run( &
real(kind=kind_phys), dimension(:,:), intent(inout) :: qr, qs, qg
! qci_conv only allocated if GF is used
real(kind=kind_phys), dimension(:,:), intent(inout) :: qci_conv
real(kind=kind_phys), dimension(:,:), intent(in) :: ud_mf
real(kind=kind_phys), dimension(:,:), intent(in) ::ud_mf
real(kind=kind_phys), dimension(:,:), intent(in) ::ud_mf_timeave, qci_conv_timeave
Copy link
Collaborator

Choose a reason for hiding this comment

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

minor detail but can we preserve the space between :: and the first variable?

Copy link
Author

Choose a reason for hiding this comment

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

Fixed

@@ -328,12 +331,13 @@ subroutine sgscloud_radpre_run( &
else ! scaling function (CB2005)
f = 1.0
endif
sigq = 1.5E-3 * ud_mf(i,k)/dt * f
sigq = ud_mf_timeave(i,k)/dt * f
!sigq = 1.5E-3 * ud_mf(i,k)/dt * f
Copy link
Collaborator

Choose a reason for hiding this comment

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

This is a bit shocking. Are you sure sigq will be ~O(1e-4) without this constant? If sigq is too big, you might always get cloud fractions of around 50%,

Copy link
Author

Choose a reason for hiding this comment

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

I can check. I have output arrays of just GF cloud fraction. Where there are clouds, cloud fraction values tend to be between 30-70%.

Copy link
Collaborator

@joeolson42 joeolson42 left a comment

Choose a reason for hiding this comment

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

I'm OK with the changes so long as the GF developers are sure they are desired. I left a couple questions just to make sure.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

5 participants