603 FIRST1=FIRSTK 

 C WRITE(I0UT,'(8HFIRST1= ,E11.4) ' )FIRST1 



X=AL0G10(FIRST1/X0RG)*3.0125 



AMP1=B0*FIRST1**B1 

 C WRITE(I0UT,'(6HAMP1= ,E11.4) ' )AMP1 



TPP=AMP1/Y0RG 



IF(TPP.LT.l.E-20) TPP=ABS(TPP) 



Y=-0. 5+ALOGlO (TPP )*1 . 35714 



CALL PL0T(X,Y,+3) 

 C CALL SYMBOL < X,Y, . 1,32,0.0,-1) 



FIRST1=((IA)*DELK) 

 C WRITE(I0UT,'(8HFIRST1= ,E11,4) ' )FIRST1 



X=AL0G10(FIRST1/X0RG)*3.0125 



AMP1=B0*FIRST1**B1 

 C WRITE(I0UT,'(5HAMP1= ,E11.4) ' )AMP1 



TPP=AMP1/Y0RG 



IF(TPP.LT.l.E-20) TPP=ABS(TPP) 



Y=-0.5+AL0G10(TPP)*l. 35714 



CALL PL0T(X,Y,+2) 

 C CALL SYMBOL ( X,Y, . 1,32,0.0,-2) 

 C CALL NEWPEN(l) 

 C PLOT SPECTRAL ESTIMATES 



X=AL0G10 (FIRSTK/X0RG)*3.0125 



TPP=ADATA(2)/Y0RG 



IF(TPP.EQ.0.0)TPP=1.E-15 



IF(TPP.LT.l.E-20) TPP=ABS(TPP) 



Y=-0.5+AL0G10(TPP)*l. 35714 



CALL PL0T(X,Y,+3) 

 C CALL SYMBOL ( X,Y. .05,32,0.0,-1) 



00 238 J=3,IA-1 



XJ=J-1+N(2) 



X=AL0G10((XJ*0ELK)/X0RG)*3.0125 



TPP=ADATA{J)/YORG 



IF(TPP.EQ.0.0)TPP=1.E-15 



IF(TPP.LT.l.E-20) TPP=ABS(TPP) 



Y=-0.5+AL0G10(TPP)*l. 35714 

 238 CALL PL0T(X,Y,+2) 

 C 238 CALL SYMBOL ( X,Y,.05,32,0.0,-2) 



CALL PLOT (0.0, 0.0, +3) 



IF(IPLPSS.EQ.l) GO TO 330 



IF(LPASS.EQ.l) GO TO 330 



IPLPSS=1 



BO=X INTER 



B1=SL0PEX 

 C INCLUDE NEXT LINE IF A PLOT OF SPATIAL DOMAIN ESTIMATE IS DESIRED 

 C GO TO 603 



C PLOT AX.IS 



330 IF(ANORM.EQ.l.O) GO TO 333 



CALL AXES(0.,-. 5, 'NORMALIZED AMPLITUDEdN KM)' ,27,9.5,90.0 

 *, 1.3514, 10. 0,0. 0,-1) 



GO TO 332 

 333 CALL AXES (0.,-. 5, 'AMPLITUDE (IN KILOMETERS )' ,25,9.5,90.0 

 *, 1.3514, 10. 0,0. 0,-1) 



192 



