| 查看: 399 | 回复: 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之间 程序的语法好像没有错误，估计是逻辑不对，请虫友帮忙看看错在哪里了？我真挺着急的，万分感谢了！ |
回复此楼» 猜你喜欢
献血感触
已经有4人回复
-
评审有感
已经有16人回复
-
反应很差，大量原料没有反应
已经有3人回复
-
Sci. Bull. 悲剧经验
已经有4人回复
-
找博士生导师
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有11人回复
-
26/27申博自荐-锂/钠电池方向
已经有4人回复
-
同样的基金本子，换个专家直接从C变A！
已经有3人回复
-
别被青基扩招骗了！26年科研内卷才刚刚开始
已经有4人回复
-
26/27博士推荐
已经有4人回复
老虎大王
木虫 (著名写手)
- 应助: 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
5