From 2f6b91c5309ed5b8a27700b7e768c47af58430ef Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Tue, 28 Sep 2021 12:43:44 +0000 Subject: [PATCH 1/5] Add fixes for 1000 mb HGT corrupted at RRFS 3KM parallel run. --- sorc/ncep_post.fd/MDL2P.f | 5 +++++ sorc/ncep_post.fd/SLP_new.f | 15 +++++---------- 2 files changed, 10 insertions(+), 10 deletions(-) diff --git a/sorc/ncep_post.fd/MDL2P.f b/sorc/ncep_post.fd/MDL2P.f index 8f1d27d79..90baa39fe 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -3838,6 +3838,7 @@ SUBROUTINE MDL2P(iostatusD3D) IF(ABS(SPL(LP)-1.0E5) <= 1.0E-5)THEN IF(LVLS(LP,IGET(012)) > 0)THEN ALPTH = LOG(1.E5) + !IF(MODELNAME == 'GFS'.or.MODELNAME=='FV3R')THEN IF(MODELNAME == 'GFS')THEN ! GFS does not want to adjust 1000 mb H to membrane SLP ! because MOS can't adjust to the much lower H @@ -3846,6 +3847,7 @@ SUBROUTINE MDL2P(iostatusD3D) DO I=1,IM IF(FSL(I,J)1000.)print*,"CCC ", GRID1(I,J) ELSE GRID1(I,J) = SPVAL ENDIF @@ -3866,6 +3868,9 @@ SUBROUTINE MDL2P(iostatusD3D) ENDIF Z1000(I,J) = GRID1(I,J)*GI GRID1(I,J) = Z1000(I,J) + if(abs(GRID1(I,J))>1000.)print*,"AAA ", GRID1(I,J) + ELSE + GRID1(I,J) = SPVAL END IF ENDDO ENDDO diff --git a/sorc/ncep_post.fd/SLP_new.f b/sorc/ncep_post.fd/SLP_new.f index 7815e81aa..e2aa20c0c 100644 --- a/sorc/ncep_post.fd/SLP_new.f +++ b/sorc/ncep_post.fd/SLP_new.f @@ -27,6 +27,8 @@ SUBROUTINE MEMSLP(TPRES,QPRES,FIPRES) ! ARE COMMENTED OUT FOR NOW ! 19-10-30 Bo CUI - REMOVE "GOTO" STATEMENT ! 21-07-26 W Meng - Restrict computation from undefined grids +! 21-09-25 W Meng - Further modification for restricting computation +! from undefined grids. ! ! USAGE: CALL SLPSIG FROM SUBROUITNE ETA2P ! @@ -109,7 +111,7 @@ SUBROUTINE MEMSLP(TPRES,QPRES,FIPRES) TNEW(I,J) = spval - LMHO(I,J) = 0 + LMHO(I,J) = LSM DONE(I,J) = .FALSE. ENDDO ENDDO @@ -127,7 +129,7 @@ SUBROUTINE MEMSLP(TPRES,QPRES,FIPRES) DO J=JSTA,JEND DO I=1,IM - HTMO(I,J,L)=spval + HTMO(I,J,L)=1. if(PSLP(I,J)0.5.AND. !HC 1 HTM2D(I+IHW(J),J-1,L)*HTM2D(I+IHE(J),J-1,L) @@ -278,8 +275,6 @@ SUBROUTINE MEMSLP(TPRES,QPRES,FIPRES) IF(HTM2D(I,J) > 0.5 .AND. tem < 0.5) then TTV(I,J) = TPRES(I,J,L)*(1.+0.608*QPRES(I,J,L)) ENDIF - else - TTV(I,J) = spval ! if(i==ii.and.j==jj)print*,'Debug:L,TTV B SMOO= ',l,TTV(I,J) end if ! spval ENDDO From fbc2c7875f7e62dfd226989091e2ac2611654bb6 Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Tue, 28 Sep 2021 12:58:03 +0000 Subject: [PATCH 2/5] Remove debugging 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 90baa39fe..db84ea521 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -3868,7 +3868,6 @@ SUBROUTINE MDL2P(iostatusD3D) ENDIF Z1000(I,J) = GRID1(I,J)*GI GRID1(I,J) = Z1000(I,J) - if(abs(GRID1(I,J))>1000.)print*,"AAA ", GRID1(I,J) ELSE GRID1(I,J) = SPVAL END IF From 6a3e2c3c9d2b3fed5c455e20d6a8d937981078a5 Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Tue, 28 Sep 2021 13:00:22 +0000 Subject: [PATCH 3/5] Remove debugging 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 db84ea521..9f00878bc 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -3847,7 +3847,6 @@ SUBROUTINE MDL2P(iostatusD3D) DO I=1,IM IF(FSL(I,J)1000.)print*,"CCC ", GRID1(I,J) ELSE GRID1(I,J) = SPVAL ENDIF From 62a187cfb8f161049ee30a543f0324c83b4ab5aa Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Tue, 28 Sep 2021 13:01:50 +0000 Subject: [PATCH 4/5] Remove debugging 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 9f00878bc..9b65497bf 100644 --- a/sorc/ncep_post.fd/MDL2P.f +++ b/sorc/ncep_post.fd/MDL2P.f @@ -3838,7 +3838,6 @@ SUBROUTINE MDL2P(iostatusD3D) IF(ABS(SPL(LP)-1.0E5) <= 1.0E-5)THEN IF(LVLS(LP,IGET(012)) > 0)THEN ALPTH = LOG(1.E5) - !IF(MODELNAME == 'GFS'.or.MODELNAME=='FV3R')THEN IF(MODELNAME == 'GFS')THEN ! GFS does not want to adjust 1000 mb H to membrane SLP ! because MOS can't adjust to the much lower H From df07eaa74bda4022425051104717caa041d0904b Mon Sep 17 00:00:00 2001 From: Wen Meng Date: Tue, 28 Sep 2021 22:04:22 +0000 Subject: [PATCH 5/5] Update VERSION with 10.0.10 --- VERSION | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/VERSION b/VERSION index 82f3d338c..6d6c7c85f 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -10.0.9 +10.0.10