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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
|
*DECK DEVDRV
SUBROUTINE DEVDRV(IPDEV,IPMTX,IGEO,NMIX,NTOT,LIMIT)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Read specifications for the rod-devices from the input file.
*
*Copyright:
* Copyright (C) 2007 Ecole Polytechnique de Montreal.
*
*Author(s):
* D. Sekki and A. Hebert
*
*Parameters: input
* IPDEV pointer to device information.
* IPMTX pointer to matex information.
* IGEO index related to the reactor geometry.
* NMIX old maximum number of material mixtures.
* NTOT old total number of all mixtures.
* LIMIT core limiting coordinates.
*
*-----------------------------------------------------------------------
*
USE GANLIB
*----
* SUBROUTINE ARGUMENTS
*----
TYPE(C_PTR) IPDEV,IPMTX
INTEGER IGEO,NMIX,NTOT
REAL LIMIT(6)
*----
* LOCAL VARIABLES
*----
PARAMETER(NSTATE=40,IOUT=6,MAXPRT=10)
CHARACTER TEXT*12,HSMG*131
TYPE(C_PTR) JPDEV,KPDEV
INTEGER ISTATE(NSTATE),NRGRP,DMIX(2,MAXPRT)
DOUBLE PRECISION DFLOT
INTEGER, ALLOCATABLE, DIMENSION(:) :: MIX
*----
* CORE LIMITS
*----
CALL LCMPUT(IPDEV,'CORE-LIMITS',6,2,LIMIT)
*----
* READ INPUT DATA
*----
IMPX=1
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(ITYP.NE.3)CALL XABORT('@DEVDRV: CHARACTER DATA EXPECTED(1).')
IF(TEXT.NE.'EDIT')GOTO 10
* PRINTING INDEX
CALL REDGET(ITYP,IMPX,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@DEVDRV: INTEGER FOR EDIT EXPECTED.')
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(ITYP.NE.3)CALL XABORT('@DEVDRV: CHARACTER DATA EXPECTED(2).')
10 IF(TEXT.NE.'NUM-ROD')CALL XABORT('@DEVDRV: KEYWORD NUM-ROD EX'
1 //'PECTED.')
* TOTAL NUMBER OF RODS
CALL REDGET(ITYP,NROD,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@DEVDRV: INTEGER TOTAL NUMBER OF ROD'
1 //'S EXPECTED.')
IF(NROD.LT.1)CALL XABORT('@DEVDRV: WRONG TOTAL NUMBER OF RODS <1')
IF(IMPX.GT.1)WRITE(IOUT,1003) LIMIT(1),LIMIT(3),LIMIT(5),LIMIT(2),
1 LIMIT(4),LIMIT(6)
IF(IMPX.GT.0)WRITE(IOUT,1000) NROD
*
MAXTOT=NTOT+NROD*2*MAXPRT
ALLOCATE(MIX(MAXTOT))
MIX(:MAXTOT)=0
CALL LCMGET(IPMTX,'MAT',MIX)
*----
* READ OPTION
*----
NRGRP=0
IMODE=1
JPDEV=LCMLID(IPDEV,'DEV_ROD',NROD)
30 CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(TEXT.EQ.'ROD')THEN
* READ INDIVIDUAL ROD DATA
CALL DEVGET(JPDEV,NROD,LIMIT,IMODE,IMPX)
ELSE IF(TEXT.EQ.'CREATE')THEN
* CREATE ROD-GROUPS
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(TEXT.NE.'ROD-GR') CALL XABORT('@DEVDRV: KEYWORD ROD-GR EX'
1 //'PECTED.')
CALL REDGET(ITYP,NRGRP,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1) CALL XABORT('@DEVDRV: INTEGER NUMBER OF ROD-GR'
1 //'OUPS EXPECTED.')
IF(NRGRP.LT.1) CALL XABORT('@DEVDRV: WRONG NUMBER OF GROUPS <1')
CALL DEVDGD(IPDEV,NROD,NRGRP,IMPX)
GO TO 40
ELSE IF(TEXT.EQ.'FADE')THEN
IMODE=1
ELSE IF(TEXT.EQ.'MOVE')THEN
IMODE=2
ELSE IF(TEXT.EQ.';') THEN
GOTO 40
ELSE
WRITE(HSMG,'(26H@DEVDRV: INVALID KEYWORD (,A,2H).)') TEXT
CALL XABORT(HSMG)
ENDIF
GOTO 30
*----
* VALIDATE ROD DATA AND SET MIXTURE INDICES
*----
40 IOFSET=0
DO 60 ID=1,NROD
CALL LCMLEL(JPDEV,ID,LENGT,ITYLCM)
IF(LENGT.EQ.0) THEN
WRITE(HSMG,'(18H@DEVDRV: ROD INDEX,I5,16H IS NOT DEFINED.)') ID
CALL XABORT(HSMG)
ENDIF
KPDEV=LCMGIL(JPDEV,ID)
CALL LCMGET(KPDEV,'ROD-PARTS',NPART)
IF(NPART.GT.MAXPRT) CALL XABORT('@DEVDRV: MAXPRT OVERFLOW.')
CALL LCMGET(KPDEV,'ROD-MIX',DMIX)
DO 55 IPART=1,NPART
DO 50 I=1,2
IOFSET=IOFSET+1
IF(IOFSET.GT.MAXTOT) CALL XABORT('@DEVDRV: MAXTOT OVERFLOW.')
MIX(NTOT+IOFSET)=DMIX(I,IPART)
DMIX(I,IPART)=NMIX+IOFSET
50 CONTINUE
55 CONTINUE
CALL LCMPUT(KPDEV,'ROD-MIX',2*NPART,1,DMIX)
60 CONTINUE
*----
* STATE-VECTORS
*----
ISTATE(:NSTATE)=0
ISTATE(1)=IGEO
ISTATE(2)=NROD
ISTATE(3)=NRGRP
ISTATE(6)=IMODE
CALL LCMPUT(IPDEV,'STATE-VECTOR',NSTATE,1,ISTATE)
IF(IMPX.GT.1)CALL LCMLIB(IPDEV)
* UPDATE MATEX
ISTATE(:NSTATE)=0
CALL LCMGET(IPMTX,'STATE-VECTOR',ISTATE)
ISTATE(2)=NMIX+IOFSET
ISTATE(5)=NTOT+IOFSET
CALL LCMPUT(IPMTX,'MAT',NTOT+IOFSET,1,MIX)
CALL LCMPUT(IPMTX,'STATE-VECTOR',NSTATE,1,ISTATE)
DEALLOCATE(MIX)
IF(IMPX.GT.4) CALL LCMLIB(IPMTX)
RETURN
*
1000 FORMAT(/1X,'DEVDRV: GIVEN TOTAL NUMBER OF ROD-DEVICES:',
1 I5//' ** READING INPUT DATA FOR RODS **')
1003 FORMAT(//5X,'--- REACTOR CORE LIMITS ---'//
1 1X,'Xmin',F10.4,5X,'Ymin',F10.4,5X,'Zmin',F10.4/
2 1X,'Xmax',F10.4,5X,'Ymax',F10.4,5X,'Zmax',F10.4/)
END
|
