5 * Revision 1.1.1.1 1995/10/24 10:20:53 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.29 by S.Giani
14 C. ******************************************************************
16 C. * Routine to extend the user supplied parameters P to *
17 C. * internal parameters for the trapezoidal shape TRAP *
19 C. * Called by : GSPOSP, GSVOLU , GGDPAR *
20 C. * Author R.Nierhaus ********* *
22 C. ******************************************************************
24 * GNOTR1 computes the coordinates for the 8 vertices
25 * X(8), Y(8) and Z(8).
26 * Then GNOTR2 is called for each of the six boundary surfaces.
28 DIMENSION P(35),X(8),Y(8),Z(8)
41 IF (HY1.LT.1.E-6) HY1=1.E-6
42 IF (HX11.LT.1.E-6) HX11=1.E-6
43 IF (HX12.LT.1.E-6) HX12=1.E-6
44 IF (HY2.LT.1.E-6) HY2=1.E-6
45 IF (HX21.LT.1.E-6) HX21=1.E-6
46 IF (HX22.LT.1.E-6) HX22=1.E-6
47 X(1)=-HZ*SX-HY1*SX1-HX11
48 X(2)=-HZ*SX-HY1*SX1+HX11
49 X(3)=-HZ*SX+HY1*SX1-HX12
50 X(4)=-HZ*SX+HY1*SX1+HX12
51 X(5)=+HZ*SX-HY2*SX2-HX21
52 X(6)=+HZ*SX-HY2*SX2+HX21
53 X(7)=+HZ*SX+HY2*SX2-HX22
54 X(8)=+HZ*SX+HY2*SX2+HX22
71 CALL GNOTR2(X,Y,Z,1,2,4,3,P(12))
72 CALL GNOTR2(X,Y,Z,5,7,8,6,P(16))
73 CALL GNOTR2(X,Y,Z,1,5,6,2,P(20))
74 CALL GNOTR2(X,Y,Z,3,4,8,7,P(24))
75 CALL GNOTR2(X,Y,Z,1,3,7,5,P(28))
76 CALL GNOTR2(X,Y,Z,2,6,8,4,P(32))