PROGRAM REZEST_FDC
      IMPLICIT NONE
      INCLUDE "lun.par"
      REAL*4 B, EPSILON, N_RL, N_RL_FRONT, LAMBDA, L, N_M, N_M_POLAR, M,
     X     P, R_M, STEREO
      REAL*4 DK

      B = 1.74
      epsilon = 0.00015
      STEREO = 6.0*3.14159/180.0
      n_rl = 0.0317
      n_rl_front = 0.0212
      lambda = 0*3.14159/180.0
      L = 0.58600 - 0.16895
      R_M = (0.58600 + 0.16895)/2.0
      n_m = 23
      n_m_polar = 10

      m = 0.139

      open(UNIT=LUNOUT, FILE='cdc_pi.txt', STATUS='UNKNOWN')
      open (UNIT=LUNOUTA, FILE='cdc_pi_angle.txt', STATUS='UNKNOWN')
      open (UNIT=LUNOUTB, FILE='cdc_pi_polar.txt', STATUS='UNKNOWN')
      P = 0.01
      DO WHILE (p < 2.0)
C         CALL calc_rez
         call calc_rez(LAMBDA, L, EPSILON, N_M, P, M, N_RL, B, DK)
C         CALL calc_angle
         call calc_azimuth_angle(LAMBDA, L, EPSILON, N_M, P, M, N_RL,
     X        N_RL_FRONT, B, R_M, DK)
         CALL CALC_POLAR_ANGLE(LAMBDA, L, EPSILON, STEREO, N_M_POLAR, P,
     X        M, N_RL, N_RL_FRONT, B)
         P = P*1.02
      END DO
      close(LUNOUT)
      close(LUNOUTA)
      close(LUNOUTB)

      m = 0.938

      open(UNIT=LUNOUT, FILE='cdc_p.txt', STATUS='UNKNOWN')
      open (UNIT=LUNOUTA, FILE='cdc_p_angle.txt', STATUS='UNKNOWN')
      open (UNIT=LUNOUTB, FILE='cdc_p_polar.txt', STATUS='UNKNOWN')
      P = 0.01
      DO WHILE (p < 2.0)
C         CALL calc_rez
         call calc_rez(LAMBDA, L, EPSILON, N_M, P, M, N_RL, B, DK)
C         CALL calc_angle
         call calc_azimuth_angle(LAMBDA, L, EPSILON, N_M, P, M, N_RL,
     X        N_RL_FRONT, B, R_M, DK)
         CALL CALC_POLAR_ANGLE(LAMBDA, L, EPSILON, STEREO, N_M_POLAR, P,
     X        M, N_RL, N_RL_FRONT, B)
         P = P*1.02
      END DO
      close(LUNOUT)
      close(LUNOUTA)
      close(LUNOUTB)

      STOP
      END