C ----------------------------------------------------------------------
      SUBROUTINE SBNONU(AB,BB)
      PARAMETER (N=32)
      COMMON /DX/DX /LSTM/M1(N),M2(N),M3(N) /LSTP/L1(N),L2(N),L3(N)
      COMMON /UF/UF(N,N,N) /BU/BU(N,N,N)
      COMMON /U/U(N,N,N) /W1/CX(N,N,N)
      COMMON /V/V(N,N,N) /W2/CY(N,N,N)
      COMMON /W/W(N,N,N) /W3/CZ(N,N,N)
      TX=24.*DX*16.**2

      DO 10 K=1,N
      DO 10 J=1,N
      DO 10 I=1,N
      CX(I,J,K)=
     &  -(-U(M2(I),J,K)+ 9.*U(M1(I),J,K)+ 9.*U(I,J,K)-U(L1(I),J,K))
     &  *( U(M2(I),J,K)-27.*U(M1(I),J,K)+27.*U(I,J,K)-U(L1(I),J,K))
      CY(I,J,K)=
     &  -(-V(M1(I),J,K)+ 9.*V(I,J,K)+ 9.*V(L1(I),J,K)-V(L2(I),J,K))
     &  *( U(I,M1(J),K)-27.*U(I,J,K)+27.*U(I,L1(J),K)-U(I,L2(J),K))
      CZ(I,J,K)=
     &  -(-W(M1(I),J,K)+ 9.*W(I,J,K)+ 9.*W(L1(I),J,K)-W(L2(I),J,K))
     &  *( U(I,J,M1(K))-27.*U(I,J,K)+27.*U(I,J,L1(K))-U(I,J,L2(K)))
   10 CONTINUE

      DO 20 K=1,N
      DO 20 J=1,N
      DO 20 I=1,N
      F=(-CX(M1(I),J,K)+9.*CX(I,J,K)+9.*CX(L1(I),J,K)-CX(L2(I),J,K)
     &   -CY(I,M2(J),K)+9.*CY(I,M1(J),K)+9.*CY(I,J,K)-CY(I,L1(J),K)
     &   -CZ(I,J,M2(K))+9.*CZ(I,J,M1(K))+9.*CZ(I,J,K)-CZ(I,J,L1(K))
     &   )/TX
      UF(I,J,K)=UF(I,J,K)+AB*F+BB*BU(I,J,K)
      BU(I,J,K)=F
   20 CONTINUE
      RETURN
      END
C ----------------------------------------------------------------------
      SUBROUTINE SBNONV(AB,BB)
      PARAMETER (N=32)
      COMMON /DX/DX /LSTM/M1(N),M2(N),M3(N) /LSTP/L1(N),L2(N),L3(N)
      COMMON /VF/VF(N,N,N) /BV/BV(N,N,N)
      COMMON /U/U(N,N,N) /W1/CX(N,N,N)
      COMMON /V/V(N,N,N) /W2/CY(N,N,N)
      COMMON /W/W(N,N,N) /W3/CZ(N,N,N)
      TX=24.*DX*16.**2

      DO 10 K=1,N
      DO 10 J=1,N
      DO 10 I=1,N
      CX(I,J,K)=
     &  -(-U(I,M1(J),K)+ 9.*U(I,J,K)+ 9.*U(I,L1(J),K)-U(I,L2(J),K))
     &  *( V(M1(I),J,K)-27.*V(I,J,K)+27.*V(L1(I),J,K)-V(L2(I),J,K))
      CY(I,J,K)=
     &  -(-V(I,M2(J),K)+ 9.*V(I,M1(J),K)+ 9.*V(I,J,K)-V(I,L1(J),K))
     &  *( V(I,M2(J),K)-27.*V(I,M1(J),K)+27.*V(I,J,K)-V(I,L1(J),K))
      CZ(I,J,K)=
     &  -(-W(I,M1(J),K)+ 9.*W(I,J,K)+ 9.*W(I,L1(J),K)-W(I,L2(J),K))
     &  *( V(I,J,M1(K))-27.*V(I,J,K)+27.*V(I,J,L1(K))-V(I,J,L2(K)))
   10 CONTINUE

      DO 20 K=1,N
      DO 20 J=1,N
      DO 20 I=1,N
      F=(-CX(M2(I),J,K)+9.*CX(M1(I),J,K)+9.*CX(I,J,K)-CX(L1(I),J,K)
     &   -CY(I,M1(J),K)+9.*CY(I,J,K)+9.*CY(I,L1(J),K)-CY(I,L2(J),K)
     &   -CZ(I,J,M2(K))+9.*CZ(I,J,M1(K))+9.*CZ(I,J,K)-CZ(I,J,L1(K))
     &   )/TX
      VF(I,J,K)=VF(I,J,K)+AB*F+BB*BV(I,J,K)
      BV(I,J,K)=F
   20 CONTINUE
      RETURN
      END
C ----------------------------------------------------------------------
      SUBROUTINE SBNONW(AB,BB)
      PARAMETER (N=32)
      COMMON /DX/DX /LSTM/M1(N),M2(N),M3(N) /LSTP/L1(N),L2(N),L3(N)
      COMMON /WF/WF(N,N,N) /BW/BW(N,N,N)
      COMMON /U/U(N,N,N) /W1/CX(N,N,N)
      COMMON /V/V(N,N,N) /W2/CY(N,N,N)
      COMMON /W/W(N,N,N) /W3/CZ(N,N,N)
      TX=24.*DX*16.**2

      DO 10 K=1,N
      DO 10 J=1,N
      DO 10 I=1,N
      CX(I,J,K)=
     &  -(-U(I,J,M1(K))+ 9.*U(I,J,K)+ 9.*U(I,J,L1(K))-U(I,J,L2(K)))
     &  *( W(M1(I),J,K)-27.*W(I,J,K)+27.*W(L1(I),J,K)-W(L2(I),J,K))
      CY(I,J,K)=
     &  -(-V(I,J,M1(K))+ 9.*V(I,J,K)+ 9.*V(I,J,L1(K))-V(I,J,L2(K)))
     &  *( W(I,M1(J),K)-27.*W(I,J,K)+27.*W(I,L1(J),K)-W(I,L2(J),K))
      CZ(I,J,K)=
     &  -(-W(I,J,M2(K))+ 9.*W(I,J,M1(K))+ 9.*W(I,J,K)-W(I,J,L1(K)))
     &  *( W(I,J,M2(K))-27.*W(I,J,M1(K))+27.*W(I,J,K)-W(I,J,L1(K)))
   10 CONTINUE

      DO 20 K=1,N
      DO 20 J=1,N
      DO 20 I=1,N
      F=(-CX(M2(I),J,K)+9.*CX(M1(I),J,K)+9.*CX(I,J,K)-CX(L1(I),J,K)
     &   -CY(I,M2(J),K)+9.*CY(I,M1(J),K)+9.*CY(I,J,K)-CY(I,L1(J),K)
     &   -CZ(I,J,M1(K))+9.*CZ(I,J,K)+9.*CZ(I,J,L1(K))-CZ(I,J,L2(K))
     &   )/TX
      WF(I,J,K)=WF(I,J,K)+AB*F+BB*BW(I,J,K)
      BW(I,J,K)=F
   20 CONTINUE
      RETURN
      END