342 YI = -Z32R 



343 CALL CFR (YR, YI, G2R, G2 I , D4 , D5 , ML) 



344 ST0R11 = K2*(ZRT4R*G2R-ZRT4I*G2I) 



345 ST0R12 = K2* ( ZRT4R*G2I+ZRT4I *G2R ) 



346 H22PR=STOR11*STHR+STOR12*CTHR 



347 H22PI=STOR1 1 *(-CTHR) +5T0R12+STHR 



348 390 IF (FLOUAD .LT. 0) GO TO 410 



349 400 H2PR = H22PR 



350 H2PI = H22PI 



351 GO TO 414 



352 410 YR = -Z32I 



353 YI = Z32R 



354 CALL CFR (YR, YI. G1R, Gil, D4 , D5 , ML) 



355 ST0R13 = K1 * { ZRT4R *G1 R-ZRT4I *G1 I ) 



356 STORM = K1 * ( ZRT4R + G1 I+ZRT4I *G1 R ) 



357 H21PR=ST0R13*(-5TP) -ST0R14+CTP 



358 H21PI=ST0R14*(-STP) +ST0R13*CTP 



359 H2PR = H21PR+H22PR 



360 H2PI = H21PI+H22PI 



361 414 IF (IH .EO. 2) GO TO 999 



362 100 IF (Firs .ME. 1 ) GO TO 420 



363 110 H1PR = GPR+SPR 



364 H1PI = GPI+SPI 



365 GO TO 999 



366 420 IF (FLOUAD .LT. 0) GO TO 440 



367 430 YR = -Z32I 



368 YI = Z32R 



369 CALL CFR (YR, YI, G1R, Gil, D4 , D5 , ML) 



370 440 ST0R15 = K1 * ( ZRT4R *G1 R -ZRT4I*G1I) 



371 ST0R16 = K1*(ZRT4R*G1 I +ZRT4I*G1R) 



372 H11PR = ST0R15*STHR -ST0R16+CTHR 



373 H11PI = ST0R15*CTHR +STQR16+STHR 



374 450 IF (FLOUAD .GT. 0) GO TO 470 



375 460 H1PR = HI 1PR 



376 HI PI = H11PI 



377 GO TO 999 



378 470 ST0R17 = K2* ( ZRT4R*G2R -ZRT4I*G2I) 



379 ST0R18 = K2* ( ZRT4R*G2I +ZRT4I*G2R) 



380 H12PR = ST0R17*(-STP) +ST0R18*CTP 



381 H12PI = ST0R17*(-CTP) -ST0R18*STP 



382 H1PR = H12PR+H11PR 



383 H1PI = H12PI+H11PI 



384 999 CONTINUE 



385 RETURN 



386 END 



•PRT.S J. CFR 



71 



