25 IF(AN0RM1.GT.1.)AN0RM=1./FL0AT(NP) 

 38 IF(CUT(1).NE.-1.0) 60 TO 29 



FUNDAMENTAL FREQUENCY,H(1)=.2,N(1)=3 



CUT(1)=.8/(FL0AT(NP)) 



H(l)=.2 



N(l)=3 



IF(ANORMl.GT.l.) AN0RM=1./(FL0AT(NP)-(2*(N(1)-1)) ) 

 29 IF(IADJ.NE.l) GO TO 383 



COMPUTE INDEX OF FIRST NON-ZERO POINT 

 28 DO 385 J=1,NP 



K1=J 



IF(DATA(J).NE.O.O) GO TO 386 



385 CONTINUE 



WRITE(I0UT,82) XORG,YORG,ANORM,ITG 

 82 F0RMAT(15H PLOT ORIGIN X=,E11,4,4H Y= ,E11.4,15HN0RMALIZATI0N =, 

 1E11.4,30A1) 

 COMPUTE INDEX OF LAST NON-ZERO POINT 



386 DO 391 J=NP,K1,-1 

 K2=J 

 IF(DATA(J).NE.O.O) GO TO 392 



391 CONTINUE 



SHIFT DATA TO LEFT AND RECOMPUTE NO. OF INPUT POINTS 



392 NP=K2-K1+1 



DO 393 J=1,NP 

 K=K1+J-1 



393 DATA(J)=DATA(K) 



DEMEAN DATA AND APPLY COSINE TAPER 

 383 SUM=0.0 

 XNP=NP 

 DO 388 J=1,NP 



388 SUM=SUM+DATA(J) 

 AVE=SUM/XNP 



IF(N(l).GE.O.AND.IMEAN.EQ.O) GO TO 381 

 DO 389 J=1,NP 



389 DATA(J)=DATA(J)-AVE 



CALL BATHXS(DATA,NP,GRIDKM) 



IF(N(l).GE.O) 60 TO 381 



DO 382 J=1,NP 



AM=J-1 



TNP=NP-1 



TNP=TNP/2.0 



AM=3.14159*(AM-TNP)/TNP 

 382 DATA(J)=0ATA(J)*0.5*(1.0+C0S(AM)) 



RAISE NO. OF PTS TO A POWER OF 2(M)AN0 STORE IN INP 

 381 K0UNT=1 



DO 5 M=l,12 



I=2**M 



IF(I-NP)5,6,7 



5 CONTINUE 



6 INP=I 

 GO TO 9 



7 INP=I 

 NPP=NP+1 



187 



