gpamela/gptof/gptofv.F

*
* $Id: gptofv.F,v 3.3 2006/04/09 23:29:22 cafagna Exp $
*
* $Log: gptofv.F,v $
* Revision 3.3  2006/04/09 23:29:22  cafagna
* Several new things, among this: ND and CARD
*
* Revision 3.2  2005/06/21 02:42:29  cafagna
* Major modification to the geometry and to the random number chain
*
* Revision 3.1.1.1  2002/07/11 16:02:01  cafagna
* First GPAMELA release on CVS
*
*
*CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
*CMZ :  2.00/00 03/03/2000  15.39.05  by  Francesco Cafagna
*CMZ :  1.02/00 21/10/99  16.30.25  by  Francesco Cafagna
*CMZ :  1.00/02 18/03/96  17.18.03  by  Francesco Cafagna
*-- Author :    Francesco Cafagna   05/12/95
      SUBROUTINE GPTOFV
************************************************************************
*                                                                      *
* Volumes definition for the scintillators counters                    *
* Called by: GPGEO                                                     *
* Author: Francesco Cafagna, 05/12/95 13.39.41                         *
*                                                                      *
************************************************************************
#include "gpgeo.inc"
#include "gpmed.inc"
*
      INTEGER IROT,IVOLU,N,NMED,NUM
      REAL X,Y,Z
*
      NMED=MAL
*
*  Define the S1     volume
*
      CALL GSVOLU('S1  ','BOX ',NMED,S1, 3,IVOLU)
*
*  Define the S2     volume
*
      CALL GSVOLU('S2  ','BOX ',NMED,S2, 3,IVOLU)
*
      NMED=MN2
*
*  Define the S1A     volume
*
      CALL GSVOLU('S1A ','BOX ',NMED,S1A, 3,IVOLU)
*
*  Define the S2A     volume
*
      CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)
*
*  Define the S3A     volume
*
      CALL GSVOLU('S3A ','BOX ',NMED,S3A, 3,IVOLU)
*
      NMED=MMYL
*
*  Define the S11M     volume
*
      CALL GSVOLU('S11M','BOX ',NMED,S11M, 3,IVOLU)
*
*  Define the S12M     volume
*
      CALL GSVOLU('S12M','BOX ',NMED,S12M, 3,IVOLU)
*
*  Define the S21M     volume
*
      CALL GSVOLU('S21M','BOX ',NMED,S21M, 3,IVOLU)
*
*  Define the S22M     volume
*
      CALL GSVOLU('S22M','BOX ',NMED,S22M, 3,IVOLU)
*
*  Define the S31M     volume
*
      CALL GSVOLU('S31M','BOX ',NMED,S31M, 3,IVOLU)
*
*  Define the S32M     volume
*
      CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)
*
*
*  Change medium just for S4.
*
      NMED=MSCIN
*
*  Define the S4     volume
*
      CALL GSVOLU('S4  ','BOX ',NMED,S4, 3,IVOLU)
*
*  Start subvolume definition.
*
*
*  Define the S11Y     volume
*
      CALL GSVOLU('S11Y','BOX ',NMED,S11Y, 3,IVOLU)
*
*  Define the S12X     volume
*
      CALL GSVOLU('S12X','BOX ',NMED,S12X, 3,IVOLU)
*
*  Define the S21X     volume
*
      CALL GSVOLU('S21X','BOX ',NMED,S21X, 3,IVOLU)
*
*  Define the S22Y     volume
*
      CALL GSVOLU('S22Y','BOX ',NMED,S22Y, 3,IVOLU)
*
*  Define the S31Y     volume
*
      CALL GSVOLU('S31Y','BOX ',NMED,S31Y, 3,IVOLU)
*
*  Define the S32X     volume
*
      CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)

