/[PAMELA software]/gpamela/gpspe/gpspev.F
ViewVC logotype

Diff of /gpamela/gpspe/gpspev.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.2 by cafagna, Tue Jun 21 02:42:27 2005 UTC revision 3.7 by cafagna, Fri Dec 16 10:20:23 2005 UTC
# Line 1  Line 1 
1    
2  *  *
3  * $Id: gpspev.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpspev.F,v 3.6 2005/12/16 09:13:59 cafagna Exp $
4  *  *
5  * $Log: gpspev.F,v $  * $Log: gpspev.F,v $
6  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna  * Revision 3.6  2005/12/16 09:13:59  cafagna
7  * First GPAMELA release on CVS  * New small valume added to the tracker frame
8  *  *
9  *  * Revision 3.5  2005/12/13 10:31:29  pam-ba
10  *CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola  * added a new volume, TPGI, in the spectrometer
11  *CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna  *
12  *CMZ :  1.02/00 07/03/97  14.48.41  by  Francesco Cafagna  * Revision 3.4  2005/12/05 12:15:21  pam-ba
13  *CMZ :  1.00/02 06/02/96  17.34.50  by  Francesco Cafagna  * new spectrometer geometry and internal magnetic field
14  *-- Author :    Francesco Cafagna   09/12/95  *
15        SUBROUTINE GPSPEV  * Revision 3.3  2005/07/25 11:53:21  cafagna
16  ************************************************************************  * Several updates. See history for details
17  *                                                                      *  *
18  * Volumes definition for the silicon tracking system and magnet        *  * Revision 3.2  2005/06/21 02:42:27  cafagna
19  * Called by: GPGEM                                                     *  * Major modification to the geometry and to the random number chain
20  * Author: Francesco Cafagna, 05/12/95 13.39.41                         *  *
21  *                                                                      *  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
22  ************************************************************************  * First GPAMELA release on CVS
23  #include "gpgeo.inc"  *
24  #include "gpmed.inc"  *
25  *  *CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
26        INTEGER IROT,IVOLU,N,NMED,NUM  *CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna
27        REAL X,Y,Z  *CMZ :  1.02/00 07/03/97  14.48.41  by  Francesco Cafagna
28  *  *CMZ :  1.00/02 06/02/96  17.34.50  by  Francesco Cafagna
29  *  Define the SPEB     volume  *-- Author :    Francesco Cafagna   09/12/95
30  *        SUBROUTINE GPSPEV
31        NMED= MN2  ************************************************************************
32        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)  *                                                                      *
33  *  * Volumes definition for the silicon tracking system and magnet        *
34  *  Define the MGPL     volume  * Called by: GPGEM                                                     *
35  *  * Author: Francesco Cafagna, 05/12/95 13.39.41                         *
36        NMED= MMAG  *                                                                      *
37        CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU)  ************************************************************************
38  *  #include "gpgeo.inc"
39  *  Define the MGPI     volume  #include "gpmed.inc"
40  *  *
41        NMED= MN2        INTEGER IROT,IVOLU,N,NMED,NUM
42        CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU)        REAL X,Y,Z
43  *  *
44  *  Define the TRPB     volume  *  Define the SPEB     volume
45  *  *
46        NMED= MAL        NMED= MN2
47        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
48  *****  *
49  cC ml 20/04/05:  *  Define the MGFR     volume
50  *  *
51  cc*  Define the TRPL     volume        NMED= MAL
52  *  Define the TPAS and TPAI volume        CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
53  *  *
54        NMED= MN2  *  Define the MGPL     volume
55  CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)  *
56        CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)        NMED= MMAG
57        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)        CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU)
58  CC end ml.  *
59  ****************  *  Define the MGPI     volume
60  *  *
61  *  Define the TRSL     volume        NMED= MN2
62  *        CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU)
63  c ml: 26/11/04:  *
64  c      NMED= MN2  *  Define the TRPB     volume
65        NMED=MSIT  *
66        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)        NMED= MAL
67  *        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
68  * Define the THBP volume  *****
69  *  cC ml 20/04/05:
70        NMED=MAL  *
71        CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)  cc*  Define the TRPL     volume
72  C end ml.  *  Define the TPAS and TPAI volume
73  *  *
74  *  Define the TSPA     volume        NMED= MN2
75  *  CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
76        NMED= MSIT        CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
77        CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU)        NMED= MAL
78  *        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
79  *  Define the TRCP     volume  CC end ml.
80  *  ****************
81        NMED= MCF  *
82        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)  *  Define the TRSL     volume
83  *  *
84  * Positioning volume TSPA     into the mother TRSL  c ml: 26/11/04:
85  *  c      NMED= MN2
86        N= 1        NMED=MSIT
87        X= 0.        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
88  c ml: 26/11/66:  *
89  c      Y= TRSL(2) - TSPA(2)  * Define the THBP volume
90        Y=+TRSL(2)/2  *
91        Z= 0.        NMED=MG10C
92        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')        CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
93        N= 2  C end ml.
94  C      Y= Y - 2.*TSPA(2)  *
95        Y=-TRSL(2)/2  *  Define the TSPA     volume
96        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  *
97  c ml: 26/11/04:        NMED= MSIT
98  *        CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU)
99  c* Positioning volume TRCP     into the mother TRSL  *
100  CC ml 20/04/05:  *  Define the TRCP     volume
101  CC* Positioning volume TRCP into mother TRPL  *
102  * Positioning volume TRCP into mother TRPB        NMED= MCF
103  *        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
104  c      N= 1  *
105  CC ml 20/04/05:  *  Define the TBAL     volume
106  cc      Z= 0.  *
107  cc      Y= 0.        NMED= MAL
108        Z=-TRPB(3)+0.11+TRCP(3)        CALL GSVOLU('TBAL','BOX ',NMED,TBAL, 3,IVOLU)
109        Y=-TRPB(2)+1.5-0.225+TRCP(2)  *
110  CC end ml.  *  Define the MGPA     volume
111        DO I=1,4  *
112  CC ml 20/04/05:                NMED= MAL
113  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)        CALL GSVOLU('MGPA','BOX ',NMED,MGPA, 3,IVOLU)
114           X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  *
115  CC end ml  *  Define the TPGA     volume
116  c         CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')  *
117           CALL GSPOS('TRCP',I,'TRPB',X,Y,Z,0,'ONLY')        NMED= MN2
118        ENDDO          CALL GSVOLU('TPGA','BOX ',NMED,TPGA, 3,IVOLU)
119  c      X= -TRSL(1) + TRCP(1)  *
120  c      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')  *  Define the TPGI     volume
121  C      N= 2  *
122  C      X= -X        NMED= MN2
123  C      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')        CALL GSVOLU('TPGI','BOX ',NMED,TPGI, 3,IVOLU)
124  *****  *
125  cc ml 20/04/05:  *  Define the TPGU     volume
126  *  *
127  *** Positioning volume TRSL     into the mother TRPL        NMED= MN2
128  * Positioning volume TRSL     into the mother TRPB        CALL GSVOLU('TPGU','BOX ',NMED,TPGU, 3,IVOLU)
129  *  *
130  c ml: 26/11/04:  *  Define the TPGD     volume
131  C      Y= 0.  *
132  cc      Y=THBP(2)        NMED= MN2
133        Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)        CALL GSVOLU('TPGD','BOX ',NMED,TPGD, 3,IVOLU)
134  C      Y=2.75  
135  CC      Z= 0.  *
136        Z=-TRPB(3)+0.11+TRCP(3)  * Positioning volume TSPA     into the mother TRSL
137        DO I=1, NTRSL  *
138  c         X= -TRPL(1) + TRSL(1) +2.*(I-1)*TRSL(1)        N= 1
139  cc         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*TRSL(1)        X= 0.
140           X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*TRSL(1)        Y= 0.
141  cc         CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')        Z= 0.
142           CALL GSPOS('TRSL',I,'TRPB',X,Y,Z,0,'ONLY')        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
143        ENDDO  C #       N= 2
144  **********  C # C      Y= Y - 2.*TSPA(2)
145  cc ml: 20/04/05  C #       Y=-TRSL(2)/2
146  *  C #       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
147  ********* Positioning volume THBP     into the mother TRPL  c ml: 26/11/04:
148  * Positioning volume THBP     into the mother TRPB  * Positioning volume TRCP into mother TPAS
149  *  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
150  cc      Y=-TRSL(2)  C #       Z=-TRPB(3)+0.11+TRCP(3)
151        Y=-TRPB(2)+1.5+THBP(2)  C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
152  cc      Z= 0.  CC end ml.
153        Z=-TRPB(3)+0.11+TRCP(3)  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
154        DO I=1, NTHBP        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
155  CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)        y=-0.225-TPAS(2)+TRCP(2)
156           X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)  C # end CAF mod
157  cc         CALL GSPOS('THBP',I,'TRPL',X,Y,Z,0,'ONLY')        TPASSPACE=TPAS(1)-(3*TRSL(1)+4*TRCP(1)+3*XGLUE)
158           CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')        DO I=1,4
159        ENDDO  CC ml 20/04/05:        
160  c end ml.  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
161  cc ml: 20/04/05  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
162  *  C #          X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
163  * Positioning volume TPAI into the mother TRPB  c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.
164  *  c              for now the glue is a space of nitrogen of 0.015 mm
165        X=-TRPB(1)+2.29+TPAI(1)  c         X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
166        Y=-TRPB(2)+1.5+2.5+TPAI(2)  C Ml. 15/12/05: introducing a space (along x) between the inner part of the
167        Z=-TRPB(3)+TPAI(3)  c  aluminum frame and respectively the first and the last carbon fiber barrel.
168        CALL GSPOS('TPAI',I,'TRPB',X,Y,Z,0,'ONLY')  c It is has been done increasing the dimension of TPAS along x:
169  *  C         X=-TPAS(1)+2*(I-1)*XGLUE+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
170  * Positioning volume TPAS into the mother TRPB           X=-TPAS(1)+TPASSPACE+2*(I-1)*XGLUE+(2*I-1)*TRCP(1)+
171  *       +      2*(I-1)*TRSL(1)
172        X=-TRPB(1)+2.29+TPAS(1)  CC end ml
173        Y=-TRPB(2)+1.5+TPAS(2)           CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
174        Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)  C # end CAF mod
175        CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')        ENDDO  
176  c end ml.  cc ml 20/04/05:
177    * Positioning volume TRSL     into the mother TPAS
178  *  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
179  * Positioning volume MGPI     into the mother MGPL  C #       Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
180  *  C      Y=2.75
181        N= 1  CC      Z= 0.
182        X= 0.  C      Z=-0.005
183        Y= 0.  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
184        Z= 0.        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
185        CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')        DO I=1, NTRSL
186  C Ml: 20/04/05:  c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.
187  c*  c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
188  c* Positioning volume TRPL     into the mother TRPB  c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
189  c*           X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
190  c      N= 1       +   (2*(MOD(I-1,3)+1)-1)*XGLUE+(2*(MOD(I-1,3)+1)-1)*TRSL(1)        
191  c      Z= 0.           IF(I.LE.3) THEN
192  c      CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')  C ml: 22/9/05:
193  c end ml.  c            Y=TPAS(2)-TRSL(2)
194  *              Y=TPAS(2)-0.1485-TRSL(2)
195  * Positioning volumes TRPB & MGPL into the mother SPEB           ELSE
196  *  C            Y=TPAS(2)-3*TRSL(2)
197        N= 1              Y=TPAS(2)-0.1485-XGLUE-3*TRSL(2)
198        X= 0.           ENDIF
199        Y= 0.           CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
200        Z= SPEB(3) - TRPB(3)        ENDDO
201        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  C # end CAF mod
202        DO I=1, 5  **********
203           Z= Z - TRPB(3) - MGPL(3)  cc ml: 20/04/05
204           N= I  * Positioning volume THBP     into the mother TPAS
205           CALL GSPOS('MGPL',N,'SPEB',X,Y,Z,0,'ONLY')  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
206           Z= Z - MGPL(3) - TRPB(3)  C #       Y=-TRPB(2)+1.5+THBP(2)
207           N= (I+1)  C # cc      Z= 0.
208           CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  C #       Z=-TRPB(3)+0.11+TRCP(3)
209        ENDDO  c ml: 22/9/05:
210  *  c      Y=-TPAS(2)+THBP(2)
211        RETURN        Y=-TPAS(2)+0.1485+THBP(2)
212        END  C      Z=-0.005
213    c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
214          Z=TRPB(3)-2*TPAI(3)-TRCP(3)
215          DO I=1, NTHBP
216    CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)
217    C #          X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)
218    c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon
219    c       ladders
220    c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)
221             X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)+
222         +      (2*(MOD(I-1,3)+1)-1)*XGLUE
223             CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
224          ENDDO
225    c end ml.
226    cc ml: 20/04/05
227    *
228    * Positioning volume TPAI into the mother TPAS
229    *
230    C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
231          X=0
232          Y=-TPAS(2)+TPAI(2)
233    c      Z=-TPAS(3)+TPAI(3)
234          Z=TRPB(3)-TPAI(3)
235    C #       X=-TRPB(1)+2.29+TPAI(1)
236    C #       Y=-TRPB(2)+1.5+2.5+TPAI(2)
237    C #       Z=-TRPB(3)+TPAI(3)
238          CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')
239    C # end CAF mod
240    *
241    * Positioning volume TPAS into the mother TRPB
242    *
243    C # Mod By Caf on 15th Jul 2005
244          X=0
245          Y=-TRPB(2)+1.5+TPAS(2)
246          Z=0
247    C #       X=-TRPB(1)+2.29+TPAS(1)
248    C #       Y=-TRPB(2)+1.5+TPAS(2)
249    C #       Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
250          CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
251    c end ml.
252    *
253    * Positioning volume MGPA     into the mother MGPL
254    *
255          N= 1
256          X= 0.
257          Y= 0.
258          Z= 0.
259          CALL GSPOS('MGPA',N,'MGPL',X,Y,Z,0,'ONLY')
260    *
261    * Positioning volume TPGD     into the mother MGFR
262    *
263          X=0.
264          Z=-MGFR(3)+TPGD(3)
265          Y=0.
266          CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
267    *
268    * Positioning volume TPGI     into the mother MGFR under MGPL
269    *
270          N=1
271          X=0.
272          Z=-MGFR(3)+2*TPGD(3)+TPGI(3)
273          Y=0.
274          CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')
275    *
276    * Positioning volume MGPI     into the mother MGPL
277    *
278          N= 1
279          X= 0.
280          Y= 0.
281    CC ML 10/11/05:  Positioning volume MGPI     into the mother MGPA
282          Z= 0.
283          CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')
284    C      Z= -MGFR(3)+2*TPGD(3)+MGPI(3)
285    C      CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY')
286    *
287    * Positioning volume MGPL     into the mother MGFR
288    *
289    c ml: 10/11/05
290          N=1
291    c      Z=-MGFR(3)+MGPL(3)
292          Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3)
293          CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
294    *
295    * Positioning volume TPGI     into the mother MGFR above MGPL
296    *
297          N=2
298          X=0.
299          Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3))+TPGI(3)
300          Y=0.
301          CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')
302    *
303    * Positioning volume TPGA     into the mother MGFR
304    *
305    C ML: 10/11/05
306          N=1
307          X=0.
308          Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)
309    C      Z=-MGFR(3)+2*MGPL(3)+TPGA(3)
310    C      Z=MGFR(3)-2*TPGU(3)-2*TRPB(3)-TPGA(3)
311          Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3))+TPGA(3)
312    C      Y=-0.75
313          CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')
314    *
315    * Positioning volumes TRPB & MGPL into the mother SPEB
316    *
317          N= 1
318          X= 0.
319    C      Y= SPEB(2) - TRPB(2)
320          Y= -MGFR(2) + TRPB(2)
321    C      Z= SPEB(3) - TRPB(3)
322    c      Z=-MGFR(3)+2*MGPL(3)+2*TPGA(3)+TRPB(3)
323    C      Z=MGFR(3)-2*TPGU(3)-TRPB(3)
324          Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3))+TRPB(3)
325    C      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
326          CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
327    c      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY')
328          Z=SPEB(3)
329          DO I=1, 5
330             X= 0.
331             Y= 0.
332    C         Z= Z - TRPB(3) - MGFR(3)
333             Z=Z-MGFR(3)
334             N= I
335             CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
336    C         Y= SPEB(2) - TRPB(2)
337    C         Z= Z - MGFR(3) - TRPB(3)
338    C         N= (I+1)
339    C         IF(I.EQ.5) THEN
340    CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
341    C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,10,'ONLY')
342    C         ELSE
343    CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
344    C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
345    C         ENDIF
346             Z=Z-MGFR(3)
347          ENDDO
348          N=2
349          X=0.
350          Y=MGFR(2) - TRPB(2)
351          Z=Z-TRPB(3)
352          ZTRPB=Z
353    c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')
354          CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
355    *
356    * Positioning volume TBAL into the mother SPEB
357    *
358          DO I=1,2
359    C DTBAL is the distance along x of the volumes TBAL from the end of TRPB
360             DTBAL=3.6
361             N=I
362             IF(I.EQ.1) THEN
363                X=TRPB(1)-3.6-TBAL(1)
364             ELSE
365                X=-(TRPB(1)-3.6-TBAL(1))
366             ENDIF
367             Y=MGFR(2) - 2*TRPB(2)-TBAL(2)
368             Z=ZTRPB
369             CALL GSPOS('TBAL',N,'SPEB',X,Y,Z,0,'ONLY')
370          ENDDO
371    *
372    * Positioning volume TPGU     into the mother MGFR
373    *
374          X=0.
375    C      Z=MGFR(3)-TPGU(3)
376          Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3)+TRPB(3))+
377         +   TPGU(3)      
378          Y=-MGFR(2)+TRPB(2)
379          CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
380    *
381          RETURN
382          END
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.7
  ViewVC Help
Powered by ViewVC 1.1.23