Skip to content

RRFS cloud base/top: diagnose using cloud fraction #791

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.

Sign up for GitHub

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

Merged
WenMeng-NOAA merged 7 commits into from
Sep 28, 2023
Merged

RRFS cloud base/top: diagnose using cloud fraction #791

Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
18ffd4b
Initial modifications toward a RRFS cloud-base diagnostic that uses c…
jaymes-kenyon Sep 22, 2023
4354d44
For RRFS, assign the cloud-base height at the interface level beneath…
jaymes-kenyon Sep 25, 2023
19de115
Merge branch 'develop' of https://github.com/jaymes-kenyon/UPP into r…
jaymes-kenyon Sep 25, 2023
8a7bf12
Merge branch 'NOAA-EMC:develop' into rrfs_cloudbase
jaymes-kenyon Sep 25, 2023
30ba647
Update change log and comments
jaymes-kenyon Sep 26, 2023
1f2278d
Merge branch 'NOAA-EMC:develop' into rrfs_cloudbase
jaymes-kenyon Sep 26, 2023
ee26910
Merge branch 'rrfs_cloudbase' of https://github.com/jaymes-kenyon/UPP…
jaymes-kenyon Sep 26, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
27 changes: 26 additions & 1 deletion sorc/ncep_post.fd/CLDRAD.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -74,6 +74,7 @@
!> 2023-02-23 | Eric James | Adding coarse PM from RRFS, and using AOD from FV3 for RRFS
!> 2023-04-04 | Li(Kate Zhang) | Add namelist optoin for CCPP-Chem (UFS-Chem) model.
!> 2023-04-17 | Eric James | Getting rid of special treatment for RRFS AOD (use RAP/HRRR approach)
!> 2023-09-26 | Jaymes Kenyon | For RRFS, use cloud fraction to diagnose cloud base/top (height and pressure)
!>
!> @author Russ Treadon W/NP2 @date 1993-08-30
!---------------------------------------------------------------------------------
Expand Down Expand Up @@ -1650,6 +1651,7 @@ SUBROUTINE CLDRAD
!
!--- Grid-scale cloud occurs when the mixing ratio exceeds QCLDmin
! or in the presence of snow when RH>=95% or at/above the PBL top.
! However, for RRFS (FV3R), simply use a threshold cloud fraction.
Comment thread
WenMeng-NOAA marked this conversation as resolved.
!
if(MODELNAME == 'RAPR') then
IBOTGr(I,J)=0
Expand All @@ -1668,6 +1670,21 @@ SUBROUTINE CLDRAD
EXIT
ENDIF
ENDDO !--- End L loop
else if (MODELNAME == 'FV3R') then ! RRFS: use cloud fraction to assign cloud base and cloud top
IBOTGr(I,J)=0
DO L=NINT(LMH(I,J)),1,-1
IF (CFR(I,J,L) >= 0.02) THEN
IBOTGr(I,J)=L
EXIT
ENDIF
ENDDO !--- End L loop
ITOPGr(I,J)=100
DO L=1,NINT(LMH(I,J))
IF (CFR(I,J,L) >= 0.02) THEN
ITOPGr(I,J)=L
EXIT
ENDIF
ENDDO !--- End L loop
else
IBOTGr(I,J) = 0
ZPBLtop = PBLH(I,J)+ZINT(I,J,NINT(LMH(I,J))+1)
Expand Down Expand Up @@ -1704,7 +1721,7 @@ SUBROUTINE CLDRAD
endif
!
!--- Combined (convective & grid-scale) cloud base & cloud top levels
IF(MODELNAME == 'NCAR' .OR. MODELNAME == 'RAPR')THEN
IF(MODELNAME == 'NCAR' .OR. MODELNAME == 'RAPR' .OR. MODELNAME == 'FV3R')THEN
IBOTT(I,J) = IBOTGr(I,J)
ITOPT(I,J) = ITOPGr(I,J)
ELSE
Expand Down Expand Up @@ -1757,6 +1774,14 @@ SUBROUTINE CLDRAD
+ZINT(I,J,IBOT+1)
ENDIF !--- End IF (IBOT == LM) ...
ENDIF !--- End IF (IBOT <= 0) ...
ELSE IF(MODELNAME == 'FV3R') then
IF (IBOT>0 .AND. IBOT<=NINT(LMH(I,J))) THEN
CLDP(I,J) = PINT(I,J,IBOT+1) ! Since IBOT corresponds to a mid-layer location, consider
CLDZ(I,J) = ZINT(I,J,IBOT+1) ! the underlying interfacial level as the cloud base
ELSE
CLDP(I,J) = -50000.
CLDZ(I,J) = -5000.
ENDIF
ELSE
IF (IBOT>0 .AND. IBOT<=NINT(LMH(I,J))) THEN
CLDP(I,J) = PMID(I,J,IBOT)
Expand Down