| 查看: 164 | 回复: 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之间 程序的语法好像没有错误，估计是逻辑不对，请虫友帮忙看看错在哪里了？我真挺着急的，万分感谢了！ |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有5人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有5人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复