24小时热门版块排行榜    

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

wanghao098

金虫 (正式写手)

[求助] 请教:挑错误

CODE:
open (16,file="x.txt")
        nstep= int ( wk /nsec)
        print *, "每一段元素的个数 nstep=" , nstep
        do i=1,nsec
                b(i)=0.d0
                do j=(i-1)*nstep+1,i*nstep
                        b(i)=b(i)+clshuzu(j) / nstep
                enddo
        enddo
        write(16,"(F15.9)")(b(i),i=1,nsec)




为什么我这一段程序输出来的结果都是0呢?
clshuzu(wk) 是wk个元素的实数数组 0.2857143   
  0.1111111   
  0.1111111   
  0.3333333   
  0.2857143   
  0.2222222   
  0.2000000   
  9.0909094E-02
  0.1250000   
  0.1111111   
  0.2222222   
  0.1818182   
  0.1111111   
  0.0000000E+00
  0.0000000E+00
  0.2000000   
  0.1818182   
  0.2222222   
  9.0909094E-02
  0.0000000E+00
  0.1666667   
  9.0909094E-02
  0.1111111   
  0.2222222

[ Last edited by jjdg on 2013-7-22 at 00:48 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

做一个坚强的弱女子
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanghao098

金虫 (正式写手)

引用回帖:
2楼: Originally posted by pippi6 at 2013-07-22 08:16:07
你首先需要给出 wk 和 nsec,看看 nstep到底多大。如果 nsec=0,那你的 do j=(i-1)*nstep+1,i*nstep 也许根本就没有执行,那 b(i) 当然就是零了。让别人查错,需要给全信息。

已解决,是我前面的代码就有问题!谢谢
做一个坚强的弱女子
7楼2013-07-22 11:58:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
你首先需要给出 wk 和 nsec,看看 nstep到底多大。如果 nsec=0,那你的 do j=(i-1)*nstep+1,i*nstep 也许根本就没有执行,那 b(i) 当然就是零了。让别人查错,需要给全信息。
2楼2013-07-22 08:16:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanghao098

金虫 (正式写手)

nsec 是一个整形的常数
integer,parameter:: nsec=10
wk 也是一个整形的数
做一个坚强的弱女子
3楼2013-07-22 08:47:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
3楼: Originally posted by wanghao098 at 2013-07-22 08:47:09
nsec 是一个整形的常数
integer,parameter:: nsec=10
wk 也是一个整形的数

wk=?你还是没说啊
再者b 怎么定义的?
4楼2013-07-22 09:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:25 by Aj1rhIDL5ixY
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 16:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见