| 查看: 2631 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xk6891至尊木虫 (著名写手)
|
[交流]
【求助】读文本数据出现段错误[fortran] 已有2人参与
|
||
|
写的一段程序出现“段错误”，还请高手指点 code： program main implicit none real a,b,c,a1,b1,c1,alpha,beta,gamma character(len=10) today,now integer i,j character ele(768) real pos1(256,3) real pos2(256,3) real ala(256,3) character(len=20) rub1,rub2,rub3,rub4 pos1=1 pos2=1 ala=1 open(17,file='inname.dat') read(17,"(A6,3(3XF6.3),3(2XF5.2),A16)"  rub1,a1,b1,c1,alpha,beta,gamma,rub2do i=1,768 read(17,"(A30,3(F8.3),A23,A1)" rub3,(pos1(i,j),j=1,3),rub4,ele(i)ala(i,1)=pos1(i,1)/a1 ala(i,2)=pos1(i,2)/b1 ala(i,3)=pos1(i,3)/c1 end do close(17) open(23,file='OH_diople.car') write(23,"(A17)" "!BIOSYM archive 3"write(23,"(A6)" "PBC=ON"write(23,"(A35)" "Materials Studio Generated CAR File"CALL DATE_AND_TIME(date=today,time=now) write(23,"(A5,2XA8,2XA10)" "!DATE",today,nowwrite(*,*)"Input the value of lattice parameters:" write(*,*)"a=" read(*,*)a write(*,*)"b=" read(*,*)b write(*,*)"c=" read(*,*)c write(23,"(A3,6(2XF8.4),1XA4)" "PBC",a,b,c,alpha,beta,gamma,"(P1)"do i=1,768 if(mod(i,3)/=0)then pos2(i,1)=ala(i,1)*a pos2(i,2)=ala(i,2)*b pos2(i,3)=ala(i,3)*c if(mod(i,2)/=0)then write(23,"(A1,3(F15.9),1XA4,1XI1,7XA1,A1)" ele(i),pos2(i,j),"XXXX",1,ele(i),"+"else write(23,"(A1,3(F15.9),1XA4,1XI1,7XA1,A1)" ele(i),pos2(i,j),"XXXX",1,ele(i),"-"end if end if end do write(23,"(A3)" "end"write(23,"(A3)" "end"close(23) stop end program 要读的文件： CRYST1 20.700 20.700 20.700 90.00 90.00 90.00 P 1 1 ATOM 1 OH2 TIP3W 5 3.668 10.082 15.904 1.00 0.00 WT1 O ATOM 2 H1 TIP3W 5 3.224 10.451 15.101 1.00 0.00 WT1 H ATOM 3 H2 TIP3W 5 3.092 10.379 16.627 1.00 0.00 WT1 H ATOM 4 OH2 TIP3W 21 5.186 16.696 12.072 1.00 0.00 WT1 O ATOM 5 H1 TIP3W 21 6.083 16.979 11.816 1.00 0.00 WT1 H ATOM 6 H2 TIP3W 21 5.337 15.832 12.490 1.00 0.00 WT1 H ATOM 7 OH2 TIP3W 44 3.786 16.318 1.310 1.00 0.00 WT1 O ATOM 8 H1 TIP3W 44 3.781 15.388 1.196 1.00 0.00 WT1 H …………………………………………………………………………………………………………………… …………………………………………………………………………………………………………………… …………………………………………………………………………………………………………………… 只有一个提示“段错误”，一天了没找到错误原因 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 本文档是我阅读过的最全的英语“葵花宝典”，里面包含各种英语书籍资料~~~
已经有563人回复
- 读取data文件时候出现如下错误：Error: Error: in read_thread_svar!
已经有11人回复
- 【求助】FORTRAN读取文本文件
已经有9人回复
xk6891
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 13735.8
- 散金: 379
- 红花: 9
- 帖子: 1514
- 在线: 546小时
- 虫号: 890342
- 注册: 2009-11-01
- 专业: 理论和计算化学
3楼2011-04-24 09:20:32
whl2dxl
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2115.6
- 红花: 3
- 帖子: 151
- 在线: 115小时
- 虫号: 1267846
- 注册: 2011-04-17
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 谢谢回复！ 2011-04-24 18:00:46
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 谢谢回复！ 2011-04-24 18:00:46
|
您好，你的这几行是给变量付初值吗？ pos1=1 pos2=1 ala=1 如果是的话，这个程序就有点问题了 因为上面几个变量你定义的都是数组，要给数组付初值的话，就要把所有的都付一遍 ， 例如： do i=1,256 do j=1,3 pos1(i,j)=1 pos2(i,j)=1 ala(i,j)=1 end do end do 还有， 你读入文件时的格式都正确吗？ [ Last edited by whl2dxl on 2011-4-24 at 01:10 ] |
2楼2011-04-24 01:09:04
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2011-04-24 11:08:37
xk6891
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 13735.8
- 散金: 379
- 红花: 9
- 帖子: 1514
- 在线: 546小时
- 虫号: 890342
- 注册: 2009-11-01
- 专业: 理论和计算化学
5楼2011-04-24 12:42:40
28