function CL(ALPHA,BETA)
    A=[   0.     0.     0.     0.     0.     0.     0.     0.     0.     0.     0.     0.
       -.001  -.004  -.008  -.012  -.016  -.019  -.020  -.020  -.015  -.008  -.013  -.015
       -.003  -.009  -.017  -.024  -.030  -.034  -.040  -.037  -.016  -.002  -.010  -.019
       -.001  -.010  -.020  -.030  -.039  -.044  -.050  -.049  -.023  -.006  -.014  -.027
        .000  -.010  -.022  -.034  -.047  -.046  -.059  -.061  -.033  -.036  -.035  -.035
        .007  -.010  -.023  -.034  -.049  -.046  -.068  -.071  -.060  -.058  -.062  -.059
        .009  -.011  -.023  -.037  -.050  -.047  -.074  -.079  -.091  -.076  -.077  -.076]'
      #
      S= 0.2 * ALPHA
      K= floor(S)+3
      if(K <= 1)
          K= 2
      end
      if(K >=  12)
          K=  11
      end
      DA= S - K +3
      L = K + fix( 1.1*sign(DA) )
      #
      S= .2* abs(BETA)
      M= floor(S)+1
      if(M <=  1)
          M= 2
      end
      if(M >=  7)
          M= 6
      end
      DB= S - M +1
      N= M + fix( 1.1*sign(DB) )
      #
      K=convert(Int64, K)
      L=convert(Int64, L)
      M=convert(Int64, M)
      N=convert(Int64, N)
      T= A[K,M]
      U= A[K,N]
      V= T + abs(DA) * (A[L,M] - T)
      W= U + abs(DA) * (A[L,N] - U)
      DUM= V + (W-V)  * abs(DB)
      return clv = DUM * sign(BETA)
end
