Skip to content

Add optional scaling to RRTMGP flux adjustment#668

Merged
grantfirl merged 14 commits into
NCAR:mainfrom
dustinswales:hotfix_AddPLimitToLWJacAdj
Jul 9, 2021
Merged

Add optional scaling to RRTMGP flux adjustment#668
grantfirl merged 14 commits into
NCAR:mainfrom
dustinswales:hotfix_AddPLimitToLWJacAdj

Conversation

@dustinswales
Copy link
Copy Markdown
Member

This PR contains improvements to how the LW heating-rate profile is updated between calls to the RRTMGP radiation scheme.
The adjusted net LW flux is computed following Hogan and Bozzo 2015 (10.1002/2015MS000455)
Optionally, the flux adjustment can be damped with height using a logistic function. This was found to be necessary when using configurations with very high/thin layers (e.g. L127). The shape of the damping function can be controlled by namelist parameters.

Copy link
Copy Markdown
Collaborator

@climbfuji climbfuji left a comment

Choose a reason for hiding this comment

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

Looks pretty straightforward, just one question for my understanding.

Comment thread physics/dcyc2.f
! Vertical ordering?
if (p_lev(1,1) .lt. p_lev(1, levs)) then
iSFC = levs
iSFC = levs + 1
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

Was that an actual bug? Is iSFC for interfaces (which is always levs+1)?

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

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

Yes^2.

@climbfuji
Copy link
Copy Markdown
Collaborator

This PR has been pulled into #680 and will be merged as part of it.

grantfirl added a commit that referenced this pull request Jul 9, 2021
Add extended diagnostic output from Thompson MP, includes "Add optional scaling to RRTMGP flux adjustment" (#668), fix index bug in m_micro.F90
@grantfirl grantfirl merged commit 3542434 into NCAR:main Jul 9, 2021
@dustinswales dustinswales deleted the hotfix_AddPLimitToLWJacAdj branch September 17, 2021 17:34
JohanaRomeroAlvarez pushed a commit to JohanaRomeroAlvarez/ccpp-physics that referenced this pull request Sep 8, 2025
…ys that are too small (NCAR#668)

* move utility routines, data structures, and name generation out of FV3GFS_io.F90

* move RRFS SD state restart io to separate file and rename clm_lake_io.F90

* move rrfs sd emissions input code to FV3GFS_rrfs_sd_io.F90

* register sfc fields in FV3GFS_sfc_io.F90

* move oro reading to FV3GFS_oro_io.F90

* gwd I/O in FV3GFS_oro_io.F90

* store quilt restart data structures in FV3GFS_sfc_io types

* rrfs sd restart in quilt (still under testing)

* bug fixes for RRFS 32-bit physics
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.

3 participants