GEANT321/geocad/gpaset.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:46  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 GPASET(PAR)
  13 C-
  14 C-   Created  12-JUN-1991   Nils Joar Hoimyr
  15 C-   Modified 21.02.1992  Jouko Vuoskoski
  16 C-
  17 C-  Calculates the coordinates of the corner points of a GEANT PARA  shape
  18 C-  from the shape parameters.   The end faces are generated by drawing
  19 C-  connected polylines between the corner points at each Z-end of the shape.
  20 C-  A ruled solid is then defined between the 2 faces, and written out in SET.
  21 C---------------------------------------------------------
  22 C
  23 #include "geant321/gcsetf.inc"
  24
  25 C
  26       DIMENSION PAR(50)
  27       REAL P1X,P1Y,P1Z,P2X,P2Y,P2Z,P3X,P3Y,P3Z,P4X,P4Y,P4Z
  28       REAL P5X,P5Y,P5Z,P6X,P6Y,P6Z,P7X,P7Y,P7Z,P8X,P8Y,P8Z
  29
  30 C
  31 C----------------------------------------------------------------------
  32 C
  33       DX=PAR(1)
  34       DY=PAR(2)
  35       DZ=PAR(3)
  36       TXY=PAR(4)
  37       TXZ=PAR(5)
  38       TYZ=PAR(6)
  39 C       Calculation of corner points (verify PHI/THETA coordinates,
  40 C       I.e. TXZ AND TYZ should possibly change place)
  41 C
  42 C       FACE 1:
  43       P1X= -DX*COS(TYZ)-DY*SIN(TXY)
  44       P1Y= -DY*COS(TXY)
  45       P1Z= -DZ*COS(TXZ)-DY*SIN(TXY)
  46 C
  47       P2X= DX*COS(TYZ)-DY*SIN(TXY)
  48       P2Y= -DY*COS(TXY)
  49       P2Z= -DZ*COS(TXZ)-DY*SIN(TXY)
  50 C
  51       P3X= DX*COS(TYZ)-DY*SIN(TXY)
  52       P3Y= DY*COS(TXY)
  53       P3Z= -DZ*COS(TXZ)+DY*SIN(TXY)
  54 C
  55       P4X= DX*COS(TYZ)-DY*SIN(TXY)
  56       P4Y= DY*COS(TXY)
  57       P4Z= -DZ*COS(TXZ)+DY*SIN(TXY)
  58 C
  59 C      FACE 2:
  60 C
  61       P5X= -DX*COS(TYZ)+DY*SIN(TXY)
  62       P5Y= -DY*COS(TXY)
  63       P5Z= DZ*COS(TXZ)-DY*SIN(TXY)
  64 C
  65       P6X= DX*COS(TYZ)+DY*SIN(TXY)
  66       P6Y= -DY*COS(TXY)
  67       P6Z= DZ*COS(TXZ)-DY*SIN(TXY)
  68 C
  69       P7X= DX*COS(TYZ)+DY*SIN(TXY)
  70       P7Y= DY*COS(TXY)
  71       P7Z= DZ*COS(TXZ)+DY*SIN(TXY)
  72 C
  73       P8X= -DX*COS(TYZ)+DY*SIN(TXY)
  74       P8Y= DY*COS(TXY)
  75       P8Z= DZ*COS(TXZ)+DY*SIN(TXY)
  76 C
  77 C       DESCRIBED BY 2 FACES represented by polylines
  78 C       connected by the corners of the shape.
  79 C *WRITE SET@103,F1#3,3,2,P1X,P1Y,P1Z,P2X,P2Y,P2Z,P3X,P3Y,P3Z,P4X,P4Y,P4Z
  80 C *WRITE SET@103,F2#3,3,2,P5X,P5Y,P5Z,P6X,P6Y,P6Z,P7X,P7Y,P7Z,P8X,P8Y,P8Z
  81 C *WRITE SET @100,N#145,!F1,!F2
  82 C
  83 C----------------------------------------------------------------------
  84 C
  85       WRITE(BLKSTR,10000)N1,P1X,P1Y,P1Z,P2X,P2Y,P2Z, P3X,P3Y,P3Z,P4X,
  86      +P4Y,P4Z
  87       CALL GJWRIT
  88       N1=N1+1
  89       WRITE(BLKSTR,10000)N1,P5X,P5Y,P5Z,P6X,P6Y,P6Z, P7X,P7Y,P7Z,P8X,
  90      +P8Y,P8Z
  91       CALL GJWRIT
  92       N1=N1+1
  93       WRITE(BLKSTR,10100)N1,N1-2,N1-1
  94       CALL GJWRIT
  95 C
  96 10000   FORMAT('@103,',I10,',:5,2#3,3,2,',G14.7,',',G14.7,','
  97      +  ,G14.7,',',G14.7,',',G14.7,',',G14.7,',',G14.7,','
  98      +  ,G14.7,',',G14.7,',',G14.7,',',G14.7,',',G14.7)
  99 10100   FORMAT('@100,',I10,',:5,2#145,!',I10,',!',I10)
 100 C
 101       END