24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2437  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 emanlee at 2012-02-10 23:31:23:
先把Ng的值改得比较小,比如改为5。然后,单步执行程序,观察各个变量的值。找到问题之后,再改回Ng的原值。

感谢你的提醒,这是一个思路
4楼2012-02-11 08:43:59
已阅   回复此楼   关注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 snoopyzhao at 2012-02-10 16:39:18:
你把 an(i,j,1) 和 q 打印出来看看不就知道了……

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

应该是我程序中还有一些条件没考虑全,谢谢
5楼2012-02-11 08:44:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农微生物,288分,三年实验经历 +9 代fish 2026-04-09 9/450 2026-04-09 21:46 by bljnqdcc
[考研] 一志愿中南大学物理学,英一66,求调剂 +3 长烟旖旎 2026-04-08 4/200 2026-04-09 21:27 by wutongshun
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +4 Ezra_Zhang 2026-04-09 4/200 2026-04-09 20:59 by boxking200
[考研] 材料295 +15 小英11 2026-04-03 16/800 2026-04-09 15:44 by only周
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 292求调剂 +4 lilllllxccc 2026-04-05 5/250 2026-04-07 09:29 by 纺大杨老师
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
信息提示
请填处理意见