24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 278  |  回复: 1

xk6891

至尊木虫 (著名写手)

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

CODE:
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)

用上面的一段代码(部分)去读下面的内容,希望能将坐标扩展一下,结果出现最后的结果,好像第一个do while做完之后就跳出了,外面的几个do循环都没有做。是因为二维数组与三维数组之间的转化么,搞不明白?求大侠指点
CODE:
          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
………………………………………………………………

CODE:
          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
………………………………………………………………

[ Last edited by 余泽成 on 2011-6-26 at 00:23 ]
回复此楼
人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)

jjdg: 恭喜 2011-06-25 11:19:08
知道哪里错了
每一次附初值应当放在每一个循环的开头,而不是放在总的开头
CODE:
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

人生自古多险阻,勤者自得天酬助。试问否泰何所依,枯藤老枝待新抽。临渊踌躇终迈步,振翅鹏起云霄冲。似是前程甚堪忧,他日振臂揽苍穹。
2楼2011-06-24 21:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xk6891 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010,总分320求调剂 +6 困困困困坤坤 2026-04-02 6/300 2026-04-02 21:28 by dongzh2009
[考研] 321求调剂 +12 y-yh 2026-04-01 14/700 2026-04-02 19:34 by 帕尔马拉特
[考研] 413求调剂 +3 柯某某 2026-03-31 3/150 2026-04-02 16:59 by zzsw+
[考研] 290求调剂085701 +3 1314捧花 2026-04-02 3/150 2026-04-02 13:34 by Ecowxq666!
[考研] 08开头看过来!!! +4 wwwwffffff 2026-03-31 6/300 2026-04-02 11:42 by 均值回归
[考研] 322求调剂 +5 熹僖XX 2026-03-31 6/300 2026-04-02 10:08 by 求调剂zz
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 英一数一408,总分284,二战真诚求调剂 +12 12.27 2026-03-30 14/700 2026-04-02 00:18 by 欣喜777
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-01 21:48 by chyhaha
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 350求调剂 +7 阿佳~ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 英一数二生物信息学287分,本科生物科学,求调剂 +3 碧水xyz 2026-03-29 4/200 2026-04-01 15:35 by plum
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-31 7/350 2026-04-01 14:50 by ZXlzxl0425
[考研] 085600,321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 329求调剂,一志愿西北工业大学,材料工程(085601) +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
[考研] 求调剂 +10 家佳佳佳佳佳 2026-03-29 10/500 2026-03-30 18:34 by 544594351
信息提示
请填处理意见