24小时热门版块排行榜    

查看: 2042  |  回复: 5
【悬赏金币】回答本帖问题,作者xingxiao将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xingxiao

铜虫 (小有名气)

[求助] 运行结果为0,怎么回事?

我编了下面一段子程序,运行后结果输出0,不知道是怎么回事,存储空间肯定能存储这么大的数,而且没有越界:
     implicit none
     real,intent(in)::q
      integer,intent(in)::Ng
      integer,dimension(1:Ng),intent(in)::ink
      real,dimension(1:Ng,1:3000,1:7),intent(in)::AN
      integer,intent(inout)::Nele,Nion
      integer i,j
      nele=0.0;nion=0.0
      do i=1,Ng
             do j=1,ink(i)
            if(AN(i,j,1)<=0.0.and.q<0.0) then   !AN(i,j,1)是路程     !  
                Nele=nele+1
            else if(AN(i,j,1)<=0.0.and.q>0.0) then
              Nion=nion+1
          end if
        end do
      end do


我想实现的目的是统计粒子在一个空间边界处的数目,粒子在电场力的作用下向前走,每走一步(如果走到了边界处)就把粒子数目加1

[ Last edited by xingxiao on 2012-2-10 at 13:22 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingxiao

铜虫 (小有名气)

送鲜花一朵
引用回帖:
: Originally posted by snoopyzhao at 2012-02-10 16:39:18:
你把 an(i,j,1) 和 q 打印出来看看不就知道了……

这个程序的判断句中就有可能跳过一些条件的啊

应该是我程序中还有一些条件没考虑全,谢谢
5楼2012-02-11 08:44:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
jjdg(金币+1): 感谢参与 2012-02-11 09:03:01
你把 an(i,j,1) 和 q 打印出来看看不就知道了……

这个程序的判断句中就有可能跳过一些条件的啊

» 本帖已获得的红花(最新10朵)

2楼2012-02-10 16:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emanlee

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
jjdg(金币+1): 感谢参与 2012-02-11 09:03:07
先把Ng的值改得比较小,比如改为5。然后,单步执行程序,观察各个变量的值。找到问题之后,再改回Ng的原值。

» 本帖已获得的红花(最新10朵)

3楼2012-02-10 23:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingxiao

铜虫 (小有名气)

送鲜花一朵
引用回帖:
: Originally posted by emanlee at 2012-02-10 23:31:23:
先把Ng的值改得比较小,比如改为5。然后,单步执行程序,观察各个变量的值。找到问题之后,再改回Ng的原值。

感谢你的提醒,这是一个思路
4楼2012-02-11 08:43:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 变性人申请基金咋算性别? +22 leicahfs 2024-11-11 22/1100 2024-11-16 15:31 by 竹子宝
[论文投稿] 找不到审稿人 +7 吉祥如意141913 2024-11-13 7/350 2024-11-16 15:09 by 羊驼22
[论文投稿] 期刊选择 +3 上善若水 cc 2024-11-15 3/150 2024-11-16 14:23 by lx7991
[论文投稿] 小论文影响大论文查重吗? +4 阿道夫深度 2024-11-14 5/250 2024-11-16 14:20 by lx7991
[教师之家] 我感觉当老师好累啊。。 +19 fairy1122 2024-11-15 20/1000 2024-11-16 14:16 by lx7991
[论文投稿] 国内期刊审稿人数量 +4 新时代核动力驴 2024-11-13 6/300 2024-11-16 08:24 by pityzhou
[教师之家] 北大教授何怀宏曾如此描述他的同行 +13 zju2000 2024-11-09 13/650 2024-11-16 07:45 by 080520dxf
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 14/700 2024-11-15 22:15 by 仇恨毁灭财富
[论文投稿] chemical science和advanced science哪个好一点 +6 yly150 2024-11-12 6/300 2024-11-15 21:06 by hcj9909053
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +7 中南大学-金属材 2024-11-14 18/900 2024-11-15 20:22 by 中南大学-金属材
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[教师之家] 大学老师 +11 考研一路顺风 2024-11-13 14/700 2024-11-15 15:24 by LNP@mRNA
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样,文章目前被录用了? +4 babybabygo 2024-11-12 5/250 2024-11-14 19:13 by 走了002
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[找工作] 咨询一下江西的高校待遇,人文氛围怎么样? +5 akslis2024 2024-11-09 5/250 2024-11-14 13:53 by 啄木鸟、
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[硕博家园] 同步辐射 +4 小张要加油努力 2024-11-10 7/350 2024-11-11 13:47 by mpdfwxgui
信息提示
请填处理意见