C                     
C	
C   TEYSEC.FOR	
C*  **********************************************
C
	COMMON //NX,NY,NZ,IXMIN,IYMIN,IZMIN,IXMAX,IYMAX,IZMAX
	DIMENSION NXYZ(3),MXYZ(3),NXYZST(3)
	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,string*80
C	character*1 aline(28),bline(131072),CLINE(131072)
	integer*2 inta(14),intb(65536),ZLP
	integer*4 linta(7),n,ii,ylayer,zlayer
	REAL*4 REALB(32768),REALC(32768)
	pi=4*ATAN(1.)
	R64K=256*256
 	R32K=128*256
2	format(a80)
3	format(32i3)
4	format(15i5) 
44 	format (15f5.0)
	WRITE(6,1000)
1000	FORMAT(//' TEY DENSITY MAP CONSTRUCTION PROGRAM Version 940401')
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')

	read(2,2) title
	write(6,*)title
	do 1201 i=1,3
	read(2,2)string
1201	write(6,*)string
	read(2,3)ixmin,ixmax,iymin,iymax,izmin,izmax
	write(6,*)ixmin,ixmax,iymin,iymax,izmin,izmax
	read(2,3)nx,ny,nz
	n=nx*ny
	write(6,*)nx,ny,nz,n
	do 1202 i=1,4
	read(2,2)string
1202	write(6,*)string

	OPEN (UNIT=3,FILE='map.TEY',STATUS='NEW',RECL=2042,form='unformatted')

	WRITE (3) TITLE
61	format(a80)
	write(6,*)'infile=',infile
	write(6,*)'izmin,max=',izmin,izmax
	DO 95 ZLP=IZMIN,IZmax
	WRITE(6,*)'ZLP=',ZLP
120	read(2,4,end=5000)(linta(i),i=2,6)
        if (zlp.eq.izmin) zoff=linta(2)
	zlayer=linta(2)-zoff
	IXMAX=linta(4)
 	write(6,*) 'zlayer=',zlayer,' ixmax=',ixmax
	do 90 iy=0,ny-1
	m=iy*nx
	read(2,44) (realb(m+ii),ii=1,nx)
 90	write(6,3) (int(realb(m+ii)),ii=1,nx)
	WRITE(3) zlayer,(LINTA(I),I=3,6)
	WRITE(3)(REALB(II),II=1,N)
95	continue
	CLOSE (UNIT=2)

5000	do 5010 i=2,6
5010	LINTA(I)=-1
	WRITE(3) (LINTA(I),I=2,6)
C  (NEGATIVE NUMBERS INDICATE NO MORE SECTIONS TO READ)
	close (unit=3)
	END