本帖最后由 15235312522 于 2021-3-8 10:59 编辑
我在现有程序里,加了如下代码,通过冒泡法,把主程序里的数据从小到大排列,求得最大和最小值的距离并输出到“RDF”文件夹里,“nDpEnd-nWallAtom”、“r(n,i)”是主程序里的
[Fortran] 纯文本查看 复制代码 subroutine suan_chang_du()
implicit none
integer :: N
integer N,I,J,TEMP
real chang_du,A(N)
open(0306, file = './data/RDF/pu_zhan.plt',status = 'unknown')
N = nDpEnd-nWallAtom
A(N) = r(n,i)
do I=N-1,1,-1
do J=1,I
if(A(J)>A(J+1))then
TEMP=A(J)
A(J)=A(J+1)
A(J+1)=TEMP
end if
end do
end do
chang_du = abs(A(N)-A(1))
write(0306,'(''Variables= "chang_du"'')') chang_du
return
end subroutine
|