diff --git a/sorc/ncep_post.fd/CALLCL.f b/sorc/ncep_post.fd/CALLCL.f
index c0cd037b2..222a969a9 100644
--- a/sorc/ncep_post.fd/CALLCL.f
+++ b/sorc/ncep_post.fd/CALLCL.f
@@ -25,7 +25,8 @@
 !> 2002-04-24 | Mike Baldwin | WRF Version            
 !> 2019-10-30 | Bo Cui       | Remove "GOTO" Statement
 !> 2021-07-28 | W Meng       | Restriction compuatation from undefined grids
-!> 2021-09-02 | Bo Cui       | Decompose UPP in X direction          
+!> 2021-09-02 | Bo Cui       | Decompose UPP in X direction     
+!> 2026-03-27 | Alyson Stahl | Remove shared DO termination labels     
 !>
 !> @author Russ Treadon W/NP2 @date 1993-03-15
 !-----------------------------------------------------------------------
@@ -80,8 +81,8 @@ SUBROUTINE CALLCL(P1D,T1D,Q1D,PLCL,ZLCL)
 !
 ! Bo Cui 10/30/2019, remove "GOTO" statement
 
-      DO 30 J=JSTA_M,JEND_M
-      DO 30 I=ISTA_M,IEND_M
+      DO J=JSTA_M,JEND_M
+      DO I=ISTA_M,IEND_M
       IF(P1D(I,J)<spval.and.Q1D(I,J)<spval)THEN
       EVP       = P1D(I,J)*Q1D(I,J)/(EPS+ONEPS*Q1D(I,J))
       RMX       = EPS*EVP/(P1D(I,J)-EVP)
@@ -103,7 +104,8 @@ SUBROUTINE CALLCL(P1D,T1D,Q1D,PLCL,ZLCL)
       ENDIF
  20   CONTINUE
       ENDIF
- 30   CONTINUE
+      ENDDO
+      ENDDO
 !     
 !     END OF ROUTINE.
 !     
diff --git a/sorc/ncep_post.fd/CALWXT.f b/sorc/ncep_post.fd/CALWXT.f
index 48aac1332..d05f47cf5 100644
--- a/sorc/ncep_post.fd/CALWXT.f
+++ b/sorc/ncep_post.fd/CALWXT.f
@@ -11,6 +11,7 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
 !     19-10-30  Bo CUI - REMOVE "GOTO" STATEMENT
 !     21-07-26  Wen Meng - Restrict computation from undefined grids
 !     21-10-31  JESSE MENG - 2D DECOMPOSITION
+!     26-03-27  Alyson Stahl - Remove shared DO termination labels
 !                              
 !
 !     ROUTINE TO COMPUTE PRECIPITATION TYPE USING A DECISION TREE
@@ -99,8 +100,8 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
 
 !
 !!$omp  parallel do private(a,lmhk,pkl,psfck,qkl,tdchk,tdkl,tdpre,tkl)
-      DO 800 J=JSTA,JEND
-      DO 800 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       LMHK=NINT(LMH(I,J))
 !
 !   SKIP THIS POINT IF NO PRECIP THIS TIME STEP 
@@ -152,12 +153,13 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
         jcontinue=.false.
       ENDIF
       enddo     ! enddo jcontinue
-  800 CONTINUE
+      ENDDO
+      ENDDO
 !
 !    LOWEST LAYER T
 !
-      DO 850 J=JSTA,JEND
-      DO 850 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       KARR(I,J)=0
       IF (PREC(I,J)<=PTHRESH) cycle    
       LMHK=NINT(LMH(I,J))
@@ -185,7 +187,8 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
           ENDIF
       ENDIF
       KARR(I,J)=1
-  850 CONTINUE
+      ENDDO
+      ENDDO
 !
 !   COMPUTE WET BULB ONLY AT POINTS THAT NEED IT
 !
@@ -196,8 +199,8 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
 !    & private(area1,areap4,areas8,dzkl,ifrzl,iwrml,lice,          &
 !    &         lmhk,pintk1,pintk2,pm150,psfck,surfc,surfw,         &
 !    &         tlmhk,twrmk)
