]> git.uio.no Git - u/mrichter/AliRoot.git/blame - GEANT321/ggeom/gnotr1.F
Cleaned up some more in AliL3MemHandler, by removing the AliL3Transform in
[u/mrichter/AliRoot.git] / GEANT321 / ggeom / gnotr1.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:20:53 cernlib
6* Geant
7*
8*
9#include "geant321/pilot.h"
10*CMZ : 3.21/02 29/03/94 15.41.29 by S.Giani
11*-- Author :
12 SUBROUTINE GNOTR1(P)
13C.
14C. ******************************************************************
15C. * *
16C. * Routine to extend the user supplied parameters P to *
17C. * internal parameters for the trapezoidal shape TRAP *
18C. * *
19C. * Called by : GSPOSP, GSVOLU , GGDPAR *
20C. * Author R.Nierhaus ********* *
21C. * *
22C. ******************************************************************
23C.
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.
27*
28 DIMENSION P(35),X(8),Y(8),Z(8)
29*
30 HZ =P(1)
31 SX =P(2)
32 SY =P(3)
33 HY1 =P(4)
34 HX11=P(5)
35 HX12=P(6)
36 SX1 =P(7)
37 HY2 =P(8)
38 HX21=P(9)
39 HX22=P(10)
40 SX2 =P(11)
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
55 Y(1)=-HZ*SY-HY1
56 Y(2)=-HZ*SY-HY1
57 Y(3)=-HZ*SY+HY1
58 Y(4)=-HZ*SY+HY1
59 Y(5)=+HZ*SY-HY2
60 Y(6)=+HZ*SY-HY2
61 Y(7)=+HZ*SY+HY2
62 Y(8)=+HZ*SY+HY2
63 Z(1)=-HZ
64 Z(2)=-HZ
65 Z(3)=-HZ
66 Z(4)=-HZ
67 Z(5)=+HZ
68 Z(6)=+HZ
69 Z(7)=+HZ
70 Z(8)=+HZ
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))
77*
78 END