c$$$<<<<<<< gptofv.F
c$$$=======
c$$$c ml: 21/12/04:
c$$$      NMED=MAL
c$$$*
c$$$*  Define the SC11     volume
c$$$*
c$$$      CALL GSVOLU('SC11','BOX ',NMED,SC11, 3,IVOLU)
c$$$*
c$$$*  Define the SC12     volume
c$$$*
c$$$      CALL GSVOLU('SC12','BOX ',NMED,SC12, 3,IVOLU)
c$$$*
c$$$*  Define the SC21     volume
c$$$*
c$$$      CALL GSVOLU('SC21','BOX ',NMED,SC21, 3,IVOLU)
c$$$*
c$$$*  Define the SC22     volume
c$$$*
c$$$      CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
c$$$
c$$$>>>>>>> 3.3
c*
c*  Define the STOF volume
c*      
c      NMED=MAL
c      CALL GSVOLU('STOF','BOX ',NMED,STOF, 3,IVOLU)
c*
c*  Define the SGLU volume
c*
cC COLLA ARALDITICA (DA VERIFICARNE LA COMPOSIZIONE)
c      NMED=MSIC
c      CALL GSVOLU('SGLU','BOX ',NMED,SGLU, 3,IVOLU)
c*
c*  Define the SPV1 volume
c*
c      NMED=? (PVC espanso)
c      CALL GSVOLU('SPV1','BOX ',NMED,SPV1, 3,IVOLU)
c*
c*  Define the SMYL volume
c*
c      NMED=? (MYLAR)
c      CALL GSVOLU('SMYL','BOX ',NMED,SMYL, 3,IVOLU)
c*
c*  Define the SPV2 volume
c*
c      NMED=? (PVC espanso)
c      CALL GSVOLU('SPV2','BOX ',NMED,SPV2, 3,IVOLU)
cc end ml.
*
* Positioning volumes
*
      X=0.
      Y=0.
      Z=0.
      CALL GSPOS('S1A ',N,'S1  ',X,Y,Z,0,'ONLY')
*
      X=0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
      DO I=1,NS12X
         N= I
         Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)
         CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
      ENDDO
      Y= 0.
      Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
      DO I=1,NS11Y
         N= I
         X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)
         CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
      ENDDO
C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')    
C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')    
*      
      X=0.
      Y=0.
      Z=0.
      CALL GSPOS('S2A ',N,'S2  ',X,Y,Z,0,'ONLY')    
*
      Y= 0.
      Z=-S2A(3)+HGLUTOF+S22M(3)
      DO I=1,NS22Y
         N= I
         X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1)        
         CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      X= 0.
      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
      DO I=1,NS21X
         N= I
         Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)
         CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
      ENDDO
C      Z=-S2A(3)+HGLUTOF+S22M(3)
C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')    
C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')    
*
      X= 0.
      Z=-S3A(3)+HAIRS3+S32M(3)
      DO I=1,NS32X
         N= I
         Y= -S3A(2)+ (2*I-1)*S32M(2)
         CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')
      ENDDO
*
      Y= 0.
      Z= -S3A(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
      DO I=1,NS31Y
         N= I
         X= -S3A(1)+(2*I-1)*S31M(1)
         CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')
      ENDDO
*
C      Z=-S3A(3)+HAIRS3+S32M(3)
C      CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')    
C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
C      CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')    
****************************************
c ml: 21/12/04:
      X=0.
      Y= 0.
      Z= 0.
      N=1
c ml: 27/10/04:
c      DO I=1,NS11Y
c         N= I
c         X= -S11M(1) + HMYLTOF + (2*I-1)*S11Y(1)
      CALL GSPOS('S11Y',N,'S11M',X,Y,Z,0,'ONLY')
c      ENDDO
c end ml.
*
      X= 0.
      Y=0.
      Z= 0.
      N=1
c ml: 27/10/04
C      DO I=1,NS12X
C         N= I
C         Y= -S11M(2) + HMYLTOF + (2*I-1)*S12X(2)
      CALL GSPOS('S12X',N,'S12M',X,Y,Z,0,'ONLY')
C      ENDDO
c end ml.
*
      X=0.
      y=0.
      Z= 0.
      N=1
C      DO I=1,NS21X
C         N= I
c      Y= -S21M(2) + HMYLTOF+ (2*I-1)*S21X(2)
      CALL GSPOS('S21X',N,'S21M',X,Y,Z,0,'ONLY')
C      ENDDO
*
      X=0.
      Y= 0.
      Z= 0.
      N=1
C      DO I=1,NS22Y
C         N= I
C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)        
      CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