-      DO 1900 J=JSTA,JEND
-      DO 1900 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
 !       IF (I == 324 .AND. J == 390) THEN
 !          LMHK=NINT(LMH(I,J))
 !          DO L=LMHK,1,-1          
@@ -313,7 +316,8 @@ SUBROUTINE CALWXT_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX,ZWET)
           IWX(I,J)=IWX(I,J)+8
         ENDIF
       ENDIF
- 1900 CONTINUE
+      ENDDO
+      ENDDO
 !---------------------------------------------------------
       DEALLOCATE (TWET)
 
diff --git a/sorc/ncep_post.fd/CALWXT_RAMER.f b/sorc/ncep_post.fd/CALWXT_RAMER.f
index 1985c67e1..64cbf8d3e 100644
--- a/sorc/ncep_post.fd/CALWXT_RAMER.f
+++ b/sorc/ncep_post.fd/CALWXT_RAMER.f
@@ -78,8 +78,8 @@ SUBROUTINE CALWXT_RAMER_POST(T,Q,PMID,PINT,LMH,PREC,PTYP)
       enddo
 
 !  BIG LOOP
-      DO 800 J=JSTA,JEND
-      DO 800 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
 !
 !   SKIP THIS POINT IF NO PRECIP THIS TIME STEP
 !
@@ -373,8 +373,8 @@ SUBROUTINE CALWXT_RAMER_POST(T,Q,PMID,PINT,LMH,PREC,PTYP)
       END IF
       IF (trace) WRITE (*,*) "Returned ptyp is:ptyp,lll ", ptyp, lll,'me=',me
       IF (trace) WRITE (*,*) "Returned icefrac is: ", icefrac,'me=',me
- 800  CONTINUE 
-
+      ENDDO
+      ENDDO
       RETURN
 !
       END
diff --git a/sorc/ncep_post.fd/CALWXT_REVISED.f b/sorc/ncep_post.fd/CALWXT_REVISED.f
index 792680d09..5b22fde21 100644
--- a/sorc/ncep_post.fd/CALWXT_REVISED.f
+++ b/sorc/ncep_post.fd/CALWXT_REVISED.f
@@ -12,6 +12,7 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
 !                TO MAKE AN ALTERNATE ALGORITHM 
 !     19-10-30  Bo CUI - REMOVE "GOTO" STATEMENT
 !     21-10-31  JESSE MENG - 2D DECOMPOSITION
+!     26-03-27  Alyson Stahl - Remove shared DO termination labels
 !                              
 !
 !     ROUTINE TO COMPUTE PRECIPITATION TYPE USING A DECISION TREE
@@ -88,8 +89,8 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
 !
 !!$omp  parallel do
 !!$omp& private(a,lmhk,pkl,psfck,qkl,tdchk,tdkl,tdpre,tkl)
-      DO 800 J=JSTA,JEND
-      DO 800 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       LMHK=NINT(LMH(I,J))
 !
 !   SKIP THIS POINT IF NO PRECIP THIS TIME STEP 
@@ -141,12 +142,13 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
         jcontinue=.false.
       ENDIF
       enddo     ! enddo jcontinue
-  800 CONTINUE
+      ENDDO
+      ENDDO
 !
 !    LOWEST LAYER T
 !
-      DO 850 J=JSTA,JEND
-      DO 850 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       KARR(I,J)=0
       IF (PREC(I,J)<=PTHRESH) cycle    
       LMHK=NINT(LMH(I,J))
@@ -174,7 +176,8 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
           ENDIF
       ENDIF
       KARR(I,J)=1
-  850 CONTINUE
+      ENDDO
+      ENDDO
 !
 !   COMPUTE WET BULB ONLY AT POINTS THAT NEED IT
 !
@@ -184,8 +187,8 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
 !!$omp& private(area1,areap4,areap0,areas8,dzkl,ifrzl,iwrml,lice,
 !!$omp&         lmhk,pintk1,pintk2,pm150,psfck,surfc,surfw,
 !!$omp&         tlmhk,twrmk)
