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
|
*DECK RESIND
SUBROUTINE RESIND(IPMAP,IPMTX,NX,NY,NZ,LX,LY,LZ,MIX,NFUEL,IMPX)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Update material index, it will store the negative fuel mixtures.
*
*Copyright:
* Copyright (C) 2007 Ecole Polytechnique de Montreal.
*
*Author(s):
* E. Varin, D. Sekki
*
*Parameters: input/output
* IPMAP pointer to fuel-map information.
* IPMTX pointer to matex information.
* NX number of elements along x-axis in fuel map.
* NY number of elements along y-axis in fuel map.
* NZ number of elements along z-axis in fuel map.
* LX number of elements along x-axis in geometry.
* LY number of elements along y-axis in geometry.
* LZ number of elements along z-axis in geometry.
* MIX renumbered index over the fuel-map geometry.
* NFUEL number of fuel types.
* IMPX printing index (=0 for no print).
*
*-----------------------------------------------------------------------
*
USE GANLIB
*----
* SUBROUTINE ARGUMENTS
*----
TYPE(C_PTR) IPMAP,IPMTX
INTEGER NX,NY,NZ,LX,LY,LZ,MIX(NX*NY*NZ),NFUEL,IMPX
*----
* LOCAL VARIABLES
*----
PARAMETER(IOUT=6)
INTEGER INDX(LX*LY*LZ),FTOT(NFUEL)
REAL MAPXX(NX+1),MAPYY(NY+1),MAPZZ(NZ+1),
1 GEOXX(LX+1),GEOYY(LY+1),GEOZZ(LZ+1)
TYPE(C_PTR) JPMAP
*----
* ALLOCATABLE ARRAYS
*----
INTEGER, ALLOCATABLE, DIMENSION(:) :: IMAT
*----
* UPDATE MATERIAL INDEX
*----
CALL LCMGET(IPMTX,'MESHX',GEOXX)
CALL LCMGET(IPMTX,'MESHY',GEOYY)
CALL LCMGET(IPMTX,'MESHZ',GEOZZ)
JPMAP=LCMGID(IPMAP,'GEOMAP')
CALL LCMGET(JPMAP,'MESHX',MAPXX)
CALL LCMGET(JPMAP,'MESHY',MAPYY)
CALL LCMGET(JPMAP,'MESHZ',MAPZZ)
INDX(:LX*LY*LZ)=0
IF(IMPX.GT.2)WRITE(IOUT,*)'UPDATING MATERIAL INDEX'
I1=0
I2=0
J1=0
J2=0
K1=0
K2=0
DO 52 KM=1,NZ
DO 51 JM=1,NY
DO 50 IM=1,NX
DO IG=1,LX
IF(MAPXX(IM).EQ.GEOXX(IG)) I1=IG
IF(MAPXX(IM+1).EQ.GEOXX(IG+1))THEN
I2=IG
GOTO 10
ENDIF
ENDDO
10 DO JG=1,LY
IF(MAPYY(JM).EQ.GEOYY(JG)) J1=JG
IF(MAPYY(JM+1).EQ.GEOYY(JG+1))THEN
J2=JG
GOTO 20
ENDIF
ENDDO
20 DO KG=1,LZ
IF(MAPZZ(KM).EQ.GEOZZ(KG)) K1=KG
IF(MAPZZ(KM+1).EQ.GEOZZ(KG+1))THEN
K2=KG
GOTO 30
ENDIF
ENDDO
30 IELM=(KM-1)*NX*NY+(JM-1)*NX +IM
DO 42 KG=K1,K2
DO 41 JG=J1,J2
DO 40 IG=I1,I2
IELG=(KG-1)*LX*LY+(JG-1)*LX+IG
INDX(IELG)=MIX(IELM)
40 CONTINUE
41 CONTINUE
42 CONTINUE
50 CONTINUE
51 CONTINUE
52 CONTINUE
* CHECK TOTAL NUMBER
ITOT=0
DO 60 IEL=1,LX*LY*LZ
IF(INDX(IEL).NE.0)ITOT=ITOT+1
60 CONTINUE
NTOT=0
CALL LCMGET(IPMTX,'FTOT',FTOT)
DO 70 IFUEL=1,NFUEL
NTOT=NTOT+FTOT(IFUEL)
70 CONTINUE
IF(ITOT.NE.NTOT) THEN
WRITE(IOUT,'(/15H @RESIND: ITOT=,I8,6H NTOT=,I8)') ITOT,NTOT
CALL XABORT('@RESIND: FOUND DIFFERENT TOTAL NUMBER OF FUEL MI'
1 //'XTURES IN FUEL-MAP AND MATEX.')
ENDIF
* STORE NEGATIVE FUEL MIXTURES
CALL LCMLEN(IPMTX,'MAT',LENGT,ITYP)
ALLOCATE(IMAT(LENGT))
IMAT(:LENGT)=0
CALL LCMGET(IPMTX,'MAT',IMAT)
DO 100 IEL=1,LX*LY*LZ
IF(INDX(IEL).NE.0)IMAT(IEL)=-INDX(IEL)
100 CONTINUE
CALL LCMPUT(IPMTX,'MAT',LENGT,1,IMAT)
DEALLOCATE(IMAT)
RETURN
END
|
