************************************************************************
      PROGRAM RAYTRACING
************************************************************************
      IMPLICIT REAL*8(A-H,O-Z),INTEGER*4(I-N)      
c10================================================================ 
c10     RAY TRACING PROGRAM FOR VLF WAVE (X-MODE) IN DE-MODEL
c10     Dipole model is selected in REF by Call dipole   
c10     IGRF model is selected in REF by Call IGRF
c10     kizami H is shown in Adam_update 
c10     Result of ray tracing (OUT1) is written in fort.21 
c10     Instantaneous H is written in fort.25                      
c10========  initial conditions ==================================

c	fq00: wave frequency in kHz
c	rr00: starting radius distance from the earth center in km
c	th00: geomagnetic colatitude in deg
c	ph00: geomagnetic longitude in deg
c	es00: angle between the initial wave normal and geomagnetic 
c            meiridian plane, in deg (Epsilon)          
c	dl00: angle between the projected direction of the initial 
c            wave normal onto the geomagnetic meridial plane and 
c            the radial direction, in deg (Delta)

      rr00 = 6370.0+91.0d0
      ph00 = 0.0d0
      ameg = 1000.0d0
      fq00 = 0.003d0*ameg
      dl00 = 0.0d0
      es00 = 0.0d0
      print*,'VLF_DE_dipole_F.f',', rr00=', rr00,', fq00=',fq00, 
     &', ph00=', ph00,', es00=', es00,', dl00=',dl00

      DO 100 I=1,5
      WRITE(21,*)
      write(25,*)
      th00 = dble(i-1)*1.0d0+40.0d0
      print*, 'th00=', th00
      CALL INIT(fq00,rr00,th00,ph00,es00,dl00)
      CALL OUT0
      CALL ADAMS
  100 CONTINUE

      STOP
      END