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
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
|
BEGIN NEW DATA CASE
C BENCHMARK DC-26
C Test of Type-59 S.M. dynamics, with load flow ("FIX SOURCE" request)
C determination of initial conditions that observe given power constraints.
C 1st of 5 subcases involving load flow solutions; 1st of 2 machine cases.
C Answers change slightly on 10 February 1999 following the massive
C changes from TEPCO (Tokyo Electric Power Company) in Japan. See April
C newsletter. Most extrema agree to 5 or 6 decimal digits.
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
POWER FREQUENCY, 50.0, { Avoid warning message about suspicious XOPT, COPT.
.000200 .500 50. 50.
1 1 0 { test } 1 -1
5 5 20 20 100 100
A1 A2 0.1 1
B1 B2 0.1 1
C1 C2 0.1 1
51A2 B2 .00528 571.5428
52A3 -.029 3125.811 .1582 17115.12
51B2 C2 A2 B2
52B3
51C2 A2 A2 B2
52C3
A2 31.416
B2 31.416
C2 31.416
1A3 A4 43.342
2B3 B4 12.546 43.342
3C3 C4 12.546 12.546 43.342
BLANK card ending branch cards.
C The next card has T-open = 80 msec minus DELTAT/2 in order to avoid
C delayed opening (1/2 cycle) for PRIME and Burroughs. WSM. 27 FEB 1982.
A3 .01990 .0799 1
BLANK card ending switch cards.
14A4 112.059 50. -20. -1.
14B4 112.059 50. -140. -1.
14C4 112.059 50. 100. -1.
C 59A1 11.3901 50.0 - 30.0 (before FIX SOURCE use)
$DEPOSIT, NOSM59=-1 { Special value will change Type-58 S.M. to Type-59 S.M.
C Note about preceding and following. Beginning 27 March 2000, this new
C variable NOSM59 of STARTUP can be used either to convert or to trap
C Type-58 or Type-59 S.M. type codes. Five values are recognized:
C NOSM59 = 0 ===> no change. Preserve the original meaning.
C 1 convert any Type-59 S.M. to Type-58 S.M.
C 2 trap any Type-59 S.M. (halt after all are counted).
C -1 convert any Type-58 S.M. to Type-59 S.M.
C -2 trap any Type-58 S.M. (halt after all are counted).
C So, the following machine really is Type-59. But 58 is used because ATP
C automatically will convert this to 59 as a result of preceding NOSM59 = -1
C This preserves the old answers, which will remain unchanged.
58A1 11.0 50.0 -40.0
58B1
58C1
C TOLERANCES 200. { Columns 11-20 EPSUPA = 200 } 20 { 51-60: NIOMAX=20
C 14 October 2001, define NIOMIN of columns 41-50 so Watcom agrees with Mingw32
C and Salford. This _does_ help. Value 8 provides better FC agreement than 5.
C Previously, Watcom was iterating a different number of times, & was producing
C a slightly different answer. So, we slow execution artificially:
TOLERANCES 200. 8 20 { 51-60: NIOMAX=20
PARAMETER FITTING 1.0
1 1 2 1.0 1.0 150.0 13.8 600.0 600.0 720.0
BLANK card for quadrature axis of machine.
0.0014 0.175 1.85 1.76 0.2575 1.76 0.18 0.18
5.74757 0.051142 0.382609 0.197985
1 1.0 50.0 1.0
BLANK card ending all (here, just one) mass cards.
12
21
3 1
BLANK card terminating all (here, three) S.M. output requests.
FINISH
BLANK card terminating all EMTP source cards.
A1 39.98755 -4.507399 10.0 12.0 -60.
B1 39.98755 -4.507399 10.0 12.0 -180. -120.
C1 39.98755 -4.507399 10.0 12.0 60.0 120.
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 1000 20 1 0.00001 0.1 2.0 2
C Max del-V: .7E-5 .7E-5 .6E-5 .6E-5 .5E-5 .5E-5 .5E-5 .4E-5 .4E-5 .4E-5 .3E-5
C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3
C Exit the load flow iteration loop with counter NEKITE = 133. If no warning
C Row Node Name Voltage magnit Degrees Real power P Reactive power
C 4 8 A1 1.13901111E+01 -30.00022 3.99871502E+01 -4.50736252E+00
C 5 9 B1 1.13901111E+01 -150.00022 3.99871502E+01 -4.50736252E+00
C 6 10 C1 1.13901111E+01 89.99978 3.99871502E+01 -4.50736252E+00
C Total network loss P-loss by summing injections = 8.018000511887E+00
C Step Time C2 B2 A2 A3 C3
C
C C1 A4 B4 C4 A3
C TERRA
C
C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1
C ID IQ I0 IF IKD
C
C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1
C IB IC EFD MFORCE MANG
C 0 0.0 .0791881332 -9.29565074 9.216462609 102.5472108 -37.8516127
C .4577064E-4 105.3010154 -85.8421743 -19.4588411 0.0
C -6.82107489 5.325616419 .391002E-14 1.014337449 -.76501E-15
C -5.68498411 -.791423626 -.176005017 .60792369 .887825386
C 1 .2E-3 -.592823298 -8.9455915 9.538414794 103.3173764 -43.8432677
C -.715886965 107.4997647 -81.1499746 -26.3497901 0.0
C -6.82114583 5.325512789 -.4718E-14 1.014337757 .6231167E-5
C -5.41035524 -1.23063667 -.176005017 .6079236529 .8878253793
1 { Request for the output of all node voltages
C Last step begin: 2500 0.5 .1013093411 -9.26844034 9.167130999 16.94173877
C Last step begin: .0149766932 105.3010154 -85.8421743 -19.4588411 0.0
C Last step begin: -6.73361875 5.375255555 -.34086E-14 1.004763786 -.001205227
C Last step begin: -5.61469909 -.863326479 -.176005017 .6045428409 .889016932
C Variable max : 11.6152661 11.51967968 12.3018062 204.4814649 119.2607541
C 12.40458274 112.0153099 112.0562692 112.0480769 6.3761474
C -2.00131682 9.101894661 .13086E-13 1.315246828 1.316347492
C 12.08207638 7.893166343 -.176005017 .60792369 .8935781536
C Times of max : .0938 .1282 .0806 .0802 .0934
C .0938 .4812 .4878 .2344 .0252
C .0782 .0708 .1114 .074 .0332
C .033 .0938 0.0 0.0 .0226
PRINTER PLOT
19410. 0. 120. MACH 1ID MACH 1IQ { Plot limits: (-2.301, 0.910)
C Note: Former IG (next plot) now (after Prof. Juan Martinez's changes of
C October, 1989) becomes IKQ. The eddy-current coil of the Q axis
C formerly was coil "KQ", now is "G". WSM and JM, 13 October 90
19410. 0. 120. MACH 1IF MACH 1IKQ { Plot limits: (-0.425, 1.636)
19410. 0. 120. BRANCH { Plot limits: (-1.874, 1.755)
A1 A2 B1 B2 C1 C2
$DEPOSIT, NOSM59=0 { Cancel special value that was set earlier in this subcase
BLANK card terminating all plot cards.
BEGIN NEW DATA CASE
C 2nd of 6 data subcases.
C Same as preceding case, only with significant imbalance in the steady
C state due to change of generator at "A4" from -20 to -10 degrees.
C The load flow does not converge, and case is only run 5 steps. The
C idea is to force use of negative-sequence correction logic of SMINIT
C of overlay 11. No other test cases do this. Eventually, replace by
C a meaningful case. 2nd of 4 subcases (2nd of two machine cases).
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
POWER FREQUENCY, 50.0, { Avoid warning message about suspicious XOPT, COPT.
.000200 .001 50. 50.
1 1 0 { test } 0 -1
5 5 20 20 100 100
A1 A2 0.1 1
B1 B2 0.1 1
C1 C2 0.1 1
51A2 B2 .00528 571.5428
52A3 -.029 3125.811 .1582 17115.12
51B2 C2 A2 B2
52B3
51C2 A2 A2 B2
52C3
A2 31.416
B2 31.416
C2 31.416
1A3 A4 43.342
2B3 B4 12.546 43.342
3C3 C4 12.546 12.546 43.342
BLANK card ending branch cards.
C The next card has T-open = 80 msec minus DELTAT/2 in order to avoid
C delayed opening (1/2 cycle) for PRIME and Burroughs. WSM. 27 FEB 1982.
A3 .01990 .0799 1
BLANK card ending switch cards.
14A4 112.059 50. -10. { Change -20 to -10 } -1.
14B4 112.059 50. -140. -1.
14C4 112.059 50. 100. -1.
C 59A1 11.3901 50.0 - 30.0 (before FIX SOURCE use)
59A1 11.0 50.0 -40.0
59B1
59C1
TOLERANCES 200. { Columns 11-20 EPSUPA = 200 } 20 { 51-60: NIOMAX=20
PARAMETER FITTING 1.0
1 1 2 1.0 1.0 150.0 13.8 600.0 600.0 720.0
BLANK card for quadrature axis of machine.
0.0014 0.175 1.85 1.76 0.2575 1.76 0.18 0.18
5.74757 0.051142 0.382609 0.197985
1 1.0 50.0 1.0
BLANK card ending all (here, just one) mass cards.
12
21
3 1
BLANK card terminating all (here, three) S.M. output requests.
FINISH
BLANK card terminating all EMTP source cards.
A1 39.98755 -4.507399 10.0 12.0 -60.
B1 39.98755 -4.507399 10.0 12.0 -180. -120.
C1 39.98755 -4.507399 10.0 12.0 60.0 120.
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 100 20 1 0.00001 0.1 2.0 2
C Max del-V: .0016 .0015 .0013 .0011 .001 .9E-3 .8E-3 .6E-3 .5E-3 .4E-3 .4E-3
C //. Warning! Divergence of the EMTP load flow has occurred. In the hope
C initialization continues, however.
C Exit the load flow iteration loop with counter NEKITE = 100. If no warning
C Row Node Name Voltage magnit Degrees Real power P Reactive power
C 4 8 A1 1.00000000E+01 -20.91363 3.45426843E+01 -2.24391750E+00
C 5 9 B1 1.20000000E+01 -150.14225 4.22174603E+01 -6.71602291E+00
C 6 10 C1 1.20000000E+01 90.73965 3.92028723E+01 -5.90634496E+00
C Total network loss P-loss by summing injections = 7.625614104022E+00
1 { Request for the output of all node voltages
C Step Time C2 B2 A2 A3 C3
C
C C1 A4 B4 C4 A3
C TERRA
C
C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1
C ID IQ I0 IF IKD
C
C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1
C IB IC EFD MFORCE MANG
C 0 0.0 -.048041917 -9.85266713 8.679836062 104.3486527 -36.7261188
C -.154907255 110.356572 -85.8421743 -19.4588411 0.0
C -6.92581041 5.259477223 -.445329E-4 .9859361614 .0320876279
C -5.54502117 -1.06865338 -.170793823 .5995965889 .8894202603
C 1 .2E-3 -1.16677225 -7.85747584 10.21232132 103.2738496 -48.1646884
C -1.32143162 111.3606393 -81.1499746 -26.3497901 0.0
C -6.87231021 5.289979997 .0014353419 .9858026654 .0272832017
C -5.20835729 -1.54659375 -.170793823 .5996092398 .8894186749
C Last step follows:
C 5 .001 -3.69592471 -6.2063594 10.61819433 99.27645368 -67.8903465
C -4.03363622 110.9684495 -59.3822228 -52.6085139 0.0
C -6.76271705 5.471445353 .0060095523 .9856033893 .017362355
C -3.71587997 -3.37711507 -.170793823 .5996729546 .8894772231
PRINTER PLOT
BLANK card terminating all plot cards.
BEGIN NEW DATA CASE
C 5-bus, 1-phase loadflow example from "Elements of Power System Analysis"
C by William D. Stevenson. McGraw-Hill Book Company, 2nd Edition (1962).
C The Gauss-Seidel solution is contained in Example 10.1, page 219 onward.
C Data is given in per unit, and we can use this directly (pretend that
C the base voltage, current, and power are all unity). We need XOPT = 60
C since impedances are in ohms. Voltages could be multiplied by 1.414 to
C give the peak values required by the EMTP, but instead we use the VSCALE
C option that takes RMS (the square root of 2 is supplied internally) as
C input and produces RMS output (at least for the load flow solution table
C if not for EMTP branch flows). Bus number one, WHITE, is taken as the
C slack bus (the only bus without a power constraint). This is the 3rd of
C 5 subcases. It has no machine dynamics at all (unlike the preceding 2).
C This data case was 1st considered by Profs. Saul Goldberg & Bill Horton
C of Cal Poly in San Luis Obispo, California, for June, 1988, EMTP course.
C Answers change 10 February 1999 following the massive changes from
C TEPCO (Tokyo Electric Power Company) in Japan. See April newsletter.
C Convergence has been slowed substantially (when time, consider this).
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
$ERASE { Erase punched cards, created by preceding subcases, from punch buffer
.000200 .000 60. { T-max = 0 means that no transient solution follows
1 1 1 0 1
WHITE RED .10 .40
C Stevenson's line between WHITE and BLUE appears on the final comment
C card of this comment block. For more generality, I have introduced an
C intermediate node DUMMY that is unconstrained. The solution of the
C network should be unchanged, of course. It is (original is ok, too).
C WHITE BLUE .15 .60
WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X
DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X
WHITE YELLOW .05 .20
RED GREEN .05 .20
RED BLUE .10 .40
GREEN YELLOW .05 .20
BLANK card ending branch cards.
BLANK card ending switch cards.
C The following 2 cards would be used if peak rather than RMS input data:
C 14WHITE 1.4425 60. 0.0 -1.
C 14GREEN 1.4708 60. 2.0 -1.
14WHITE 1.02 60. 0.0 -1.
14GREEN 1.04 60. 2.0 -1.
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C WHITE |RED *DUMMY *YELLOW*
C RED |WHITE *BLUE *GREEN *
C DUMMY |WHITE *BLUE *
C BLUE |RED *DUMMY *
C YELLOW |WHITE *GREEN *
C GREEN |RED *YELLOW*
C --------------+------------------------------
BLANK card terminating all EMTP source cards.
C Next come power constraints of the load flow. There will be one
C for each non-slack generator. So, 3 of them will apply to nodes
C that have no Type-14 source as required by the algorithm. The
C program will define these internally. Yet, this is only possible
C for TMAX non-positive (no transient continuation). If data is
C modified to make TMAX > 0, a KILL = 40 error termination will
C will result. Yet, the transient simulation is possible as a
C 2nd simulation that replaces the power constraints at load nodes
C by constant-impedance loads. Branch cards for these will be
C punched by the present subcase, and the transient continuation
C will be illustrated by the following (4th of 4) subcase.
C The following 4 cards would be used if peak rather than RMS input data.
C Note that average of Vmin & Vmax = 0.5 ( 1.0 + 1.828428 ) = 1.414212
C (power constraints RED, BLUE, and YELLOW), and 1.4708 = 1.414 * 1.02
C RED -0.6 -0.3 1.01.828428
C 1GREEN 1.0 1.4708 -20. 20.
C BLUE -0.4 -0.1 1.01.828428
C YELLOW -0.6 -0.2 1.01.828428
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
RED -0.6 -0.3 .85 1.15
1GREEN 1.0 1.02 -20. 20.
BLUE -0.4 -0.1 .85 1.15
YELLOW -0.6 -0.2 .85 1.15
C Max del-V: .0293 .0153 .0141 .0138 .0118 .0096 .0085 .0076 .0068 .0061 .0055
C Source No. -2 -3 -4 -1 -3 -3 -3 -3 -3 -3 -3
C Max del-V: .0023 .0021 .002 .0018 .0017 .0016 .0015 .0014 .0013 .0012 .0011
C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3
C Max del-V: .6E-3 .5E-3 .5E-3 .5E-3 .4E-3 .4E-3 .4E-3 .4E-3 .4E-3 .3E-3 .3E-3
C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3
C The following load-flow miscellaneous data card has two peculiarities. The
C use of VSCALE = 1.414 is the special flag requesting RMS rather than peak
C voltages. The use of KTAPER = 0 ensures constant acceleration factors
C (this works well for this problem whether RMS or peak values are used).
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 500 20 1 .00001 0.0 0.0 1.414 0
C Max del-V: .6E-6 .6E-6 .6E-6 .5E-6
C Source No. -3 -3 -3 -3
C Exit the load flow iteration loop with counter NEKITE = 164. If no warning
C Row Node Name Voltage magnit Degrees Real power P Reactive power
C 3 3 RED 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01
C 2 7 GREEN 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01
C 4 5 BLUE 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02
C 5 6 YELLOW 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01
C Total network loss P-loss by summing injections = 5.282912652566E-02
C Next, show the punched cards that $PUNCH (next) will display:
C $UNITS card should give X, C units XUNITS, CUNITS = 2.6526E-03 1.0000E-06
C Following are branch cards for either loads or generator impedances:
C $VINTAGE, 1,
C Polar V, P, Q = 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01
C RED .591308166 .295659302
C Polar V, P, Q = 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01
C GREEN -.45374803 15276.6334
C Polar V, P, Q = 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02
C BLUE .985230493 .246313607
C Polar V, P, Q = 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01
C YELLOW .724153874 .241387449
C $VINTAGE, -1, { Last of punched cardss that the following $PUNCH will display
$PUNCH { Flush punched cards: R-L branches for equivalent impedance loads
C Note about preceding card. Since no images anyway, remove. WSM. 15 Mar 98
BLANK card ending requests for output variables
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 6 data subcases is a continuation of the 3rd into the time-step
C loop. The network differs in that constant-impedance loads (using the
C branch cards that were punched by execution of the 3rd subcase) replace
C the original power constraints at nodes without Type-14 sources. The
C time-step size DELTAT and end-time TMAX have been arbitrarily to
C selected to show that the transient continuation is smooth for just over
C a cycle. Note that this continuation involves a load flow, too, since
C otherwise angles of the Type-14 generators would not correspond to the
C network solution. The load flow could be omitted if the user were
C willing to re-key generator angles, but it generally is easier to let
C the computer solve the load flow a second time (this time with constant-
C impedance loads).
C Branch voltage output of the time-step loop, if declared by a 2-punch
C in column 80, was incompatible with load flow prior to 21 July 2002.
C So, to prove it now works, add 2-punch to 1st branch on this date. An
C older version of the program will simply ignore the column-80 punch.
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
.000200 .020 60. { T-max = 0 means that no transient solution follows
1 1 1 0 1 -1
5 5 20 20
WHITE RED .10 .40 2
WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X
DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X
WHITE YELLOW .05 .20
RED GREEN .05 .20
RED BLUE .10 .40
GREEN YELLOW .05 .20
C <++++++> Cards punched by support routine on 20-Oct-90 05.01.52 <++++++>
$UNITS, 0.0, 0.0 { 1st of cards punched by preceding 3rd subcase: mHenry, microF
$VINTAGE, 1, { 2nd of cards punched by preceding 3rd subcase: use wide formt
C Polar V, P, Q = 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01
RED 1.18261633E+00 1.56852131E+00
C Polar V, P, Q = 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02
BLUE 1.97046099E+00 1.30673427E+00
C Polar V, P, Q = 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01
YELLOW 1.44830775E+00 1.28060019E+00
$UNITS, -1, -1 { Restore whatever XOPT and COPT existed before preceding $UNITS
$VINTAGE, 0, { Last of cards punched by preceding 3rd subcase
BLANK card ending branch cards.
BLANK card ending switch cards.
14WHITE 1.02 60. 0.0 -1.
14GREEN 1.04 60. 2.0 -1.
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C WHITE |RED *DUMMY *YELLOW*
C RED |TERRA *WHITE *BLUE *GREEN *
C DUMMY |WHITE *BLUE *
C BLUE |TERRA *RED *DUMMY *
C YELLOW |TERRA *WHITE *GREEN *
C GREEN |RED *YELLOW*
C TERRA |RED *BLUE *YELLOW*
C --------------+------------------------------
BLANK card terminating all EMTP source cards.
C Max del-V: .0282 .0015 .0013 .0011 1.E-3 .9E-3 .8E-3 .7E-3 .6E-3 .5E-3 .4E-3
C Source No. -1 1 1 1 1 1 1 1 1 1 1
C Max del-V: .1E-3 .9E-4 .8E-4 .7E-4 .6E-4 .5E-4 .5E-4 .4E-4 .4E-4 .3E-4 .3E-4
C Source No. 1 1 1 1 1 1 1 1 1 1 1
C Max del-V: .7E-5 .6E-5 .5E-5 .4E-5 .4E-5 .3E-5 .3E-5 .3E-5 .2E-5 .2E-5 .2E-5
C Source No. 1 1 1 1 1 1 1 1 1 1 1
C Exit the load flow iteration loop with counter NEKITE = 58. If no warning o
C Row Node Name Voltage magnit Degrees Real power P Reactive power
C 2 7 GREEN 1.02000000E+00 2.43435 9.99985664E-01 3.82679671E-01
C Total network loss P-loss by summing injections = 1.652795946964E+00
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
1GREEN 1.0 1.02 -20. 20.
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 500 20 1 .00001 0.0 0.0 1.414 0
C Step Time YELLOW BLUE DUMMY RED WHITE
C
C 0 0.0 1.38883765 1.28149564 1.336017508 1.328924768 1.442497834
C 1 .2E-3 1.388413609 1.291450977 1.337993527 1.331824678 1.438399553
C 2 .4E-3 1.380101179 1.294070033 1.332367678 1.32715794 1.426128
C 3 .6E-3 1.363947438 1.289337599 1.319171785 1.314950922 1.405752903
1 { Request for the output of all node voltages
C 100 .02 .4736844387 .567774976 .485772514 .4950043848 .4457563449
C Variable max: 1.389411211 1.294070033 1.33815194 1.331824678 1.442497834
C Times of max: .0168 .4E-3 .0168 .2E-3 0.0
C Variable min: -1.38955118 -1.29326058 -1.33749084 -1.33106894 -1.44204228
C Times of min: .0084 .0088 .0084 .0086 .0084
PRINTER PLOT
144 2. 0.0 20. RED { Axis limits: (-1.331, 1.332)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 5th of 6 data subcases is a variation of the 3rd and 4th. One line
C was omitted from the 5-bus network in order to make a radial system.
C Specifically, the connection (WHITE, BLUE) was removed, leaving just
C 4 lines connecting the 5 busses. I.e., there is no loop. The solution
C should not be much changed since loads of the two connecting nodes were
C increased or decreased by the known amount of the branch flow. Yet, a
C 2nd change must be mentioned: conversion of bus GREEN from a generator
C to a load. This was to make the present case more like NINEBUS from
C Jian Jiang of Clemson University in South Carolina. See EEUG list server
C mail dated October 7th. WSM, 10 October 2001.
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
.000200 .000 60. { T-max = 0 means that no transient solution follows
1 1 1 0 1
WHITE YELLOW .05 .20
YELLOWGREEN .05 .20
GREEN RED .05 .20
RED BLUE .10 .40
BLANK card ending branch cards.
BLANK card ending switch cards.
14WHITE 1.02 60. 0.0 -1.
BLANK card terminating EMTP source cards.
C Next come power constraints of the load flow. There will be one
C for each non-slack generator:
RED -0.4 -.14 .85 1.15
GREEN 1.0 0.3 .85 1.15
BLUE -.15 .025 .85 1.15
YELLOW -0.6 -0.2 .85 1.15
C The following load-flow miscellaneous data card has two peculiarities. The
C use of VSCALE = 1.414 is the special flag requesting RMS rather than peak
C voltages. The use of KTAPER = 0 ensures constant acceleration factors
C (this works for this data).
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 500 20 1 .00001 0.1 7.0 1.414 0
C Max del-V: .0249 .0254 .0259 .0216 .0158 .0092 .0059 .006 .0039 .004 .0059 .0038 .0037 .0031 .002 .0025 .0022 .0028 .0018 .0015
C Source No. 2 2 2 -3 2 -3 4 2 -2 -1 -2 -1 -2 -2 -1 -2 -1 -2 -1 -1
C Max del-V: .0016 .0012 .0013 .0011 .0012 .001 1.E-3 .9E-3 .9E-3 .9E-3 .8E-3 .8E-3 .8E-3 .7E-3 .7E-3 .7E-3 .7E-3 .6E-3 .6E-3 .6E-3
C Source No. -2 -1 -2 -1 -2 -3 -1 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 ... (etc.)
C Max del-V: .2E-5 .2E-5 .2E-5 .2E-5 .2E-5 .1E-5 .1E-5 .1E-5 .1E-5
C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3
C Exit the load flow iteration loop with counter NEKITE = 169. If no warning on the preceding line, convergence was attained.
C Row Node Name Voltage magnit Degrees Real power P Reactive power
C 2 5 RED 9.61553517E-01 -3.15303 -4.00001772E-01 -1.39992436E-01
C 3 4 GREEN 1.02208195E+00 2.93622 9.99997725E-01 3.00004959E-01
C 4 6 BLUE 9.54082493E-01 -7.05955 -1.50001447E-01 2.50097905E-02
C 5 3 YELLOW 9.83589266E-01 -1.67083 -6.00000908E-01 -1.99997137E-01
BLANK card ending requests for output variables
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 6th of 6 data subcases is a variation of the 3rd. The solution is
C identical. But punched cards differ because T-stop of Type-14 source
C GREEN is nonzero. Previously, this data field was left blank. Nonzero
C T-stop is a flag that affects only the punched output. A nonzero value
C in columns 71-80 is taken to be internal source reactance. Instead of an
C equivalent impedance load, ATP will create an equivalent Type-14 voltage
C source behind this reactance. The A6 bus name will have "G" appended
C to it on the right to produce the name of the internal node. Thus GREEN
C yields GREENG as the internal node. Should the original name involve
C all 6 bytes, the last will be overlaid (e.g., PURPLE would be changed
C to PURPLG). This is for every Type-14 source except the one that will
C be used as the slack bus. The slack bus (here, WHITE) is the exception,
C and no output is created for it. This is why T-stop of WHITE has
C been left blank. Adding reactance to WHITE will not change the output.
C Let's summarize the change. Using X = 0.5 will change the following
C two 2 punched cards that follow the ruler:
C 12345678901234567890123456789012345678901234567890123456789012345678901234
C Polar V, P, Q = 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01
C GREEN -.45374803 15276.6334
C (one comment card followed by one wide-format R-C branch card) to:
C 14GREENG 1.84312934 60.24.5286405 -1.0
C GREENGGREEN 0.0 0.5
C (one Type-14 source card followed by one R-L branch card).
C Note that reactance X normally will be positive, and the R-L branch
C card carries this positive value in columns 43-58. This differs from
C the preceding equivalent impedance which might be either R-L or R-C
C depending on the power factor. Finally, note the raised voltage. Most
C of the increase from 1.02 to 1.84 is due to the need for peak rather
C than RMS on the Type-14 source card. As the reactance approaches zero,
C the voltage should approach 1.02 * 1.414 (square root of 2). This 6th
C subcase is being added 8 August 2007 following an inquiry by Prof.
C Mustafa Kizilcay of the University of Siegen in Germany. The use of
C T-stop > 0 was never documented in the ATP Rule Book. Neither was it
C illustrated in standard test case. That is why this subcase now is
C being added. WSM.
C Previous comment about the slack bus being missing is correct for normal
C data. But another undocumented option allows the user to name his slack
C bus manually at the end of power constraints. If this is done, there will
C be output for the slack generator, too. This option depends on the request
C "last is slack bus" which can appear in any case and any location to the
C right of the name field in columns 3-8. WSM. 9 August 2007
C The preceding continues to be honored, but is believed to be clumsy.
C The more modern way to name the slack bus omits one of the 2 associated
C data cards by use of a special, in-line comment. The interpretation is
C improved at the same time. WSM. 15 August 2007
C Final improvement was suggested by Prof. Kizilcay, who wondered why the
C user should be forced to sort the punched output manually. Prof. Kizilcay
C suggested that /SOURCE might precede each Type-14 source card & /BRANCH
C might follow it. Done. WSM. 15 August 2007.
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns
FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections.
$ERASE { Erase punched cards, created by preceding subcases, from punch buffer
.000200 .000 60. { T-max = 0 means that no transient solution follows
1 1 1 0 1
WHITE RED .10 .40
WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X
DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X
WHITE YELLOW .05 .20
RED GREEN .05 .20
RED BLUE .10 .40
GREEN YELLOW .05 .20
BLANK card ending branch cards.
BLANK card ending switch cards.
14WHITE 1.02 60. 0.0 -1. 0.4
14GREEN 1.04 60. 2.0 -1. 0.5
BLANK card terminating all EMTP source cards.
RED -0.6 -0.3 .85 1.15
1GREEN 1.0 1.02 -20. 20.
BLUE -0.4 -0.1 .85 1.15
YELLOW -0.6 -0.2 .85 1.15
WHITE { Slack bus for punched output
C The preceding special in-line comment identifies the name of the slack bus.
C The older way, which continues to be honored, is to follow the slack bus name
C by a 2nd card which carries the special text (any case) "last is slack bus"
C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER
1 500 20 1 .00001 0.0 0.0 1.414 0
$PUNCH { Flush punched cards: R-L branches for equivalent impedance loads
BLANK card ending requests for output variables
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
EOF
|
