Line data Source code
1 : MODULE m_points 2 : CONTAINS 3 0 : SUBROUTINE points(x,n) 4 : c ********************************************************* 5 : c generate random points, in internal coordinates, 6 : c within the unit cell omega-tilda 7 : c ********************************************************* 8 : USE m_qranf 9 : IMPLICIT NONE 10 : C .. Scalar Arguments .. 11 : INTEGER n 12 : C .. 13 : C .. Array Arguments .. 14 : REAL x(3,n) 15 : C .. 16 : C .. Local Scalars .. 17 : REAL r 18 : INTEGER i,j 19 : C .. 20 : 21 : C .. 22 : C .. Intrinsic Functions .. 23 : INTRINSIC sqrt 24 : C .. 25 0 : r = sqrt(13.) 26 0 : j = 1 27 0 : DO i = 1,n 28 0 : x(1,i) = qranf(r,j) 29 0 : x(2,i) = qranf(r,j) 30 0 : x(3,i) = qranf(r,j) - 0.5 31 : ENDDO 32 0 : RETURN 33 : END SUBROUTINE 34 : END