SUBROUTINE VTRI0 ( A,B,C,F,Y,K,irun)

      _RL A(irun,K),B(irun,K),C(irun,K),Y(irun,K+1)

        A(I,L) = 1. / ( A(I,L) - B(I,LM1) * C(I,L) )

       Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L)

      SUBROUTINE VTRI1 ( A,B,Y,K,irun)

      _RL A(irun,K),B(irun,K),Y(irun,K+1)

        Y(I,L) = B(I,L) * Y(I,L+1) * A(I,L)

      SUBROUTINE VTRI2 ( A,B,C,F,Y,K,irun)

      _RL A(irun,K),B(irun,K),C(irun,K),F(irun,K)

        Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) *  A(I,L)

        B(I,J) = A(I+OFFSET-1,J)

        B(I,J) = A(I+OFFSET-1,J)

        B(Lena+I,J) = A(Lena+OFFSET-1,J)

      SUBROUTINE PASTE(B,A,IB,IA,L,K)

      _RL A(IA,L), B(IB,L)

      A(I+OFFSET-1,J) = B(I,J)

      SUBROUTINE PSTBMP(B,A,IB,IA,L,K)

      _RL A(IA,L), B(IB,L)

      A(I+OFFSET-1,J) = A(I+OFFSET-1,J) + B(I,J)

      SUBROUTINE STRINT(A,B,IA,IB,L,K)

      INTEGER A(IA,L), B(IB,L)

        B(I,J) = A(I+OFFSET-1,J)

        B(I,J) = A(I+OFFSET-1,J)

        B(Lena+I,J) = A(Lena+OFFSET-1,J)

      SUBROUTINE STRIP(A,B,IA,IB,L,K)

      _RL A(IA,L), B(IB,L)

      SUBROUTINE GRD2MSC(A,IM,JM,IGRD,B,MXCHPS,NCHP)

      _RL A(IM*JM), B(MXCHPS)

         B(I) = A(IGRD(I))

      SUBROUTINE MSC2GRD(IGRD,CHFR,B,MXCHPS,NCHP,FRACG,A,IM,JM)

      _RL A(IM*JM), B(MXCHPS), CHFR(MXCHPS), FRACG(IM*JM)

         VT1(IGRD(I)) = VT1(IGRD(I)) + B(I)*CHFR(I)

        B(I,J) = A(I+OFFSET-1,J)

        B(I,J) = A(I+OFFSET-1,J)

        B(Lena+I,J) = A(Lena+OFFSET-1,J)

      SUBROUTINE STRIPINT(A,B,IA,IB,L,K)

      INTEGER A(IA,L), B(IB,L)

         B = 4. _d 0*d2p1 - d2

         dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

         B = 4. _d 0*d2 - d2m1

         dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

      _RL d2,d2p1,d2m1,A,B,C,D

         B = 4. _d 0*d2p1 - d2

         dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

         B = 4. _d 0*d2 - d2m1

         dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

      _RL d2,d2p1,d2m1,A,B,C,D

         B = 4. _d 0*d2p1 - d2

         dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

         B = 4. _d 0*d2 - d2m1

         dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)

      _RL d2,d2p1,d2m1,A,B,C,D

      B  = (5.03217E-5*P+3.359406)*P+21582.27

      x = x*(1.0 - P/B) + (V350P+x)*P*SR/(B*(B+SR))

      SR = V350P*(1.0 - P/B)

      SBO_RHO = 1028.106331 + P/B/SR - x / (SR*(SR+x))

      _RL PLAT,D,C1,P,PR,Q,X,SR,V350P,B

        B(I,J)=QS1*6.11 _d +00*EXP(17.2694 _d +00

        A2(I,J)=-D1*UG(I,J)*tsurfLoc(I,J)-D1W*UG(I,J)*B(I,J)

      _RL B          (1:sNx,1:sNy)

        B(I,J)=QS1*6.11 _d +00*EXP(17.2694 _d +00

        A2(I,J)=-D1*UG(I,J)*tsurfLoc(I,J)-D1W*UG(I,J)*B(I,J)

      _RL B          (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

      REAL    A, B, E, F, Q, DA, DB, XLAT

         B    = 0.0

            B    = B+DB

            B    = B+DB

            B    = B+DB

            HS(M+M+1) = B