-      DO 1900 J=JSTA,JEND
-      DO 1900 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       IF(KARR(I,J)>0)THEN
         LMHK=NINT(LMH(I,J))
         LICE=LICEE(I,J)
@@ -307,7 +310,8 @@ SUBROUTINE CALWXT_REVISED_POST(T,Q,PMID,PINT,HTM,LMH,PREC,ZINT,IWX)
           IWX(I,J)=IWX(I,J)+8
         ENDIF
       ENDIF
- 1900 CONTINUE
+      ENDDO
+      ENDDO
 !      print *, 'revised check ', IWX(500,800)
 
 !---------------------------------------------------------
diff --git a/sorc/ncep_post.fd/FRZLVL.f b/sorc/ncep_post.fd/FRZLVL.f
index 6fd1b3c43..ac471a756 100644
--- a/sorc/ncep_post.fd/FRZLVL.f
+++ b/sorc/ncep_post.fd/FRZLVL.f
@@ -36,6 +36,7 @@
 !> 2019-10-30 | Bo Cui       | Remove "GOTO" statement
 !> 2020-11-10 | Jesse Meng   | Use UPP_PHYSICS module
 !> 2021-10-15 |JESSE MENG    | 2D DECOMPOSITION
+!> 2026-03-27 | Alyson Stahl | Remove shared DO termination labels
 !>
 !> @author Russ Treadon W/NP2 @date 1992-12-22
 !-----------------------------------------------------------------------------
@@ -80,8 +81,8 @@ SUBROUTINE FRZLVL(ZFRZ,RHFRZ,PFRZL)
 !    &         qsat,qsfc,qsfrz,rhsfc,rhz,tsfc,                        &
 !    &         zl,zu)
 
-       DO 20 J=JSTA,JEND
-       DO 20 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
          HTSFC    = FIS(I,J)*GI
          LLMH     = NINT(LMH(I,J))
          RHFRZ(I,J) = D00
@@ -225,7 +226,8 @@ SUBROUTINE FRZLVL(ZFRZ,RHFRZ,PFRZL)
                EXIT             
             ENDIF
  10      CONTINUE
-20   CONTINUE
+      ENDDO
+      ENDDO
 !     
 !     END OF ROUTINE.
 !     
diff --git a/sorc/ncep_post.fd/FRZLVL2.f b/sorc/ncep_post.fd/FRZLVL2.f
index e1537a57a..112656eb8 100644
--- a/sorc/ncep_post.fd/FRZLVL2.f
+++ b/sorc/ncep_post.fd/FRZLVL2.f
@@ -36,6 +36,7 @@
 !> 2020-11-10 | Jesse Meng   | Use UPP_PHYSICS module
 !> 2021-10-15 | JESSE MENG   | 2D DECOMPOSITION
 !> 2021-07-28 | W. Meng      | Restrict compuatation from undefined grids
+!> 2026-03-27 | Alyson Stahl | Remove shared DO termination labels
 !>
 !> @author Russ Treadon W/NP2 @date 1992-12-22
 !-------------------------------------------------------------------------------
@@ -80,8 +81,8 @@ SUBROUTINE FRZLVL2(ISOTHERM,ZFRZ,RHFRZ,PFRZL)
 !     LOOP OVER HORIZONTAL GRID.
 !     
 
-      DO 20 J=JSTA,JEND
-      DO 20 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
          IF(FIS(I,J)<spval)THEN
          HTSFC    = FIS(I,J)*GI
          LLMH     = NINT(LMH(I,J))
@@ -230,7 +231,8 @@ SUBROUTINE FRZLVL2(ISOTHERM,ZFRZ,RHFRZ,PFRZL)
                RHFRZ(I,J) = spval
                ZFRZ(I,J)  = spval
          ENDIF
- 20   CONTINUE
+      ENDDO
+      ENDDO
 !     
 !     END OF ROUTINE.
 !     
