171 260 DELTI=DELTI-VI 



172 270 V1=1.E-18 



173 C** SET UP DETERMINANT FOR PHASE VELOCITY V + VI 



174 280 CALL SETUP 



175 CALL DETNT (N.DET.DETI) 

 173 IF (K6.NE.1) GO TO 300 



177 PRINT 1330, V, VI, DET, DETI, SLR, SLI 



178 300 TEMNR = DELTA 



179 TEMNI = DELTI 



180 TEMDR=VEL-DET 



181 TEMDI=VELI-DETI 



182 TEIV1DEN = TEMDR*TEMDR+TEMDI*TEMDI 



183 IF (TEMDEN.EO.O. ) GO TO 320 



184 TEmr;;,o-T[;mNR + TEIVIDR+TEIVINI*TEMDI 



185 TEMINU=TEMNI*TEMDR-TEMNR*TEMDI 



186 SLR = TEMRNU / TEMDEN 



187 SLI = TEMINU / TEMDEN 



188 IF (J .GT. 3) GO TO 125 



189 * SR(4-NUMBER) = SLR 



190 SI{4-NUMBER) = SLI 



191 125 DELTA = DET * SLR - DETI * SLI 



192 DELTI = DET * SLI + DETI * SLR 



193 SIZE=DELTA*DELTA+DELTI*DELTI 



194 C** DISCONTINUE ITERATION AFTER 2ND STEP IF CORRECTION STEP IS MORE THAN 



195 C** PREVIOUS STEP. 



19C IF ( (SIZE.GT.SIZE2) .AND. (d.GT.2)) GO TO 320 



197 SIZE2=SIZE*2. 



198 VEL=DET 



199 VELI=DETI 



200 310 CONTINUE 



201 320 CONTINUE 



202 51 PV(4-NUMBER) = V 



203 PVI(4-NUMBER)= VI 



204 NUMBER = NUMBER + 1 



205 GO TO 107 



206 C** START STANDARD STEP, EXTRAPOLATE PHASE VELOCITY AND SLOPE. 



207 111 INCA = DP 



208 INCB = 



209 INCC = 



210 122 INCD = 



211 T(1 ) = 



212 T(2) = 



213 T(3) = 



214 DO 112 



215 W(IS+4) = -SR(IS) / T(IS) 



216 WI(IS + 4) = -SKIS) / T(IS) 



217 W(IS) = -PV(IS) / T(IS) 



218 112 WI(IS) = -PVI(IS) / T(IS) 



219 113 INCD = INCA + INCB 



220 INCE = INCD + INCC 



221 T(1 ) = INCD * INCE 



222 T(2) = INCA * INCE 



223 T(3) = INCA * INCD 



224 SLOP = 0. 



225 SLOPI = 0. 



226 SUM = 0. 



227 SUMI = 0. 



DP 





DP 





-INCB • 



- INCC 



INCB * 



INCD 



INCB * 



INCC 



INCD * 



INCC 



IS = 1 



,3 



97 



