Fixes for 1000 mb HGT corrupted at RRFS 3KM parallel run. (#396)

* Add fixes for 1000 mb HGT corrupted at RRFS 3KM parallel run.

* Remove debugging line.

* Remove debugging line.

* Remove debugging line.

* Update VERSION with 10.0.10

VERSION

@@ -1 +1 @@
-10.0.9
+10.0.10

sorc/ncep_post.fd/MDL2P.f

@@ -3866,6 +3866,8 @@ SUBROUTINE MDL2P(iostatusD3D)
                      ENDIF
                      Z1000(I,J) = GRID1(I,J)*GI
                      GRID1(I,J) = Z1000(I,J)
+                    ELSE
+                     GRID1(I,J) = SPVAL
                     END IF
                    ENDDO
                  ENDDO    

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)<spval) then
 
             PSFC = PSLP(I,J)
@@ -259,12 +261,7 @@ SUBROUTINE MEMSLP(TPRES,QPRES,FIPRES)
         DO J=JSTA_M,JEND_M
           DO I=2,IM-1
 
-           !if(PSLP(I,J)<spval) then
-           if(HTM2D(I-1,J)<spval.and.HTM2D(I+1,J)<spval.and. &
-              HTM2D(I,J-1)<spval.and.HTM2D(I,J+1)<spval.and. &
-              HTM2D(I-1,J-1)<spval.and.HTM2D(I+1,J-1)<spval.and. &
-              HTM2D(I-1,J+1)<spval.and.HTM2D(I+1,J+1)<spval.and. &
-              TPRES(I,J,L)<spval.and.QPRES(I,J,L)<spval)then
+           if(PSLP(I,J)<spval) then
 
 !HC        IF(HTM2D(I,J,L)>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