[Fortran] 纯文本查看 复制代码
PROGRAM CEA2
IMPLICIT NONE
INCLUDE 'cea.inc'
C LOCAL VARIABLES
CHARACTER*15 ensert(20)
CHARACTER*200 infile,ofile
CHARACTER*196 prefix
LOGICAL caseok,ex,readok
INTEGER i,inc,iof,j,ln,n
INTEGER INDEX
REAL*8 xi,xln
REAL*8 DLOG
SAVE caseok,ensert,ex,i,inc,infile,iof,j,ln,n,ofile,prefix,readok,
& xi,xln
C
WRITE (*,99001)
READ (*,99002) prefix
ln = INDEX(prefix,' ') - 1
infile = prefix(1:ln)//'.inp'
ofile = prefix(1:ln)//'.out'
Pfile = prefix(1:ln)//'.plt'
INQUIRE (FILE=infile,EXIST=ex)
IF ( .NOT.ex ) THEN
PRINT *,infile,' DOES NOT EXIST'
GOTO 400
ENDIF
OPEN (IOINP,FILE=infile,STATUS='old',FORM='formatted')
OPEN (IOOUT,FILE=ofile,STATUS='unknown',FORM='formatted')
OPEN (IOSCH,STATUS='scratch',FORM='unformatted')
OPEN (IOTHM,FILE='thermo.lib',FORM='unformatted')
OPEN (IOTRN,FILE='trans.lib',FORM='unformatted')
WRITE (IOOUT,99006)
WRITE (IOOUT,99007)
WRITE (IOOUT,99006)
readok = .TRUE.
Newr = .FALSE.
100 Iplt = 0
Nplt = 0
CALL INPUT(readok,caseok,ensert)