MINICERN/mathlib/gen/v/ran3d.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:56  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE RAN3D (X,Y,Z,XLONG)
  11 C         SUBROUTINE TO GENERATE RANDOM THREE-DIMENSIONAL VECTOR
  12 C         OF LENGTH XLONG.       F. JAMES, MAY, 1976
  13       DIMENSION RBUF(30)
  14       SAVE IBUF,RBUF
  15       DATA IBUF/30/
  16     1 CONTINUE
  17       IF (IBUF .LT. 30)  GO TO 10
  18       CALL NRAN(RBUF,30)
  19       IBUF = 1
  20    10 CONTINUE
  21       A = RBUF(IBUF) - 0.5
  22       B = RBUF(IBUF+1) - 0.5
  23       C = RBUF(IBUF+2) - 0.5
  24       IBUF = IBUF + 3
  25       R2=A**2+B**2+C**2
  26       IF(R2.GT.0.25) GO TO 1
  27       RINV = XLONG/SQRT(R2)
  28       X=A*RINV
  29       Y=B*RINV
  30       Z=C*RINV
  31       RETURN
  32       END