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))
	

	CHARACTER*12 AQDATE2, AQTIME2
	INTEGER*4 ICCOUNTS(2)
	EQUIVALENCE (ICCOUNTS,BLINE2(25))
	INTEGER*4 EXPTIME
	EQUIVALENCE (EXPTIME,BLINE2(33))
	REAL*4 LAMBDA2,DISTANCE2,PHI_START2,PHI_END2
	EQUIVALENCE (LAMBDA2,BLINE2(77)),(DISTANCE2,BLINE2(81)),
     -  (PHI_START2,BLINE2(85)),(PHI_END2,BLINE2(89)),
     -  (AQDATE2,BLINE2(101)),(AQTIME2,BLINE2(113))
	CHARACTER*12 IPSYSTEM
	EQUIVALENCE (IPSYSTEM,BLINE2(125))
	I64K=256*256
	I32K=I64K/2

	DATA NXYZST/3*0/, CNV/57.29578/
C
 40	WRITE(6,1000)
1000	FORMAT(//' R-AXIS TO MARS IMAGE FILE CONVERSION PROGRAM',
     .		'  V  940105'/)
	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	MODE=1
	  WRITE(6,1600)
1600	  FORMAT('$Output filename:  ')
	  READ(5,1200) OUTFILE
       OPEN (UNIT=3,FILE=OUTFILE,STATUS='NEW',RECL=2400,
     &  CARRIAGECONTROL='NONE',RECORDTYPE='FIXED')


	
C	FIND MIN, MAX, MEAN DENSITY
	DMEAN=0
	DMAX=-1E32
	DMIN=1E32
	  NX=1200
	  NY=1200
	  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        
	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)


	AQDATE2=DATOPERTARG(1:12)
	AQTIME2=DATOPERTARG(13:24)
	EXPTIME=EXPOSURE*60
	ICCOUNTS(1)=0
	ICCOUNTS(2)=0
	LAMBDA2=WAVELENGTH
	DISTANCE2=CAMLENGTH
	PHI_START2=PHI(2)
	PHI_END2=PHI(3)
	IPSYSTEM='R-AXIS'
   	write(3,155) (BLINE2(II),II=1,2400)
	  DO 88 IY=1,350
 88          READ(1,155) BLINE(1)
	  DO 95 IY = 1,NY
	    READ(1,155) (BLINE(II),II=1,700),(BLINE2(II),II=1,2400)
C skip over 350 values (700 bytes), take next 1200
c 	    READ(1,155) (BLINE(II),II=1,4096)
c	    DO 90 IX=1,NX
c 90	    ILINE2(IX)=ILINE(IX+350)           
 95   	     write(3,155) (BLINE2(II),II=1,2400)
	CLOSE (UNIT=3)
995	CLOSE (UNIT=1)
	GOTO 40
	END