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

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

        C(I,L) = C(I,L) * A(I,LM1)

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

        F(I,L) = F(I,L) + F(I,LM1) * C(I,L)

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

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

        F(I,L) = F(I,L) + F(I,L-1) * C(I,L)

         C = d2

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

         C = 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

         C = d2

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

         C = 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

         C = d2

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

         C = 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

      REAL    SINA, COSA, C, S

               C = SHC(K+M)

               SHC(K+M) = COSA*C + SINA*S

               SHC(K-M) =-SINA*C + COSA*S

      REAL    C, CM, AK

      C  = SNGL(DSQRT(1.0D0-DBLE(S)*DBLE(S)))

            CM   =    C * CM

      SUBROUTINE R8FFTB1 (N,C,CH,WA,IFAC)

      DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)

         CALL R8ADB4 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3))

  101    CALL R8ADB4 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3))

         CALL R8ADB2 (IDO,L1,C,CH,WA(IW))

  104    CALL R8ADB2 (IDO,L1,CH,C,WA(IW))

         CALL R8ADB3 (IDO,L1,C,CH,WA(IW),WA(IX2))

  107    CALL R8ADB3 (IDO,L1,CH,C,WA(IW),WA(IX2))

         CALL R8ADB5 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3),WA(IX4))

  110    CALL R8ADB5 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3),WA(IX4))

         CALL R8ADBG (IDO,IP,L1,IDL1,C,C,C,CH,CH,WA(IW))

  113    CALL R8ADBG (IDO,IP,L1,IDL1,CH,CH,CH,C,C,WA(IW))

         C(I) = CH(I)

      SUBROUTINE R8FFTF1 (N,C,CH,WA,IFAC)

      DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)

         CALL R8ADF4 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3))

  101    CALL R8ADF4 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3))

         CALL R8ADF2 (IDO,L1,C,CH,WA(IW))

  103    CALL R8ADF2 (IDO,L1,CH,C,WA(IW))

         CALL R8ADF3 (IDO,L1,C,CH,WA(IW),WA(IX2))

  105    CALL R8ADF3 (IDO,L1,CH,C,WA(IW),WA(IX2))

         CALL R8ADF5 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3),WA(IX4))

  107    CALL R8ADF5 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3),WA(IX4))

         CALL R8ADFG (IDO,IP,L1,IDL1,C,C,C,CH,CH,WA(IW))

  109    CALL R8ADFG (IDO,IP,L1,IDL1,CH,CH,CH,C,C,WA(IW))

         C(I) = CH(I)

      SUBROUTINE RFFTB1 (N,C,CH,WA,IFAC)

      DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)

         CALL RADB4 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3))

  101    CALL RADB4 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3))

         CALL RADB2 (IDO,L1,C,CH,WA(IW))

  104    CALL RADB2 (IDO,L1,CH,C,WA(IW))

         CALL RADB3 (IDO,L1,C,CH,WA(IW),WA(IX2))

  107    CALL RADB3 (IDO,L1,CH,C,WA(IW),WA(IX2))

         CALL RADB5 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3),WA(IX4))

  110    CALL RADB5 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3),WA(IX4))

         CALL RADBG (IDO,IP,L1,IDL1,C,C,C,CH,CH,WA(IW))

  113    CALL RADBG (IDO,IP,L1,IDL1,CH,CH,CH,C,C,WA(IW))

         C(I) = CH(I)

      SUBROUTINE RFFTF1 (N,C,CH,WA,IFAC)

      DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)

         CALL RADF4 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3))

  101    CALL RADF4 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3))

         CALL RADF2 (IDO,L1,C,CH,WA(IW))

  103    CALL RADF2 (IDO,L1,CH,C,WA(IW))

         CALL RADF3 (IDO,L1,C,CH,WA(IW),WA(IX2))

  105    CALL RADF3 (IDO,L1,CH,C,WA(IW),WA(IX2))

         CALL RADF5 (IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3),WA(IX4))

  107    CALL RADF5 (IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3),WA(IX4))

         CALL RADFG (IDO,IP,L1,IDL1,C,C,C,CH,CH,WA(IW))

  109    CALL RADFG (IDO,IP,L1,IDL1,CH,CH,CH,C,C,WA(IW))

         C(I) = CH(I)

          C = 0.6 _d 0 *gravity*kappa*kappa*Nw/

          Cs = rho_air*cp*C*

          Cf = htil*2.97 _d 12*C/(T35(no_so)**3)*(

      _RL Cl, Cf, Cs, C