Appendix D 



NOTE: Many of the subroutines required by the following 

 routines are common to the province selection 

 software listed in Appendix B.2. 



C GENERATE AMPLITUDE SPECTRUM OF DEPTHS OUTPUT FROM PROVINCE PICKER. 



C PROGRAM IS A SLIGHT MODIFICATION OF THE ORIGINAL CREATED BY 



C T.M. DAVIS, IN WHICH THE INPUTS HAVE BEEN SET WITHIN THE CODE 



C AND ONLY THOSE PARAMETERS NECESSARY FOR DAILY USE HAVE BEEN LEFT 



C FOR INPUT BY THE USER. INPUT IS EXPECTED FROM UNIT 13, WHICH IS 



C THE UNIT NUMBER USED BY PROVPICK. THE PROGRAM ADOS AN INTERACTIVE 



C GRAPHICS SECTION WHICH ALLOWS MODIFICATION OF THE INPUT SERIES, 



C DELETION OF SECTIONS, AND CONTROL OF THE INDEPENDENT VARIABLE FOR 



C REGRESSION ANALYSIS. ROUGHNESS MODEL PARAMETERS A & B ALONG WITH 



C LATITUDE AND LONGITUDE INFORMATION FOR EACH PROFILE ARE OUTPUT 



C TO UNIT 14 FOR LATER USE IN PROVCHART PROGRAM. 



C 



C PROGRAMMMED BY C.G.FOX,ADVANCED TECHNOLOGY STAFF, NAVOCEANO, 5/17/84 



C 



C THE FOLLOWING COMMENTS ARE FROM THE ORIGINAL PROGRAM BY T.M.DAVIS 



C 



C 



C PROGRAM TO COMPUTE PREWHITENED, CORRECTED AND SMOOTHED AMPLITUDE 



C SPECTRUM CUT,H,N=PARAMETERS TO CONTROL PREWHITENING AND SMOOTHING 



C FILTERS ,ANORM=SPECTRUM NORMALIZATION IN TERMS OF DATA INTERVALS 



C XORG,YORG=LOG-LOG PLOT ORIGIN IN POWERS OF 10,ITG=HEADING 



C IUNIT=5 IF CONTROL DATA ON CARDS = TAPE UNIT IF ON TAPE 



C SAME FOR JUNIT FOR INPUT DATA 



C TAPE UNIT 2=0UTPUT 



C CODE IFE0F=1 IF EOF FOLLOWS EACH SET 



C FIRST SET OF FILTER PARAMETERS ARE FOR HIGH PASS PREWHITENING 



C SECOND SET ARE FOR LOW PASS SMOOTHING 



C IF N(l)=-1 COSINE TAPER IS APPLIED, =0 NO TAPER OR PREWHITENING 



C LEAVE HIGH OR LOW PASS PARAMETERS BLANK IF NO FILTER DESIRED 



C INPUT DATA LIMIT IS 2048 PTS,NP=NO.OF INPUT PTS 



C IF ANORM IS BLANK SET TO 1 IN PROGRAM 



C USES SUBROTINE NLOGN FROM ROBINSON 



C IN FIRST CONTROL CARD NSETS=NO.OF SETS OF DATA THIS RUN 



C SET KPHA=1 IF PHASE SPECTRUM IS DESIRED 



C IORGN=NO.OF INPUT DATA PT.TO USE AS ORIGIN FOR PHASE SPECTRUM 



C IF PHASE IS DESIRED CODE IPHA=0 IF PHASE IN DEGREES OR CODE IPHA 



C =NO.OF DATA INTERVALS/INCH FOR PLOT IF YOU WANT PHASE IN DATA INT. 



C SET JC0DE=1 IF ONLY PLOT OF FINAL SMOOTHED SPECTRUM DESIRED 



C SET JC0DE=-1 FOR NO PLOT AT ALL 



C CODE ILIST=1 IF YOU DESIRE ONLY LISTING OF SMOOTHED SPECTRUM 



C **IMPORTANT-IF YOU HAVE ALREADY ADDED ZEROES TO THE BEGINNING OR END 



C OF YOUR INPUT DATA AND YOU REQUEST PREWHITENING OR COSINE TAPER 



C YOU MUST CODE IADJ=1 FOR PROPER EXECUTION 



C IMEAN=1 IF YOU WANT THIS MEAN REMOVED 



DIMENSION DATA(4500),ST0(2500),CUT(2),H(2),N(2),ADATA(4100) 



185 



