Skip to content

Bug fix in tcdc on isobaric levels #540

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 5 commits into from
Aug 4, 2022
Merged

Bug fix in tcdc on isobaric levels #540

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
80d257c
Restrict undefined points involved in rounding cloud faction uppper a…
WenMeng-NOAA Aug 3, 2022
cd7184e
Merge remote-tracking branch 'upstream/develop' into post_tcdc
WenMeng-NOAA Aug 4, 2022
8c8120f
tweaking changes
WenMeng-NOAA Aug 4, 2022
5cd9e1e
add a change for underground points
WenMeng-NOAA Aug 4, 2022
7a726ea
Remove comment line
WenMeng-NOAA Aug 4, 2022
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
10 changes: 6 additions & 4 deletions sorc/ncep_post.fd/MDL2P.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,6 +30,7 @@
!> 2021-03-11 | B Cui | Change local arrays to dimension (im,jsta:jend)
!> 2021-04-01 | J Meng | Computation on defined points only
!> 2021-07-07 | J MENG | 2D DECOMPOSITION
!> 2022-08-03 | W Meng | Modify total cloud fraction(331)
!>
!> @author T Black W/NP2 @date 1999-09-23
SUBROUTINE MDL2P(iostatusD3D)
Expand Down Expand Up @@ -718,7 +719,7 @@ SUBROUTINE MDL2P(iostatusD3D)
FRIME(I,J) = 1.
RAD(I,J) = 0.
O3SL(I,J) = O3(I,J,LLMH)
CFRSL(I,J) = 0.
IF(CFR(I,J,1)<SPVAL)CFRSL(I,J) = 0.
END IF
! Compute heights by interpolating from heights on interface for NAM but
! hydrostaticJ integration for GFS
Expand Down Expand Up @@ -1316,9 +1317,10 @@ SUBROUTINE MDL2P(iostatusD3D)
DO J=JSTA,JEND
DO I=ISTA,IEND
GRID1(I,J) = SPVAL
CFRSL(I,J) = MIN(MAX(0.0,CFRSL(I,J)),1.0)
IF(abs(CFRSL(I,J)-SPVAL) > SMALL) &
GRID1(I,J) = CFRSL(I,J)*H100
IF(abs(CFRSL(I,J)-SPVAL) > SMALL) THEN
CFRSL(I,J) = MIN(MAX(0.0,CFRSL(I,J)),1.0)
GRID1(I,J) = CFRSL(I,J)*H100
ENDIF
ENDDO
ENDDO
if(grib == 'grib2')then
Expand Down