C USE FOLLOWING CALL FOR CALCOMP PLOTTER OUTPUT 



C CALL PL0TS(0,0,10) 



C USE FOLLOWING CALL FOR ZETA PLOTTER OUTPUT 



CALL PLOTS (53, 0,-10) 



CALL FACTOR (.5) 



CALL PLOTd. 0,1. 0,-3) 

 C READ NAME OF FILE CONTAINING LIST OF DATA FILE NAMES 



24 READ(5,'(A20)') INFILE 



25 0PEN(UNIT=18,FILE=INFILE,STATUS='0LD' ) 



C READ LIST OF DATA FILE NAMES FROM FILE INFILE 



DO 5 IK=1,JT0TL 



READ(18,'(A20)') FILE(IK) 

 5 CONTINUE 



CL0SE(UNIT=18) 

 C READ X,Y,Z, VALUES FOR THIS TRACK. FIRST POINT IS 

 C ASSUMED TO BE ORIGIN, PUT END OF FILE AT THE END OF EACH TRACK. 



440 IK0NT=IK0NT+1 



C ZERO OUT TEMPORARY ARRAYS 

 DO 441 J=l,10 

 DO 441 I=1,ISIZ 



441 TZ(I,J)= 0.0 



DO 442 I=1,ISIZ 



AX(I)=0.0 



AY(I)=0.0 



AZ(I)=0.0 



Z(I)=0.0 



Y(I)=0.0 



442 X(I)=0.0 

 IAVPLT=0 



C INPUT X(I)=LONGITUDE(DEC. DEG. ),Y(I)=LAT,Z(i)=DEPTH 

 C N=# OF POINTS - 1 

 IKNT=IKONT 

 JPTS=KPTS 



IF(ILIST.EQ.l) WRITE(6,45) FILE(IKONT) 

 45 FORMATC OPENING FILE ',A20) 

 C ROUTINE TO READ LAT,LON, DEPTH FROM FILE # IKNT 

 CALL PROVRD(JPTS,Y,X,Z, IKNT, FILE) 

 4 N=JPTS-1 



IF(ILIST.EQ.l) WRITE(6,'(24H NUMBER OF INPUT POINTS ,15)' )N 

 C INTERPOLATE DATA ON STRAIGHT LINE 



CALL MAPCTN(GRID,X,Y,Z,N,AX,AY,AZ,JCT) 



IF(ILIST.EQ.1)WRITE(6,'(31H NUMBER OF INTERPOLATED POINTS ,15)') 

 IJCT 

 839 C=JCT 



AZAVE=0.0 

 C PLOT INTERPOLATED DATA 

 DO 261 I=1,JCT 

 261 AZAVE=AZAVE+(AZ(I)/C) 

 IF(IPLOT.NE.l) GO TO 26 

 XL=(C*.02)+0.5 

 DV=50.*GRID 

 CALL NEWPEN(l) 



151 



