24小时热门版块排行榜    

查看: 2585  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中,共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见