24小时热门版块排行榜    

查看: 2619  |  回复: 24

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by snoopyzhao at 2012-12-16 22:40:57
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:


program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b ...

您好,我又来找您帮忙了:我写了下面一个fortran的程序。也用您原来教我的write(*,*)查看了,查看了里面的ispin的确是2,因为在91行的循环中我用了i=1,ispin。可是怎么每次i只能等于1,而不能i=2啊,也就是91行的循环怎么不起作用啊?
多次麻烦您,谢谢您了!


dos-procar1.f90(5.78KB)
http://kuai.xunlei.com/d/LnU4DYxJXv3aUAQA6ad?p=130497


DOSCAR(810.18KB)
http://kuai.xunlei.com/d/LnU4DZBJsv3aUAQA876?p=130497


EIGENVAL(27.2KB)
http://kuai.xunlei.com/d/LnU4DZZJ-P3aUAQA9ec?p=130497


PROCAR(1.69MB)
http://kuai.xunlei.com/d/LnU4DZhJHv7aUAQA1d3?p=130497
21楼2012-12-26 21:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by snoopyzhao at 2012-12-16 22:40:57
查一下 Fortran 中有关“内部文件”的用法:

这里举一个小例子:


program test
character(len = 27) :: a = '1.0 2.0 3.0 4.0 5.0 6.0 7.0'
real, dimension(7) :: b

read(a, *) b(1:7)
write(*,*) b ...

哎,问题终于找到了。
我定义了phase这个矩阵,但是这个矩阵有的时候要用到,有的时候用不到。
定了allocatable的phase矩阵后,还必须得allocate phase矩阵。
allocate phase矩阵后可以不用,但如果不用,就不能用deallocate phase矩阵。
deallocate的矩阵必须是用过的矩阵。

呵呵,以上是我摸索的结果。谢谢您啦,因为有您的帮助,才让我有信心去开始编程!
22楼2012-12-26 22:45:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
22楼: Originally posted by Gina88 at 2012-12-26 22:45:27
哎,问题终于找到了。
我定义了phase这个矩阵,但是这个矩阵有的时候要用到,有的时候用不到。
定了allocatable的phase矩阵后,还必须得allocate phase矩阵。
allocate phase矩阵后可以不用,但如果不用,就不能 ...

如果不能 deallocate 会有错误提示的。根据这个去修改应该不是很难。

不必客气,编程的进步很多时候就是从 debug 开始的。加油……
23楼2012-12-27 07:08:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
23楼: Originally posted by snoopyzhao at 2012-12-27 07:08:53
如果不能 deallocate 会有错误提示的。根据这个去修改应该不是很难。

不必客气,编程的进步很多时候就是从 debug 开始的。加油……...

呵呵,编译的时候没有任何错误。
可是运行的时候总有fault的提示,但是并不指明是哪行。
然后我用您以前教我的
write(*,*)一步一步的去试的。
最后试出来上面的错误。
前面进行了好多这样的试,但是找不到错误,就上来问了。
后来我干脆一行一行的试,呵呵,后来就成功。
看来调试的时候还得仔细,不能浮躁啊。
24楼2012-12-27 10:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
24楼: Originally posted by Gina88 at 2012-12-27 10:06:45
呵呵,编译的时候没有任何错误。
可是运行的时候总有fault的提示,但是并不指明是哪行。
然后我用您以前教我的
write(*,*)一步一步的去试的。
最后试出来上面的错误。
前面进行了好多这样的试,但是找不到 ...

你所提及的问题的确只是在运行时才会出现。

如果你会用一些 debugger 的话,会方便一些,呵呵……

是的, 调试程序比写程序还累人……
25楼2012-12-27 13:38:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 22:16 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-16 21:39 by MHf40yA0Uhyb
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-16 21:26 by MHf40yA0Uhyb
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 yaCh4X3wd045 2026-09-14 3/150 2026-09-16 20:39 by MHf40yA0Uhyb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-16 18:39 by rRvxEj3zorCy
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-16 18:00 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:51 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 15:12 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
信息提示
请填处理意见