Skip to content

Update surface net heat fluxes for AQM coupling#953

Closed
rmontuoro wants to merge 1 commit into
NCAR:mainfrom
rmontuoro:bugfix/aqm-heat-fluxes
Closed

Update surface net heat fluxes for AQM coupling#953
rmontuoro wants to merge 1 commit into
NCAR:mainfrom
rmontuoro:bugfix/aqm-heat-fluxes

Conversation

@rmontuoro
Copy link
Copy Markdown
Contributor

This PR implements proper computation of surface net heat fluxes provided to the coupled NOAA Air Quality Model (AQM).

The update uses PBL-computed fluxes both for fractional and non-fractional grids unless cplflx is set to .true..

Fixes #945.

dtsfci_cpl(i) = cp * rho * hflx_wat(i) ! sensible heat flux over open ocean
dqsfci_cpl(i) = hvap * rho * evap_wat(i) ! latent heat flux over open ocean
else ! use results from PBL scheme for 100% open ocean
if (cplaqm) then
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

When running AQM, it only couples to ATM and ATM doesn't couple to Ocean or Ice model. I would suggest to make it simple by changing the code to:
if (cplaqm .and. .not. cplflx) then
do i = 1, im
dtsfci_cpl(i) = dtsfc1(i)*hffac(i)
dqsfci_cpl(i) = dqsfc1(i)
end do
end if

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

I would prefer to leave in place the ability to use AQM in fully coupled mode for future applications.

Copy link
Copy Markdown
Collaborator

@grantfirl grantfirl 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 OK to me if @ChunxiZhang-NOAA is happy with this.

@ChunxiZhang-NOAA
Copy link
Copy Markdown
Contributor

@grantfirl I am ok with it.

@grantfirl
Copy link
Copy Markdown
Collaborator

Closed in favor of #962

@grantfirl grantfirl closed this Sep 6, 2022
grantfirl added a commit that referenced this pull request Sep 9, 2022
Update surface net heat fluxes for AQM coupling (Was #953)
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.

Surface upward heat fluxes for coupling are not properly calculated when (cplaqm .and. .not.cplflx) is ture in GFS_PBL_generic_post.F90

3 participants