C      ENDDO
c*
c*  Positioning volumes into S11:
c+
c      X=0.
c      Y=0.
c      Z=-S11(3)+SMYL(3)
c      N=1
c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
c      Z=Z+SMYL(3)+S11Y(3)
c      DO I=1,NS11Y
c         N= I
c         X= -S11(1) + (2*I-1)*S11Y(1)
c         CALL GSPOS('S11Y',N,'S11 ',X,Y,Z,0,'ONLY')
c      ENDDO
c      N=2
c      Z=Z+S11Y(3)+SMYL(3)
c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
c      Z=Z+SMYL(3)+SPV1(3)
c      N=1
c      CALL GSPOS('SPV1',N,'SPV1',X,Y,Z,0,'ONLY')
c      Z=Z+SPV1(3)+SGLU(3)
c      N=1
c      CALL GSPOS('SGLU',N,'SGLU',X,,Y,Z,0,'ONLY')
c      Z=Z+SGLU(3)+STOF(3)
c      N=1
c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
c*
c* Positioning volumes into S12:
c*
c      Z=-S12(3)+STOF(3)
c      N=2
c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
c      Z=Z+STOF(3)+SGLU(3)
c      CALL GSPOS('SGLU',N,'SGLU',X,Y,Z,0,'ONLY')
c      Z=Z+SGLU(3)+SPV2
      Y= 0.
      Z= -S31M(3)+HMYLTOF+S31Y(3)
      N=1
C      DO I=1,NS31Y
C         N= I
C         X= -S31M(1) + HMYLTOF+(2*I-1)*S31Y(1)
      X= -S31M(1) + 2*(HMYLTOF/2)+S31Y(1)
      CALL GSPOS('S31Y',N,'S31M',X,Y,Z,0,'ONLY')
C      ENDDO
*
      X= 0.
      Z= -S32M(3)+HMYLTOF+S32X(3)
      N=1
C      DO I=1,NS32X
C         N= I
C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S32X(2)
      Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
      CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
C      ENDDO
*
      RETURN
      END
