主题:格式问题第二弹 依然求解答!
想读入的文件是NLIST.DAT 的全部四行
--------------------------------------------------
NODE X Y Z
122 -5.00000000000 -3.00000000000 0.00000000000
123 5.00000000000 -3.00000000000 0.00000000000
124 -4.80000000000 -3.00000000000 0.00000000000
125 -4.60000000000 -3.00000000000 0.00000000000
126 -4.40000000000 -3.00000000000 0.00000000000
127 -4.20000000000 -3.00000000000 0.00000000000
128 -4.00000000000 -3.00000000000 0.00000000000
129 -3.80000000000 -3.00000000000 0.00000000000
130 -3.60000000000 -3.00000000000 0.00000000000
131 -3.40000000000 -3.00000000000 0.00000000000
132 -3.20000000000 -3.00000000000 0.00000000000
133 -3.00000000000 -3.00000000000 0.00000000000
134 -2.80000000000 -3.00000000000 0.00000000000
135 -2.60000000000 -3.00000000000 0.00000000000
136 -2.40000000000 -3.00000000000 0.00000000000
137 -2.20000000000 -3.00000000000 0.00000000000
138 -2.00000000000 -3.00000000000 0.00000000000
139 -1.80000000000 -3.00000000000 0.00000000000
140 -1.60000000000 -3.00000000000 0.00000000000
141 -1.40000000000 -3.00000000000 0.00000000000
NODE X Y Z
142 -1.20000000000 -3.00000000000 0.00000000000
143 -1.00000000000 -3.00000000000 0.00000000000
.................
类似这样的数据块 ( 空两行 + 一行英文标识 +20行数据)一共有69块
--------------------
程序部分
REAL MATNO(2000),LNODS(2000,9),TITLE(30)
REAL COORD(2000,3)
NELEM=1280
NNODE=4
OPEN(6,FILE='ELIST.DAT',STATUS='OLD')
OPEN(7,FILE='NLIST.DAT',STATUS='OLD')
OPEN(8,FILE='ELEMENT.DAT',STATUS='UNKNOWN')
OPEN(9,FILE='NODE COORDINATES.DAT',STATUS='UNKNOWN')
DO 20 I=1,69
READ(7,*)
READ(7,*)
DO 20 J=1,20
READ(7,200) IPOIN,(COORD(IPOIN,IDIME),IDIME=1,3)
WRITE(9,200) IPOIN,(COORD(IPOIN,IDIME),IDIME=1,3)
200 FORMAT(3X,I5,3X,F14.12,3X,F14.12)
20 CONTINUE
输出结果
-----------------------------------------------
122 ************** **************
123
124 ************** **************
125
126 ************** **************
127
128 ************** **************
129
130 ************** **************
131
132 ************** **************
133
134 ************** **************
135
136 ************** **************
137
138 ************** **************
139
140 ************** **************
141
数据只有这么多 而且显示数组COORD(2000,3)越界 我调整数组为双精度时 依然显示越界
求解答!!!!
--------------------------------------------------
NODE X Y Z
122 -5.00000000000 -3.00000000000 0.00000000000
123 5.00000000000 -3.00000000000 0.00000000000
124 -4.80000000000 -3.00000000000 0.00000000000
125 -4.60000000000 -3.00000000000 0.00000000000
126 -4.40000000000 -3.00000000000 0.00000000000
127 -4.20000000000 -3.00000000000 0.00000000000
128 -4.00000000000 -3.00000000000 0.00000000000
129 -3.80000000000 -3.00000000000 0.00000000000
130 -3.60000000000 -3.00000000000 0.00000000000
131 -3.40000000000 -3.00000000000 0.00000000000
132 -3.20000000000 -3.00000000000 0.00000000000
133 -3.00000000000 -3.00000000000 0.00000000000
134 -2.80000000000 -3.00000000000 0.00000000000
135 -2.60000000000 -3.00000000000 0.00000000000
136 -2.40000000000 -3.00000000000 0.00000000000
137 -2.20000000000 -3.00000000000 0.00000000000
138 -2.00000000000 -3.00000000000 0.00000000000
139 -1.80000000000 -3.00000000000 0.00000000000
140 -1.60000000000 -3.00000000000 0.00000000000
141 -1.40000000000 -3.00000000000 0.00000000000
NODE X Y Z
142 -1.20000000000 -3.00000000000 0.00000000000
143 -1.00000000000 -3.00000000000 0.00000000000
.................
类似这样的数据块 ( 空两行 + 一行英文标识 +20行数据)一共有69块
--------------------
程序部分
REAL MATNO(2000),LNODS(2000,9),TITLE(30)
REAL COORD(2000,3)
NELEM=1280
NNODE=4
OPEN(6,FILE='ELIST.DAT',STATUS='OLD')
OPEN(7,FILE='NLIST.DAT',STATUS='OLD')
OPEN(8,FILE='ELEMENT.DAT',STATUS='UNKNOWN')
OPEN(9,FILE='NODE COORDINATES.DAT',STATUS='UNKNOWN')
DO 20 I=1,69
READ(7,*)
READ(7,*)
DO 20 J=1,20
READ(7,200) IPOIN,(COORD(IPOIN,IDIME),IDIME=1,3)
WRITE(9,200) IPOIN,(COORD(IPOIN,IDIME),IDIME=1,3)
200 FORMAT(3X,I5,3X,F14.12,3X,F14.12)
20 CONTINUE
输出结果
-----------------------------------------------
122 ************** **************
123
124 ************** **************
125
126 ************** **************
127
128 ************** **************
129
130 ************** **************
131
132 ************** **************
133
134 ************** **************
135
136 ************** **************
137
138 ************** **************
139
140 ************** **************
141
数据只有这么多 而且显示数组COORD(2000,3)越界 我调整数组为双精度时 依然显示越界
求解答!!!!