1 SUBROUTINE SETUP 



2 IMPLICIT DOUBLE PRECISION (A-H,0-Z) 



3 DOUBLE PRECISION LAMBDA , LAMBDI 



4 COMMON /HAN/ H2R , H2 I , HI R , HI I , H2PR , H2PI ,H1 PR , HI PI , EXPONT 



5 COMMON /EXPO/ EXSUM, CNTR, RATI0(25) 



6 COMMON/DETMNT/ A ( 25 , 4 ) ,0( 25 , 4 ) 



7 COMMON/INPUT/ Z{12), N, OMEGA, V, VI, C0N(12), GS0(12), 



8 1 CAY(12), LAMBDA, LAMBDI, G(12) 



9 2,RH0(12), GI(12), G S0I(12), CAYI(12) 



10 COMMON /LIMIT/ TLIM, EXPON , SLIM 



11 COMMON/PARTS/ ZT ( 1 2 ) , ZTI ( 1 2 ) , ZB( 1 2 ) , ZBI ( 1 2 ) 



12 DENOM = V * V + VI * VI 



13 LAMBDA = OMEGA * V / DENOM 



14 LAMBDI ^ -OMEGA * VI/ DENOM 



15 M = N - 1 



16 DO 10 I = 0,M 



17 IF ( I .EO. 0) GO TO 35 



18 IF (ZR .GT. -7.4) GO TO 25 



19 IF (G( I ) .LT. 0. ) GO TO 25 



20 ZE = G(I) * (Z(I+1) - Z(I)) + ZE 



21 IF (ZE . LT. -7.5) ZE = -7.5 



22 GO TO 26 



23 25 CONTINUE 



24 ZE = G(I) * (Z(I+1) - Z(I)) + ZR 



25 IF (ZE . LT. SLIM) ZE = SLIM 



26 26 CONTINUE 



27 ZO = GI(I) * (Z(I+1) - Z(I)) + ZI 



28 30 ZB(I) = ZE 



29 ZBI(I) = ZO 



30 CALL HANKEL(ZE,ZO,0) 



31 ZB(I) = ZE 



32 ZBI(I) - ZO 



33 RATI0(2*I) = EXPONT 



34 A(2*I , 1 ) = H2R * RHO(I) 



35 0(2*1 , 1 ) = H2I * RHO(I) 



36 A(2*I ,2) = H1R * RHG(I) 



37 0(2*1 ,2) = HI I * RHO(I) 



38 A(2*I+1,1) = H2PR * G( I ) - H2PI * GI(I) 



39 0(2+1+1,1) = H2PI * G(I) + H2PR * GI ( I ) 



40 A(2*I+1,2) = H1PR * G(I) - H1PI ♦ GI(I) 



41 Q(2*l+1,2) = H1PI * G(I) + H1PR * GI( I) 



42 35 CONTINUE 



43 GSABS = G S0(I+1)**2 + G S0I(I+1)**2 



44 XI = CAY( 1+1 ) - LAMBDA 



45 X2 = CAY(I+1) + LAMBDA 



46 X3 =CAYI(I+1) - LAMBDI 



47 X4 =CAYI(I+1) + LAMBDI 



48 X = X1 * X2 - X3 * X4 



49 Y = XI * X4 + X3 * X2 



50 ZT(I + 1) = (X ♦ G S0(I + 1) + Y * G SQId + D) / GSABS 



51 2TI(I + 1) = (Y ♦ G S0(I + 1) - X ♦ G SQId + l)) / GSABS 



52 ZR = ZT(I+1) 



53 21 = ZTI(I+1 ) 

 54 



55 ZE = ZR 



56 ZQ = ZI 



59 



