利用FORTRAN 读取文件行数的小代码

zsy011210804

fortran小白想问下，这个是有什么问题呢，
program getn
    integer n,k,i
   
    real*8 a(100),b(100),cc(100)   
     OPEN(1,FILE='inputabcc.txt')   
     do k=1,9
         read(1,*)
     enddo
     do while (.not.eof(1))
         read(1,*)i,a(k),b(k),cc(k)
     enddo
     close(1)
     n=i
     print*,n
     end program
txt文件如图，前九行不要，从第十行起读数据，i 就是行数，
fortran运行问题是出现断点
txt名字inputabcc.txt
数据如下
ndv
      3
     idv        xl        xu
      1     0.00000E+00     0.20000E+01
      2     0.20000E+01     0.80000E+01
      3     0.60000E+01     0.24000E+02
     nsp
     32
    isp              x1              x2              x3
      1     0.71875E+00     0.26562E+01     0.22594E+02
      2     0.11562E+01     0.73438E+01     0.22031E+02
      3     0.21875E+00     0.24688E+01     0.15281E+02
      4     0.84375E+00     0.43438E+01     0.79688E+01
      5     0.93750E-01     0.39688E+01     0.19781E+02
      6     0.15312E+01     0.22812E+01     0.90938E+01
      7     0.78125E+00     0.67812E+01     0.14719E+02
      8     0.18438E+01     0.69688E+01     0.20344E+02
      9     0.19062E+01     0.52812E+01     0.74062E+01
     10     0.34375E+00     0.30312E+01     0.68438E+01
     11     0.46875E+00     0.79062E+01     0.18094E+02
     12     0.19688E+01     0.35938E+01     0.17531E+02
     13     0.13438E+01     0.32188E+01     0.21469E+02
     14     0.14062E+01     0.41562E+01     0.14156E+02
     15     0.12812E+01     0.65938E+01     0.62812E+01
     16     0.40625E+00     0.47188E+01     0.16406E+02
     17     0.96875E+00     0.20938E+01     0.18656E+02
     18     0.10938E+01     0.37812E+01     0.10219E+02
     19     0.12188E+01     0.49062E+01     0.11906E+02
     20     0.28125E+00     0.64062E+01     0.23156E+02
     21     0.90625E+00     0.75312E+01     0.13594E+02
     22     0.59375E+00     0.50938E+01     0.20906E+02
     23     0.14688E+01     0.58438E+01     0.19219E+02
     24     0.16562E+01     0.45312E+01     0.23719E+02
     25     0.53125E+00     0.34062E+01     0.11344E+02
     26     0.65625E+00     0.60312E+01     0.96562E+01
     27     0.17188E+01     0.77188E+01     0.10781E+02
     28     0.15938E+01     0.62188E+01     0.15844E+02
     29     0.10312E+01     0.54688E+01     0.16969E+02
     30     0.31250E-01     0.56562E+01     0.12469E+02
     31     0.17812E+01     0.28438E+01     0.13031E+02
     32     0.15625E+00     0.71562E+01     0.85312E+01
请大大们指点一下

狼牙李

感谢楼主分享

贫道敬孔

多谢楼主分享{:5_213:}

polestar7

谢谢楼主分享，很有用，已经运行出来了~

tzclf

dddddddddddddddddddddddd

violetgyd

idl中直接有现成的函

转音2

感谢分享~~~~~~~~~非常感谢

1530718856

谢谢分享！！！！

qyx527161

感谢大大们的回复 用到啦

qyx527161

感谢大大们的回复 用到啦