CHARACTER DAT*12,TIM*12,filedate*23 CHARACTER*80 TITLE CHARACTER*50 INFILE,OUTFILE CHARACTER*50 FILENAME CHARACTER*40 FNAME,ASTRING,outfilmx CHARACTER*1 OUTSTRING(3800),HSTRING(16),BOXFLAG LOGICAL EFLG integer*2 X,Y,Z WRITE(6,1000) 1000 FORMAT(//' Read x, y, z from tffc outp (act fft line) ') 5001 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 (EFLG) GOTO 5 write(6,11001) goto 5001 11001 format(/' FILE DOES NOT EXIST ') 5 OPEN (UNIT=2,FILE=INFILE,STATUS='OLD') OPEN(UNIT=3,FILE='RTOMOD.COM',STATUS='NEW') WRITE(3,*)' $cd b:[berry.p21]' WRITE(3,*)' $r prg:rtomod' WRITE(3,*)' 80 120 140' WRITE(3,*)' 118.0 178.0 199.8 ' WRITE(3,*)' 106.0 ' !this is beta angle- only for second [.p21] version of rtomod! 370 READ(2,101,END=500)Y, X, Z 101 FORMAT (53X,I4,7X,I4,4X,I4) WRITE(3,*)' rot1.o' WRITE(3,*) X,Y,Z WRITE(3,*)' tran1.o' GOTO 370 500 CLOSE (UNIT=2) CLOSE (UNIT=3) end C 1 2 3 4 5 6 7 8 C 12345678901234567890123456789012345678901234567890123456789012345678901234567890 C Overall maximum rho = 0.0141 on section Y = 3 at X = 11 Z = 54