diff --git a/sorc/ncep_post.fd/LFMFLD.f b/sorc/ncep_post.fd/LFMFLD.f
index f44c8c5b9..53e86ac78 100644
--- a/sorc/ncep_post.fd/LFMFLD.f
+++ b/sorc/ncep_post.fd/LFMFLD.f
@@ -39,6 +39,7 @@
 !> 2019-10-30 | Bo Cui       | Remove "GOTO" statement
 !> 2020-11-10 | Jesse Meng   | Use UPP_PHYSICS Module
 !> 2021-10-14 | JESSE MENG   | 2D DECOMPOSITION
+!> 2026-03-27 | Alyson Stahl | Remove shared DO termination labels
 !>
 !> @author Russ Treadon W/NP2 @date 1992-12-22
 !--------------------------------------------------------------------------------------
@@ -80,8 +81,8 @@ SUBROUTINE LFMFLD(RH3310,RH6610,RH3366,PW3310)
 !
 !     LOOP OVER HORIZONTAL GRID.
 !     
-      DO 30 J=JSTA,JEND
-      DO 30 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
 !     
 !        ZERO VARIABLES.
          RH3310(I,J) = D00
@@ -172,7 +173,8 @@ SUBROUTINE LFMFLD(RH3310,RH6610,RH3366,PW3310)
          ELSE
             RH3366(I,J) = SPVAL
          ENDIF
- 30   CONTINUE
+      ENDDO
+      ENDDO
 !     
 !     
 !     END OF ROUTINE.
diff --git a/sorc/ncep_post.fd/LFMFLD_GFS.f b/sorc/ncep_post.fd/LFMFLD_GFS.f
index e23729898..59efb8705 100644
--- a/sorc/ncep_post.fd/LFMFLD_GFS.f
+++ b/sorc/ncep_post.fd/LFMFLD_GFS.f
@@ -40,6 +40,7 @@
 !> 2019-10-30 | Bo Cui       | Remove "GOTO" statement
 !> 2020-11-10 | Jesse Meng   | Use UPP_PHYSICS Module
 !> 2021-10-14 | JESSE MENG   | 2D DECOMPOSITION
+!> 2026-03-27 | Alyson Stahl | Remove shared DO termination labels
 !>
 !> @author Russ Treadon W/NP2 @date 1992-12-22
 !---------------------------------------------------------------------------
@@ -87,8 +88,8 @@ SUBROUTINE LFMFLD_GFS(RH4410,RH7294,RH4472,RH3310)
 !
 !     LOOP OVER HORIZONTAL GRID.
 !     
-      DO 30 J=JSTA,JEND
-      DO 30 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
 !     
 !        ZERO VARIABLES.
          RH4410(I,J) = D00
@@ -194,7 +195,8 @@ SUBROUTINE LFMFLD_GFS(RH4410,RH7294,RH4472,RH3310)
          ELSE
             RH7294(I,J) = SPVAL
          ENDIF
- 30   CONTINUE
+      ENDDO
+      ENDDO
 !     
 !     END OF ROUTINE.
 !     
diff --git a/sorc/ncep_post.fd/MDL2AGL.f b/sorc/ncep_post.fd/MDL2AGL.f
index 19bfedb20..bc4a6ed8e 100644
--- a/sorc/ncep_post.fd/MDL2AGL.f
+++ b/sorc/ncep_post.fd/MDL2AGL.f
@@ -20,6 +20,7 @@
 !! 2023-03-02  S TRAHAN - copy lightning threat index 3 element-by-element
 !! 2023-10-23  J Kenyon - HAILCAST output enabled in RRFS
 !! 2025-04-01  W Meng - Bug fix in HAILCAST
+!! 2026-03-27  A Stahl - Remove shared DO termination labels
 !!     
 !! USAGE:    CALL MDL2P
 !!   INPUT ARGUMENT LIST:
@@ -1037,8 +1038,8 @@ SUBROUTINE MDL2AGL
 	   call exch(VH(ISTA_2L:IEND_2U,JSTA_2L:JEND_2U,LL))
 	  END DO
 	 END IF   