1	cafagna	3.1	*
2	pam-ba	3.4	* $Id: gptofv.F,v 3.3 2006/04/09 23:29:22 cafagna Exp $
3	cafagna	3.2	*
4			* $Log: gptofv.F,v $
5	pam-ba	3.4	* Revision 3.3 2006/04/09 23:29:22 cafagna
6			* Several new things, among this: ND and CARD
7			*
8	cafagna	3.3	* Revision 3.2 2005/06/21 02:42:29 cafagna
9			* Major modification to the geometry and to the random number chain
10			*
11	cafagna	3.2	* Revision 3.1.1.1 2002/07/11 16:02:01 cafagna
12			* First GPAMELA release on CVS
13	cafagna	3.1	*
14			*
15			*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
16			*CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna
17			*CMZ : 1.02/00 21/10/99 16.30.25 by Francesco Cafagna
18			*CMZ : 1.00/02 18/03/96 17.18.03 by Francesco Cafagna
19			*-- Author : Francesco Cafagna 05/12/95
20			SUBROUTINE GPTOFV
21			************************************************************************
22			* *
23			* Volumes definition for the scintillators counters *
24			* Called by: GPGEO *
25			* Author: Francesco Cafagna, 05/12/95 13.39.41 *
26			* *
27			************************************************************************
28			#include "gpgeo.inc"
29			#include "gpmed.inc"
30			*
31			INTEGER IROT,IVOLU,N,NMED,NUM
32			REAL X,Y,Z
33			*
34	pam-ba	3.4	NMED=MAL
35			*
36			* Define the S1 volume
37			*
38			CALL GSVOLU('S1 ','BOX ',NMED,S1, 3,IVOLU)
39			*
40			* Define the S2 volume
41			*
42			CALL GSVOLU('S2 ','BOX ',NMED,S2, 3,IVOLU)
43			*
44	cafagna	3.1	NMED=MN2
45			*
46	pam-ba	3.4	* Define the S1A volume
47			*
48			CALL GSVOLU('S1A ','BOX ',NMED,S1A, 3,IVOLU)
49			*
50			* Define the S2A volume
51			*
52			CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)
53			*
54			* Define the S3A volume
55	cafagna	3.1	*
56	pam-ba	3.4	CALL GSVOLU('S3A ','BOX ',NMED,S3A, 3,IVOLU)
57	cafagna	3.1	*
58	pam-ba	3.4	NMED=MMYL
59	cafagna	3.1	*
60	pam-ba	3.4	* Define the S11M volume
61	cafagna	3.1	*
62	pam-ba	3.4	CALL GSVOLU('S11M','BOX ',NMED,S11M, 3,IVOLU)
63	cafagna	3.1	*
64	pam-ba	3.4	* Define the S12M volume
65	cafagna	3.1	*
66	pam-ba	3.4	CALL GSVOLU('S12M','BOX ',NMED,S12M, 3,IVOLU)
67	cafagna	3.1	*
68	pam-ba	3.4	* Define the S21M volume
69	cafagna	3.1	*
70	pam-ba	3.4	CALL GSVOLU('S21M','BOX ',NMED,S21M, 3,IVOLU)
71	cafagna	3.1	*
72	pam-ba	3.4	* Define the S22M volume
73	cafagna	3.1	*
74	pam-ba	3.4	CALL GSVOLU('S22M','BOX ',NMED,S22M, 3,IVOLU)
75			*
76			* Define the S31M volume
77			*
78			CALL GSVOLU('S31M','BOX ',NMED,S31M, 3,IVOLU)
79			*
80			* Define the S32M volume
81			*
82			CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)
83	cafagna	3.1	*
84			*
85			* Change medium just for S4.
86			*
87			NMED=MSCIN
88			*
89			* Define the S4 volume
90			*
91			CALL GSVOLU('S4 ','BOX ',NMED,S4, 3,IVOLU)
92			*
93			* Start subvolume definition.
94			*
95			*
96	cafagna	3.2	* Define the S11Y volume
97	cafagna	3.1	*
98	cafagna	3.2	CALL GSVOLU('S11Y','BOX ',NMED,S11Y, 3,IVOLU)
99	cafagna	3.1	*
100	cafagna	3.2	* Define the S12X volume
101	cafagna	3.1	*
102	cafagna	3.2	CALL GSVOLU('S12X','BOX ',NMED,S12X, 3,IVOLU)
103	cafagna	3.1	*
104			* Define the S21X volume
105			*
106			CALL GSVOLU('S21X','BOX ',NMED,S21X, 3,IVOLU)
107			*
108			* Define the S22Y volume
109			*
110			CALL GSVOLU('S22Y','BOX ',NMED,S22Y, 3,IVOLU)
111			*
112			* Define the S31Y volume
113			*
114			CALL GSVOLU('S31Y','BOX ',NMED,S31Y, 3,IVOLU)
115			*
116			* Define the S32X volume
117			*
118			CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)
119
120	pam-ba	3.4	c$$$<<<<<<< gptofv.F
121			c$$$=======
122			c$$$c ml: 21/12/04:
123			c$$$ NMED=MAL
124			c$$$*
125			c$$$* Define the SC11 volume
126			c$$$*
127			c$$$ CALL GSVOLU('SC11','BOX ',NMED,SC11, 3,IVOLU)
128			c$$$*
129			c$$$* Define the SC12 volume
130			c$$$*
131			c$$$ CALL GSVOLU('SC12','BOX ',NMED,SC12, 3,IVOLU)
132			c$$$*
133			c$$$* Define the SC21 volume
134			c$$$*
135			c$$$ CALL GSVOLU('SC21','BOX ',NMED,SC21, 3,IVOLU)
136			c$$$*
137			c$$$* Define the SC22 volume
138			c$$$*
139			c$$$ CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
140			c$$$
141			c$$$>>>>>>> 3.3
142	cafagna	3.2	c*
143			c* Define the STOF volume
144			c*
145			c NMED=MAL
146			c CALL GSVOLU('STOF','BOX ',NMED,STOF, 3,IVOLU)
147			c*
148			c* Define the SGLU volume
149			c*
150			cC COLLA ARALDITICA (DA VERIFICARNE LA COMPOSIZIONE)
151			c NMED=MSIC
152			c CALL GSVOLU('SGLU','BOX ',NMED,SGLU, 3,IVOLU)
153			c*
154			c* Define the SPV1 volume
155			c*
156			c NMED=? (PVC espanso)
157			c CALL GSVOLU('SPV1','BOX ',NMED,SPV1, 3,IVOLU)
158			c*
159			c* Define the SMYL volume
160			c*
161			c NMED=? (MYLAR)
162			c CALL GSVOLU('SMYL','BOX ',NMED,SMYL, 3,IVOLU)
163			c*
164			c* Define the SPV2 volume
165			c*
166			c NMED=? (PVC espanso)
167			c CALL GSVOLU('SPV2','BOX ',NMED,SPV2, 3,IVOLU)
168			cc end ml.
169	cafagna	3.1	*
170	pam-ba	3.4	* Positioning volumes
171	cafagna	3.1	*
172	pam-ba	3.4	X=0.
173			Y=0.
174			Z=0.
175			CALL GSPOS('S1A ',N,'S1 ',X,Y,Z,0,'ONLY')
176	cafagna	3.1	*
177	pam-ba	3.4	X=0.
178			Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
179			DO I=1,NS12X
180	cafagna	3.2	N= I
181	pam-ba	3.4	Y= -S1A(2) + HMYLTOF + (2I-1)S12M(2)
182			CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
183	cafagna	3.1	ENDDO
184			Y= 0.
185	pam-ba	3.4	Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
186			DO I=1,NS11Y
187	cafagna	3.1	N= I
188	pam-ba	3.4	X= -S1A(1) + HMYLTOF + (2I-1)S11M(1)
189			CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
190	cafagna	3.1	ENDDO
191	pam-ba	3.4	C Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
192			C CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
193			C Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
194			C CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
195			*
196			X=0.
197			Y=0.
198			Z=0.
199			CALL GSPOS('S2A ',N,'S2 ',X,Y,Z,0,'ONLY')
200	cafagna	3.1	*
201			Y= 0.
202	pam-ba	3.4	Z=-S2A(3)+HGLUTOF+S22M(3)
203			DO I=1,NS22Y
204	cafagna	3.1	N= I
205	pam-ba	3.4	X= -S2A(1) + HMYLTOF + (2I-1)S22M(1)
206			CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
207	cafagna	3.1	ENDDO
208			*
209	pam-ba	3.4	X= 0.
210			Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
211	cafagna	3.1	DO I=1,NS21X
212			N= I
213	pam-ba	3.4	Y= -S2A(2) + HMYLTOF+ (2I-1)S21M(2)
214			CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
215	cafagna	3.1	ENDDO
216	pam-ba	3.4	C Z=-S2A(3)+HGLUTOF+S22M(3)
217			C CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
218			C Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
219			C CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
220			*
221	cafagna	3.1	X= 0.
222	pam-ba	3.4	Z=-S3A(3)+HAIRS3+S32M(3)
223			DO I=1,NS32X
224	cafagna	3.1	N= I
225	pam-ba	3.4	Y= -S3A(2)+ (2I-1)S32M(2)
226			CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')
227	cafagna	3.1	ENDDO
228			*
229	pam-ba	3.4	Y= 0.
230			Z= -S3A(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
231			DO I=1,NS31Y
232	cafagna	3.1	N= I
233	pam-ba	3.4	X= -S3A(1)+(2I-1)S31M(1)
234			CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')
235	cafagna	3.1	ENDDO
236	pam-ba	3.4	*
237			C Z=-S3A(3)+HAIRS3+S32M(3)
238			C CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')
239			C Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
240			C CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')
241			****************************************
242			c ml: 21/12/04:
243			X=0.
244			Y= 0.
245			Z= 0.
246			N=1
247			c ml: 27/10/04:
248			c DO I=1,NS11Y
249			c N= I
250			c X= -S11M(1) + HMYLTOF + (2I-1)S11Y(1)
251			CALL GSPOS('S11Y',N,'S11M',X,Y,Z,0,'ONLY')
252			c ENDDO
253			c end ml.
254			*
255	cafagna	3.1	X= 0.
256	pam-ba	3.4	Y=0.
257			Z= 0.
258			N=1
259			c ml: 27/10/04
260			C DO I=1,NS12X
261			C N= I
262			C Y= -S11M(2) + HMYLTOF + (2I-1)S12X(2)
263			CALL GSPOS('S12X',N,'S12M',X,Y,Z,0,'ONLY')
264			C ENDDO
265			c end ml.
266			*
267			X=0.
268			y=0.
269			Z= 0.
270			N=1
271			C DO I=1,NS21X
272			C N= I
273			c Y= -S21M(2) + HMYLTOF+ (2I-1)S21X(2)
274			CALL GSPOS('S21X',N,'S21M',X,Y,Z,0,'ONLY')
275			C ENDDO
276			*
277			X=0.
278	cafagna	3.1	Y= 0.
279	pam-ba	3.4	Z= 0.
280			N=1
281			C DO I=1,NS22Y
282			C N= I
283			C X= -S22M(1) + HMYLTOF + (2I-1)S22Y(1)
284			CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
285			C ENDDO
286	cafagna	3.2	c*
287			c* Positioning volumes into S11:
288			c+
289			c X=0.
290			c Y=0.
291			c Z=-S11(3)+SMYL(3)
292			c N=1
293			c CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
294			c Z=Z+SMYL(3)+S11Y(3)
295			c DO I=1,NS11Y
296			c N= I
297			c X= -S11(1) + (2I-1)S11Y(1)
298			c CALL GSPOS('S11Y',N,'S11 ',X,Y,Z,0,'ONLY')
299			c ENDDO
300			c N=2
301			c Z=Z+S11Y(3)+SMYL(3)
302			c CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
303			c Z=Z+SMYL(3)+SPV1(3)
304			c N=1
305			c CALL GSPOS('SPV1',N,'SPV1',X,Y,Z,0,'ONLY')
306			c Z=Z+SPV1(3)+SGLU(3)
307			c N=1
308			c CALL GSPOS('SGLU',N,'SGLU',X,,Y,Z,0,'ONLY')
309			c Z=Z+SGLU(3)+STOF(3)
310			c N=1
311			c CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
312			c*
313			c* Positioning volumes into S12:
314			c*
315			c Z=-S12(3)+STOF(3)
316			c N=2
317			c CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
318			c Z=Z+STOF(3)+SGLU(3)
319			c CALL GSPOS('SGLU',N,'SGLU',X,Y,Z,0,'ONLY')
320			c Z=Z+SGLU(3)+SPV2
321	cafagna	3.1	Y= 0.
322	pam-ba	3.4	Z= -S31M(3)+HMYLTOF+S31Y(3)
323			N=1
324			C DO I=1,NS31Y
325			C N= I
326			C X= -S31M(1) + HMYLTOF+(2I-1)S31Y(1)
327			X= -S31M(1) + 2*(HMYLTOF/2)+S31Y(1)
328			CALL GSPOS('S31Y',N,'S31M',X,Y,Z,0,'ONLY')
329			C ENDDO
330	cafagna	3.1	*
331			X= 0.
332	pam-ba	3.4	Z= -S32M(3)+HMYLTOF+S32X(3)
333			N=1
334			C DO I=1,NS32X
335			C N= I
336			C Y= -S32M(2) + HMYLTOF+ (2I-1)S32X(2)
337			Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
338			CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
339			C ENDDO
340	cafagna	3.1	*
341			RETURN
342			END