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