C                     
C	
C   TEYSEC.FOR	
C*  **********************************************
C
	COMMON //NX,NY,NZ,IXMIN,IYMIN,IZMIN,IXMAX,IYMAX,IZMAX
	DIMENSION NXYZ(3),MXYZ(3),NXYZST(3)
c	CHARACTER*80 TITLE
	DIMENSION IXYZMIN(3),IXYZMAX(3)
	CHARACTER*50 INFILE,OUTFILE,FILENAME
        CHARACTER*40 FNAME,ASTRING
	CHARACTER*1 OUTSTRING(1200),HSTRING(16),BOXFLAG
	LOGICAL EFLG
	character LABEL*80,title*80,junk*2
C	character*1 aline(28),bline(131072),CLINE(131072)
	integer*2 inta(14),intb(65536),ZLP
	integer*4 linta(7),n,ii,ylayer
	REAL*4 REALB(32768),REALC(32768)
	pi=4*ATAN(1.)
	R64K=256*256
 	R32K=128*256
 
	WRITE(6,1000)
1000	FORMAT(//' FIND MAX DENSITY IN TEY DENSITY MAP')
5	WRITE(6,1100)
1100	FORMAT(/'$Input filename (CR to quit):  ')
	READ(5,1200) INFILE
1200	FORMAT(A)
	IF (INFILE(:5).EQ.'     ') STOP
	  INQUIRE(FILE=INFILE, NAME=FILENAME,EXIST=EFLG)
	WRITE(6,*)'INFILE=',INFILE,'EFLG=',EFLG,' FILENAME=',FILENAME
	  IF (.NOT.EFLG) GOTO 5
        OPEN (UNIT=2,FILE=INFILE,STATUS='OLD',form='unformatted')

	read(2) title
11	format (A2,a80)
	write(6,*)title
	DMEAN=0
	DMAX=-1E32
	DMIN=1E32
	NY=0
	IYMIN=10000
	IYMAX=-10000
20	read(2) (linta(i),i=2,6)
21	format(A2,24a1)
c	write(6,31) (inta(i),i=1,12)
31	format(' ',14i6)
	ylayer=linta(2)
	write(6,*) 'ylayer=',ylayer
	if(ylayer.LT.0) goto 985
	IF (ylayer.lt.iymin)iymin=ylayer	
	IF (ylayer.gt.iymax)iymax=ylayer	
	IXMIN=linta(3)
	IXMAX=linta(4)
	IZMIN=linta(5)
	IZMAX=linta(6)
	write(6,33) ylayer,IXMIN,IXMAX,IZMIN,IZMAX
33	FORMAT(' ylayer:',I5,' IXMIN:',I5,' IXMAX:',I5,' IZMIN:',I5,
     .  ' IZMAX',I5)

	nx=(1+IXMAX-IXMIN)
	nz=(1+IZMAX-IZMIN)
	ny=NY+1
	n=nx*nz
	read(2) (realb(ii),ii=1,n)
41	format(2x,2042a1)
51	FORMAT('AA',2042A1)

	do 80 i=1,nz
	DO 80 J=1+(I-1)*NX,I*NX
	IF (REALB(J).GT.DMAX) THEN
	 DMAX=REALB(J)
	 XM=J-(i-1)*nx
	 YM=NY
	 ZM=I
	 ENDIF
	IF (REALB(J).LT.DMIN) THEN
	 DMIN=REALB(J)
	 XN=J-(i-1)*nx
	 YN=NY
	 ZN=I
	 ENDIF
80	CONTINUE
    	goto 20

985	CLOSE(UNIT=2)

	WRITE(6,*)'MAX DENSITY AT X,Y,Z=',XM,YM,ZM,' VALUE IS',DMAX
	WRITE(6,*)'MIN DENSITY AT X,Y,Z=',XN,YN,ZN,' VALUE IS',DMIN
	write(6,*)'That is Assuming Y-sections, Z is medium axis!'
	END