| 查看: 127 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】程序出错，寻求高手指点
|
|||
|
我编了一个子程序，在一系列坐标点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之间 程序的语法好像没有错误，估计是逻辑不对，请虫友帮忙看看错在哪里了？我真挺着急的，万分感谢了！ |
回复此楼» 猜你喜欢
一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂!
已经有4人回复
-
298求调剂
已经有4人回复
-
材料学硕333求调剂
已经有11人回复
-
0856求调剂
已经有7人回复
-
356求调剂
已经有3人回复
-
332求调剂
已经有7人回复
-
348求调剂
已经有4人回复
-
305求调剂
已经有8人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有13人回复
-
调剂310
已经有9人回复
10