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
|
*DECK CVR
SUBROUTINE CVR(NENTRY,HENTRY,IENTRY,JENTRY,KENTRY)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Perform reordering of fuel-regions properties in the reactor core,
* according to the specified voiding pattern.
*
*Copyright:
* Copyright (C) 2007 Ecole Polytechnique de Montreal.
*
*Author(s):
* D. Sekki
*
*Parameters: input
* NENTRY number of data structures transfered to this module.
* HENTRY name of the data structures.
* IENTRY data structure type where:
* IENTRY=1 for LCM memory object;
* IENTRY=2 for XSM file;
* IENTRY=3 for sequential binary file;
* IENTRY=4 for sequential ASCII file.
* JENTRY access permission for the data structure where:
* JENTRY=0 for a data structure in creation mode;
* JENTRY=1 for a data structure in modifications mode;
* JENTRY=2 for a data structure in read-only mode.
* KENTRY data structure pointer.
*
*Comments:
* The CVR: module specification is:
* FMAPV := CVR: FMAP :: (descrcvr) ;
* where
* FMAP : name of a read-only \emph{fmap} object,
* created in the RESINI: module. This object must contain the non-perturbed
* fuel-cell properties.
* FMAPV : name of a new \emph{fmap} object,
* that will contain the modified fuel-type indices and reordered coolant
* densities according to the specified core-voiding pattern.
* (descrcvr) : structure describing the input data to the CVR: module.
*
*-----------------------------------------------------------------------
*
USE GANLIB
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER NENTRY,IENTRY(NENTRY),JENTRY(NENTRY)
TYPE(C_PTR) KENTRY(NENTRY)
CHARACTER HENTRY(NENTRY)*12
*----
* LOCAL VARIABLES
*----
PARAMETER(NSTATE=40)
CHARACTER TEXT*12,HSIGN*12
INTEGER ISTATE(NSTATE),IGST(NSTATE)
DOUBLE PRECISION DFLOT
TYPE(C_PTR) IPMAP,JPMAP
*----
* PARAMETER VALIDATION
*----
IF(NENTRY.NE.2)CALL XABORT('@CVR: TWO PARAMETERS EXPECTED.')
TEXT=HENTRY(1)
IF((IENTRY(1).NE.1).AND.(IENTRY(1).NE.2))CALL XABORT('@CVR:'
1 //' LCM OBJECT EXPECTED AT LHS ('//TEXT//').')
IF(JENTRY(1).NE.0)CALL XABORT('@CVR: FUEL MAP OBJECT IN CRE'
1 //'ATE MODE EXPECTED AT LHS ('//TEXT//').')
TEXT=HENTRY(2)
IF((IENTRY(2).NE.1).AND.(IENTRY(2).NE.2))CALL XABORT('@CVR:'
1 //' LCM OBJECT EXPECTED AT RHS ('//TEXT//').')
IF(JENTRY(2).NE.2)CALL XABORT('@CVR: FUEL MAP OBJECT IN REA'
1 //'D-ONLY MODE EXPECTED AT RHS ('//TEXT//').')
CALL LCMGTC(KENTRY(2),'SIGNATURE',12,HSIGN)
IF(HSIGN.NE.'L_MAP')THEN
TEXT=HENTRY(2)
CALL XABORT('@CVR: SIGNATURE OF '//TEXT//' IS '//HSIGN//
1 '. L_MAP EXPECTED.')
ENDIF
IPMAP=KENTRY(1)
CALL LCMEQU(KENTRY(2),IPMAP)
*----
* RECOVER INFORMATION
*----
ISTATE(:NSTATE)=0
CALL LCMGET(IPMAP,'STATE-VECTOR',ISTATE)
NB=ISTATE(1)
NCH=ISTATE(2)
NFUEL=ISTATE(7)
NPARM=ISTATE(8)
* FUEL-MAP GEOMETRY
JPMAP=LCMGID(IPMAP,'GEOMAP')
IGST(:NSTATE)=0
CALL LCMGET(JPMAP,'STATE-VECTOR',IGST)
IF(IGST(1).NE.7)CALL XABORT('@CVR: ONLY 3-D CART'
1 //'ESIAN GEOMETRY ALLOWED.')
NX=IGST(3)
NY=IGST(4)
NZ=IGST(5)
* PRINTING INDEX
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(ITYP.NE.3)CALL XABORT('@CVR: CHARACTER DATA EXPECTED.')
IF(TEXT.NE.'EDIT')CALL XABORT('@CVR: KEYWORD EDIT EXPECTED.')
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@CVR: INTEGER FOR EDIT EXPECTED.')
IMPX=MAX(0,NITMA)
* READ INPUT DATA
CALL CVRDRV(IPMAP,NCH,NB,NFUEL,NPARM,NX,NY,NZ,NVD,IVD,IMPX)
* UPDATE STATE-VECTOR
ISTATE(10)=NVD
ISTATE(11)=IVD
CALL LCMPUT(IPMAP,'STATE-VECTOR',NSTATE,1,ISTATE)
IF(IMPX.GT.1) CALL LCMLIB(IPMAP)
RETURN
END
|