-         DO 230 J=JSTART,JSTOP
-         DO 230 I=ISTART,ISTOP
+        DO J=JSTART,JSTOP
+        DO I=ISTART,ISTOP
           LL=NL1X(I,J)
 !---------------------------------------------------------------------
 !***  VERTICAL INTERPOLATION OF GEOPOTENTIAL, TEMPERATURE, SPECIFIC
@@ -1128,7 +1129,8 @@ SUBROUTINE MDL2AGL
      &	       VH(IE,JS,NINT(LMV(IE,JS)))+VH(IW,JS,NINT(LMV(IW,JS))))/4.0
            END IF
           END IF
-  230    CONTINUE
+        ENDDO
+        ENDDO
 !
 !     
 !---------------------------------------------------------------------
@@ -1238,8 +1240,8 @@ SUBROUTINE MDL2AGL
 !chc        I=IHOLD(NN)
 !chc        J=JHOLD(NN)
 !        DO 220 J=JSTA,JEND
-            DO 240 J=JSTA_2L,JEND_2U
-              DO 240 I=ISTA_2L,IEND_2U
+            DO J=JSTA_2L,JEND_2U
+              DO I=ISTA_2L,IEND_2U
                 LL = NL1X(I,J)
 !---------------------------------------------------------------------
 !***  VERTICAL INTERPOLATION OF GEOPOTENTIAL, TEMPERATURE, SPECIFIC
@@ -1295,7 +1297,8 @@ SUBROUTINE MDL2AGL
                 UAGL(I,J) = UH(I,J,NINT(LMV(I,J)))
                 VAGL(I,J) = VH(I,J,NINT(LMV(I,J)))
               END IF
-  240 CONTINUE
+      ENDDO
+      ENDDO
 !
 !
 !---------------------------------------------------------------------
diff --git a/sorc/ncep_post.fd/MDL2SIGMA.f b/sorc/ncep_post.fd/MDL2SIGMA.f
index f34e6b861..afe234c95 100644
--- a/sorc/ncep_post.fd/MDL2SIGMA.f
+++ b/sorc/ncep_post.fd/MDL2SIGMA.f
@@ -22,6 +22,7 @@
 !!   21-03-11  B Cui - change local arrays to dimension (im,jsta:jend)
 !!   21-10-14  J MENG - 2D DECOMPOSITION
 !!   2022-09-01 S Trahan - fixed bugs where extreme atmospheric conditions can cause out-of-bounds access
+!!   26-03-27  Alyson Stahl - Remove shared DO termination labels
 !!  
 !! USAGE:    CALL MDL2P
 !!   INPUT ARGUMENT LIST:
@@ -210,8 +211,8 @@ SUBROUTINE MDL2SIGMA
         ENDDO
        END DO	
        END DO
-       DO 167 J=JSTA,JEND
-        DO 167 I=ISTA_2L,IEND_2U
+       DO J=JSTA,JEND
+        DO I=ISTA_2L,IEND_2U
 	 DONEFSL1=.FALSE.
          PFSIGO=PTSIGO
          APFSIGO=LOG(PFSIGO)
@@ -307,7 +308,8 @@ SUBROUTINE MDL2SIGMA
      &      AKH(I,J)=EXCH_H(I,J,LL-1)+(EXCH_H(I,J,LL-1)                 &
      &	    -EXCH_H(I,J,LL-2))*FACT 
 	 END IF    
- 167    CONTINUE
+        ENDDO
+        ENDDO
 ! OUTPUT FIRST LAYER GEOPOTENTIAL
 !       GEOPOTENTIAL (SCALE BY GI)
         IF (IGET(205)>0) THEN
@@ -407,12 +409,12 @@ SUBROUTINE MDL2SIGMA
 !$omp  parallel do private(i,j,ll,llmh,psigo,apsigo,fact,dum,pl,     &
 !$omp &         zl,tl,ql,ai,bi,qsat,rhl,tvrl,tvrblo,tblo,tmt0,       &
 !$omp &         qblo,pnl1,fac,ahf)
