1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
*DECK MCGDTV
SUBROUTINE MCGDTV(NDIM,NFI,NREG,NSOU,NSEG,NMU,LMCU,LMXMCU,NZONA,
1 NRSEG,MCUW,MCUI,WEI2D,SEGLEN,WZMU,SURFD,CYCLIC,
2 ACFLAG,ZMU,XSIXYZ,CAZ)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Compute the contribution of a track to the numerical surfaces and
* connection matrices for an EXCELT tracking.
*
*Copyright:
* Copyright (C) 2002 Ecole Polytechnique de Montreal
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version
*
*Author(s): R. Le Tellier
*
*Parameters: input
* NDIM number of dimensions for the geometry.
* NFI total number of volumes and surfaces.
* NREG number of regions.
* NSOU number of external surfaces.
* NSEG number of segments for this track.
* NMU number of polar angles.
* LMXMCU maximum dimension for the connection matrix.
* NZONA index-number of the mixture/albedo type assigned to
* each volume/surface.
* NRSEG vector containing the region number of the different segments
* of this track.
* WEI2D weight for this track.
* SEGLEN vector containing the length of the different segments of this
* track.
* ZMU polar quadrature points.
* WZMU polar quadrature weights.
* CYCLIC cyclic tracking flag.
* ACFLAG preconditioning techniques flag.
* CAZ directional cosines.
*
*Parameters: input/output
* LMCU number of elements in the connection matrix.
* MCUW temporary connection matrix.
* MCUI temporary connection matrix.
* SURFD numerical surfaces.
* XSIXYZ XSI for B1 leakage.
*
*-----------------------------------------------------------------------
*
IMPLICIT NONE
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER NDIM,NFI,NREG,NSOU,NSEG,NMU,LMCU,LMXMCU,
1 NZONA(NFI),NRSEG(NSEG),MCUW(LMXMCU),MCUI(LMXMCU)
REAL WZMU(NMU)
DOUBLE PRECISION WEI2D,SEGLEN(NSEG),SURFD(NSOU)
LOGICAL CYCLIC,ACFLAG
REAL ZMU(NMU)
DOUBLE PRECISION CAZ(NDIM),XSIXYZ(NSOU,3)
DOUBLE PRECISION OMEGA2(3)
INTEGER IDIR
*---
* LOCAL VARIABLES
*---
INTEGER II,NOMCEL,IMU,ITEMP
DOUBLE PRECISION WEIGHT
IF(.NOT.CYCLIC) THEN
* non cyclic tracking: calculate numerical surfaces
DO II=1,NSEG,NSEG-1
NOMCEL=-NRSEG(II)
IF (NOMCEL.GT.0) THEN
IF (NDIM.EQ.2) THEN
DO IMU=1,NMU
OMEGA2(1)=(CAZ(1)/DBLE(ZMU(IMU)))**2
OMEGA2(2)=(CAZ(2)/DBLE(ZMU(IMU)))**2
OMEGA2(3)=1.0D0-1.0D0/DBLE(ZMU(IMU)**2)
WEIGHT=WEI2D*DBLE(WZMU(IMU))
SURFD(NOMCEL)=SURFD(NOMCEL)+WEIGHT
DO IDIR=1,3
XSIXYZ(NOMCEL,IDIR)=XSIXYZ(NOMCEL,IDIR)+
1 3.0D0*OMEGA2(IDIR)*WEIGHT
ENDDO
ENDDO
ELSE
WEIGHT=WEI2D
SURFD(NOMCEL)=SURFD(NOMCEL)+WEIGHT
DO IDIR=1,3
XSIXYZ(NOMCEL,IDIR)=XSIXYZ(NOMCEL,IDIR)+
1 3.0D0*CAZ(IDIR)*CAZ(IDIR)*WEIGHT
ENDDO
ENDIF
ENDIF
ENDDO
ENDIF
*
IF(ACFLAG) THEN
* SCR or ACA acceleration required
DO II=1,NSEG
IF(NRSEG(II).LT.0) THEN
NRSEG(II)=NREG-NRSEG(II)
ELSE IF(NRSEG(II).EQ.0) THEN
NRSEG(II)=NREG+1
ENDIF
ENDDO
IF (CYCLIC) THEN
* cyclic tracking: "unfold" the tracking line
* calculate connection matrices
CALL MCGTRK(NFI,NZONA,NSEG,NRSEG,SEGLEN)
CALL MOCCAL(NSEG,NRSEG,NREG,MCUW,MCUI,LMCU,LMXMCU)
DO II=1,NSEG/2
ITEMP=NRSEG(II)
NRSEG(II)=NRSEG(NSEG+1-II)
NRSEG(NSEG+1-II)=ITEMP
ENDDO
CALL MOCCAL(NSEG,NRSEG,NREG,MCUW,MCUI,LMCU,LMXMCU)
ELSE
* non-cyclic tracking: calculate connection matrices
CALL MCGCAL(NSEG,NRSEG,NREG,MCUW,MCUI,LMCU,LMXMCU)
DO II=1,NSEG/2
ITEMP=NRSEG(II)
NRSEG(II)=NRSEG(NSEG+1-II)
NRSEG(NSEG+1-II)=ITEMP
ENDDO
CALL MCGCAL(NSEG,NRSEG,NREG,MCUW,MCUI,LMCU,LMXMCU)
ENDIF
ENDIF
*
RETURN
END
|
