From 80d257c402c96da4dd7546416fa88d01cc7551d3 Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Wed, 3 Aug 2022 13:55:38 +0000 Subject: [PATCH 1/4] Restrict undefined points involved in rounding cloud faction uppper and lower bounds --- sorc/ncep_post.fd/MDL2P.f | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index ccb2b5c43..cf613be70 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -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) @@ -1316,9 +1317,9 @@ 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 + CFRSL(I,J) = MIN(MAX(0.0,CFRSL(I,J)),1.0) + GRID1(I,J) = CFRSL(I,J)*H100 ENDDO ENDDO if(grib == 'grib2')then From 8c8120f9cd66f6da4fdfc3a858d6379ad7630991 Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Thu, 4 Aug 2022 12:41:51 +0000 Subject: [PATCH 2/4] tweaking changes --- sorc/ncep_post.fd/MDL2P.f | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index cf613be70..5ca7ae64a 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -1317,9 +1317,10 @@ SUBROUTINE MDL2P(iostatusD3D) DO J=JSTA,JEND DO I=ISTA,IEND GRID1(I,J) = SPVAL - IF(abs(CFRSL(I,J)-SPVAL) > SMALL) & + 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 From 5cd9e1e918eb307c39ff12593980a7fb8d38eec2 Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Thu, 4 Aug 2022 15:47:00 +0000 Subject: [PATCH 3/4] add a change for underground points --- sorc/ncep_post.fd/MDL2P.f | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index 5ca7ae64a..9777a7abb 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -719,7 +719,8 @@ SUBROUTINE MDL2P(iostatusD3D) FRIME(I,J) = 1. RAD(I,J) = 0. O3SL(I,J) = O3(I,J,LLMH) - CFRSL(I,J) = 0. + !wmCFRSL(I,J) = 0. + IF(CFR(I,J,1) Date: Thu, 4 Aug 2022 15:48:09 +0000 Subject: [PATCH 4/4] Remove comment line --- sorc/ncep_post.fd/MDL2P.f | 1 - 1 file changed, 1 deletion(-) diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index 9777a7abb..ce3269951 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -719,7 +719,6 @@ SUBROUTINE MDL2P(iostatusD3D) FRIME(I,J) = 1. RAD(I,J) = 0. O3SL(I,J) = O3(I,J,LLMH) - !wmCFRSL(I,J) = 0. IF(CFR(I,J,1)