GEANT321/geocad/gckpar.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:47  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.32  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GCKPAR (ISH,NPAR,PAR)
  13 *
  14 #include "geant321/gconsp.inc"
  15 *
  16       DIMENSION PAR(50)
  17 *
  18       IF (ISH.EQ.4) THEN
  19 * ***    trap
  20          NPAR = 11
  21          PH = 0.
  22          IF (PAR(2).NE.0.) PH = ATAN2(PAR(3),PAR(2))*RADDEG
  23          TT = SQRT(PAR(2)**2+PAR(3)**2)
  24          PAR(2) = ATAN(TT)*RADDEG
  25          IF (PH.LT.0.0) PH = PH + 360.0
  26          PAR(3) = PH
  27          PAR(7) = ATAN(PAR(7))*RADDEG
  28          IF (PAR(7).GT.90.0) PAR(7) = PAR(7)-180.0
  29          PAR(11)= ATAN(PAR(11))*RADDEG
  30          IF (PAR(11).GT.90.0) PAR(11) = PAR(11)-180.0
  31 *
  32       ELSE IF (ISH.EQ.10) THEN
  33 * ***    para
  34          PH = 0.
  35          IF (PAR(5).NE.0.) PH = ATAN2(PAR(6),PAR(5))*RADDEG
  36          TT = SQRT(PAR(5)**2+PAR(6)**2)
  37          PAR(4) = ATAN(PAR(4))*RADDEG
  38          IF (PAR(4).GT.90.0) PAR(4) = PAR(4)-180.0
  39          PAR(5) = ATAN(TT)*RADDEG
  40          IF (PH.LT.0.0) PH = PH + 360.0
  41          PAR(6) = PH
  42 *
  43       ELSE IF (ISH.EQ.28) THEN
  44 * ***    gtra
  45          NPAR = 12
  46       ENDIF
  47 *
  48       END