subroutine openf77() OPEN(UNIT=10,FILE='f77speed.unf',FORM='UNFORMATTED') RETURN END subroutine writef77(array) REAL*8 array(128,128,128) WRITE(10) array RETURN END subroutine closef77() CLOSE(10) RETURN END