-
Notifications
You must be signed in to change notification settings - Fork 83
/
pglib_opf_case39_epri.m
307 lines (302 loc) · 18.3 KB
/
pglib_opf_case39_epri.m
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
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% %%%%%
%%%% IEEE PES Power Grid Library - Optimal Power Flow - v23.07 %%%%%
%%%% (https://github.com/power-grid-lib/pglib-opf) %%%%%
%%%% Benchmark Group - Typical Operations %%%%%
%%%% 23 - July - 2023 %%%%%
%%%% %%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% CASE39 Power flow data for 39 bus New England system.
%
% Data taken from [1] with the following modifications/additions:
%
% - renumbered gen buses consecutively (as in [2] and [4])
% - added Pmin = 0 for all gens
% - added Qmin, Qmax for gens at 31 & 39 (copied from gen at 35)
% - added Vg based on V in bus data (missing for bus 39)
% - added Vg, Pg, Pd, Qd at bus 39 from [2] (same in [4])
% - added Pmax at bus 39: Pmax = Pg + 100
% - added line flow limits and area data from [4]
% - added voltage limits, Vmax = 1.06, Vmin = 0.94
% - added identical quadratic generator costs
% - increased Pmax for gen at bus 34 from 308 to 508
% (assumed typo in [1], makes initial solved case feasible)
% - re-solved power flow
%
% Notes:
% - Bus 39, its generator and 2 connecting lines were added
% (by authors of [1]) to represent the interconnection with
% the rest of the eastern interconnect, and did not include
% Vg, Pg, Qg, Pd, Qd, Pmin, Pmax, Qmin or Qmax.
% - As the swing bus, bus 31 did not include and Q limits.
% - The voltages, etc in [1] appear to be quite close to the
% power flow solution of the case before adding bus 39 with
% it's generator and connecting branches, though the solution
% is not exact.
% - Explicit voltage setpoints for gen buses are not given, so
% they are taken from the bus data, however this results in two
% binding Q limits at buses 34 & 37, so the corresponding
% voltages have probably deviated from their original setpoints.
% - The generator locations and types are as follows:
% 1 30 hydro
% 2 31 nuke01
% 3 32 nuke02
% 4 33 fossil02
% 5 34 fossil01
% 6 35 nuke03
% 7 36 fossil04
% 8 37 nuke04
% 9 38 nuke05
% 10 39 interconnection to rest of US/Canada
%
% This is a solved power flow case, but it includes the following
% violations:
% - Pmax violated at bus 31: Pg = 677.87, Pmax = 646
% - Qmin violated at bus 37: Qg = -1.37, Qmin = 0
%
% References:
% [1] G. W. Bills, et.al., "On-Line Stability Analysis Study"
% RP90-1 Report for the Edison Electric Institute, October 12, 1970,
% pp. 1-20 - 1-35.
% prepared by E. M. Gulachenski - New England Electric System
% J. M. Undrill - General Electric Co.
% "generally representative of the New England 345 KV system, but is
% not an exact or complete model of any past, present or projected
% configuration of the actual New England 345 KV system.
% [2] M. A. Pai, Energy Function Analysis for Power System Stability,
% Kluwer Academic Publishers, Boston, 1989.
% (references [3] as source of data)
% [3] Athay, T.; Podmore, R.; Virmani, S., "A Practical Method for the
% Direct Analysis of Transient Stability," IEEE Transactions on Power
% Apparatus and Systems , vol.PAS-98, no.2, pp.573-584, March 1979.
% URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4113518&isnumber=4113486
% (references [1] as source of data)
% [4] Data included with TC Calculator at http://www.pserc.cornell.edu/tcc/
% for 39-bus system.
%
% Copyright (c) 1989 by The Institute of Electrical and Electronics Engineers (IEEE)
% Licensed under the Creative Commons Attribution 4.0
% International license, http://creativecommons.org/licenses/by/4.0/
%
% Contact M.E. Brennan (me.brennan@ieee.org) for inquries on further reuse of
% this dataset.
%
function mpc = pglib_opf_case39_epri
mpc.version = '2';
mpc.baseMVA = 100.0;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 1 97.6 44.2 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
2 1 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
3 1 322.0 2.4 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
4 1 500.0 184.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
5 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
6 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
7 1 233.8 84.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
8 1 522.0 176.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
9 1 6.5 -66.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
10 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
11 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
12 1 8.53 88.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
13 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
14 1 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
15 1 320.0 153.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
16 1 329.0 32.3 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
17 1 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
18 1 158.0 30.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
19 1 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
20 1 680.0 103.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
21 1 274.0 115.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
22 1 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
23 1 247.5 84.6 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
24 1 308.6 -92.2 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
25 1 224.0 47.2 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
26 1 139.0 17.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
27 1 281.0 75.5 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
28 1 206.0 27.6 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
29 1 283.5 26.9 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
30 2 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
31 3 9.2 4.6 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
32 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
33 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
34 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
35 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
36 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
37 2 0.0 0.0 0.0 0.0 2 1.00000 0.00000 345.0 1 1.06000 0.94000;
38 2 0.0 0.0 0.0 0.0 3 1.00000 0.00000 345.0 1 1.06000 0.94000;
39 2 1104.0 250.0 0.0 0.0 1 1.00000 0.00000 345.0 1 1.06000 0.94000;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
mpc.gen = [
30 520.0 270.0 400.0 140.0 1.0 100.0 1 1040.0 0.0; % NUC
31 323.0 100.0 300.0 -100.0 1.0 100.0 1 646.0 0.0; % COW
32 362.5 225.0 300.0 150.0 1.0 100.0 1 725.0 0.0; % COW
33 326.0 125.0 250.0 0.0 1.0 100.0 1 652.0 0.0; % COW
34 254.0 83.5 167.0 0.0 1.0 100.0 1 508.0 0.0; % COW
35 343.5 100.0 300.0 -100.0 1.0 100.0 1 687.0 0.0; % COW
36 290.0 120.0 240.0 0.0 1.0 100.0 1 580.0 0.0; % COW
37 282.0 125.0 250.0 0.0 1.0 100.0 1 564.0 0.0; % COW
38 432.5 75.0 300.0 -150.0 1.0 100.0 1 865.0 0.0; % COW
39 550.0 100.0 300.0 -100.0 1.0 100.0 1 1100.0 0.0; % COW
];
%% generator cost data
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0.0 0.0 3 0.000000 6.724778 0.000000; % NUC
2 0.0 0.0 3 0.000000 14.707625 0.000000; % COW
2 0.0 0.0 3 0.000000 24.804734 0.000000; % COW
2 0.0 0.0 3 0.000000 34.844643 0.000000; % COW
2 0.0 0.0 3 0.000000 24.652994 0.000000; % COW
2 0.0 0.0 3 0.000000 32.306483 0.000000; % COW
2 0.0 0.0 3 0.000000 18.157477 0.000000; % COW
2 0.0 0.0 3 0.000000 31.550181 0.000000; % COW
2 0.0 0.0 3 0.000000 22.503168 0.000000; % COW
2 0.0 0.0 3 0.000000 27.434444 0.000000; % COW
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
1 2 0.0035 0.0411 0.6987 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
1 39 0.001 0.025 0.75 1000.0 1000.0 1000.0 0.0 0.0 1 -30.0 30.0;
2 3 0.0013 0.0151 0.2572 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0;
2 25 0.007 0.0086 0.146 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0;
2 30 0.0 0.0181 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0;
3 4 0.0013 0.0213 0.2214 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0;
3 18 0.0011 0.0133 0.2138 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0;
4 5 0.0008 0.0128 0.1342 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
4 14 0.0008 0.0129 0.1382 500.0 500.0 500.0 0.0 0.0 1 -30.0 30.0;
5 6 0.0002 0.0026 0.0434 1200.0 1200.0 1200.0 0.0 0.0 1 -30.0 30.0;
5 8 0.0008 0.0112 0.1476 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
6 7 0.0006 0.0092 0.113 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
6 11 0.0007 0.0082 0.1389 480.0 480.0 480.0 0.0 0.0 1 -30.0 30.0;
6 31 0.0 0.025 0.0 1800.0 1800.0 1800.0 1.07 0.0 1 -30.0 30.0;
7 8 0.0004 0.0046 0.078 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
8 9 0.0023 0.0363 0.3804 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
9 39 0.001 0.025 1.2 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
10 11 0.0004 0.0043 0.0729 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
10 13 0.0004 0.0043 0.0729 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
10 32 0.0 0.02 0.0 900.0 900.0 2500.0 1.07 0.0 1 -30.0 30.0;
12 11 0.0016 0.0435 0.0 500.0 500.0 500.0 1.006 0.0 1 -30.0 30.0;
12 13 0.0016 0.0435 0.0 500.0 500.0 500.0 1.006 0.0 1 -30.0 30.0;
13 14 0.0009 0.0101 0.1723 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
14 15 0.0018 0.0217 0.366 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
15 16 0.0009 0.0094 0.171 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
16 17 0.0007 0.0089 0.1342 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
16 19 0.0016 0.0195 0.304 600.0 600.0 2500.0 0.0 0.0 1 -30.0 30.0;
16 21 0.0008 0.0135 0.2548 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
16 24 0.0003 0.0059 0.068 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
17 18 0.0007 0.0082 0.1319 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
17 27 0.0013 0.0173 0.3216 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
19 20 0.0007 0.0138 0.0 900.0 900.0 2500.0 1.06 0.0 1 -30.0 30.0;
19 33 0.0007 0.0142 0.0 900.0 900.0 2500.0 1.07 0.0 1 -30.0 30.0;
20 34 0.0009 0.018 0.0 900.0 900.0 2500.0 1.009 0.0 1 -30.0 30.0;
21 22 0.0008 0.014 0.2565 900.0 900.0 900.0 0.0 0.0 1 -30.0 30.0;
22 23 0.0006 0.0096 0.1846 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
22 35 0.0 0.0143 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0;
23 24 0.0022 0.035 0.361 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
23 36 0.0005 0.0272 0.0 900.0 900.0 2500.0 0.0 0.0 1 -30.0 30.0;
25 26 0.0032 0.0323 0.531 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
25 37 0.0006 0.0232 0.0 900.0 900.0 2500.0 1.025 0.0 1 -30.0 30.0;
26 27 0.0014 0.0147 0.2396 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
26 28 0.0043 0.0474 0.7802 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
26 29 0.0057 0.0625 1.029 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
28 29 0.0014 0.0151 0.249 600.0 600.0 600.0 0.0 0.0 1 -30.0 30.0;
29 38 0.0008 0.0156 0.0 1200.0 1200.0 2500.0 1.025 0.0 1 -30.0 30.0;
];
% INFO : === Translation Options ===
% INFO : Phase Angle Bound: 30.0 (deg.)
% INFO : Gen Active Cost Model: stat
% INFO : Setting Flat Start
% INFO :
% INFO : === Generator Classification Notes ===
% INFO : NUC 1 - 3.97
% INFO : COW 9 - 96.03
% INFO :
% INFO : === Generator Active Cost Stat Model Notes ===
% INFO : Updated Generator Cost: NUC - 0.2 0.3 0.01 -> 0 6.72477811338 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 14.7076252515 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 24.8047338223 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 34.84464286 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 24.6529938713 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 32.306483114 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 18.1574766212 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 31.5501806971 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 22.5031675377 0
% INFO : Updated Generator Cost: COW - 0.2 0.3 0.01 -> 0 27.4344440931 0
% INFO :
% INFO : === Generator Bounds Update Notes ===
% INFO :
% INFO : === Base KV Replacement Notes ===
% INFO :
% INFO : === Transformer Setting Replacement Notes ===
% WARNING : Transformer 23-36 connects the same voltage levels (345.0, 345.0) and has no phase shift, changing tap ratio 1.0 => 0.0
% INFO :
% INFO : === Line Capacity Monotonicity Notes ===
% INFO :
% INFO : === Voltage Setpoint Replacement Notes ===
% INFO : Bus 1 : V=1.0393836, theta=-13.536602 -> V=1.0, theta=0.0
% INFO : Bus 2 : V=1.0484941, theta=-9.7852666 -> V=1.0, theta=0.0
% INFO : Bus 3 : V=1.0307077, theta=-12.276384 -> V=1.0, theta=0.0
% INFO : Bus 4 : V=1.00446, theta=-12.626734 -> V=1.0, theta=0.0
% INFO : Bus 5 : V=1.0060063, theta=-11.192339 -> V=1.0, theta=0.0
% INFO : Bus 6 : V=1.0082256, theta=-10.40833 -> V=1.0, theta=0.0
% INFO : Bus 7 : V=0.99839728, theta=-12.755626 -> V=1.0, theta=0.0
% INFO : Bus 8 : V=0.99787232, theta=-13.335844 -> V=1.0, theta=0.0
% INFO : Bus 9 : V=1.038332, theta=-14.178442 -> V=1.0, theta=0.0
% INFO : Bus 10 : V=1.0178431, theta=-8.170875 -> V=1.0, theta=0.0
% INFO : Bus 11 : V=1.0133858, theta=-8.9369663 -> V=1.0, theta=0.0
% INFO : Bus 12 : V=1.000815, theta=-8.9988236 -> V=1.0, theta=0.0
% INFO : Bus 13 : V=1.014923, theta=-8.9299272 -> V=1.0, theta=0.0
% INFO : Bus 14 : V=1.012319, theta=-10.715295 -> V=1.0, theta=0.0
% INFO : Bus 15 : V=1.0161854, theta=-11.345399 -> V=1.0, theta=0.0
% INFO : Bus 16 : V=1.0325203, theta=-10.033348 -> V=1.0, theta=0.0
% INFO : Bus 17 : V=1.0342365, theta=-11.116436 -> V=1.0, theta=0.0
% INFO : Bus 18 : V=1.0315726, theta=-11.986168 -> V=1.0, theta=0.0
% INFO : Bus 19 : V=1.0501068, theta=-5.4100729 -> V=1.0, theta=0.0
% INFO : Bus 20 : V=0.99101054, theta=-6.8211783 -> V=1.0, theta=0.0
% INFO : Bus 21 : V=1.0323192, theta=-7.6287461 -> V=1.0, theta=0.0
% INFO : Bus 22 : V=1.0501427, theta=-3.1831199 -> V=1.0, theta=0.0
% INFO : Bus 23 : V=1.0451451, theta=-3.3812763 -> V=1.0, theta=0.0
% INFO : Bus 24 : V=1.038001, theta=-9.9137585 -> V=1.0, theta=0.0
% INFO : Bus 25 : V=1.0576827, theta=-8.3692354 -> V=1.0, theta=0.0
% INFO : Bus 26 : V=1.0525613, theta=-9.4387696 -> V=1.0, theta=0.0
% INFO : Bus 27 : V=1.0383449, theta=-11.362152 -> V=1.0, theta=0.0
% INFO : Bus 28 : V=1.0503737, theta=-5.9283592 -> V=1.0, theta=0.0
% INFO : Bus 29 : V=1.0501149, theta=-3.1698741 -> V=1.0, theta=0.0
% INFO : Bus 30 : V=1.0499, theta=-7.3704746 -> V=1.0, theta=0.0
% INFO : Bus 31 : V=0.982, theta=0.0 -> V=1.0, theta=0.0
% INFO : Bus 32 : V=0.9841, theta=-0.1884374 -> V=1.0, theta=0.0
% INFO : Bus 33 : V=0.9972, theta=-0.19317445 -> V=1.0, theta=0.0
% INFO : Bus 34 : V=1.0123, theta=-1.631119 -> V=1.0, theta=0.0
% INFO : Bus 35 : V=1.0494, theta=1.7765069 -> V=1.0, theta=0.0
% INFO : Bus 36 : V=1.0636, theta=4.4684374 -> V=1.0, theta=0.0
% INFO : Bus 37 : V=1.0275, theta=-1.5828988 -> V=1.0, theta=0.0
% INFO : Bus 38 : V=1.0265, theta=3.8928177 -> V=1.0, theta=0.0
% INFO : Bus 39 : V=1.03, theta=-14.535256 -> V=1.0, theta=0.0
% INFO :
% INFO : === Generator Setpoint Replacement Notes ===
% INFO : Gen at bus 30 : Pg=250.0, Qg=161.762 -> Pg=520.0, Qg=270.0
% INFO : Gen at bus 30 : Vg=1.0499 -> Vg=1.0
% INFO : Gen at bus 31 : Pg=677.871, Qg=221.574 -> Pg=323.0, Qg=100.0
% INFO : Gen at bus 31 : Vg=0.982 -> Vg=1.0
% INFO : Gen at bus 32 : Pg=650.0, Qg=206.965 -> Pg=362.5, Qg=225.0
% INFO : Gen at bus 32 : Vg=0.9841 -> Vg=1.0
% INFO : Gen at bus 33 : Pg=632.0, Qg=108.293 -> Pg=326.0, Qg=125.0
% INFO : Gen at bus 33 : Vg=0.9972 -> Vg=1.0
% INFO : Gen at bus 34 : Pg=508.0, Qg=166.688 -> Pg=254.0, Qg=83.5
% INFO : Gen at bus 34 : Vg=1.0123 -> Vg=1.0
% INFO : Gen at bus 35 : Pg=650.0, Qg=210.661 -> Pg=343.5, Qg=100.0
% INFO : Gen at bus 35 : Vg=1.0494 -> Vg=1.0
% INFO : Gen at bus 36 : Pg=560.0, Qg=100.165 -> Pg=290.0, Qg=120.0
% INFO : Gen at bus 36 : Vg=1.0636 -> Vg=1.0
% INFO : Gen at bus 37 : Pg=540.0, Qg=-1.36945 -> Pg=282.0, Qg=125.0
% INFO : Gen at bus 37 : Vg=1.0275 -> Vg=1.0
% INFO : Gen at bus 38 : Pg=830.0, Qg=21.7327 -> Pg=432.5, Qg=75.0
% INFO : Gen at bus 38 : Vg=1.0265 -> Vg=1.0
% INFO : Gen at bus 39 : Pg=1000.0, Qg=78.4674 -> Pg=550.0, Qg=100.0
% INFO : Gen at bus 39 : Vg=1.03 -> Vg=1.0
% INFO :
% INFO : === Writing Matpower Case File Notes ===