-!hc        DO 220 NN=1,NHOLD
+!hc        DO NN=1,NHOLD
 !hc        I=IHOLD(NN)
 !hc        J=JHOLD(NN)
-        DO 220 J=JSTA,JEND      ! Moorthi on Nov 26 2014
-!       DO 220 J=JSTA_2L,JEND_2U
-        DO 220 I=ISTA,IEND
+        DO J=JSTA,JEND      ! Moorthi on Nov 26 2014
+!       DO J=JSTA_2L,JEND_2U
+        DO I=ISTA,IEND
         LL=NL1X(I,J)
 !---------------------------------------------------------------------
 !***  VERTICAL INTERPOLATION OF GEOPOTENTIAL, TEMPERATURE, SPECIFIC
@@ -556,7 +558,8 @@ SUBROUTINE MDL2SIGMA
 	  QG1(I,J)=0.
           CFRSIG(I,J)=0.
         END IF
-  220   CONTINUE
+        ENDDO
+        ENDDO
 !
 ! OBTAIN GEOPOTENTIAL AND KH ON INTERFACES 
        DO J=JSTA_2L,JEND_2U
@@ -792,9 +795,9 @@ SUBROUTINE MDL2SIGMA
         ENDDO
         ENDDO
 !
-        DO 230 J=JSTA,JEND
-!        DO 230 I=1,IM-MOD(j,2)
-        DO 230 I=ISTA,IEND-MOD(j,2) !Jesse 20211014
+        DO J=JSTA,JEND
+!        DO I=1,IM-MOD(j,2)
+        DO I=ISTA,IEND-MOD(j,2) !Jesse 20211014
 
          LLMH = NINT(LMH(I,J))
 
@@ -865,7 +868,8 @@ SUBROUTINE MDL2SIGMA
           IF(UH(I,J,LLMH)<SPVAL)USL(I,J)=UH(I,J,LLMH)
 	  IF(VH(I,J,LLMH)<SPVAL)VSL(I,J)=VH(I,J,LLMH)
         END IF
-  230   CONTINUE
+        ENDDO
+        ENDDO
         JJB=JSTA 
         IF(MOD(JSTA,2)==0)JJB=JSTA+1
         JJE=JEND
@@ -905,8 +909,8 @@ SUBROUTINE MDL2SIGMA
          ENDDO
          ENDDO
 !
-         DO 231 J=JSTA,JEND_M
-         DO 231 I=ISTA,IEND_M
+         DO J=JSTA,JEND_M
+         DO I=ISTA,IEND_M
 	  PDV=0.25*(PINT(I,J,LP1)+PINT(I+1,J,LP1)                       &
              +PINT(I,J+1,LP1)+PINT(I+1,J+1,LP1))
           PSIGO=PTSIGO+ASIGO(LP)*(PDV-PTSIGO)
@@ -943,7 +947,8 @@ SUBROUTINE MDL2SIGMA
            IF(UH(I,J,LLMH)<SPVAL)USL(I,J)=UH(I,J,LLMH)
 	   IF(VH(I,J,LLMH)<SPVAL)VSL(I,J)=VH(I,J,LLMH)
           END IF
-  231   CONTINUE
+        ENDDO
+        ENDDO
 	
 	
 	
diff --git a/sorc/ncep_post.fd/MDL2SIGMA2.f b/sorc/ncep_post.fd/MDL2SIGMA2.f
index a02107e10..5bbbfe88a 100644
--- a/sorc/ncep_post.fd/MDL2SIGMA2.f
+++ b/sorc/ncep_post.fd/MDL2SIGMA2.f
@@ -21,6 +21,7 @@
 !!   21-03-11  B Cui - change local arrays to dimension (im,jsta:jend)
 !!   21-07-26  W Meng - Restrict compuatation from undefined grids
 !!   21-10-14  J MENG - 2D DECOMPOSITION
