PROGRAM RAY TO WAV
CALCULATES REFLECTION COEFFICIENTS FROM RAY TRACING IN COMPLEX SPACE
      COMPLEX I,TEE,TEM,TME,TMM,TOO,TOX,TXO,TXX 
      DATA (ZERO=0.)
      DATA(I=0.,1.),(PI=3.141592654)
      DATA (C=.2997925) 
      DATA(HITE=50.)
      COMPLEX BETA,AZI,H,X1,X2,TERAYV,TERAYH,RERAYV,RERAYH,EL,AZ,EV,EH
      REAL LN10 
      DATA (LN10=2.3025851),(EV=2(0.)),(EH=2(0.)) 
      DIMENSION BUFF(20),W(26)
      DIMENSION TOO(4)
      EQUIVALENCE (TOX,TOO(2)),(TXO,TOO(3)),(TXX,TOO(4))
      EQUIVALENCE (ID1,W(1)),(ID2,W(3)),   (F,W(4)),(AZI,W(5))
     1,(BETA,W(7)),(DELTA,W(9)),(H,W(10)),(X1,W(12)),(X2,W(14)) 
     2,(AMP,W(16)),(PHASE,W(17)),(TIME,W(18)),(TERAYV,W(19))
     3,(TERAYH,W(21)),(RERAYV,W(23)),(RERAYH,W(25)),(HI,W(11))
     4,(NGRONK,W(2))
      CALL Q9EXUN 
  1   BUFFER IN(60,1)(BUFF(1),BUFF(20)) 
  6   IF(UNIT,60) 6,7,10,11 
  7   CONTINUE
      DECODE (144,2,BUFF) W 
  2   FORMAT(A1,R1,A2,23A6) 
      IF(F.EQ.F0) GO TO 15
      IF(CABS(TXO).EQ.0.) TXO= TOX
      TEE= .5*(  -TOO  -TOX  -TXO  -TXX)
      TEM= .5*(-I*TOO+I*TOX-I*TXO+I*TXX)
      TME= .5*(-I*TOO-I*TOX+I*TXO+I*TXX)
      TMM= .5*(   TOO  -TOX  -TXO  +TXX)
      TEEA=20.*ALOG10(CABS(TEE))
      TEMA=20.*ALOG10(CABS(TEM))
      TMEA=20.*ALOG10(CABS(TME))
      TMMA=20.*ALOG10(CABS(TMM))
      TEEP=180./PI*CANG(TEE)
      TEMP=180./PI*CANG(TEM)
      TMEP=180./PI*CANG(TME)
      TMMP=180./PI*CANG(TMM)
      TEEP=TEEP+DPHASE
      TEMP=TEMP+DPHASE
      TMEP=TMEP+DPHASE
      TMMP=TMMP+DPHASE
      TEEP=FMOD(-180.,TEEP, 180.) 
      TEMP=FMOD(-180.,TEMP, 180.) 
      TMEP=FMOD(-180.,TMEP, 180.) 
      TMMP=FMOD(-180.,TMMP, 180.) 
      PRINT 16,        TEEA,TEEP,TEMA,TEMP,TMEA,TMEP,TMMA,TMMP
  16  FORMAT(//1X,8F12.4) 
      TOO=TOX=TXO=TXX=0.
  15  CONTINUE
      DPHASE=2.*HITE*360.    *F/C 
      EL=180./PI*BETA 
      AZ=180./PI*AZI
      DEL=180./PI*DELTA 
      NPUNTYP=NGRONK/8
      TOO(NPUNTYP)=       CEXP(CMPLX(LN10/20.*AMP,PI/180.*PHASE)) 
C     TAMPV=20.*ALOG10(CABS(TERAYV))
C     RAMPV=20.*ALOG10(CABS(RERAYV))
C     RAMPH=20.*ALOG10(CABS(RERAYH))
C     TAMPH=20.*ALOG10(CABS(TERAYH))
C     TPHASEV=180./PI*CANG(TERAYV)
C     RPHASEV=180./PI*CANG(RERAYV)
C     RPHASEH=180./PI*CANG(RERAYH)
C     TPHASEH=180./PI*CANG(TERAYH)
      IF(F.NE.F0) PRINT 5 
  5   FORMAT(1H1,*IID*,3X,*TYPE*,5X,*F*,3X,*AZIMUTH*,6X,*BETA*,5X,
     1 *DELTA*,5X,*HIMAG*,5X,*XREAL*, 5X,*XIMAG*,7X,*AMP*,5X,*PHASE*, 
     2 6X,*TIME*) 
      PRINT 3,ID1,ID2,NPUNTYP,F,AZ,EL,DEL,HI,X1,AMP,PHASE,TIME
  3   FORMAT(1X,A1,A5,I2,5F10.4,C(F10.4,F10.4),3F10.4)
      F0=F
      GO TO 1 
  10  CALL EXIT 
  11  PRINT 12
  12  FORMAT(* PARITY ERROR IN READING*)
      GO TO 7 
      END 
      FUNCTION MOD(A,X,B) 
      N=(X-A)/(B-A) 
      IF(X.LT.A) N=N-1
      FMOD=X-N*(B-A)
      RETURN
      END