C* RAXIS2MRC.FOR   
C*  DERIVED FROM LABEL.FOR**********************************************
C
	COMMON //NX,NY,NZ,IXMIN,IYMIN,IZMIN,IXMAX,IYMAX,IZMAX
	DIMENSION ALINE(9998),TITLE(20),NXYZ(3),MXYZ(3),NXYZST(3)
	DIMENSION IXYZMIN(3),IXYZMAX(3),OUT(9998)
	DIMENSION LABELS(20,10),CELL(6)
	COMPLEX CLINE(4999),COUT(4999)
	CHARACTER*50 INFILE,OUTFILE
      	REAL*8 DOUBLMEAN
	CHARACTER*1 BLINE(4096),BLINE2(2400)
	INTEGER*2 ILINE(2048),ILINE2(1200)
 	EQUIVALENCE (BLINE,ILINE),(BLINE2,ILINE2)
	EQUIVALENCE (NX,NXYZ), (ALINE,CLINE), (OUT,COUT), (BLINE,ILINE)
	EQUIVALENCE (IXYZMIN, IXMIN), (IXYZMAX, IXMAX)
	CHARACTER*20 DEVNAMVERS
       	EQUIVALENCE (DEVNAMVERS,BLINE(1))
	CHARACTER*32 CRYSTAL
       	EQUIVALENCE (CRYSTAL,BLINE(21))
	REAL*4 CELLPAR(6)
	EQUIVALENCE (CELLPAR,BLINE(53))
	CHARACTER*12 SPACEGROUP
	EQUIVALENCE (SPACEGROUP,BLINE(77))
	REAL*4 MOSAICITY
	EQUIVALENCE (MOSAICITY,BLINE(89))
	CHARACTER*80 MEMO
	EQUIVALENCE (MEMO,BLINE(93))
	CHARACTER*36 DATOPERTARG
	EQUIVALENCE (DATOPERTARG,BLINE(257))
	REAL*4 WAVELENGTH
	EQUIVALENCE (WAVELENGTH,BLINE(293))
	CHARACTER*20 MONOCHROM
	EQUIVALENCE (MONOCHROM,BLINE(297))
	REAL*4 TWOTHETA
	EQUIVALENCE (TWOTHETA,BLINE(317))
	CHARACTER*24 COLLFILT
	EQUIVALENCE (COLLFILT,BLINE(321))
	REAL*4 CAMLENGTH
	EQUIVALENCE (CAMLENGTH,BLINE(345))
	REAL*4 XRAYPRM(2)
	EQUIVALENCE (XRAYPRM,BLINE(349))
	CHARACTER*8  AXES
	EQUIVALENCE (AXES,BLINE(513))
	REAL*4 PHI(3)
		EQUIVALENCE (PHI,BLINE(521))
	INTEGER*4  NOSCIL
	EQUIVALENCE (NOSCIL,BLINE(533))
	REAL*4 EXPOSURE
	EQUIVALENCE (EXPOSURE,BLINE(537))
	REAL*4 BEAMCENTER(2)
	EQUIVALENCE (BEAMCENTER,BLINE(541))
	

	R64K=256*256
	I32K=I64K/2

	DATA NXYZST/3*0/, CNV/57.29578/
C
 40	WRITE(6,1000)
1000	FORMAT(//' R-AXIS HEADER DUMP PROGRAM',
     .		'  V  941228'/)
50	WRITE(6,1100)
1100	FORMAT('$Input filename (CR to quit):  ')
	READ(5,1200) INFILE
1200	FORMAT(A)
	IF (INFILE(:5).EQ.'     ') STOP
        OPEN (UNIT=1,FILE=INFILE,STATUS='OLD')
	
C	FIND MIN, MAX, MEAN DENSITY
	DMEAN=0
	DMAX=-1E32
	DMIN=1E32
	  NX=1900
	  NY=1900
C	  type *,nx,ny,nz

	SCALE = 1.0
	xmin = 0.


C	read HEADER into bline, encode into titles:
	READ(1,155) (BLINE(II),II=1,4096)
155	FORMAT (4096A1)
	TYPE *,'DEVNAM: ',DEVNAMVERS
	TYPE *,'CRYSTAL: ',CRYSTAL
	TYPE *,'CELLPAR: ',CELLPAR
	TYPE *,'SPACE GROUP: ',SPACEGROUP
	TYPE *,'COMMENT: ',MEMO(:70)        
	TYPE *,'DATE,OPERATOR, TARGET: ',DATOPERTARG
	TYPE *,'WAVELENGTH: ',WAVELENGTH
	TYPE *,'MONOCHROM: ',MONOCHROM
	TYPE *,'COLLIMATOR, FILTER: ',COLLFILT
	TYPE *,'CAMERA LENGTH: ',CAMLENGTH
	TYPE *,'MOUNT, BEAM AXES: ',AXES
	TYPE *,'BEAMCENTER: ',(BEAMCENTER(I),I=1,2)
	TYPE *,'EXPOSURE, PHI(1-3):',EXPOSURE, (PHI(I),I=1,3)
	TYPE *,'READING IMAGE.  FINDING MIN, MAX DENSITY'

	  DO 85 IY = 1,NY

	READ(1,155) (BLINE(II),II=1,3800)
	    DO 80 IX=1,NX
C		TYPE *, ILINE(IX)
 	     X=ILINE(IX)           
		IF (X.LT.0)X=X+R64K
		DMEAN=DMEAN+X
		IF (X.GT.DMAX) DMAX=X
 80		IF (X.LT.DMIN) DMIN=X
 85	CONTINUE

	DMEAN=DMEAN/(NX*NY)
C  	ENCODE(80,6750,TITLE) DMIN,DMAX,DMEAN
	WRITE(6,6750)DMIN,DMAX,DMEAN
6750	FORMAT('Min/Max/mean density = ',2f8.0,f8.2)          

 	CLOSE (UNIT=1)
	GOTO 50
	END