171 XRE = (T * TEMP + TI * CI(I)) / TEM 



172 XIM = -(T * CI(I) - TI * TEMP) / TEM 



173 TEM = XRE**2 + XIM**2 



174 CB(I) = SORT((XRE + SC5RT(TEM)) ♦ ,5) 



175 CBI(I) = .5 * XIM / CB(I) 



176 GAMMAK I ) = 0. 



177 GO TO 33 



178 29 TEMP = C(I) - 2. ♦ GAMMA( I ) * (2(1+1) - Z( I ) ) 



179 CB( I) = SORT (T / TEMP) 



180 GAMMAI(I) = .5 * (CI(I) - TI / CB(I)**2) / (Z(I+1) - 2(1)) 



181 GO TO 33 



182 27 N = N - 1 



183 33 CONTINUE 



184 C COMPUTE USEFULL QUANTITIES 



185 PRINT 58 



186 58 FORMAT (7X,6H RE M ,8X,6H IM M ,9X,5H L/KM,8X,6H RE C ,8X, 



187 * 6K IM C ,5X,12H RE C BOT TOM , 4X , 1 2H IM C BOTTOM, 1 OX , 9H GRADIENT ) 



188 OMEGA = 6.283185307DO ♦ FREO 



189 DO 40 I = 1 ,N 



190 TEMP = C(I)**2 + CI(I)**2 



191 CAY(I) ^ OP/IEGA * C(I) / TEMP 



192 CAYI(I) = -OMEGA * CI(I) / TEMP 



193 CAY SQ(I) = CAY(I)**2 - CAYI(I)**2 



194 CAY SOI(I) = 2. DO * CAY(I) * CAYI(I) 



195 TEMDR = -2. * (GAMMA(I) * CAY SO(I) " GAMMAI(I) * CAY SOI(I)) 



196 TEMDI = -2. * (GAMMA(I) * CAY SOI(I) + GAMMAI(I) » CAY S(3(I)) 



197 G CU(I) = (TEMDR * C(I) + TEMDI ♦ CI(I)) / TEMP 



198 G CUI(I) = (TEMDI * C(I) - TEMDR ♦ CI(I)) / TEMP 



199 TEM1 = DCBRT(-DSQRT( GAMMA(I)**2 + GAMMAI ( I ) **2 ) • 2.*0MEGA**2) 



200 TEM1I = DATAN ( ABS ( GAMMAI ( I ) / GAMMA(I)))/ 3. 



201 CRTG = TEM1 * DCOS(TEMII) 



202 CRTGI = TEM1 * DSIN(TEMII) 



203 IF (GAMMA(I) .LT. 0.) CRTG = -CRTG 



204 IF (GAMMAId ) . LT. 0.) CRTGI = -CRTGI 



205 G(I) = (C(I) * CRTG + CI(I) * CRTGI) / TEMP 



206 GI(I) = (C(I) * CRTGI - CI(I) * CRTG) / TEMP 



207 CON(I) = G(I) * C(I) - GI(I) * CI(I) 



208 CON(I) = 0MEGA**2 / C0N(1)**2 



209 XMI = -GI(I) ♦ (Z(I+1) - Z(I)) 



210 XM = -G(I) * (Z(I + 1 ) - Z(I)) 



211 DPK(I) = -8686. DO * CAYI(I) 



212 PRINT 30, XM, XMI, DPK(I), C(I), CI(I), CB(I), CBI(I) 



213 * ,GAMMA(I). GAMMAI (I) 



214 G SQI(I) = 2. * G( I) * GI(I) 



215 40 G SO(I) = G(I)**2 - GI(I)**2 



216 C FIND MODES 



217 NXTRA=0 



218 IJ FLAG=0 



219 NN = NN + 1 



220 IF (K5 .EO. 1 ) GO TO 15 



221 DO 50 NN = 1 ,350 



222 15 IF (IJ FLAG .EQ. 1) GO TO 53 



223 52 IF (NXTRA .GT. 0) GO TO 44 



224 READ 60, V , VI , STEP , STEPI , NXTRA 



225 60 FORMAT (4D10.4,I10) 



226 IF (NXTRA .GE.O) GO TO 62 



227 V = V + VI * 1 .0-10 



47 



