CHARACTER*1 BLINE(1024) CHARACTER*36 RECLINE(30) EQUIVALENCE (RECLINE(1),BLINE(1)) real*4 rline(256),TCELL(6),CELL(6) REAL*4 RMX(3,3),RMY(3,3),RMZ(3,3),TEMP(3,3),PROD(3,3),X1(3),X2(3) EQUIVALENCE (RLINE(1),BLINE(1)),(TCELL(1),RLINE(11)) INTEGER*4 LINT(256),ITRANS(3) EQUIVALENCE (LINT(1),BLINE(1)) CHARACTER*80 TTITLE,TITLE,TSYMTRY,SYMTRY EQUIVALENCE (TTITLE,RLINE(56)) EQUIVALENCE (TSYMTRY,RLINE(1)) CHARACTER*50 INFILE,OUTFILE,FILENAME LOGICAL*1 EFLG DIMENSION RMAP(100,100,300) WRITE(6,1001) 1001 FORMAT(//' PROGRAM TO READ 2-FOLD SECTION WRITTEN BY CCP4 ROTATE') 5 WRITE(6,1002) 1002 FORMAT(/'$Input filename (SECTION.DAT): ') READ(5,1003) INFILE 1003 FORMAT(A) IF (INFILE(:5).EQ.' ') INFILE='SECTION.DAT' INQUIRE(FILE=INFILE, NAME=FILENAME,EXIST=EFLG) IF (.NOT.EFLG) GOTO 5 OPEN (UNIT=2,FILE=INFILE,STATUS='OLD',FORM='UNFORMATTED') WRITE(6,*)'ENTER RECORD LENGTH AND # RECORDS:' READ(5,*) NWORDS,NRECORDS DO 290 K=1,NRECORDS READ(2)(RLINE(I),I=1,NWORDS) DO 280 I=1,NWORDS 280 IF (RLINE(I).LT.-9.) RLINE(I)=-9. 290 WRITE(6,291)(RLINE(I),I=1,NWORDS) 291 FORMAT(' ',50F3.0) END