程序如下：
subroutine generate_start_config(Np,box,T,x,v)
        implicit none
        integer maxNp
        parameter(maxNp=500000)
        integer Np
        real*8  x(maxNp,3),v(maxNp,3),T,box(3)
        real*8  invbox(3),vec(3),dist2,mindist2,mdtime
        common /mymdtime/ mdtime
        integer i,j,k,isbound(maxNp,2),N(maxNp)
        logical overlap
        mindist2 = 0.85**2                                        ! EDIT SETTING
        mdtime   = 0.D0
        do k=1,3; invbox(k)=1.D0/box(k); enddo
        do k=1,3; x(1,k) = (rand(100000)-0.5D0)*box(k); enddo
        !do k=1,3; x(1,k) = (rand()-0.5D0)*box(k); enddo
        do i=2,Np
         if (mod(i,max(1,Np/10)).eq.1) then
          print *,'inserted particle ',i,' of ',Np
         end if
         overlap = .true.
         do while (overlap)
          overlap = .false.
          do k=1,3; x(i,k) = (rand()-0.5D0)*box(k); enddo
          do j=1,i-1
           do k=1,3; vec(k) = x(i,k)-x(j,k); enddo
           call boundary(box,invbox,vec)
           dist2 = vec(1)**2+vec(2)**2+vec(3)**2
           if (dist2.lt.mindist2) then
            overlap = .true.
            goto 1
           end if
          enddo
1          continue
         enddo               
        enddo
        ! create velocities
        do i=1,Np
        do k=1,3
         v(i,k)=-6.D0
         do j=1,12; v(i,k) = v(i,k) + rand(); enddo
        enddo
        enddo
        call Tcontrol(Np,x,v,T,0.D0)
        ! create empty isbound and N
        do i=1,Np
         N(i)         = 1
         isbound(i,1) = 0
         isbound(i,2) = 0
        enddo
        ! writing files
        call writing_config(Np,Np,N,box,x,v)
        call writing_isbound(Np,isbound)
100     format(A50,10(I10,1x))
101     format(A50,10(F10.4,1x))
       
        return
        end

编译出错：
Error: This name does not have a type, and must have an explicit type.   [RAND]
Rand()表示的是生成0~1之间的随机数。不知道如何改进，请帮忙，谢谢 返回小木虫查看更多