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
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
|
BEGIN NEW DATA CASE
C BENCHMARK DCNEW-27
C 1st of 10 subcases that verify correct use of manually-specified initial
C conditions for problems involving various distributed line models. Prior
C to beginning on 22 November 1998, there was no such illustration. Manual
C specification of initial conditions can, however, be found in many test
C cases (e.g., DC-34 that involves Prof. Hian Lauw's U.M.). But none used
C a distributed line. Begin with 1-phase, constant-parameter distributed.
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
.000050 .002
1 1 1 1
C Create two parallel, independent, identical circuits of line plus resistor.
C The (SEND, REC) alternative will involve no phasor solution whereas the
C alternative (GEN, END) will have a phasor solution. Answers agree:
-1SEND REC .306 5.82 .012 100. { Half the length of original
-1GEN END .306 5.82 .012 100. { Half the length of original
C 19 July 2009, we add branch name "RESIST" for use with initial conditions.
C 34567890123 Name RESIST --- Equivalent, alternative request is left-adjusted
REC Name:RESIST 100. { Name branch using right-adjusted " Name:" } 1
END 100. 1
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60. { Excitation of (SEND, REC) is transient only
14GEN 100. 60. { Note phasor presence } -1.
BLANK card ending source data
C (SEND, REC) had no phasor solution, initialize it using values of (GEN, END):
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0 60.
2REC 19.956813332748-34.69967901024
C Distributed: I-k real I-k imag I-m real I-m imag
3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239
C 19 July 2009, convert the following linear branch current from the original
C use of terminal node names (REC, TERRA) to the use of branch name RESIST.
C Just today, such alternative use of branch names is being added to ATP. WSM
C 3REC .19956813332748 { Unnecessary, since purely resistive branch
3Name: RESIST.19956813332748 { Unnecessary, since purely resistive branch
SEND GEN REC END { List of nodes for node-voltage output
C Travel time is less than 17 steps, so document the solution this far:
C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Next 2 output variables are branch currents (flowing from the upper node to the lower node);
C Step Time SEND GEN REC END REC END
C TERRA TERRA
C 0 0.0 100. 100. 19.95681 19.95681 .1995681 .1995681
C 1 .5E-4 99.98224 99.98224 20.60655 20.60655 .2060655 .2060655
C 2 .1E-3 99.92895 99.92895 21.2497 21.2497 .212497 .212497
C 3 .15E-3 99.84016 99.84016 21.88529 21.88529 .2188529 .2188529
C 4 .2E-3 99.71589 99.71589 22.51311 22.51311 .2251311 .2251311
C 5 .25E-3 99.5562 99.5562 23.13293 23.13293 .2313293 .2313293
C 6 .3E-3 99.36113 99.36113 23.74453 23.74453 .2374453 .2374453
C 7 .35E-3 99.13076 99.13076 24.3477 24.3477 .243477 .243477
C 8 .4E-3 98.86517 98.86517 24.94221 24.94221 .2494221 .2494221
C 9 .45E-3 98.56446 98.56446 25.52786 25.52786 .2552786 .2552786
C 10 .5E-3 98.22873 98.22873 26.10444 26.10444 .2610444 .2610444
C 11 .55E-3 97.85809 97.85809 26.67175 26.67175 .2667175 .2667175
C 12 .6E-3 97.45269 97.45269 27.22958 27.22958 .2722958 .2722958
C 13 .65E-3 97.01266 97.01266 27.77774 27.77774 .2777774 .2777774
C 14 .7E-3 96.53816 96.53816 28.31603 28.31603 .2831603 .2831603
C 15 .75E-3 96.02937 96.02937 28.84425 28.84425 .2884425 .2884425
C 16 .8E-3 95.48645 95.48645 29.36223 29.36223 .2936223 .2936223
C 17 .85E-3 94.90961 94.90961 29.86972 29.86972 .2986972 .2986972
BLANK card ending node voltage outputs
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd of 10 subcases is like 1st, but more general. The line is 3-phase, and
C excitation and connection are unbalanced in order that all line modes have
C strong phasor history. The line is half the length originaly used in DC-38
C and nonlinear elements at the far end have been replaced by linear, 200-ohm
C resistors. Phase c of excitation is omitted, with the voltage source is
C replaced by a 4th 200-ohm resistor to ground. All voltages and currents
C are substantial.
PRINTED NUMBER WIDTH, 11, 2, { Reduced precision
POWER FREQUENCY, 60., { Definition here allows omision on initial condition card
.000050 .00100 { TMAX = 1 msec slightly exceeds the travel time of the line
1 -1 1 1
-1SENDA RECA .305515.8187.01210 100. 0 { 100-mile, constant-
-2SENDB RECB .031991.5559.01937 100. 0 { parameter, 3-phase
-3SENDC RECC { transmission line.
RECA 200.
RECB 200.
RECC 200.
SENDC 200. { Replace phase c source by this resistor
BLANK card follows the last branch card
BLANK line terminates the last (here, nonexistent) switch
14SENDA 100. 60. 0.0
14SENDB 100. 60. -120.
BLANK card ending source data
C The following initial-condition cards were constructed from branch flow
C output with sources present for a phasor solution. Full precision of
C the ATP output is being used, note. Purely resistive branches are not
C supposed to require initialization, so these have been omitted.
C 2SENDA 100. 0.0 60.
C Note omission of frequency on right of preceding card. Since equal to
C the power frequency, its appearance on the 1st node voltage initial
C condition card is not required. The power frequency will be used.
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SENDA 100.0 0.0
2RECA 77.881089324669-16.76951438169
2SENDB -50.-86.60254037844
2RECB -84.53015167415-56.11022779303
2SENDC 8.9524884789137-4.001192729873
2RECC -7.4834609042614.3043134435609
C Distributed: I-k real I-k imag I-m real I-m imag
3SENDA RECA .38833508724138-.020723270343 -.3894054466233.08384757190847
3SENDB RECB -.3776059443784-.3319155592773.42265075837076.28055113896513
3SENDC RECC -.0447624423946.02000596364937 .0374173045213-.0215215672178
C Step Time RECC RECB RECA SENDC SENDA SENDB
C 0 0.0 -7.483461 -84.53015 77.881089 8.9524885 100. -50.
C 1 .5E-4 -7.563122 -83.45714 78.183167 9.0261691 99.982235 -48.35879
C 2 .1E-3 -7.640246 -82.35489 78.457621 9.0967927 99.928947 -46.70041
C 3 .15E-3 -7.714656 -81.22339 78.704198 9.1641842 99.840155 -45.02543
C 4 .2E-3 -7.786324 -80.06303 78.922812 9.2283197 99.71589 -43.33445
C 5 .25E-3 -7.855227 -78.87422 79.113386 9.2891764 99.556196 -41.62808
C 6 .3E-3 -7.921338 -77.65739 79.27585 9.3467328 99.361131 -39.90692
C 7 .35E-3 -7.984635 -76.41297 79.410149 9.4009683 99.130763 -38.17157
C 8 .4E-3 -8.045095 -75.14139 79.516233 9.4518636 98.865174 -36.42267
C 9 .45E-3 -8.102697 -73.84312 79.594066 9.4994008 98.56446 -34.66082
C 10 .5E-3 -8.15742 -72.51862 79.643619 9.5435629 98.228725 -32.88666
C 11 .55E-3 -8.209246 -71.16835 79.664875 9.5843341 97.85809 -31.10082
C 12 .6E-3 -8.258233 -69.7927 79.657818 9.6216977 97.452687 -29.30393
C 13 .65E-3 -8.304204 -68.39235 79.622466 9.6556448 97.01266 -27.49662
C 14 .7E-3 -8.347225 -66.9677 79.558825 9.6861614 96.538164 -25.67954
C 15 .75E-3 -8.387279 -65.51925 79.466917 9.7132365 96.029369 -23.85335
C 16 .8E-3 -8.424354 -64.04753 79.346775 9.7368605 95.486454 -22.01867
C 17 .85E-3 -8.458403 -62.55302 79.198474 9.7570132 94.909614 -20.17617
C 18 .9E-3 -8.489362 -61.03619 79.022118 9.7736699 94.299054 -18.32651
C 19 .95E-3 -8.517428 -59.49781 78.817564 9.786901 93.654989 -16.47033
C 20 .1E-2 -8.542467 -57.93828 78.585006 9.7966549 92.977649 -14.6083
C Note about preceding solution. This not only is smooth, it was verified
C by adding the phasor solution and omitting initial conditions. Agreement
C is exact.
1 { Request node voltage output of all nodes
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 3rd of 10 subcases is like 1st, but without the built-in phasor solution,
C and with the constant-parameter line replaced by JMarti frequency-dependent
C line modeling. Line length has changed, unfortunately. Rather than build
C data to derive Marti branch cards, we borrowed those already punched by
C the 2nd subcase of DCNEW-5. Line length is 108.7 km (7 time steps)
C Reduce 2 in col. 31 to 0 so all Marti data card input is seen.
PRINTED NUMBER WIDTH, 10, 2,
.000050 .001
1 1 1 1
C <++++++> Cards punched by support routine on 11-Nov-18 11.00.00 <++++++>
C ***** TRANSPOSED JMARTI line segment ******
C JMARTI SETUP
C BRANCH JDA LMA JDB LMB JDC LMC
C LINE CONSTANTS
C METRIC
C 1.3161 .1151 4 2.355 -6.5 22. 8.0
C 0.5 4.189 4 0.914 4.5 27.3 15.7
C BLANK card ending conductor cards of imbedded "LINE CONSTANTS" data
C 1500. 60. 1 108.7
C 1500. .01 1 108.7 9 10
C BLANK card ending frequency cards of inbedded "LINE CONSTANTS" data
C BLANK card ending "LINE CONSTANTS" cases (back to "JMARTI SETUP")
C 1 { Note use of default conductance G to ground }
C .30 30 0 1 3 0 0
-1SEND REC 0. 0.00 -2 1
25 4.6006674814379022100E+02
-3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02
1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03
2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03
2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04
6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06
5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07
1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07
1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07
3.72990420164290220E+07
2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01
1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00
9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01
5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03
1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04
4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06
6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07
1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07
2.89354778735701405E+07
14 3.8038812947515656600E-04
1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00
9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01
2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04
-3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10
-7.93497908986478425E+10 3.97029750605793992E+10
3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02
1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03
1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04
5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05
1.82589462762646523E+05 1.82772052225409076E+05
0.58951119
0.00000000
REC 100. 1
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60.
BLANK card ending source data
C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678
C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865
C
C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945
C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15
C The following initial condition cards are built from preceding phasor flow.
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0 60.
2REC 14.563049373597-25.89164864838
C Distributed: I-k real I-k imag I-m real I-m imag
3SEND REC .14695409358625-.2529237293961-.145630493736 .25891648648378
C First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Next 1 output variables are branch currents (flowing from the upper node to the lower node);
C Step Time REC SEND REC
C TERRA
C 0 0.0 14.56305 100. .1456305
C 1 .5E-4 15.04761 99.98224 .1504761
C 2 .1E-3 15.52757 99.92895 .1552757
C 3 .15E-3 16.00208 99.84016 .1600208
C 4 .2E-3 16.47093 99.71589 .1647093
C 5 .25E-3 16.93392 99.5562 .1693392
C 6 .3E-3 17.3909 99.36113 .173909
C 7 .35E-3 17.8417 99.13076 .178417
C 8 .4E-3 18.28609 98.86517 .1828609
C 9 .45E-3 18.72385 98.56446 .1872385
C 10 .5E-3 19.155 98.22873 .19155
C 11 .55E-3 19.57944 97.85809 .1957944
C 12 .6E-3 19.99694 97.45269 .1999694
C 13 .65E-3 20.40735 97.01266 .2040735
C 14 .7E-3 20.8105 96.53816 .208105
C 15 .75E-3 21.20625 96.02937 .2120625
C 16 .8E-3 21.59436 95.48645 .2159436
C Showing 16 steps is more than enough. Without manual initial conditions, the
C surge arrives at the receiving end on step 8. Above exceeds 2 travel times.
C Furthermore, above agrees exactly to output that includes a phasor solution.
1 { Request output of all node voltages
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 10 subcases is like 2nd, but with the constant-parameter line
C replaced by JMarti frequency-dependent line modeling. Line length has
C changed, unfortunately. Rather than build data to derive Marti branch
C cards, we borrowed those used by DCNEW-4. Line length is 138 miles.
C Reduce 2 in col. 31 to 1 so some Marti data card input is seen.
PRINTED NUMBER WIDTH, 10, 2, { Reduced precision
.000050 .001 { Note TMAX exceeds travel time of .74 msec by about 40%
1 1 1 1 { Request phasor output to prove there is none
C PUNCHED CARD OUTPUT OF "JMARTI SETUP" WHICH BEGAN AT 14.32.51 09/14/82
C 1.3636 .05215 4 1.602 -20.75 50. 50.
C 1.3636 .05215 4 1.602 -19.25 50. 50.
C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5
C 2.3636 .05215 4 1.602 0.75 77.5 77.5
C 3.3636 .05215 4 1.602 19.25 50. 50.
C 3.3636 .05215 4 1.602 20.75 50. 50.
C 0.5 2.61 4 0.386 -12.9 98.5 98.5
C 0.5 2.61 4 0.386 12.9 98.5 98.5
C
C 100. 5000. 1 138. 1 3-2
C 100. 60.00 1 138. 1 3
C 100. .01 1 138. 1 9 10 3
C
-1SENDA RECA 1. -2 3 1
18 0.47487144790587633736E+03
-0.117967342612603676E+01 -0.316671345021416945E+01 0.255662317832426673E+02
-0.601654631032706213E+02 -0.118825306534453652E+03 0.128169142098635430E+04
0.589610442538186533E+04 0.292891488173549769E+05 0.115665757459954655E+06
0.456394278488149546E+06 0.153681652821077613E+07 0.681848054466802243E+07
0.155477562610202006E+08 0.231789897204889809E+08 0.143827766961221036E+08
0.241524708800130989E+08 0.245049034784953832E+08 0.476808364145524129E+08
0.321645453133055891E+00 0.897616891499021768E+00 0.165632367032913183E+01
0.166669419493737986E+01 0.253073682990780913E+01 0.200014655383098283E+02
0.124257951002701473E+03 0.657345842950598183E+03 0.278623017775457413E+04
0.117145310573639467E+05 0.427410884649070877E+05 0.202553264108715757E+06
0.968789502295299084E+06 0.296908522006270138E+07 0.737968316276072920E+07
0.121009167337626640E+08 0.132891445947112725E+08 0.252097199908134253E+08
14 0.86849922355034181931E-03
0.580543547692016874E-01 0.269744993210885436E+00 0.358646519088149721E+00
0.775351125410302933E+00 0.512138885239557462E+01 0.838265420474258915E+01
0.700229782114550385E+02 0.393084113996385128E+03 0.138208543918962840E+04
0.706160326497018531E+03 0.527246586822120662E+04 0.528343251309079615E+05
0.178307801042203880E+08 -0.178914532139174500E+08
0.228360173022402533E+02 0.102750699821169636E+03 0.139222638455845377E+03
0.269608322018776320E+03 0.389220463867018722E+03 0.427793721751318500E+03
0.132148883486994549E+04 0.277178340638843952E+04 0.449341296469570784E+04
0.596603161966094024E+04 0.114714152641659587E+05 0.268314807408425590E+05
0.204027813729699328E+05 0.204231841543429041E+05
-2SENDB RECB 1. -2 3 1
13 0.28580875773089406522E+03
0.381761604345214016E+04 -0.269529222495194318E+04 0.445018224900628539E+03
0.214658130219351907E+03 0.897325461203551775E+02 0.130480129629523390E+03
0.683037485206990933E+02 0.746626715384045969E+02 0.133938729387751110E+04
0.249954495902361043E+04 0.187764116219525090E+05 0.391031620802855759E+06
0.562668431190421223E+07
0.358668999787749571E+01 0.376679186833120605E+01 0.706439034439901437E+01
0.109721634064781350E+02 0.140506334140286961E+02 0.234292596911152251E+02
0.407686998675834982E+02 0.797598814775078289E+02 0.129106739969348078E+04
0.242138030939135859E+04 0.182362646041542116E+05 0.379814243183306855E+06
0.548808385606378759E+07
13 0.74123751236187229426E-03
0.192189595814863363E+02 0.210970481418782463E+01 0.383406546226866283E+02
0.307999912886148377E+02 0.634013410317139090E+03 -0.253243195247046863E+03
0.210540381301789816E+05 0.755533947306502050E+05 0.229576557459444489E+07
0.607777584135737231E+13 -0.156620677810598129E+13 0.153143851607186420E+13
-0.604300997216751013E+13
0.325736126030671755E+04 0.366694913707514260E+03 0.675727697588602189E+04
0.516024976370228876E+04 0.249456113770459378E+05 0.115067772934717399E+06
0.105205466546493019E+06 0.240930799961836921E+06 0.633276538663731131E+06
0.975820481990592802E+06 0.976796302472583498E+06 0.974176721795814519E+06
0.975150898517610403E+06
-3SENDC RECC 1. -2 3 1
13 0.27248404992352924125E+03
0.928429573805075023E+03 0.205756155237142632E+03 0.538106630554037608E+03
0.105199818376817531E+03 0.821382476778671116E+02 0.133356875032521778E+03
0.654774478488469374E+02 0.612546121562179922E+02 0.905544916173387175E+02
0.973244316091135204E+03 0.267481871266201051E+04 0.138566411938941123E+05
0.646075558955914967E+06
0.297277477955690506E+01 0.413177069940672514E+01 0.739410817322770197E+01
0.103072536754267949E+02 0.135874091964585515E+02 0.234792992956042887E+02
0.384075708508238121E+02 0.692539345057236435E+02 0.997101922255402791E+02
0.941460937860312015E+03 0.259295974388396348E+04 0.134538093446956691E+05
0.628668738677182439E+06
20 0.74871165651714274120E-03
0.249456855794677017E-01 0.516637614002139312E+01 0.435227595525510835E+01
0.696293545843379624E+01 0.113102519337405665E+02 0.107077756312168120E+02
0.136526319661978881E+02 0.103057494642840792E+03 0.647491597081196815E+03
0.221213037293261823E+04 0.164010494681972204E+05 0.517776934139508076E+05
0.393233870054795443E+06 -0.101893253559274247E+05 -0.173413454816459129E+06
0.524218185885897590E+06 0.711044562955687729E+12 -0.559640645251973480E+12
0.582892127161358521E+12 -0.734296849907947113E+12
0.887741958857318814E+01 0.183177145564733675E+04 0.151655674045953921E+04
0.245910936284398730E+04 0.383273362703089043E+04 0.399285585783856465E+04
0.496339281485561617E+04 0.253464455063045348E+05 0.272156609769823021E+05
0.475215596432249440E+05 0.149791320991727040E+06 0.333884409922466803E+06
0.541583000670874084E+06 0.670984155731562481E+06 0.569769858667429973E+06
0.132622344670840449E+07 0.440534234233555838E+07 0.440974768467789434E+07
0.436787151891064062E+07 0.437223939042955148E+07
0.57154434 0.70710678 -0.41761362
0.00000000 0.00000000 0.00000000
0.58879039 0.00000000 0.80696823
0.00000000 0.00000000 0.00000000
0.57154434 -0.70710678 -0.41761362
0.00000000 0.00000000 0.00000000
RECA 200.
RECB 200.
RECC 200.
SENDC 200. { Replace phase c source by this resistor
BLANK card follows the last branch card
BLANK line terminates the last (here, nonexistent) switch
14SENDA 100. 60. 0.0
14SENDB 100. 60. -120.
BLANK card ending source data
C The following phasor solution was obtained by setting T-start = -1.0 for
C each of the two preceding sources. This provides basis for initialization.
C SENDA 100. 100. .37357344098553 .37380400389422 18.678672049276 3.7819577671008
C 0.0 0.0 -.0131269767114 -2.0124852 .65634883557087 0.6563488
C
C RECA 74.582077276654 77.192523685071 -.3729103863833 .38596261842535 -14.89671428218
C -19.90475978176 -14.9430397 .09952379890879 165.0569603 -.451244944E-14
C
C SENDB -50. 100. -.3833667890536 .47733933630699 21.899004304761 -1.819030186807
C -86.60254037844 -120.0000000 -.2843989223551 -143.4302949 9.4902958555069 9.4902959
C
C RECB -87.9019678687 97.402329523617 .43950983934349 .48701164761808 -23.71803449157
C -41.95781025551 -154.4836884 .20978905127755 25.5163116 -.213197515E-14
C
C SENDC 8.9664534487851 11.539874097595 -.0448322672439 .05769937048798 -.3329217354709 -.5808639378985
C -7.264393074381 -39.0135136 .03632196537191 140.9864864 .2136962482E-15 0.0000000
C
C RECC -7.027068223793 9.9587590075801 .03513534111896 .0497937950379 -.2479422024276
C 7.0567126304834 134.8794005 -.0352835631524 -45.1205995 .3134699531E-16
C The following initial-condition cards were constructed from branch flow
C output with sources present for a phasor solution. Full precision of
C the ATP output is being used, note. Purely resistive branches are not
C supposed to require initialization, so these were omitted.
2SENDA 100. 0.0 60.
2RECA 74.582077276654-19.90475978176
2SENDB -50.-86.60254037844
2RECB -87.9019678687-41.95781025551
2SENDC 8.9664534487851-7.264393074381
2RECC -7.0270682237937.0567126304834
C Distributed: I-k real I-k imag I-m real I-m imag
3SENDA RECA .37357344098553-.0131269767114-.3729103863833.09952379890879
3SENDB RECB -.3833667890536-.2843989223551.43950983934349.20978905127755
3SENDC RECC -.0448322672439.03632196537191.03513534111896-.0352835631524
1 { Request node voltage outputs for each network node
C First 6 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Next 3 output variables are branch currents (flowing from the upper node to the lower node);
C Step Time RECC RECB RECA SENDC SENDA SENDB SENDA SENDB SENDC
C RECA RECB RECC
C 0 0.0 -7.02707 -87.902 74.58208 8.966453 100. -50. .3735734 -.383367 -.044832
C 1 .5E-4 -7.15772 -87.0944 74.94189 9.101568 99.98224 -48.3588 .3737522 -.377936 -.045508
C 2 .1E-3 -7.28683 -86.2569 75.27724 9.233558 99.92895 -46.7004 .3738012 -.372372 -.046168
C 3 .15E-3 -7.41337 -85.3888 75.58585 9.362285 99.84016 -45.0254 .3737174 -.366676 -.046811
C 4 .2E-3 -7.5373 -84.4904 75.86758 9.487714 99.71589 -43.3345 .3735005 -.360851 -.047439
C 5 .25E-3 -7.65858 -83.5619 76.12234 9.609793 99.5562 -41.6281 .3731507 -.354897 -.048049
C 6 .3E-3 -7.77715 -82.6038 76.35003 9.728472 99.36113 -39.9069 .3726683 -.348817 -.048642
C 7 .35E-3 -7.89297 -81.6164 76.55059 9.843704 99.13076 -38.1716 .3720534 -.342613 -.049219
C 8 .4E-3 -8.00599 -80.6 76.72395 9.955443 98.86517 -36.4227 .3713062 -.336287 -.049777
C 9 .45E-3 -8.11617 -79.5549 76.87004 10.06365 98.56446 -34.6608 .3704271 -.329843 -.050318
C 10 .5E-3 -8.22347 -78.4816 76.98882 10.16828 98.22873 -32.8867 .3694163 -.32328 -.050841
C 11 .55E-3 -8.32786 -77.3804 77.08024 10.2693 97.85809 -31.1008 .3682743 -.316604 -.051347
C 12 .6E-3 -8.42928 -76.2517 77.14427 10.36668 97.45269 -29.3039 .3670014 -.309814 -.051833
C 13 .65E-3 -8.52771 -75.0959 77.1809 10.46037 97.01266 -27.4966 .3655981 -.302915 -.052302
C 14 .7E-3 -8.62311 -73.9135 77.1901 10.55035 96.53816 -25.6795 .364065 -.295908 -.052752
C 15 .75E-3 -8.71543 -72.7047 77.17185 10.63653 96.02937 -23.8533 .362402 -.288795 -.053183
C 16 .8E-3 -8.80466 -71.47 77.12599 10.71883 95.48645 -22.0187 .3606083 -.281579 -.053594
C 17 .85E-3 -8.89082 -70.21 77.05293 10.79745 94.90961 -20.1762 .3586886 -.274264 -.053987
C 18 .9E-3 -8.97382 -68.925 76.95252 10.87225 94.29905 -18.3265 .3566416 -.266852 -.054361
C 19 .95E-3 -9.05361 -67.6156 76.82477 10.94318 93.65499 -16.4703 .354468 -.259345 -.054716
C 20 .1E-2 -9.13015 -66.2822 76.66974 11.01022 92.97765 -14.6083 .3521686 -.251745 -.055051
BLANK card ending plot cards
BEGIN NEW DATA CASE
C BENCHMARK DC-9
C 5th of 10 subcases is like preceding, but with JMarti frequency
C dependence replaced by Semlyen branch cards from DC-31. This example
C comes from BPA's Dr. Tsu-huei Liu, who first made it work.
PRINTED NUMBER WIDTH, 11, 2, { Reduced precision
.000050 .00100
1 1 1 1 1
C L= 138.0MILES, RHO= 27.0, SS FREQ= 60.00, NSS=0, KFIT= 0, KPS=2, KYC=10
C 1.3636 .05215 4 1.602 -20.75 50. 50.
C 1.3636 .05215 4 1.602 -19.25 50. 50.
C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5
C 2.3636 .05215 4 1.602 0.75 77.5 77.5
C 3.3636 .05215 4 1.602 19.25 50. 50.
C 3.3636 .05215 4 1.602 20.75 50. 50.
C 0.5 2.61 4 0.386 -12.9 98.5 98.5
C 0.5 2.61 4 0.386 12.9 98.5 98.5
C
C 27. 5000. 1 138.
C 27. 60.00 1 138.
C 27. 6.00 1 138. 6 20
C
-1SENDA RECA 0.55456E-02 0.77998E-03 1 1 2 2 3
0.14938602E+02 0.90516313E+02-0.16261865E-03 0.14500324E-02 0.60000000E+02
0.00000E+00 0.15161E+05 0.75119E+00 0.00000E+00 0.17105E+04 0.24881E+00
0.00000E+00 0.59584E+03-0.11954E-02 0.00000E+00 0.39933E+05-0.74162E-03
-1SENDB RECB 0.74392E-02 0.74149E-03 2 2 2 2 3
0.83801231E+00 0.38634735E+02-0.44536501E-04 0.20637999E-02 0.60000000E+02
0.00000E+00 0.26608E+06 0.83767E+00 0.00000E+00 0.35627E+05 0.16233E+00
0.00000E+00 0.13630E+03-0.49171E-03 0.00000E+00 0.32940E+05-0.75349E-04
-1SENDC RECC 0.42823E-02 0.74017E-03 3 3 2 2 3
0.26535168E+01 0.67894100E+02-0.45225808E-04 0.11682248E-02 0.60000000E+02
0.00000E+00 0.61698E+06 0.90903E+00 0.00000E+00 0.15239E+05 0.90969E-01
0.00000E+00 0.25667E+03-0.34465E-03 0.00000E+00 0.11831E+05-0.28233E-04
0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00
0.71114E+00 0.00000E+00-0.41983E-15 0.00000E+00 0.10000E+01 0.00000E+00
0.10000E+01 0.00000E+00-0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00
0.42027E+00 0.00000E+00 0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00
0.22471E+00 0.00000E+00-0.27373E-15 0.00000E+00 0.84054E+00 0.00000E+00
0.42027E+00 0.00000E+00-0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00
BLANK card follows the last branch card
BLANK line terminates the last (here, nonexistent) switch
14SENDA 408000. 60. 0.0
14SENDB 408000. 60. -120.
14SENDC 408000. 60. 120.
BLANK card ending source data
C The following initial-condition cards were constructed from branch flow
C output with sources present for a phasor solution. Full precision of
C the ATP output is being used, note. Purely resistive branches are not
C supposed to require initialization, so these were omitted.
C
C SENDA 408000. 408000. -14.35539062039 427.20215266352 -.29284996866E7 -.29284996866E7
C 0.0 0.0 426.96089048118 91.9256883 -.87100021658E8 87100021.6581596
C
C RECA 426275.19929468 426279.30459602 .1509903313E-13 .9664560746E-12 .23142460224E-8
C 1870.8268620134 0.2514568 -.966338121E-12 -89.1048263 .20597711132E-6
C
C SENDB -204000. 408000. 390.21858493854 442.08847622165 -.30943350074E7 -.30943350074E7
C -353338.3647441 -120.0000000 -207.7779506506 -28.0337631 -.90132949314E8 90132949.3138173
C
C RECB -211353.2339786 423078.09601739 .6536993169E-12 .7963152328E-12 .14252537328E-7
C -366503.5959121 -119.9709473 -.454747351E-12 -34.8244892 -.1678477368E-6
C
C SENDC -204000. 408000. -374.9796048273 439.89286922318 -.2384468398E7 -.2384468398E7
C 353338.36474405 120.0000000 -229.9913745272 -148.4773916 -.89706460393E8 89706460.3928106
C
C RECC -210971.8157184 425738.55951138 -.170530257E-12 .190658623E-12 .33753610037E-7
C 369789.41849017 119.7055709 .8526512829E-13 153.4349488 -.2253587274E-7
2SENDA 408000. 0.0 60.
2RECA 426275.199295 1870.82686201
2SENDB -204000. -353338.36474
2RECB -211353.23398 -366503.59591
2SENDC -204000. 353338.36474
2RECC -210971.81572 369789.41849
3SENDA RECA -14.35539062039426.96089048118.1509903313E-13-.966338121E-12
3SENDB RECB 390.21858493854-207.7779506506.6536993169E-12-.454747351E-12
3SENDC RECC -374.9796048273-229.9913745272-.170530257E-12.8526512829E-13
1 { Request for the output of each node voltage
C First 6 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Step Time RECC RECB RECA SENDA SENDB SENDC
C 0 0.0 -210971.8 -211353.2 426275.2 408000. -204000. -204000.
C 1 .5E-4 -217903.1 -204391.8 426140.86 407927.52 -197303.9 -210623.6
C 2 .1E-3 -224764. -197371.2 425871.55 407710.1 -190537.7 -217172.4
C 3 .15E-3 -231544. -190283.3 425452.69 407347.83 -183703.7 -223644.1
C 4 .2E-3 -238240.9 -183128.9 424883.53 406840.83 -176804.6 -230036.3
C 5 .25E-3 -244852.7 -175910. 424163.92 406189.28 -169842.6 -236346.7
C 6 .3E-3 -251377.2 -168628.7 423293.95 405393.41 -162820.2 -242573.2
C 7 .35E-3 -257812.1 -161287.4 422273.86 404453.51 -155740. -248713.5
C 8 .4E-3 -264155.2 -153888.8 421103.95 403369.91 -148604.5 -254765.4
C 9 .45E-3 -270404.2 -146435.4 419784.61 402142.99 -141416.2 -260726.8
C 10 .5E-3 -276557. -138929.9 418316.3 400773.2 -134177.6 -266595.6
C 11 .55E-3 -282611.4 -131374.9 416699.51 399261.01 -126891.3 -272369.7
C 12 .6E-3 -288565.2 -123773.2 414934.81 397606.96 -119560. -278046.9
C 13 .65E-3 -294416.4 -116127.4 413022.82 395811.65 -112186.2 -283625.4
C 14 .7E-3 -300162.8 -108440.2 410964.19 393875.71 -104772.5 -289103.2
C 15 .75E-3 -305803.1 -100717.9 408762.02 391799.82 -97321.65 -294478.2
C 16 .8E-3 -311336.8 -92970.64 406424.38 389584.73 -89836.18 -299748.6
C 17 .85E-3 -316755.1 -85178.69 403935.32 387231.23 -82318.79 -304912.4
C 18 .9E-3 -322059.1 -77353.3 401302.31 384740.14 -74772.16 -309968.
C 19 .95E-3 -327248.4 -69499.21 398526.63 382112.35 -67198.95 -314913.4
C 20 .1E-2 -332321.6 -61619.94 395609.19 379348.81 -59601.88 -319746.9
BLANK card ending plot cards
BEGIN NEW DATA CASE
PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns)
.000050 .001
1 1 1 1
C 6th of 10 subcases illustrates multiple frequencies for initial condit.
C Create 3 parallel, independent, identical circuits of line plus resistor.
C The (SEND, REC) alternative will be initialized for 60 Hz whereas the
C alternative (GEN, END) and (GEN, MARTI) will use 50 Hz --- with the
C former being constant-parameter and the latter JMARTI. It works. There is
C perfect agreement with related solution that begins with phasors. The
C Marti line is identical to what already was used (see 3rd subcase), so
C documentation will be minimized. Comment cards have been dropped, and the
C 1 of col. 31 is restored to 2.
-1SEND REC .306 5.82 .012 100. { Half the length of original
-1GEN END .306 5.82 .012 100. { Half the length of original
-1GEN MARTI 2. 0.00 -2 1
25 4.6006674814379022100E+02
-3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02
1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03
2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03
2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04
6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06
5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07
1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07
1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07
3.72990420164290220E+07
2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01
1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00
9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01
5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03
1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04
4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06
6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07
1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07
2.89354778735701405E+07
14 3.8038812947515656600E-04
1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00
9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01
2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04
-3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10
-7.93497908986478425E+10 3.97029750605793992E+10
3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02
1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03
1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04
5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05
1.82589462762646523E+05 1.82772052225409076E+05
0.58951119
0.00000000
REC 100.
END 100.
MARTI 100.
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60.
14GEN 100. 50.
C 14GEN2 100. 50. -1.
BLANK card ending source data
C GEN 100. 100. .26567063730001 .43724454992683 13.283531865001 3.0530027479459
C 0.0 0.0 -.3472778554952 -52.5836910 17.363892774758 17.3638928
C
C END 25.859996851553 45.233901262338 -.2585999685155 .45233901262338 -10.23052911705
C -37.11288706432 -55.1314894 .37112887064316 124.8685106 -.737257477E-17
C
C GEN 100. 100. .18468803864924 .33562381569663 9.2344019324622 3.4793989521433
C 0.0 0.0 -.2802386019853 -56.6135784 14.011930099264 14.0119301
C
C MARTI 18.346330639813 33.926399692036 -.1834633063981 .33926399692036 -5.755002980319
C -28.53791772569 -57.2640259 .28537917725693 122.7359741 .7014788056E-15
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0
2REC 19.956813332748-34.69967901024
2GEN 100. 0.0
2END 25.859996851553-37.11288706432
2MARTI 18.346330639813-28.53791772569
C Distributed: I-k real I-k imag I-m real I-m imag [Hz]
3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239 60.
3GEN END .26567063730001-.3472778554952-.2585999685155.37112887064316 50.
3GEN MARTI .18468803864924-.2802386019853-.1834633063981.28537917725693
C Note that preceding branch card has blank cols. 75-80 for the frequency. By
C definition, frequency will remain unchanged at 50 Hz (2nd Type-3 card).
SEND GEN REC END MARTI { Names of nodes for voltage output
BLANK card ending node voltage outputs
C First 5 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Step Time SEND GEN REC END MARTI
C 0 0.0 100. 100. 19.956813 25.859997 18.346331
C 1 .5E-4 99.982235 99.987663 20.606552 26.439082 18.791563
C 2 .1E-3 99.928947 99.950656 21.249697 27.012298 19.232816
C 3 .15E-3 99.840155 99.888987 21.885291 27.578849 19.669386
C 4 .2E-3 99.71589 99.802673 22.51311 28.138595 20.101117
C 5 .25E-3 99.556196 99.691733 23.13293 28.691398 20.527892
C 6 .3E-3 99.361131 99.556196 23.744531 29.237122 20.949603
C 7 .35E-3 99.130763 99.396096 24.347696 29.775632 21.366146
C 8 .4E-3 98.865174 99.21147 24.94221 30.306796 21.777336
C 9 .45E-3 98.56446 99.002366 25.527863 30.830482 22.183024
C 10 .5E-3 98.228725 98.768834 26.104445 31.346561 22.583287
C 11 .55E-3 97.85809 98.510933 26.671752 31.854905 22.97808
C 12 .6E-3 97.452687 98.228725 27.229584 32.35539 23.367227
C 13 .65E-3 97.01266 97.922281 27.77774 32.847892 23.750616
C 14 .7E-3 96.538164 97.591676 28.316028 33.332289 24.128147
C 15 .75E-3 96.029369 97.236992 28.844255 33.808461 24.499706
C 16 .8E-3 95.486454 96.858316 29.362233 34.276292 24.865139
C 17 .85E-3 94.909614 96.455742 29.869718 34.735588 25.224341
C 18 .9E-3 94.299054 96.029369 30.366492 35.186194 25.577322
C 19 .95E-3 93.654989 95.579301 30.852618 35.628304 25.924072
C 20 .1E-2 92.977649 95.105652 31.327783 36.061624 26.264482
BLANK card ending plot cards
BEGIN NEW DATA CASE
C BENCHMARK DC-9
C 7th of 10 subcases extends the preceding. It has two 3-phase Semlyen
C lines with different frequencies (50 and 60 Hz). Semlyen branch cards
C are copied from the 5th subcase (from DC-31), so omit the comments.
PRINTED NUMBER WIDTH, 11, 2, { Allow output line nicely within 80 columns
.000050 .001 { Just 20 steps. Extend to 40 msec if plot card is enabled
1 -1 1 1 { Note IPLOT of cols. 9-16 disables plotting
-1SENDA RECA 0.55456E-02 0.77998E-03 1 1 2 2 3
0.14938602E+02 0.90516313E+02-0.16261865E-03 0.14500324E-02 0.60000000E+02
0.00000E+00 0.15161E+05 0.75119E+00 0.00000E+00 0.17105E+04 0.24881E+00
0.00000E+00 0.59584E+03-0.11954E-02 0.00000E+00 0.39933E+05-0.74162E-03
-1SENDB RECB 0.74392E-02 0.74149E-03 2 2 2 2 3
0.83801231E+00 0.38634735E+02-0.44536501E-04 0.20637999E-02 0.60000000E+02
0.00000E+00 0.26608E+06 0.83767E+00 0.00000E+00 0.35627E+05 0.16233E+00
0.00000E+00 0.13630E+03-0.49171E-03 0.00000E+00 0.32940E+05-0.75349E-04
-1SENDC RECC 0.42823E-02 0.74017E-03 3 3 2 2 3
0.26535168E+01 0.67894100E+02-0.45225808E-04 0.11682248E-02 0.60000000E+02
0.00000E+00 0.61698E+06 0.90903E+00 0.00000E+00 0.15239E+05 0.90969E-01
0.00000E+00 0.25667E+03-0.34465E-03 0.00000E+00 0.11831E+05-0.28233E-04
0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00
0.71114E+00 0.00000E+00-0.41983E-15 0.00000E+00 0.10000E+01 0.00000E+00
0.10000E+01 0.00000E+00-0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00
0.42027E+00 0.00000E+00 0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00
0.22471E+00 0.00000E+00-0.27373E-15 0.00000E+00 0.84054E+00 0.00000E+00
0.42027E+00 0.00000E+00-0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00
RECA 200.
RECB 200.
RECC 200.
SENDC 200. { Replace phase c source by this resistor
C Illustrate reference branch copying for Semlyen model. 2nd circuit follows.
C But because the 2nd will be excited at 50 Hz whereas the 1st has its phasor
C parameters derived at 60 Hz, there will be warning messages about 2nd. Note
C that the phasor solution of the second is substantially in error (50 Hz is
C not close to 60 Hz). Yet, the result is amazingly smooth. To see this,
C extend TMAX to a couple of cycles (40 msec), and enable plotting (now,
C IPLOT = -1 will ignore the plot card, which involves RECA, RECB, RECC).
-1GENA ENDA SENDA RECA
-1GENB ENDB SENDB RECB
-1GENC ENDC SENDC RECC
ENDA 200.
ENDB 200.
ENDC 200.
GENC 200. { Replace phase c source by this resistor
BLANK card follows the last branch card
BLANK line terminates the last (here, nonexistent) switch
14SENDA 100. 60. 0.0
14SENDB 100. 60. -120.
14GENA 100. 50. 0.0
14GENB 100. 50. -120.
BLANK card ending source data
C The following initial-condition cards were constructed from branch flow
C output with sources present for a phasor solution. Full precision of
C the ATP output is being used, note. Purely resistive branches are not
C supposed to require initialization, so these were omitted.
C The following 3 branch flows are for 60 Hz:
C SENDA 100. 100. .36998975926495 .37180820218231 18.499487963247 2.9848546654506
C 0.0 0.0 -.0367276088128 -5.6689829 1.8363804406416 1.8363804
C
C RECA 74.690024760283 78.777238585258 -.3734501238014 .39388619292629 -15.5146332978
C -25.04502985478 -18.5373124 .12522514927392 161.4626876 -.253833413E-14
C
C SENDB -50. 100. -.3685220242714 .49798629886772 23.71618850994 -1.806979904073
C -86.60254037844 -120.0000000 -.3349356228992 -137.7334962 7.5840811711722 7.5840812
C
C RECB -86.37515570642 101.04091926346 .43187577853211 .5052045963173 -25.52316841401
C -52.4270907289 -148.7435235 .26213545364451 31.2564765 -.116668827E-13
C
C SENDC 9.9402547548508 10.505728992727 -.0497012737743 .05252864496364 -.2759258541716 -.4312032208709
C -3.40024662007 -18.8841944 .01700123310035 161.1158056 -.204060401E-15 0.0000000
C
C RECC -7.010739885979 7.8810498462914 .0350536994299 .03940524923146 -.1552773666993
C 3.6000656842441 152.8191734 -.0180003284212 -27.1808266 -.154803741E-15
C The following 3 branch flows are for 50 Hz:
C GENA 100. 100. .36998975926495 .37180820218231 18.499487963247 2.9848546654506
C 0.0 0.0 -.0367276088128 -5.6689829 1.8363804406416 1.8363804
C
C ENDA 74.690024760283 78.777238585258 -.3734501238014 .39388619292629 -15.5146332978
C -25.04502985478 -18.5373124 .12522514927392 161.4626876 -.865843855E-14
C
C GENB -50. 100. -.3685220242714 .49798629886772 23.71618850994 -1.806979904073
C -86.60254037844 -120.0000000 -.3349356228992 -137.7334962 7.5840811711722 7.5840812
C
C ENDB -86.37515570642 101.04091926346 .43187577853211 .5052045963173 -25.52316841401
C -52.4270907289 -148.7435235 .26213545364451 31.2564765 .1027303242E-13
C
C GENC 9.9402547548508 10.505728992728 -.0497012737743 .05252864496364 -.2759258541716 -.4312032208709
C -3.40024662007 -18.8841944 .01700123310035 161.1158056 -.319839641E-16 0.0000000
C
C ENDC -7.010739885979 7.8810498462914 .0350536994299 .03940524923146 -.1552773666993
C 3.6000656842441 152.8191734 -.0180003284212 -27.1808266 -.56249764E-16
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SENDA 100.0 0.0
2RECA 74.690024760283-25.04502985478
2SENDB -50.-86.60254037844
2RECB -86.37515570642 -52.4270907289
2SENDC 9.9402547548508 -3.40024662007
2RECC -7.0107398859793.6000656842441
2GENA 100.0 0.0
2ENDA 74.690024760283-25.04502985478
2GENB -50.-86.60254037844
2ENDB -86.37515570642 -52.4270907289
2GENC 9.9402547548508 -3.40024662007
2ENDC -7.0107398859793.6000656842441
C Distributed: I-k real I-k imag I-m real I-m imag [Hz]
3SENDA RECA .36998975926495-.0367276088128-.3734501238014.12522514927392 60.
3SENDB RECB -.3685220242714-.3349356228992.43187577853211.26213545364451
3SENDC RECC -.0497012737743.01700123310035 .0350536994299-.0180003284212
3GENA ENDA .36998975926495-.0367276088128-.3734501238014.12522514927392 50.
3GENB ENDB -.3685220242714-.3349356228992.43187577853211.26213545364451
3GENC ENDC -.0497012737743.01700123310035 .0350536994299-.0180003284212
C About phasor frequency on the right, note that only the first of each group
C of 3 coupled conductors requires a nonzero frequency. A blank for following
C phases means that the declared value will remain in effect. As for the 1st
C frequency (60 Hz), this could be omitted provided it is equal to the power
C frequency. But different places in the world have different power frequency
C so it is prudent practice to declare this frequency explicitely.
RECA ENDA RECB ENDB RECC ENDC { Names of nodes for voltage output
BLANK card ending node names for voltage output
C Step Time RECA ENDA RECB ENDB RECC ENDC
C 0 0.0 74.690025 74.690025 -86.37516 -86.37516 -7.01074 -7.01074
C 1 .5E-4 75.163363 75.132945 -85.37169 -84.34315 -7.059456 -7.969316
C 2 .1E-3 75.595736 75.435454 -84.33743 -82.81562 -7.123114 -8.190566
C 3 .15E-3 76.001094 75.753382 -83.27324 -81.54915 -7.18451 -8.234526
C 4 .2E-3 76.379406 76.079746 -82.17945 -80.40443 -7.243416 -8.226221
C 5 .25E-3 76.730556 76.398672 -81.05645 -79.31142 -7.299777 -8.200803
C 6 .3E-3 77.054425 76.700684 -79.90465 -78.23596 -7.353565 -8.170835
C 7 .35E-3 77.350898 76.981153 -78.72446 -77.16148 -7.404759 -8.141389
C 8 .4E-3 77.619873 77.238057 -77.51631 -76.07989 -7.453337 -8.114506
C 9 .45E-3 77.861255 77.470657 -76.28063 -74.98724 -7.499282 -8.090883
C 10 .5E-3 78.07496 77.678819 -75.01785 -73.88158 -7.542576 -8.070602
C 11 .55E-3 78.260914 77.862692 -73.72842 -72.76194 -7.583202 -8.05347
C 12 .6E-3 78.419051 78.02254 -72.41281 -71.62787 -7.621145 -8.039181
C 13 .65E-3 78.549315 78.158673 -71.07147 -70.47916 -7.65639 -8.027389
C 14 .7E-3 78.651663 78.271409 -69.70489 -69.31574 -7.688924 -8.017745
C 15 .75E-3 78.726915 78.367856 -68.31324 -68.14096 -7.71979 -8.014505
C 16 .8E-3 78.777534 78.482038 -66.89676 -66.94954 -7.75238 -8.021203
C 17 .85E-3 78.796385 78.578276 -65.45863 -65.69399 -7.778268 -7.991112
C 18 .9E-3 78.78656 78.66135 -63.99679 -64.41226 -7.800598 -7.945947
C 19 .95E-3 78.749002 78.71769 -62.51169 -63.1208 -7.819641 -7.907277
C 20 .1E-2 78.683581 78.746939 -61.00415 -61.8202 -7.835692 -7.875092
C ----------------- Preceding was before Tsu-huei's OVER8 correction.
C Step Time RECA ENDA RECB ENDB RECC ENDC
C 0 0.0 74.690025 74.690025 -86.37516 -86.37516 -7.01074 -7.01074
C 1 .5E-4 75.163363 74.744884 -85.37169 -84.65358 -7.059456 -8.300563
C 2 .1E-3 75.595736 75.003289 -84.33743 -83.20507 -7.123114 -8.557454
C 3 .15E-3 76.001094 75.325776 -83.27324 -81.97004 -7.18451 -8.597739
C 4 .2E-3 76.379406 75.66166 -82.17945 -80.84254 -7.243416 -8.582616
C 5 .25E-3 76.730556 75.989351 -81.05645 -79.76068 -7.299777 -8.551501
C 6 .3E-3 77.054425 76.299083 -79.90465 -78.69329 -7.353565 -8.516962
C 7 .35E-3 77.350898 76.58653 -78.72446 -77.62521 -7.404759 -8.483715
C 8 .4E-3 77.619873 76.849919 -77.51631 -76.54911 -7.453337 -8.453539
C 9 .45E-3 77.861255 77.088671 -76.28063 -75.46146 -7.499282 -8.42697
C 10 .5E-3 78.07496 77.302754 -75.01785 -74.36051 -7.542576 -8.403987
C 11 .55E-3 78.260914 77.492377 -73.72842 -73.24544 -7.583202 -8.384334
C 12 .6E-3 78.419051 77.657849 -72.41281 -72.11583 -7.621145 -8.367661
C 13 .65E-3 78.549315 77.799509 -71.07147 -70.97152 -7.65639 -8.353591
C 14 .7E-3 78.651663 77.917697 -69.70489 -69.81243 -7.688924 -8.341754
C 15 .75E-3 78.726915 77.996637 -68.31324 -68.63631 -7.71979 -8.317196
C 16 .8E-3 78.777534 78.003092 -66.89676 -67.41624 -7.75238 -8.217882
C 17 .85E-3 78.796385 78.078973 -65.45863 -66.12507 -7.778268 -8.151926
C 18 .9E-3 78.78656 78.16892 -63.99679 -64.82237 -7.800598 -8.093319
C 19 .95E-3 78.749002 78.231147 -62.51169 -63.5241 -7.819641 -8.050547
C 20 .1E-2 78.683581 78.263977 -61.00415 -62.22435 -7.835692 -8.017481
144 4. 0.0 40. ENDA ENDB ENDC
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 8th of 10 subcases demonstrates manual initialization of a JMarti line.
C Use two different frequencies. Single phase is used for simplicity.
C One frequency will be very low (.01 Hz) to prove no need for true dc.
C This subcase will use phasor solution for .01 Hz. The following subcase
C will replace the phasor solution here by initial condition cards.
PRINTED NUMBER WIDTH, 13, 2,
.000050 .001
1 1 1 1
-1SEND REC 1. 0.00 -2 1
25 4.6006674814379022100E+02
-3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02
1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03
2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03
2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04
6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06
5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07
1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07
1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07
3.72990420164290220E+07
2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01
1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00
9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01
5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03
1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04
4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06
6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07
1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07
2.89354778735701405E+07
14 3.8038812947515656600E-04
1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00
9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01
2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04
-3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10
-7.93497908986478425E+10 3.97029750605793992E+10
3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02
1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03
1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04
5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05
1.82589462762646523E+05 1.82772052225409076E+05
0.58951119
0.00000000
REC 100.
-1GEN END SEND REC { 2nd JMarti line is disconnected from the first
END 100.
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60. { No phasor presence; I.C. cards will be used
14GEN 100. .010 { Low-freq phasor solution } -1.
BLANK card ending source data
C The following initial-condition cards were constructed from branch flow
C output with the source present for a phasor solution. Full precision of
C the ATP output is being used, note.
C The following branch flow is for 60 Hz:
C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678
C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865
C
C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945
C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15
C The following initial condition cards are built from preceding phasor flow.
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0
2REC 14.563049373597-25.89164864838
C Distributed: I-k real I-k imag I-m real I-m imag F [Hz]
3SEND REC .14695409358625-.2529237293961-.145630493736 .25891648648378 60.
C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C Step Time SEND REC GEN END
C 0 0.0 100. 14.56304937 100. 74.15408515
C 1 .5E-4 99.98223524 15.04761164 100. 74.15408201
C 2 .1E-3 99.92894726 15.52757104 100. 74.15407898
C 3 .15E-3 99.84015501 16.00208444 100. 74.15407596
C 4 .2E-3 99.71589003 16.47092745 99.99999999 74.15407293
C 5 .25E-3 99.55619646 16.93392289 99.99999999 74.1540699
C 6 .3E-3 99.36113105 17.39090307 99.99999998 74.15406686
C 7 .35E-3 99.13076311 17.84170463 99.99999998 74.15406383
C 8 .4E-3 98.86517447 18.28608879 99.99999997 74.15406078
C 9 .45E-3 98.56445951 18.72384749 99.99999996 74.1540577
C 10 .5E-3 98.22872507 19.15499653 99.99999995 74.15405464
C 11 .55E-3 97.85809043 19.57943512 99.99999994 74.1540516
C 12 .6E-3 97.45268728 19.99693921 99.99999993 74.15404857
C 13 .65E-3 97.01265965 20.40734552 99.99999992 74.15404553
C 14 .7E-3 96.53816388 20.81050376 99.9999999 74.15404249
C 15 .75E-3 96.02936857 21.20624627 99.99999989 74.15403945
C 16 .8E-3 95.48645447 21.59436163 99.99999987 74.1540364
C Showing 16 steps is more than enough. Without manual initial conditions, the
C surge arrives at the receiving end on step 8. Above exceeds 2 travel times.
C Furthermore, above agrees exactly to output that includes a phasor solution.
SEND REC GEN END { Node names for node voltage output
BLANK card ending output requests
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 9th of 10 subcases has the same two JMarti lines as the preceding.
C It replaces the low-frequency phasor solution of the preceding by
C by initial condition cards. This is manual, multi-frequency JMarti.
PRINTED NUMBER WIDTH, 13, 2,
.000050 .001
1 1 1 1
-1SEND REC 1. 0.00 -2 1
25 4.6006674814379022100E+02
-3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02
1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03
2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03
2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04
6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06
5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07
1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07
1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07
3.72990420164290220E+07
2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01
1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00
9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01
5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03
1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04
4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06
6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07
1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07
2.89354778735701405E+07
14 3.8038812947515656600E-04
1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00
9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01
2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04
-3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10
-7.93497908986478425E+10 3.97029750605793992E+10
3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02
1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03
1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04
5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05
1.82589462762646523E+05 1.82772052225409076E+05
0.58951119
0.00000000
REC 100.
-1GEN END SEND REC
C GEN END 100.
END 100.
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60.
C 14GEN 100. 60. -1.
14GEN 100. .010
BLANK card ending source data
C The following initial-condition cards were constructed from branch flow
C output with the source present for a phasor solution. Full precision of
C the ATP output is being used, note.
C The following branch flow is for 60 Hz:
C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678
C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865
C
C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945
C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15
C The following branch flow is for .01 Hz:
C GEN 100. 100. .74154511249386 .74160775929634 37.077255624693 9.5784693158927
C 0.0 0.0 .00963923155659 0.7447371 -.4819615778294 -0.4819616
C
C END 74.154085146507 74.16034831202 -.7415408514651 .7416034831202 -27.4987863088
C .96380384134885 0.7446492 -.0096380384135 -179.2553508 .5416403003E-15
C
C The following initial condition cards are built from preceding phasor flow.
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0
2GEN 100. 0.0
2REC 14.563049373597-25.89164864838
2END 74.154085146507.96380384134885
C Distributed: I-k real I-k imag I-m real I-m imag F [Hz]
3SEND REC .14695409358625-.2529237293961 -.145630493736.25891648648378 60.
3GEN END .74154511249386.00963923155659-.7415408514651-.0096380384135 .01
C Showing 16 steps is more than enough. Without manual initial conditions, the
C surge arrives at the receiving end on step 8. Above exceeds 2 travel times.
C Furthermore, above agrees exactly to output that includes a phasor solution.
SEND REC GEN END { Node names for node voltage output
C Step Time SEND REC GEN END
C 0 0.0 100. 14.56304937 100. 74.15408515
C 1 .5E-4 99.98223524 15.04761164 100. 74.15408201
C 2 .1E-3 99.92894726 15.52757104 100. 74.15407898
C 3 .15E-3 99.84015501 16.00208444 100. 74.15407596
C 4 .2E-3 99.71589003 16.47092745 99.99999999 74.15407293
C 5 .25E-3 99.55619646 16.93392289 99.99999999 74.1540699
C 6 .3E-3 99.36113105 17.39090307 99.99999998 74.15406686
C 7 .35E-3 99.13076311 17.84170463 99.99999998 74.15406383
C 8 .4E-3 98.86517447 18.28608879 99.99999997 74.15406078
C 9 .45E-3 98.56445951 18.72384749 99.99999996 74.1540577
C 10 .5E-3 98.22872507 19.15499653 99.99999995 74.15405464
C 11 .55E-3 97.85809043 19.57943512 99.99999994 74.1540516
C 12 .6E-3 97.45268728 19.99693921 99.99999993 74.15404857
C 13 .65E-3 97.01265965 20.40734552 99.99999992 74.15404553
C 14 .7E-3 96.53816388 20.81050376 99.9999999 74.15404249
C 15 .75E-3 96.02936857 21.20624627 99.99999989 74.15403945
C 16 .8E-3 95.48645447 21.59436163 99.99999987 74.1540364
BLANK card ending output requests
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 10th of 10 subcases is like the 1st except that a repeated initial
C condition has been added to illustrate the new trap for these. Use
C 1st subcase data prior to 19 July 2009 (see comments in 1st subcase).
C WSM appends this subcase 1 August 2009.
.000050 .002
1 1 0 0 { No need to see connectivity or phasor output
-1SEND REC .306 5.82 .012 100. { Half the length of original
-1GEN END .306 5.82 .012 100. { Half the length of original
REC 100. 1
END 100. 1
BLANK card terminating branch data
BLANK card terminating all (in this case, nonexistent) switches
14SEND 100. 60. { Excitation of (SEND, REC) is transient only
14GEN 100. 60. { Note phasor presence } -1.
BLANK card ending source data
C (SEND, REC) had no phasor solution, initialize it using values of (GEN, END):
C Name K V-k real V-k imag ----- Node V initial condition is type 2
2SEND 100. 0.0 60.
2REC 19.956813332748-34.69967901024
C Distributed: I-k real I-k imag I-m real I-m imag
3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239
3REC .19956813332748 { Unnecessary, since purely resistive branch
3REC .19956813332748 { This is the duplicate initial condition
SEND GEN REC END { List of nodes for node-voltage output
BLANK card ending node voltage outputs
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
|
