| 查看: 115 | 回复: 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之间 程序的语法好像没有错误，估计是逻辑不对，请虫友帮忙看看错在哪里了？我真挺着急的，万分感谢了！ |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复