-
Notifications
You must be signed in to change notification settings - Fork 153
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
Fixed a bug related to sensible/latent heat fluxes for AQM coupling (Issue#945) #951
Conversation
…n GFS_PBL_generic_post.F90
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me.
Does anyone know if @rmontuoro (the original author of this section) agrees with the solution? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'm assuming that the bug is a result of copy/paste for the case when cplflx = T above? Regardless, the fix should work to fix the reported issue.
@grantfirl Maybe we should invite Raffaele to review this code change. |
@ChunxiZhang-NOAA, @grantfirl - Thanks. I will take a look and get back to you. |
The proposed changes do not address the lack of surface net heat fluxes over land, which need to be provided to AQM. A new formulation is proposed in PR #953. |
@ChunxiZhang-NOAA @rmontuoro Should we close this (and upstream PRs) in favor of merging #953 (and upstream PRs)? |
@grantfirl We can close this PR. |
Closing this PR in favor of continuing discussion of the solution in #953 |
A bug is fixed which is described in Issue#945.
Currently, this PR only fixed the bug that caused model crashes due to the unallocated variables dtsfc_cice and dqsfc_cice which are only allocated when cplflx=.true. in GFS_typedefs.F90.
There are potential other issues. For instance, when (cplaqm .and. .not.cplflx) is ture, dtsfci_cpl and dqsfci_cpl are only assigned values over ocean (e.g., oceanfrac > 0).