[求助] 利用数组转化坐标文件出现问题（已解决）

open(17,file=in_filename,status='old')
read(17,*)atomnum
read(17,*)step,stepnum

maxn=max(na,nb,nc)

allocate(label(atomnum))
allocate(coordinate(atomnum,3))
allocate(new_coordi(atomnum,3,maxn))

  do i1=1,atomnum
     read(17,*)label(i1),(coordinate(i1,j1),j1=1,3)
  end do
close(17)

a_count=1
b_count=1
c_count=1
new_atomnum=atomnum*na*nb*nc
open(18,file=out_filename)
write(18,"(10XI3)")new_atomnum
write(18,"(1XA5,5XI6)")step,stepnum

do i1=1,atomnum
do while(a_count.le.na)
   new_coordi(i1,1,a_count)=coordinate(i1,1)+a_count*a
     do while(b_count.le.nb)
     new_coordi(i1,2,b_count)=coordinate(i1,2)+b_count*b
       do while(c_count.le.nc)
       new_coordi(i1,3,c_count)=coordinate(i1,3)+c_count*c
       write(18,"(1XA1,3(F14.6))")label(i1),new_coordi(i1,1,a_count),new_coordi(i1,2,b_count),new_coordi(i1,3,c_count)
       c_count=c_count+1
       end do
     b_count=b_count+1
     end do
   a_count=a_count+1
end do   
end do
close(18)

deallocate(label)
deallocate(coordinate)
deallocate(new_coordi)

          28
STEP:           1
C      7.273878      8.481164      4.787954
C      8.285965      2.209745      9.083592
C     10.699671      5.799165      2.790090
C      9.290149      5.432365      7.050378
H      6.812623      8.515214      5.770640
………………………………………………………………

          224
STEP:          1
C      7.273878      8.481164      4.787954
C      7.273878      8.481164     13.306053
          224
STEP:          2
C      7.273347      8.482496      4.788510
………………………………………………………………

do i1=1,atomnum
a_count=1
do while(a_count.le.na)
   new_coordi(i1,1,a_count)=coordinate(i1,1)+a_count*a
     b_count=1
     do while(b_count.le.nb)
     new_coordi(i1,2,b_count)=coordinate(i1,2)+b_count*b
       c_count=1
       do while(c_count.le.nc)
       new_coordi(i1,3,c_count)=coordinate(i1,3)+c_count*c
       write(18,"(1XA1,3(F14.6))")label(i1),new_coordi(i1,1,a_count),new_coordi(i1,2,b_count),new_coordi(i1,3,c_count)
       c_count=c_count+1
       end do
     b_count=b_count+1
     end do
   a_count=a_count+1
end do   
end do