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
|
*DECK USPMIX
SUBROUTINE USPMIX(IPMTX,NEL,NREFL,NFUEL,MAT,RMIX,FMIX,INDX,NMIX)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Recover and check the material mixtures.
*
*Copyright:
* Copyright (C) 2007 Ecole Polytechnique de Montreal.
*
*Author(s):
* D. Sekki
*
*Parameters: input
* IPMTX pointer to matex information.
* NEL total number of volumes in reactor geometry.
* NREFL total number of reflector types.
* NFUEL total number of fuel types.
* MAT material index from geometry.
* RMIX reflector-type mixtures indices.
* FMIX fuel-type mixtures indices.
*
*Parameters: output
* INDX renumbered material index.
* NMIX total number of non-virtual volumes.
*
*-----------------------------------------------------------------------
*
USE GANLIB
*----
* SUBROUTINE ARGUMENTS
*----
TYPE(C_PTR) IPMTX
INTEGER NEL,NREFL,NFUEL,MAT(NEL),RMIX(NREFL),FMIX(NFUEL),
1 INDX(NEL),NMIX
*----
* LOCAL VARIABLES
*----
INTEGER, ALLOCATABLE, DIMENSION(:) :: RTOT,FTOT
*----
* SCRATCH STORAGE ALLOCATION
*----
ALLOCATE(RTOT(NREFL),FTOT(NFUEL))
*----
* MATERIAL INDEX
*----
RTOT(:NREFL)=0
FTOT(:NFUEL)=0
INDX(:NEL)=0
NMIX=0
DO IEL=1,NEL
IF(MAT(IEL).NE.0)THEN
NMIX=NMIX+1
INDX(IEL)=NMIX
ENDIF
ENDDO
IF((NFUEL.EQ.0).AND.(NREFL.EQ.0)) GOTO 20
* CHECK MIXTURES
DO 10 IEL=1,NEL
IMIX=MAT(IEL)
IF(IMIX.EQ.0)GOTO 10
IF(NREFL.EQ.0)GOTO 5
DO IREFL=1,NREFL
IF(IMIX.EQ.RMIX(IREFL))THEN
RTOT(IREFL)=RTOT(IREFL)+1
GOTO 10
ENDIF
ENDDO
5 IF(NFUEL.EQ.0)GOTO 10
DO IFUEL=1,NFUEL
IF(IMIX.EQ.FMIX(IFUEL))THEN
FTOT(IFUEL)=FTOT(IFUEL)+1
GOTO 10
ENDIF
ENDDO
10 CONTINUE
* STORAGE
20 CALL LCMPUT(IPMTX,'MAT',NEL,1,MAT)
CALL LCMPUT(IPMTX,'INDEX',NEL,1,INDX)
IF(NREFL.NE.0) THEN
CALL LCMPUT(IPMTX,'RMIX',NREFL,1,RMIX)
CALL LCMPUT(IPMTX,'RTOT',NREFL,1,RTOT)
ENDIF
IF(NFUEL.NE.0) THEN
CALL LCMPUT(IPMTX,'FMIX',NFUEL,1,FMIX)
CALL LCMPUT(IPMTX,'FTOT',NFUEL,1,FTOT)
ENDIF
*----
* SCRATCH STORAGE DEALLOCATION
*----
DEALLOCATE(FTOT,RTOT)
RETURN
END
|
