我只修改一段代码,剩下的没动。
加了一个变量:
Character(Len=512) :: cStr !// 读取用
[Fortran] 纯文本查看 复制代码 ! 读取数据
OPEN(UNIT = 415, FILE = 'plot.out')
Do !// 一直读
READ(415 ,IOSTAT = ERROR) cStr !// 先读入字符串
if ( Error /= 0 ) exit !// 如果读错,则退出循环(文件结束)
read( cStr , * ) cTime
if ( cTime == "Time" ) THEN
IF(IT == 1) THEN ! 如果是第一个时间步数据
D(IT) = DTIME
READ(415,IOSTAT = ERROR) IWELL,WELL,BB,CC,OR,WR,GR,G,WO,WC,FB
WELLNAME(IWELL) = WELL
OILRATE(IWELL,IT) = OR
WATERRATE(IWELL,IT) = WR
GASRATE(IWELL,IT) = GR
GOR(IWELL,IT) = G
WOR(IWELL,IT) = WO
WCUT(IWELL,IT) = WC
FBHP(IWELL,IT) = FB
ELSE
我算出的结果也是:
TIME OILRATE WATERRATE GASRATE GOR WOR WCUT FBHP
0.00000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
|