是的我又来了
同样的档案 其中一行像是底下这样
然后有好几千行
想要一次读取1000行
他原来是char 想要改成 integer or real
档案其中一行是这样
1
-6.08974 -0.389376 -2.5e-05 1096.12 157.843 1 20.7615 938.272 2212 1 -2.16598
-0.777673 -566.665 0 0 0
还有几千行 下略
我现在可以读成char一次n行
不过不知道怎么转成integer 和 real
看了别人写的照写似乎没有用
还请前辈们指点迷津
Character*80 NP(50),Cg(50),PID(50),X(50),Y(50),Z(50),TE(50)
& ,E(50),WT(50),T(50),MS(50),U(50),V(50),W(50),Px(50),Py(50),
& Pz(50),MU(50),MV(50),MW(50)
open(unit=80,file='test.phsp',access='sequential')
do i=1,10
read(80,'(I2)') NP(i)
read (80,120) X(i),Y(i),Z(i),TE(i),E(i),WT(i),T(i),MS(i),
& PID(i),Cg(i),MU(i),MV(i),MW(i),Px(i),Py(i),Pz(i)
120 format ('F6.3,F6.3,E9.3,F8.3,F7.3,I2,F8.5,F8.3,I4,I2,F5.3,F5.3,
& F5.3,I1,I1,I1')
end do
do i=1,50
write(*,*) NP(i),X(i),Y(i),Z(i),TE(i),E(i),WT(i),
& T(i),MS(i),PID(i),Cg(i),MU(i),MV(i),MW(i),Px(i),Py(i),Pz(i)
end do
close(80)
stop
end
非常感谢