24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2142  |  回复: 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 的主题更新
信息提示
请填处理意见