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
	integer*2 inta(14)
	integer*4 linta(7),n,ii,ylayer
	REAL*4 REALB(500,500)
c	REAL*4 REALB(250000)
	pi=4*ATAN(1.)
	R64K=256*256
 	R32K=128*256
 
	WRITE(6,1000)
1000	FORMAT(//' EXCHANGE FAST AND MEDIUM AXES OF TEY 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')
        OPEN (UNIT=3,FILE='new.map',STATUS='NEW',form='unformatted')

	read(2) title
11	format (A2,a80)
	write(6,*)title
	WRITE(3) TITLE
	DMEAN=0
	NY=0
	IYMIN=10000
	IYMAX=-10000
20	read(2) (linta(i),i=2,6)
21	format(A2,24a1)
31	format(' ',14i6)
	ylayer=linta(2)
	write(6,*) 'ylayer=',ylayer
	IF (ylayer.lt.iymin)iymin=ylayer	
	IF (ylayer.gt.iymax)iymax=ylayer	
	IXMIN=linta(3)
	IXMAX=linta(4)
	IZMIN=linta(5)
	IZMAX=linta(6)

	linta(3)=IZMIN
	linta(4)=IZMAX
	linta(5)=IXMIN
	linta(6)=IXMAX
  	WRITE(3) (linta(i),i=2,6)
	if(ylayer.LT.0) goto 985

	nx=(1+IXMAX-IXMIN)
	nz=(1+IZMAX-IZMIN)
	ny=NY+1
	n=nx*nz
c	write(6,*)'nx, nz, n=',nx,nz,n
	n=1
c	read(2) (realb(ii),ii=1,n)
	READ(2) ((REALB(II,JJ),II=1,NX),JJ=1,NZ)
	WRITE(3) ((REALB(II,JJ),JJ=1,NZ),II=1,NX)
    	goto 20

985	CLOSE(UNIT=2)
	CLOSE(UNIT=3)

	write(6,*)'That is Assuming Y-sections, Z is medium axis!'
	END