PROGRAM TEST

USE HELBIG_ANGSIZ_OLD

REAL :: LAMBDA, OMEGA, ETA=1.0, ZMAX, Z1, Z2, D

LOGICAL ::VARETA=.FALSE., DEBUG=.TRUE., MAXZ, BOUNCE

INTEGER :: WMTYPE, ERROR, I, J

REAL, DIMENSION(18) :: A_LAMBDA = &
(/-5.0,-4.0,-3.0,-2.0,0.0,0.1,0.0,0.0,0.0,0.2,0.3,1.0,0.7,0.5,0.1,2.0,3.0,2.0/)

REAL, DIMENSION(18) :: A_OMEGA = &
(/0.0,1.0,4.0,5.0,6.0,7.0,0.0,0.5,1.0,0.0,0.4,0.0,0.3,1.0,2.52777,0.5,1.0,0.0/)

REAL, DIMENSION(10) :: A1_Z=(/0.0,1.0,0.0,2.0,0.5,2.0,0.7,0.7,3.0,1.5/)

REAL, DIMENSION(10) :: A2_Z=(/0.0,0.1,0.0,0.5,0.15,0.3,0.2,0.2,0.3,0.2/)

LOGICAL :: RISK

COMMON /ARTHUR/ RISK

RISK = .TRUE.

DO, I = 1,18,1

   LAMBDA = A_LAMBDA(I)

   OMEGA  = A_OMEGA(I)

   CALL INICOS(LAMBDA,OMEGA,ETA,VARETA,DEBUG,WMTYPE,MAXZ,ZMAX,BOUNCE,ERROR)

   IF (.NOT. MAXZ) ZMAX = 9999999.0

   PRINT*, "WMTYPE, ZMAX: ", WMTYPE, ZMAX

   IF (WMTYPE > 16) A1_Z = A2_Z

   DO, J = 1, 10, 2

      Z1 = A1_Z(J)

      Z2 = A1_Z(J+1)

      CALL ANGSIZ(Z1,Z2,D,ERROR)

      PRINT*, "D: ", D

   END DO

END DO

END PROGRAM TEST