J1=J+N(K) 



J2=JJ-(J-1) 



CDATA(J1)=BDATA(J) 

 52 CDATA(J2)=CDATA(J1) 

 C COMPUTE AMPLITUDE SPECTRUM OF PREWITENED DATA AND PLOT 



GO TO 53 

 79 CONTINUE 

 C STORE HIGH PASS WEIGHTS IN DATA AND AMP.SPECT.OF PREWHITENED 



C DATA IN CDATA 



JJ=2*N(K)+1 



DO 102 J=1,JJ 



102 DATA(J)=CDATA(J) 

 JJ=JJ+1 



DO 104 J=JJ,INP 

 104 DATA(J)=0.0 

 DO 103 J=1,JX 



103 CDATA(J)=BDATA(J) 



C COMPUTE AMPLITUDE SPECTRUM OF PREWH ITEMING FILTER 



K0UNT=3 

 GO TO 53 



236 CONTINUE 



C OUTPUT SMOOTHED SPECTRUM AND PLOT 



FIRST=(N(2)+1)*DELR 

 FIRSTK=(N(2)+1)*DELK 

 WRITE(I0UT,237) FIRSTK 



237 F0RMAT(38H FINAL SMOOTHED AND CORRECTED SPECT. ,16HFIRST FREQUENC 

 1Y=,E8.3) 



CALL PARSVL ( ADATA,NP , INP , AVE , ANO IS ,N ( 2 ) .ANORM ) 

 WRITE(I0UT,491)AVE,AN0IS 



491 FORMATC VALUE OF AMP SPECT OF NOISE = ' ,E10.4, ' RMS WHITE NOISE LE 

 1VEL = ',E10.4) 



WRITE(I0UT,81)(ADATA(J), J=2,IA) 

 4911 CALL P0WWGT(B0,B1,FIRSTK,DELK,ADATA(2),IA-1,0,CT0FF1,CT0FF2) 

 BOSAVE=BO 

 B1SAVE=B1 

 C CALCULATE RMS FOR WHITE NOISE WITH INTERCEPT(BO) 



AVE=BO 



AVE =AVE /ANORM 

 AVE=AVE**2 



ANOISE=(AVE*INP)/(NP*INP) 

 ANOISE=SQRT(ANOISE) 

 WRITE (I0UT,492)AVE,AN0ISE 



492 FORMATC AVE= ',F8.3,' RMS WHITE NOISE = •,F8.3) 

 IF(JCODE.LT.O) GO TO 486 

 ENC00E(13,601,ITH) Bl 



501 FORMATC 'SLOPE =',F6.3) 



CALL SYMB0L(7.5,8.5,.20,ITH,0.,13) 

 ENCODE ( 19,502, ITH)B0 

 602 FORMATC 'INTERCEPT =',E8.3) 



CALL SYMB0L(7.5,8.0,.20,ITH,0.,19) 

 C PLOT REGRESSION LINE FROM POWWGT 



C CALL NEWPEN(3) 



IPLPSS=0 



191 