+!!   26-03-27  Alyson Stahl - Remove shared DO termination labels
 !!  
 !! USAGE:    CALL MDL2P
 !!   INPUT ARGUMENT LIST:
@@ -176,8 +177,8 @@ SUBROUTINE MDL2SIGMA2
 !hc        J=JHOLD(NN)
 !         DO 220 J=JSTA,JEND
 !         DO 220 J=JSTA_2L,JEND_2U
-          DO 220 J=JSTA,JEND           ! Moorthi on Nov 26, 2014
-            DO 220 I=ISTA,IEND
+          DO J=JSTA,JEND           ! Moorthi on Nov 26, 2014
+            DO I=ISTA,IEND
               LL=NL1X(I,J)
 !---------------------------------------------------------------------
 !***  VERTICAL INTERPOLATION OF GEOPOTENTIAL, TEMPERATURE, SPECIFIC
@@ -252,7 +253,8 @@ SUBROUTINE MDL2SIGMA2
 !
                 TSL(I,J) = TBLO
               END IF
-  220       CONTINUE
+            ENDDO
+            ENDDO
 
 !---------------------------------------------------------------------
 !        *** PART II ***
diff --git a/sorc/ncep_post.fd/WETBULB.f b/sorc/ncep_post.fd/WETBULB.f
index f63b9c73b..e1b9bf712 100644
--- a/sorc/ncep_post.fd/WETBULB.f
+++ b/sorc/ncep_post.fd/WETBULB.f
@@ -9,6 +9,7 @@ SUBROUTINE WETBULB(T,Q,PMID,HTM,KARR,TWET)
 !     02-01-15  MIKE BALDWIN - WRF VERSION
 !     21-07-26  Wen Meng  - Restrict compuation from undefined grids
 !     21-09-13  Jesse Meng- 2D DECOMPOSITION
+!     26-03-27  Alyson Stahl - Remove shared DO termination labels
 !
 !-----------------------------------------------------------------------
 !     ROUTINE TO COMPUTE WET BULB TEMPERATURES USING THE LOOK UP TABLE
@@ -62,9 +63,9 @@ SUBROUTINE WETBULB(T,Q,PMID,HTM,KARR,TWET)
 !!$omp&         presk,qbtk,qqk,sqk,sqs00k,sqs10k,tbtk,thesp,tpspk,
 !!$omp&         tqk,tthbtk,tthk)
 !-----------------------------------------------------------------------
-      DO 300 L=1,LM
-      DO 125 J=JSTA,JEND
-      DO 125 I=ISTA,IEND
+      DO L=1,LM
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
         IF (HTM(I,J,L)<1.0) THEN
           THESP(I,J)=273.15
           cycle    
@@ -128,13 +129,15 @@ SUBROUTINE WETBULB(T,Q,PMID,HTM,KARR,TWET)
         THESP(I,J)=spval
         ENDIF !end t(i,j,l)<spval
 !      ENDIF
-  125 CONTINUE
+      ENDDO
+      ENDDO
+
 !--------------SCALING PRESSURE & TT TABLE INDEX------------------------
       KNUML=0
       KNUMH=0
 !
-      DO 280 J=JSTA,JEND
-      DO 280 I=ISTA,IEND
+      DO J=JSTA,JEND
+      DO I=ISTA,IEND
       KLRES(I,J)=0
       KHRES(I,J)=0
 !
@@ -150,7 +153,8 @@ SUBROUTINE WETBULB(T,Q,PMID,HTM,KARR,TWET)
           KHRES(I,J)=1
         ENDIF
 !      ENDIF
- 280  CONTINUE
+      ENDDO
+      ENDDO
 !***
 !***  COMPUTE PARCEL TEMPERATURE ALONG MOIST ADIABAT FOR PRESSURE<PL
 !**
@@ -169,6 +173,6 @@ SUBROUTINE WETBULB(T,Q,PMID,HTM,KARR,TWET)
       ENDIF
 !-----------------------------------------------------------------------
 !-----------------------------------------------------------------------
-  300 CONTINUE
+      ENDDO
       RETURN
       END