| 查看: 464 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】程序出错，寻求高手指点
|
|||
|
我编了一个子程序，在一系列坐标点pp（np,mp）其中np是粒子的标号1，2，3......，mp代表x,y,z方向，想删除20个任意标号的坐标点，并且删除后把标号重新排列，我编的程序如下： INTEGER ITF,i1,j1,j2,j3 DOUBLE PRECISION x1,PP_temp(MP,ND) DO 75 i1 =NL+1 , NP PP_temp(i1,X) = PP(i1,X) PP_temp(i1,Y) = PP(i1,Y) PP_temp(i1,Z) = PP(i1,Z) PP(i1,X)=0 PP(i1,Y)=0 PP(i1,Z)=0 75 continue CALL RANDOM_seed() do 610 j1=1,20 call random_number(x1) ITF=NINT(x1*(NP-NL-1))+NL+1 print*,itf NP=NP-1 do 76 j2=NL+1,ITF-1 PP(j2, x) = PP_temp(j2,x) PP(j2, y) = PP_temp(j2,y) PP(j2, z) = PP_temp(j2,z) 76 continue do 77 j3=ITF,NP PP(j3, x) = PP_temp(j3+1,x) PP(j3, y) = PP_temp(j3+1,y) PP(j3, z) = PP_temp(j3+1,z) 77 continue 610 continue ITF是产生的随机编号，在NL+1~NP之间 程序的语法好像没有错误，估计是逻辑不对，请虫友帮忙看看错在哪里了？我真挺着急的，万分感谢了！ |
回复此楼» 猜你喜欢
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
40+中年人失业后，还能找到好工作吗？
已经有5人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★
jutun(金币+5,VIP+0):恩，我明白了，能不能再指点一下，该如何修改啊？ 5-11 17:02
jutun(金币+5,VIP+0):恩，我明白了，能不能再指点一下，该如何修改啊？ 5-11 17:02
|
do 76 j2=NL+1,ITF-1 PP(j2, x) = PP_temp(j2,x) PP(j2, y) = PP_temp(j2,y) PP(j2, z) = PP_temp(j2,z) 76 continue do 77 j3=ITF,NP PP(j3, x) = PP_temp(j3+1,x) PP(j3, y) = PP_temp(j3+1,y) PP(j3, z) = PP_temp(j3+1,z) 77 continue 610 continue 问题出在这里。你每一次都从PP_temp这个数组出发来更新PP数组，但是PP_temp并没有更新。 比如，你的PP_temp数组（我们简化一下，假设它是个1维数组），比如，第一个数是1，第二个数是2，第三个数是3。。。第N个数是n。现在： 1. 第一次，j1=1时，假设随机数是5，就是删除第五个数，那么按你的程序，执行过后，在你的PP数组里，应有： 第一个数是1，第二个数是2，第三个数是3，第四个数是4，第五个数是6，第六个数是7。。。。。。但你的PP_TEMP数组没有变化。 2. 第二次，j1＝2时，假设随机数是3，就是再删除第三个数，那么按你的程序，执行过后，在你的PP数组里，应有： 第一个数是1，第二个数是2，第三个数是4，第四个数是5，第五个数是6。。。。 问题来了吧。。呵呵。 [ Last edited by 老虎大王 on 2009-5-11 at 15:14 ] |
2楼2009-05-11 15:12:33
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
3楼2009-05-11 15:17:32
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2009-05-11 20:12:16
200
