forked from etsap-TIMES/TIMES_model
-
Notifications
You must be signed in to change notification settings - Fork 0
/
dumpsol1.mod
204 lines (177 loc) · 5.8 KB
/
dumpsol1.mod
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
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
*++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
* Copyright (C) 2000-2020 Energy Technology Systems Analysis Programme (ETSAP)
* This file is part of the IEA-ETSAP TIMES model generator, licensed
* under the GNU General Public License v3.0 (see file LICENSE.txt).
*=============================================================================*
* DUMPSOLV.MOD outputs the actual solution values
* Dumps output SET, SCALARS, PARAMETER, TABLE
*
* 1 - type indicator = 'L'evel/'M'arginal
* 2-n - for rest of line = component names
*
* - Note that (declaration) is not output
* - $BATINCLUDE dumpsol.mod ITEM1 ITEM2 ... must fit on one line
*
*=============================================================================*
*GaG Questions/Comments:
*-----------------------------------------------------------------------------
file.tf = 0;
$offlisting
$set preditem ''
* hold the first parameter for later checking
$set whtype %1
*GG* later take 2nd as name of set controlling the TABLE column header to improve xeq
$shift
$label more
* only one call executes, so done when no more parameters
$if %1a == a $goto done
$if declared %1 $goto declared
* out error message if passed entity not declared
put / '*** UNKNOWN':20,'NAME':6,'%1':10 /;
$set preditem 'unknown'
$shift goto more
$label declared
file.nr=3;
$if dimension 0 %1 $goto scalar
$if not declared u1 alias(u1,u2,u3,u4,u5,u6,u7,u8,u9,*);
$if dimension 1 %1 $set r 'u1' set c '' set row row1
$if dimension 2 %1 $set r 'u1' set c '' set row row2
$if dimension 3 %1 $set r 'u1' set c 'u2' set row row3
$if dimension 4 %1 $set r 'u1,u2' set c 'u3' set row row4
$if dimension 5 %1 $set r 'u1,u2,u3' set c 'u4' set row row5
$if dimension 6 %1 $set r 'u1,u2,u3,u4' set c 'u5' set row row6
$if dimension 7 %1 $set r 'u1,u2,u3,u4,u5' set c 'u6' set row row7
$if dimension 8 %1 $set r 'u1,u2,u3,u4,u5,u6' set c 'u7' set row row8
$if dimension 9 %1 $set r 'u1,u2,u3,u4,u5,u6,u7' set c 'u8' set row row9
$set rc %c%,%r%
file.tf=3;
$if dimension 1 %1 $goto param1
$if dimension 2 %1 $goto param2
set %row%(%r%),col(%c%);
$onuni
* handle the time-dependant tables by putting out everything
$if '%whtype%' == L $GOTO SUFFIX
$if '%whtype%' == M $GOTO SUFFIX
$GOTO OTHER
$LABEL SUFFIX
*-------------------------------------------
loop((REG,%c%,%r%)$%1.%whtype%(REG,%c%,%r%),
col(%c%) = yes;
%row%(%r%) = yes);
* handle parameters with > 1 index
*$ if '%whtype%' == 'P' $goto param3
scalar colcnt;
IF(CARD(%row%) > 0,
$ set preditem 'table'
put / ' ':9,'%whtype%: Solution Reference ','%1':10," '",%1.ts:40,"'" / @24;
colcnt = 0;
loop(col$(colcnt < 19),
put col.tl:11;
colcnt = colcnt + 1
);
loop((REG,%row%)$(DUMP0 + SUM(col,%1.%whtype%(REG,col,%row%))),
put / @4 REG.TL:0,'.',%row%.te(%row%) @24;
colcnt = 0;
loop(col,
if(colcnt = 19,
colcnt = 0
put / @10 col.tl:4 @20 '+++' @24);
colcnt = colcnt + 1;
if(%1.%whtype%(REG,col,%row%),
put %1.%whtype%(REG,col,%row%):<10 ' '
else
put ' ':11 )));
put /;
);
$offuni
option clear=%row%,clear=col;
$shift goto more
*---------------------------------
$LABEL OTHER
* handle the time-dependant tables by putting out everything
$if NOT '%whtype%' == 'P' col(%whtype%) = yes;
loop((REG,%c%,%r%)$%1(REG,%c%,%r%),
col(%c%) = yes;
%row%(%r%) = yes);
* handle parameters with > 1 index
$ if '%whtype%' == 'P' $goto param3
scalar colcnt;
IF(CARD(%row%) > 0,
$ set preditem 'table'
put / ' ':9,'TABLE':6,'%1':10," '",%1.ts:40,"'" / @44;
colcnt = 0;
loop(col$(colcnt < 19),
put col.tl:11;
colcnt = colcnt + 1
);
loop((REG,%row%)$(DUMP0 + SUM(col,%1(REG,col,%row%))),
put / @4 REG.TL:0,'.',%row%.te(%row%) @44;
colcnt = 0;
loop(col,
if(colcnt = 19,
colcnt = 0
put / @10 col.tl:4 @40 '+++' @44);
colcnt = colcnt + 1;
if(%1(REG,col,%row%),
put %1(REG,col,%row%):<10 ' '
else
put ' ':11 )));
put /;
);
$goto donetab
* handle parameters with >2 index
$label param3
IF(CARD(%row%) > 0,
put / ' ':9,'PARAMETER'@20,'%1':10," '",%1.ts:40,"'";
loop((REG,col,%row%)$(DUMP0 + (%1(REG,col,%row%) NE 0)),
put / @15, REG.TL:0,'.'col.tl'.'%row%.te(%row%) @66 %1(REG,col,%row%):<10; );
put /;
);
$label donetab
$offuni
option clear=%row%,clear=col;
$shift goto more
$label scalar
* simple scalar
$if not '%preditem%' == 'scalar' put / ' ':9, 'SCALARS' /
$set preditem 'scalar'
put ' ':15,'%1':10,"/ " %1:10 " /" /;
$shift goto more
$label param2
$set preditem 'parameter'
$hidden works for any dimension
$set rc %r%
set %row%(%rc%);
$ onuni
LOOP(REG,
%row%(%rc%)$%1(REG,%rc%) = yes;
);
IF(CARD(%row%) > 0,
put / ' ':9,'PARAMETER' , @20,'%1':10,"'",%1.ts:40,"' /";
loop((REG,%row%)$(DUMP0 + (%1(REG,%row%) NE 0)),
put / @15, REG.TL:0'.'%row%.te(%row%) @44 %1(REG,%row%):<10; );
put @69 '/' /;
option clear=%row%;
put /;
);
$ offuni
$shift goto more
$label param1
$set preditem 'parameter'
$hidden works for any dimension
$set rc %r%
set %row%(%rc%);
$ onuni
%row%(%rc%)$%1(%rc%) = yes;
IF(CARD(%row%) > 0,
put / ' ':9,'PARAMETER' , @20,'%1':10,"'",%1.ts:40,"' /";
loop(%row%$(DUMP0 + (%1(%row%) NE 0)),
put / @15, %row%.te(%row%) @44 %1(%row%):<10; );
put @69 '/' /;
option clear=%row%;
put /;
);
$offuni
option clear=%row%,clear=col;
$shift goto more
$label done