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
|
*DECK MOVDEV
SUBROUTINE MOVDEV(NENTRY,HENTRY,IENTRY,JENTRY,KENTRY)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Simulate the time-dependent displacement of individual devices
* and/or of groups of devices in the reactor core.
*
*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 MOVDEV: module specification is:
* DEVICE := MOVDEV: DEVICE :: (descmove) ;
* where
* DEVICE : name of the \emph{device} object that will be modified by the
* module. The rods positions are updated according to the current time step
* of movement.
* (descmove) : structure describing the input data to the MOVDEV: 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,IOUT=6)
CHARACTER TEXT*12,HSIGN*12
INTEGER ISTATE(NSTATE),DGRP
DOUBLE PRECISION DFLOT
*----
* PARAMETER VALIDATION
*----
IF(NENTRY.GT.1)CALL XABORT('@MOVDEV: ONE PARAMETER EXPECTED.')
IF((IENTRY(1).NE.1).AND.(IENTRY(1).NE.2))CALL XABORT('@MOV'
1 //'DEV: LCM OBJECT EXPECTED AT LHS.')
CALL LCMGTC(KENTRY(1),'SIGNATURE',12,HSIGN)
IF(HSIGN.NE.'L_DEVICE')THEN
TEXT=HENTRY(1)
CALL XABORT('@MOVDEV: SIGNATURE OF '//TEXT//' IS '//HSIGN//
1 '. L_DEVICE EXPECTED.')
ENDIF
IF(JENTRY(1).NE.1)CALL XABORT('@MOVDEV: MODIFICATION MODE EX'
1 //'PECTED FOR L_DEVICE.')
*
CALL LCMGET(KENTRY(1),'STATE-VECTOR',ISTATE)
IGEO=ISTATE(1)
IF(IGEO.NE.7)CALL XABORT('@MOVDEV: ONLY 3D-CARTESIAN GEOMETR'
1 //'Y ALLOWED.')
NDEV=ISTATE(2)
DGRP=ISTATE(3)
IMODE=ISTATE(6)
IF(IMODE.EQ.0)CALL XABORT('@MOVDEV: IMODE NOT SET.')
*----
* RECOVER INFORMATION
*----
IMPX=1
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(ITYP.NE.3)CALL XABORT('@MOVDEV: CHARACTER DATA EXPECTED.')
IF(TEXT.NE.'EDIT')GOTO 10
* PRINTING INDEX
CALL REDGET(ITYP,IMPX,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@MOVDEV: INTEGER FOR EDIT EXPECTED.')
* TIME STEP INCREMENT
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
10 IF(TEXT.NE.'DELT')CALL XABORT('@MOVDEV: KEYWORD DELT EXPECTED.')
CALL REDGET(ITYP,NITMA,DELT,TEXT,DFLOT)
IF(ITYP.NE.2)CALL XABORT('@MOVDEV: REAL FOR DELT EXPECTED.')
IF(DELT.LE.0.)CALL XABORT('@MOVDEV: VALUE OF DELT SHOULD B'
1 //'E POSITIVE.')
ND=0
NG=0
20 ND=ND+1
CALL REDGET(ITYP,NITMA,FLOT,TEXT,DFLOT)
IF(TEXT.EQ.'ROD')THEN
*----
* ROD OPTION
*----
CALL REDGET(ITYP,ID,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@MOVDEV: INTEGER ROD-ID NUMB'
1 //'ER EXPECTED.')
IF((ID.GT.NDEV).OR.(ID.EQ.0))THEN
WRITE(IOUT,*)'@MOVDEV: READ CURRENT ROD-ID #',ID
CALL XABORT('@MOVDEV: WRONG ROD-ID NUMBER.')
ENDIF
IF(IMPX.GT.0)WRITE(IOUT,1000)ID
CALL MOVPOS(KENTRY(1),IMODE,ID,DELT,IMPX)
ELSEIF(TEXT.EQ.'GROUP')THEN
*----
* GROUP OPTION
*----
CALL REDGET(ITYP,IGRP,FLOT,TEXT,DFLOT)
IF(ITYP.NE.1)CALL XABORT('@MOVDEV: INTEGER GROUP-ID NUM'
1 //'BER EXPECTED.')
IF((IGRP.GT.DGRP).OR.(IGRP.EQ.0))THEN
WRITE(IOUT,*)'@MOVDEV: READ CURRENT GROUP-ID #',IGRP
CALL XABORT('@MOVDEV: WRONG GROUP-ID NUMBER.')
ENDIF
IF(IMPX.GT.0)WRITE(IOUT,1001)IGRP
CALL MOVGRP(KENTRY(1),IMODE,IGRP,NDGR,DELT,IMPX)
ND=ND+NDGR-1
NG=NG+1
*
ELSEIF(TEXT.EQ.';')THEN
GOTO 30
ELSE
WRITE(IOUT,*)'@MOVDEV: WRONG KEYWORD : ',TEXT
CALL XABORT('@MOVDEV: KEYWORD ROD OR GROUP EXPECTED.')
ENDIF
GOTO 20
30 IF(IMPX.GT.0)WRITE(IOUT,1002)NG,ND-1
IF(IMPX.GT.4)CALL LCMLIB(KENTRY(1))
RETURN
*
1000 FORMAT(/5X,'MOVING ROD #',I3.3)
1001 FORMAT(/5X,'MOVING GROUP #',I2.2)
1002 FORMAT(
1 /5X,'-------------------------------------'/
2 5X,'TOTAL NUMBER OF DISPLACED GROUPS : ',I2/
3 5X,'TOTAL NUMBER OF DISPLACED RODS : ',I3/)
END
|
