C   *********************************************                  
C	
C   ANOSCALE2HEAVY.FOR	
C*  **********************************************
C
	CHARACTER*80 TITLE
	CHARACTER*50 INFILE,OUTFILE,FILENAME
        CHARACTER*40 FNAME,ASTRING
	real int1, int2, sigmaI1, sigmaI2, SigFcut
	real fobs1, fobs2, sigmaF1,  sigmaF2 
	real  ss1, ss2, ss
	integer*2 h,k,l
	LOGICAL EFLG

	WRITE(6,1000)
1000	FORMAT(//' Enter name of anomalous sca file to convert:')
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)
	  IF (.NOT.EFLG) GOTO 5
        OPEN (UNIT=2,FILE=INFILE,readonly,STATUS='OLD')
c	[skip over first 3 lines)
	do 140 i=1,3
140	read(2,'(A)')astring

c       [input sigmaF cutoff]
c         write (6, 250)
250      format ('$Input sigmaF cutoff:' )
         read  (5, * ) SigFcut 

	  WRITE(6,1600)
1600	  FORMAT('$Output filename:  ')
	  READ(5,1200) OUTFILE
	  OPEN(UNIT=3,FILE=OUTFILE,STATUS='unknown')
c       OPEN (UNIT=3,FILE='temp.rfl',STATUS='NEW')
200	read(2,220,end=990)h,k,l,int1,sigmaI1,  int2, sigmaI2
220      format (3I4, 4f8.1)

	if (sigmaI1.ne.0.) write(3,300)  h,k,l,int1,sigmaI1
	if (sigmaI2.ne.0.) write(3,300)  -h,-k,-l, int2, sigmaI2
300	FORMAT (3I4,F10.2,F8.2,F10.2,F8.2 )
	goto 200

990	CLOSE(UNIT=2)
    	CLOSE (UNIT=3)
	END