GEANT321/gdraw/gdcirr.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:22  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.26  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GDCIRR(R,THROT)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *       Draw a circle of radius R, centered on the origin,       *
  17 C.    *       on a plane whose normal is parallel to X-Y plane         *
  18 C.    *       with an angle THROT with X-axis                          *
  19 C.    *                                                                *
  20 C.    *    ==>Called by : GDRAWS                                       *
  21 C.    *       Author : P.Zanarini   *********                          *
  22 C.    *                                                                *
  23 C.    ******************************************************************
  24 C.
  25 #include "geant321/gcdraw.inc"
  26 #include "geant321/gconsp.inc"
  27       DIMENSION X(3,41),U(41),V(41)
  28 C.
  29 C.    ------------------------------------------------------------------
  30 C.
  31       IF (R.LT.0.0001) GO TO 999
  32 C
  33       GTHROT=THROT*DEGRAD
  34       STHROT=SIN(GTHROT)
  35       CTHROT=COS(GTHROT)
  36 C
  37       DO 10 I=1,40
  38          X(1,I)=R*GCOS(I)*CTHROT
  39          X(2,I)=R*GSIN(I)
  40          X(3,I)=R*GCOS(I)*STHROT
  41    10 CONTINUE
  42 C
  43       CALL GDFR3D(X,40,U,V)
  44 C
  45 CM             Short circuit in case of GMR
  46 C
  47       IF ( IGMR .NE. 0 ) GO TO 999
  48 CM
  49       U(41)=U(1)
  50       V(41)=V(1)
  51       CALL GDRAWV(U,V,41)
  52 C
  53   999 RETURN
  54       END