subroutine openf77()
	OPEN(UNIT=10,FILE='f77speed.unf',FORM='UNFORMATTED')
	RETURN
	END

	subroutine writef77(array)
	REAL*8 array(64,64,64)
	WRITE(10) array
	RETURN
	END

	subroutine closef77()
	CLOSE(10)
	RETURN
	END