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
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
1536
1537
1538
1539
1540
1541
1542
1543
1544
1545
1546
1547
1548
1549
1550
1551
1552
1553
1554
1555
1556
1557
1558
1559
1560
1561
1562
1563
1564
1565
1566
1567
1568
1569
1570
1571
1572
1573
1574
1575
1576
1577
1578
1579
1580
1581
1582
1583
1584
1585
1586
1587
1588
1589
1590
1591
1592
1593
1594
1595
1596
1597
1598
1599
1600
1601
1602
1603
1604
1605
1606
1607
1608
1609
1610
1611
1612
1613
1614
1615
1616
1617
1618
1619
1620
1621
1622
1623
1624
1625
1626
1627
1628
1629
1630
1631
1632
1633
1634
1635
1636
1637
1638
1639
1640
1641
1642
1643
1644
1645
1646
1647
1648
1649
1650
1651
1652
1653
1654
1655
1656
1657
1658
1659
1660
1661
1662
1663
1664
1665
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
1690
1691
1692
1693
1694
|
BEGIN NEW DATA CASE
C 1st of 18 subcases illustrates Robert Meredith's Type-68 TACS device as
C first described in the April, 1998, newsletter. This is AUT5 data case.
C test120a.dat (1 HERTZ AT 120% VOLTAGE)
C without resistances
C ENERGIZED AT ZERO TIME VIA 0.100 OHM RESISTORS -steady state
C FREE FORMAT TIME STEP- 1000/CYCLE FOR 1.5 SECOND+
C .000016666666, 1.500 , 60.000 , 0.,,,,,,, { Meredith's original params
.000066666667, .100 , 60.000 , 0.,,,,,,, { Much reduced burden on computer
C PRN> PLT> NET> SS> MAX> SAVE> NENERG> >
1 5 0 0 1 -1 0
5 5 20 20 100 100 500 500
C $NEW EPSILN, 1.E-10, { Add 24 July 2008 to compensate for protection in REDUCT
$NEW EPSILN, 1.E-11, { Add 31 July 2008 to allow solution with MATFUL = 1, too.
C Note about preceding line. The default 1.E-8 is too big, and so is 1.E-9.
C This change to DCNEW-25 is made at the same time as that to the 2nd subcase
C of DC-31. See that file for what happens without a small enough EPSILN.
TACS HYBRID
C $INCLUDE aut5-tac.inc
C $INCLUDE aut5-tck.inc { Uses new Type-68 TACS device, so much smaller file
C TACS STEADY-STATE CORE NFLUX INITIALIZATION FOR TRANSFORMER AUT5
C NFLUX IS PROPORTIONAL TO THE INDUCTIVE CURRENT IN THESE SWITCHES:
C TYPE 91 TACS SOURCES DERIVED FROM NETWORK SWITCH CURRENTS:
C <NAME> <----A---><----A---><----A---> <T START ><T STOP >
91AUT5D1 -1.0 999.
91AUT5D2 -1.0 999.
91AUT5D3 -1.0 999.
91AUT5D4 -1.0 999.
91AUT5D5 -1.0 999.
C ZERO-ORDER BLOCKS TO SCALE ABOVE TO NFLUX LINKED (NREF*PHI) WITHIN CORE VOLUME
C GAIN IS BASED ON 10000-HENRY FLUX-MONITORING INDUCTANCES.
C <NAME> +<NAME> +<NAME> +<NAME> +<NAME> +<NAME> <GAIN><FXLO><FXHI><NMLO><NMHI>
0AUT5E1 +AUT5D1 1.E4
0AUT5E2 +AUT5D2 1.E4
0AUT5E3 +AUT5D3 1.E4
0AUT5E4 +AUT5D4 1.E4
0AUT5E5 +AUT5D5 1.E4
C Correction of Type-68 TACS device lines on 6 April 2000. VAX ATP error was
C traced to the following data cards (one line for each Type-68 device) having
C E-field numbers that were not right-adjusted. The first two show how data
C columns were misalligned prior to correction:
C AUT5T168AUT5E1 13.821 27.642 2.91727E-3 8.44316E-4
C AUT5T268AUT5E1 27.642 62.194 2.43012E-3 2.53199E-4
C The associated data format for 4E16.0 read in GUTS2A was not properly
C protected by the addition of BN, so exponents were scaled by a factor of
C 100. For example, the first value on the first card was 2.91E-300 instead
C of E-3 as intended. After proper allignment, here is the data:
C NEW TACS DEVICE 68 MODELS HYSTERETIC LOSS/RESIDUAL MMF COMPONENTS.
C <IRES>68<FLUX><LOWER FLUX VAL><HIGHR FLUX VAL><LINEAR COEFF. ><SQUARED COEFF.>
99AUT5T168AUT5E1 13.821 27.642 2.91727E-3 8.44316E-4
99AUT5T268AUT5E1 27.642 62.194 2.43012E-3 2.53199E-4
99AUT5T368AUT5E1 62.194 82.925 1.23361E-3 4.76044E-4
99AUT5T468AUT5E1 82.925 103.66 1.22440E-3 5.90543E-4
99AUT5T568AUT5E1 103.66 120.25 1.79633E-3 1.56980E-3
99AUT5T668AUT5E1 120.25 131.34 1.44373E-3 3.08189E-3
99AUT5T768AUT5E1 131.34 139.51 4.95355E-3 2.07092E-2
99AUT5U168AUT5E2 13.821 27.642 2.91727E-3 8.44316E-4
99AUT5U268AUT5E2 27.642 62.194 2.43012E-3 2.53199E-4
99AUT5U368AUT5E2 62.194 82.925 1.23361E-3 4.76044E-4
99AUT5U468AUT5E2 82.925 103.66 1.22440E-3 5.90543E-4
99AUT5U568AUT5E2 103.66 120.25 1.79633E-3 1.56980E-3
99AUT5U668AUT5E2 120.25 131.34 1.44373E-3 3.08189E-3
99AUT5U768AUT5E2 131.34 139.51 4.95355E-3 2.07092E-2
99AUT5V168AUT5E3 13.821 27.642 2.91727E-3 8.44316E-4
99AUT5V268AUT5E3 27.642 62.194 2.43012E-3 2.53199E-4
99AUT5V368AUT5E3 62.194 82.925 1.23361E-3 4.76044E-4
99AUT5V468AUT5E3 82.925 103.66 1.22440E-3 5.90543E-4
99AUT5V568AUT5E3 103.66 120.25 1.79633E-3 1.56980E-3
99AUT5V668AUT5E3 120.25 131.34 1.44373E-3 3.08189E-3
99AUT5V768AUT5E3 131.34 139.51 4.95355E-3 2.07092E-2
99AUT5W168AUT5E4 14.120 28.239 2.69692E-3 7.64028E-4
99AUT5W268AUT5E4 28.239 63.538 2.24656E-3 2.29122E-4
99AUT5W368AUT5E4 63.538 84.717 1.14043E-3 4.30775E-4
99AUT5W468AUT5E4 84.717 105.90 1.13191E-3 5.34391E-4
99AUT5W568AUT5E4 105.90 122.84 1.66066E-3 1.42058E-3
99AUT5W668AUT5E4 122.84 134.17 1.33473E-3 2.79024E-3
99AUT5W768AUT5E4 134.17 142.42 4.58432E-3 1.91911E-2
99AUT5X168AUT5E5 14.120 28.239 2.69692E-3 7.64028E-4
99AUT5X268AUT5E5 28.239 63.538 2.24656E-3 2.29122E-4
99AUT5X368AUT5E5 63.538 84.717 1.14043E-3 4.30775E-4
99AUT5X468AUT5E5 84.717 105.90 1.13191E-3 5.34391E-4
99AUT5X568AUT5E5 105.90 122.84 1.66066E-3 1.42058E-3
99AUT5X668AUT5E5 122.84 134.17 1.33473E-3 2.79024E-3
99AUT5X768AUT5E5 134.17 142.42 4.58432E-3 1.91911E-2
C CURRENTS INJECTED INTO NETWORK FOR HYSTERESIS MODELING AND OTHER COMPENSATION.
C PHASE LEGS HAVE CURRENTS INJECTED TO COMPENSATE COUPLING XFMR. MAGNET. AMPS.,
C FLUX-MONITORING INDUCTANCE CURRENT AND SMALL CORE GAPS:
99AUT5TP =+AUT5D1/( 1.30893E-2)+ 1*(+AUT5T1+AUT5T2+AUT5T3+AUT5T4+AUT5T5)
99AUT5A5 =+AUT5TP+ 1*(+AUT5T6+AUT5T7)
99AUT5TQ =+AUT5D2/( 1.30893E-2)+ 1*(+AUT5U1+AUT5U2+AUT5U3+AUT5U4+AUT5U5)
99AUT5B5 =+AUT5TQ+ 1*(+AUT5U6+AUT5U7)
99AUT5TR =+AUT5D3/( 1.30893E-2)+ 1*(+AUT5V1+AUT5V2+AUT5V3+AUT5V4+AUT5V5)
99AUT5C5 =+AUT5TR+ 1*(+AUT5V6+AUT5V7)
C YOKES ARE COMPENSATED FOR MONITORING INDUCTANCE CURRENTS:
99AUT5TS =+AUT5D4+ 1*(+AUT5W1+AUT5W2+AUT5W3+AUT5W4+AUT5W5)
99AUT501 =+AUT5TS+ 1*(+AUT5W6+AUT5W7)
99AUT5TT =+AUT5D5+ 1*(+AUT5X1+AUT5X2+AUT5X3+AUT5X4+AUT5X5)
99AUT502 =+AUT5TT+ 1*(+AUT5X6+AUT5X7)
C YOKE & PHASE HYSTERETIC CURRENTS ARE EXTRACTED, EXCEPT AT GROUND:
99AUT5B0 =-AUT501-AUT5B5
99AUT5C0 =-AUT502-AUT5C5
C End of $INCLUDE. File name = aut5-tck.inc
C PARTS OF HYSTERESIS COMPONENTS FOR PLOTTING
C FORTRAN STATEMENTS; 99= INPUT; 98= OUTPUT; 88= INSIDE
C <NAME> =< FREE FORMAT FORTRAN TO COL 80 ----->
C TOTAL OF ALL HYSTER INJECTIONS FOR LEGS A,B,C
99ALLINA =+AUT5T1+AUT5T2+AUT5T3+AUT5T4+AUT5T5+AUT5T6+AUT5T7
99ALLINB =+AUT5U1+AUT5U2+AUT5U3+AUT5U4+AUT5U5+AUT5U6+AUT5U7
99ALLINC =+AUT5V1+AUT5V2+AUT5V3+AUT5V4+AUT5V5+AUT5V6+AUT5V7
C
C TYPE 90 TACS SOURCES DERIVED FROM NETWORK VOLTAGES:
C <NAME> <----A---><----A---><----A---> <T START ><T STOP >
C TYPE 91 TACS SOURCES DERIVED FROM NETWORK SWITCH CURRENTS (1ST NODE)
C <NAME> <----A---><----A---><----A---> <T START ><T STOP >
91AUT5AM 0. 999.
91AUT5BM 0. 999.
91AUT5CM 0. 999.
C FORTRAN STATEMENTS; 99= INPUT; 98= OUTPUT; 88= INSIDE
C <NAME> =< FREE FORMAT FORTRAN TO COL 80 ----->
C NEXT IS SUM OF COUPLING OUTPUT AND ALL INJECTED AMPS - USED BY CORE MODELS
98LEGAMP = AUT5AM + AUT5A5
98LEGBMP = AUT5BM + AUT5B5
98LEGCMP = AUT5CM + AUT5C5
C NEXT IS APPARENT LEG AMPS FOR HYSTER PLOTTING
98NETAMP = LEGAMP - ALLINA
98NETBMP = LEGBMP - ALLINB
98NETCMP = LEGCMP - ALLINC
C S-BLOCKS OF ORDER 1 IN COL 2 INTEGRATE VOLTAGES TO GET FLUXES (TIMES 250.0).
C <NAME> +<NAME> +<NAME> +<NAME> +<NAME> +<NAME> <GAIN><FXLO><FXHI><NMLO><NMHI>
C N0 & D0><N1 & D1 ><N2 & D2 ><N3 & D3 ><N4 & D4 ><N5 & D5 ><N6 & D6 ><N7 & D7 >
C <NAME>CD+< IN1> +< IN2> +< IN3> +< IN4> +< IN5> < A >< B >< C >< D >< E >
C TACS OUTPUT REQUESTS - TYPE 33
C <NAME><NAME><NAME><NAME><NAME><NAME><NAME><NAME><NAME><NAME><NAME><NAME><TY33>
33AUT5E1 ALLINA LEGAMP NETAMP
33AUT5E2 ALLINB LEGBMP NETBMP
33AUT5E3 ALLINC LEGCMP NETCMP
C 33AUT5T1 AUT5T2 AUT5T3 AUT5T4 AUT5T5 AUT5T6 AUT5
33AUT5U1 AUT5U2 AUT5U3 AUT5U4 AUT5U5 AUT5U6 AUT5U7
BLANK CARD ENDING ALL ATP-SORTED TACS CARDS (from blank.inc)
C THIS MUST BE EDITED TO ADD AUT5A5 >> AUT5AM, ETC METERING
C $INCLUDE aut5-brn.inc
C WSM makes non-Unix by adding a card using DOS editor.
C 3-PHASE, 3-LEG SATURABLE XFMR. MODEL FROM MODELAD.FTN - 60. HZ IMPEDANCES
C CHANGED AUT5A5, B5, C5 TO AUT5AM, BM, CM ; NOW REQUIRES METERING SWITCHES
C AUT5 = IDENTIFYING BUS CODE PREFIX OF INTERNAL NODES.
C AUT5AM = NAME OF A METERING NODE ADDED AT AUT5A5
C EACH PHASE LEG HAS 3 WINDINGS WHOSE RATED VOLTAGES AND TURNS ARE:
C WINDING 1 (INNERMOST) = 13.130 KV. TURNS = 44.998
C WINDING 1 INCLUDES THE EFFECT OF 0.5670 OHMS/PHASE INTERNAL (DELTA) REACTORS.
C REACTOR RESISTANCE/PHASE MUST BE INCLUDED WITH ENTERED WINDING RESISTANCES.
C WINDING 2 = 79.674 KV. TURNS = 273.052
C WINDING 3 = 119.655 KV. TURNS = 410.072
C INTERNAL REFERENCE = 29.179 KV. TURNS = 100.000
C *********BEGIN USER COMMENTS ************
C EAST GARDEN CITY AUTOTRANSFORMER*** one *** AT 345/138 KV TAP
C ABBREVIATED HYSTERESIS MODEL; SEVEN OF POSSIBLE 11 COMPONENTS
C .95 * .7274 SQUARE METER LEG; .97 * .7274 SQUARE METER YOKE
C ACTUAL DELTA IS GROUNDED BETWEEN OUTER PHASES AT Y1.
C EACH DELTA REACTOR X =.82;OR .5467 WHEN EQUALIZED; REACTOR R = .00427 EACH
C ACTUAL DELTA WINDING R = .00510; + .002847 TO EQUALIZE = .0079467 ENTERED
C AIR-CORE COUPLINGS FROM COUPCOIL, BASED ON SMIT DIMENSIONS.
C B-H CURVE BASED ON ARMCO TRAN-COR H-0 from smit
C B-H SHIFTED TO HIGHER CURRENTS IN 1.5 - 1.85 T RANGE
C MILD OPEN HEARTH STEEL FOR TANK
C TANK HAS 21.65 M CIRCUM; FLUX SHIELD = 3.15 M HIGH
C 64-36 EDDY-HYST LOSS RATIO AT RATED;
C ********** END USER COMMENTS ************
C COUPLING, LEAKAGE IMPEDANCE & WINDING LOSS MODEL FOR TRANSFORMER AUT5:
C <NAME><NAME><NAME><NAME>< R >< X >< R >< X >< R >< X >
51Y1 Y2 .01016 10124.87377
52EGCT1A 0.000 61435.16312.04458372798.88582
53EGCHEAEGCT1A 0.000 92263.79188 0.000559873.90936.12388840864.72923
54AUT5A1 0.000 22499.42979 0.000136528.47263 0.000205039.49115
0.000 50000.88921
55AUT5A2 0.000 22499.42071 0.000136530.26291 0.000205042.81249
0.000 50000.86902 0.000 50001.52468
56AUT5A3 0.000 22499.42350 0.000136530.70112 0.000205053.04695
0.000 50000.87522 0.000 50001.68516 0.000 50004.18094
57AUT5A4 0.000 22499.25174 0.000136528.91954 0.000205046.26382
0.000 50000.49352 0.000 50001.03270 0.000 50002.52680
0.000 50002.52680
58AUT5AM 0.000 22499.02967 0.000136526.09974 0.000205035.90210
0.000 50000.00000 0.000 50000.00000 0.000 50000.00000
0.000 50000.00000 0.000 50000.00000
51Y2 Y3 .01016 10124.87377
52EGCT1B 0.000 61435.16312.04458372798.88582
53EGCHEBEGCT1B 0.000 92263.79188 0.000559873.90936.12388840864.72923
54AUT5B1AUT5B0 0.000 22499.42979 0.000136528.47263 0.000205039.49115
0.000 50000.88921
55AUT5B2AUT5B0 0.000 22499.42071 0.000136530.26291 0.000205042.81249
0.000 50000.86902 0.000 50001.52468
56AUT5B3AUT5B0 0.000 22499.42350 0.000136530.70112 0.000205053.04695
0.000 50000.87522 0.000 50001.68516 0.000 50004.18094
57AUT5B4AUT5B0 0.000 22499.25174 0.000136528.91954 0.000205046.26382
0.000 50000.49352 0.000 50001.03270 0.000 50002.52680
0.000 50002.52680
58AUT5BMAUT5B0 0.000 22499.02967 0.000136526.09974 0.000205035.90210
0.000 50000.00000 0.000 50000.00000 0.000 50000.00000
0.000 50000.00000 0.000 50000.00000
51Y3 Y1 .01016 10124.87377
52EGCT1C 0.000 61435.16312.04458372798.88582
53EGCHECEGCT1C 0.000 92263.79188 0.000559873.90936.12388840864.72923
54AUT5C1AUT5C0 0.000 22499.42979 0.000136528.47263 0.000205039.49115
0.000 50000.88921
55AUT5C2AUT5C0 0.000 22499.42071 0.000136530.26291 0.000205042.81249
0.000 50000.86902 0.000 50001.52468
56AUT5C3AUT5C0 0.000 22499.42350 0.000136530.70112 0.000205053.04695
0.000 50000.87522 0.000 50001.68516 0.000 50004.18094
57AUT5C4AUT5C0 0.000 22499.25174 0.000136528.91954 0.000205046.26382
0.000 50000.49352 0.000 50001.03270 0.000 50002.52680
0.000 50002.52680
58AUT5CMAUT5C0 0.000 22499.02967 0.000136526.09974 0.000205035.90210
0.000 50000.00000 0.000 50000.00000 0.000 50000.00000
0.000 50000.00000 0.000 50000.00000
C PHASE LEG N-L MODELS FOR TRANSFORMER: AUT5
C PHASE LEG CORES ARE 3.60 METERS LONG WITH AREAS OF 0.6910 SQUARE METERS
C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS:
C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00)
C FOR SS CALCS EACH LEG HAS REACT. OF 15740. OHMS (AT 100.0 TURNS)
C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA.
C HYSTERESIS CURRENT OFFSET OF 3.115 AMPS IS MODELED AT SATURATION.
C LAMINATIONS ARE MODELED BY 8 LAYERS; EDDY LOSSES USE ENTERED CORE RESISTIVITY.
C <NAME><NAME><NAME><NAME><I SS><P SS> (--FOR TYPE 98 N/L REACTOR--) P
C CURRENT -----><NFLUX -------->
98AUT5A5AUT5A6 4.912234.877
0.14400 3.45524
0.39528 6.91044
0.81635 15.54839
1.12175 20.73117
1.41336 24.18640
2.00198 25.91428
3.13598 28.16067
4.23399 30.06154
8.81421 31.96448
19.73322 32.83478
56.33246 33.72028
543.11535 34.87714
1083.11535 35.20278
9999
98AUT5A6AUT5A7AUT5A5AUT5A64.912234.877
98AUT5A7AUT5A8AUT5A5AUT5A64.912234.877
98AUT5A8 AUT5A5AUT5A64.912234.877
98AUT5B5AUT5B6AUT5A5AUT5A64.912234.877
98AUT5B6AUT5B7AUT5A5AUT5A64.912234.877
98AUT5B7AUT5B8AUT5A5AUT5A64.912234.877
98AUT5B8AUT5B0AUT5A5AUT5A64.912234.877
98AUT5C5AUT5C6AUT5A5AUT5A64.912234.877
98AUT5C6AUT5C7AUT5A5AUT5A64.912234.877
98AUT5C7AUT5C8AUT5A5AUT5A64.912234.877
98AUT5C8AUT5C0AUT5A5AUT5A64.912234.877
C HIGH RESISTANCE ACROSS EACH UNPARALLELED TYPE 98 ELEMENT ABOVE:
C <NAME><NAME><----><---->RRRRRRXXXXXXCCCCCC
AUT5A6AUT5A5 .1E+11
AUT5A7AUT5A6AUT5A6AUT5A5
AUT5A8AUT5A7AUT5A6AUT5A5
AUT5B6AUT5B5AUT5A6AUT5A5
AUT5B7AUT5B6AUT5A6AUT5A5
AUT5B8AUT5B7AUT5A6AUT5A5
AUT5C6AUT5C5AUT5A6AUT5A5
AUT5C7AUT5C6AUT5A6AUT5A5
AUT5C8AUT5C7AUT5A6AUT5A5
C YOKE MODELS FOR TRANSFORMER: AUT5
C EACH TOP+BOTTOM YOKE LENGTH IS 3.40 METERS WITH AREA OF 0.7060 SQUARE METERS
C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS:
C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00)
C FOR SS CALCS EACH YOKE HAS REACT. OF 17026. OHMS (AT 100.0 TURNS)
C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA.
C HYSTERESIS CURRENT OFFSET OF 2.943 AMPS IS MODELED AT SATURATION.
C LAMINATIONS ARE MODELED BY 8 LAYERS; EDDY LOSSES USE ENTERED CORE RESISTIVITY.
C <NAME><NAME><NAME><NAME><I SS><P SS> (--FOR TYPE 98 N/L REACTOR--) P
C CURRENT -----><NFLUX -------->
98AUT501AUT50A 4.636135.606
0.13600 3.52991
0.37332 7.05979
0.77100 15.88444
1.05943 21.17924
1.33484 24.70914
1.89076 26.47434
2.96176 28.76924
3.99877 30.71116
8.32457 32.65503
18.63694 33.54364
53.20871 34.44660
512.94277 35.60585
1022.94277 35.91340
9999
98AUT50AAUT50BAUT501AUT50A4.636135.606
98AUT50BAUT50CAUT501AUT50A4.636135.606
98AUT50CAUT5B0AUT501AUT50A4.636135.606
98AUT502AUT50FAUT501AUT50A4.636135.606
98AUT50FAUT50GAUT501AUT50A4.636135.606
98AUT50GAUT50HAUT501AUT50A4.636135.606
98AUT50HAUT5C0AUT501AUT50A4.636135.606
C HIGH RESISTANCE ACROSS EACH UNPARALLELED TYPE 98 ELEMENT ABOVE
C <NAME><NAME><----><---->RRRRRRXXXXXXCCCCCC
AUT50AAUT501AUT5A6AUT5A5
AUT50BAUT50AAUT5A6AUT5A5
AUT50CAUT50BAUT5A6AUT5A5
AUT50FAUT502AUT5A6AUT5A5
AUT50GAUT50FAUT5A6AUT5A5
AUT50HAUT50GAUT5A6AUT5A5
$VINTAGE, 1
C <NAME><NAME><----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P
C SHUNT RESISTANCES FOR EDDY LOSS MODELING OF 1/8 LAMINATION THICKNESS LAYERS.
C CROSS-FLUX CORE RESISTIVITY IS 12.00 E-8 OHM-METERS.
C LAMINATIONS ARE 0.2286 MM THICK; MODELED LAYERS ARE 0.0286 MM THICK.
AUT5A5 141050.02696
AUT5A6 70525.01348
AUT5A7 70525.01348
AUT5A8 70525.01348
AUT5B5AUT5B0 141050.02696
AUT5B6AUT5B0 70525.01348
AUT5B7AUT5B0 70525.01348
AUT5B8AUT5B0 70525.01348
AUT5C5AUT5C0 141050.02696
AUT5C6AUT5C0 70525.01348
AUT5C7AUT5C0 70525.01348
AUT5C8AUT5C0 70525.01348
AUT501AUT5B0 152575.09208
AUT50AAUT5B0 76287.54604
AUT50BAUT5B0 76287.54604
AUT50CAUT5B0 76287.54604
AUT502AUT5C0 152575.09208
AUT50FAUT5C0 76287.54604
AUT50GAUT5C0 76287.54604
AUT50HAUT5C0 76287.54604
C COMBINED WINDING 1-2 AND WINDING 2-3 STRAY (EDDY CUR.) LOSS MODEL:
C 3-PHASE W1-W2 STRAY LOSSES AT RATED MVA = 636.953 KW.
C WITH W1-W2 REF.VOLTS = 2.085 KV.; AND W2-W3 REF.VOLTS = 0.359 KV.
AUT5A1AUT5A2 20.47385
AUT5B1AUT5B2 20.47385
AUT5C1AUT5C2 20.47385
C 3-PHASE W2-W3 STRAY LOSSES AT RATED MVA = 127.094 KW.
C WITH W1-W2 REF.VOLTS = -0.495 KV.; AND W2-W3 REF.VOLTS = 7.203 KV.
AUT5A2AUT5A3 1706.84712
AUT5B2AUT5B3 1706.84712
AUT5C2AUT5C3 1706.84712
C TANK WALL MULTIPLE CHAINED PI MODEL RESISTANCES REFERRED TO 100.000 TURNS:
C MODELED SURFACE LAYERS ARE 0.1279 MM THICK; RHOTANK = 15.00 E-8 O-M.
C <NAME><NAME><----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P
AUT581 1.67957937
AUT582 3.35915874
AUT583 6.71831748
AUT584 13.43663496
AUT585 26.87326992
AUT586 53.74653984
AUT587 80.61980976
AUT504 161.23961952
$VINTAGE, 0
C TRANSFORMER TANK WALL MODEL FOR TRANSFORMER AUT5
C TANK WALL AND VERTICAL BRACING HAS A CROSS-SECTIONAL AREA OF 0.354 SQ. METERS.
C ITS INSIDE PERIMETER IS 21.65 METERS. ITS AVERAGE THICKNESS IS 0.016 METERS.
C EQUIVALENT HEIGHT AFFECTED BY ZERO SEQUENCE FLUX IS ASSUMED AT 3.15 METERS.
C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS:
C LINEAR:(0.20, 70.00) AND FULLY SATURATED:(2.05, 75000.)
C STEADY-STATE CALCS USE THE SLOPES AT POINTS 3 BELOW. (1.0000 TESLA)
C <NAME><NAME><NAME><NAME><I SS><P SS> (--FOR TYPE 98 N/L REACTOR--) P
C CURRENT ----->< FLUX -------->
98AUT581 21.31036.333
2.20500 3.54468
5.04000 10.63404
10.39500 17.72340
22.05000 23.04042
50.40000 26.58510
110.25000 28.35744
226.80000 30.12978
441.00000 31.90212
1134.00000 35.44680
2362.50000 36.33297
4725.00000 38.03714
9999
98AUT582AUT581 21.31018.166
2.20500 1.77234
5.04000 5.31702
10.39500 8.86170
22.05000 11.52021
50.40000 13.29255
110.25000 14.17872
226.80000 15.06489
441.00000 15.95106
1134.00000 17.72340
2362.50000 18.16648
4725.00000 19.01857
9999
98AUT583AUT582 21.3109.0832
2.20500 0.88617
5.04000 2.65851
10.39500 4.43085
22.05000 5.76011
50.40000 6.64628
110.25000 7.08936
226.80000 7.53245
441.00000 7.97553
1134.00000 8.86170
2362.50000 9.08324
4725.00000 9.50929
9999
98AUT584AUT583 21.3104.5416
2.20500 0.44309
5.04000 1.32926
10.39500 2.21543
22.05000 2.88005
50.40000 3.32314
110.25000 3.54468
226.80000 3.76622
441.00000 3.98777
1134.00000 4.43085
2362.50000 4.54162
4725.00000 4.75464
9999
98AUT585AUT584 21.3102.2708
2.20500 0.22154
5.04000 0.66463
10.39500 1.10771
22.05000 1.44003
50.40000 1.66157
110.25000 1.77234
226.80000 1.88311
441.00000 1.99388
1134.00000 2.21543
2362.50000 2.27081
4725.00000 2.37732
9999
98AUT586AUT585 21.3101.1354
2.20500 0.11077
5.04000 0.33231
10.39500 0.55386
22.05000 0.72001
50.40000 0.83078
110.25000 0.88617
226.80000 0.94156
441.00000 0.99694
1134.00000 1.10771
2362.50000 1.13541
4725.00000 1.18866
9999
98AUT587AUT586 21.310.56770
2.20500 0.05539
5.04000 0.16616
10.39500 0.27693
22.05000 0.36001
50.40000 0.41539
110.25000 0.44309
226.80000 0.47078
441.00000 0.49847
1134.00000 0.55386
2362.50000 0.56770
4725.00000 0.59433
9999
98AUT504AUT587 21.310.56770
2.20500 0.05539
5.04000 0.16616
10.39500 0.27693
22.05000 0.36001
50.40000 0.41539
110.25000 0.44309
226.80000 0.47078
441.00000 0.49847
1134.00000 0.55386
2362.50000 0.56770
4725.00000 0.59433
9999
C FLUX SHIELD N-L MODEL FOR TRANSFORMER: AUT5
C FLUX SHIELDS ARE 3.15 METERS LONG WITH AREA OF 0.3652 SQUARE METERS.
C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS:
C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00)
C FOR SS CALCS EACH F.S. HAS REACT. OF 9505. OHMS (AT 100.0 TURNS)
C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA.
C <NAME><NAME><NAME><NAME><I SS><P SS> (--FOR TYPE 98 N/L REACTOR--) P
C CURRENT -----><NFLUX -------->
98AUT503AUT504 2.896673.034
0.12600 7.30340
0.20475 14.60680
0.30870 32.86530
0.39690 43.82040
0.50400 51.12380
0.94500 54.77550
1.73250 59.52271
2.52000 63.53958
6.30000 67.55645
15.75000 69.38230
47.25000 71.20815
472.50000 73.03400
945.00000 75.06272
9999
C RESISTANCE UNDERESTIMATING FLUX SHIELD LOSSES, EQUAL TO LINEAR REACTANCE.
C <NAME><NAME><----><---->RRRRRRXXXXXXCCCCCC P
AUT503AUT504 9505.5
C TANK TOP MULTIPLE CHAINED PI MODEL RESISTANCES REFERRED TO 100.000 TURNS:
C MODELED SURFACE LAYERS ARE 0.1172 MM THICK; RHOTANK = 15.00 E-8 O-M.
$VINTAGE, 1
C <NAME><NAME><----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P
AUT591 2.00000000
AUT592 4.00000000
AUT593 8.00000000
AUT594 16.00000000
AUT595 32.00000000
AUT596 64.00000000
AUT597 96.00000000
AUT509 192.00000000
$VINTAGE, 0
C TRANSFORMER TANK TOP MODEL FOR TRANSFORMER AUT5
C TANK TOP EDDY CURRENT PATH IS MODELED WITH A CIRCUMFERENCE OF 15.00 METERS.
C AND A WIDTH OF 2.00 METERS; THICKNESS IS 0.015 METERS.
C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS:
C LINEAR:(0.20, 70.00) AND FULLY SATURATED:(2.05, 75000.)
C STEADY-STATE CALCS USE THE SLOPES AT POINTS 3 BELOW.
C <NAME><NAME><NAME><NAME><I SS><P SS> (--FOR TYPE 98 N/L REACTOR--) P
C CURRENT ----->< FLUX -------->
98AUT591 13.53023.062
1.40000 2.25000
3.20000 6.75000
6.60000 11.25000
14.00000 14.62500
32.00000 16.87500
70.00000 18.00000
144.00000 19.12500
280.00000 20.25000
720.00000 22.50000
1500.00000 23.06250
3000.00000 24.14423
9999
98AUT592AUT591 13.53011.531
1.40000 1.12500
3.20000 3.37500
6.60000 5.62500
14.00000 7.31250
32.00000 8.43750
70.00000 9.00000
144.00000 9.56250
280.00000 10.12500
720.00000 11.25000
1500.00000 11.53125
3000.00000 12.07212
9999
98AUT593AUT592 13.5305.7656
1.40000 0.56250
3.20000 1.68750
6.60000 2.81250
14.00000 3.65625
32.00000 4.21875
70.00000 4.50000
144.00000 4.78125
280.00000 5.06250
720.00000 5.62500
1500.00000 5.76562
3000.00000 6.03606
9999
98AUT594AUT593 13.5302.8828
1.40000 0.28125
3.20000 0.84375
6.60000 1.40625
14.00000 1.82812
32.00000 2.10938
70.00000 2.25000
144.00000 2.39062
280.00000 2.53125
720.00000 2.81250
1500.00000 2.88281
3000.00000 3.01803
9999
98AUT595AUT594 13.5301.4414
1.40000 0.14063
3.20000 0.42187
6.60000 0.70312
14.00000 0.91406
32.00000 1.05469
70.00000 1.12500
144.00000 1.19531
280.00000 1.26563
720.00000 1.40625
1500.00000 1.44141
3000.00000 1.50901
9999
98AUT596AUT595 13.530.72070
1.40000 0.07031
3.20000 0.21094
6.60000 0.35156
14.00000 0.45703
32.00000 0.52734
70.00000 0.56250
144.00000 0.59766
280.00000 0.63281
720.00000 0.70312
1500.00000 0.72070
3000.00000 0.75451
9999
98AUT597AUT596 13.530.36035
1.40000 0.03516
3.20000 0.10547
6.60000 0.17578
14.00000 0.22852
32.00000 0.26367
70.00000 0.28125
144.00000 0.29883
280.00000 0.31641
720.00000 0.35156
1500.00000 0.36035
3000.00000 0.37725
9999
98AUT509AUT597 13.530.36035
1.40000 0.03516
3.20000 0.10547
6.60000 0.17578
14.00000 0.22852
32.00000 0.26367
70.00000 0.28125
144.00000 0.29883
280.00000 0.31641
720.00000 0.35156
1500.00000 0.36035
3000.00000 0.37725
9999
C ZERO SEQUENCE MODELS ARE BASED ON AIR-CORE MODELS, PLUS DERIVED VALUE FOR
C X100T OF 3.29690 OHMS/PHASE TANK EFFECT, REFERRED TO 100 TURNS.
C THE (UNSHIELDED) TANK TOP IS MODELED AND HAS STRAY LOSSES.
C THE TANK WALL IS SATURABLE; BUT ANY STRAY LOSSES ARE NOT MEASURABLE.
C THE FLUX SHIELD IS SATURABLE AND LACKS SIGNIFICANT DISTRIBUTED GAPS.
C COUPLING REPRESENTING TOP AND SIDE GAP REACTANCES:
C PER-PHASE SIDE GAPS = 1.846 OHMS; TOP GAPS = 1.814 OHMS AT 100 TURNS.
51AUT503 0.000 10.978117
52AUT509 0.000 5.4405450.000 5.440545
C SHUNT FLUX PATHS BYPASSING LEGS AND COILS OF EACH PHASE FOLLOW:
C <NAME><NAME><----><---->RRRRRRXXXXXXXXXXXX
51AUT5A4AUT5B0 0.000 7.463166
51AUT5B4AUT5C0 0.000 7.463166
51AUT5C4AUT503 0.000 7.463166
C ZERO-SEQUENCE FLUXES FROM EACH PHASE TO TANK ARE EQUALIZED BY THE FOLLOWING:
C <NAME><NAME><----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX
51AUT503AUT502 0.000 3296.896640
52AUT502AUT501 0.000 3296.8933430.000 3296.896640
53AUT501 0.000 3296.8933430.000 3296.8933430.000 3296.896640
C NEXT COUPLING SUMS ZERO-SEQUENCE VOLTAGES OF OUTER WINDINGS AT NODE (0,8):
C <NAME><NAME><----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX
51AUT5A3 0.000 354922.42080
52AUT506 0.000 354922.065870.000 354922.42080
51AUT5B3AUT5B0 0.000 354922.42080
52AUT507AUT506 0.000 354922.065870.000 354922.42080
51AUT5C3AUT5C0 0.000 354922.42080
52AUT508AUT507 0.000 354922.065870.000 354922.42080
C NON-DEDICATED STRAY ZERO-SEQUENCE LOSS RESISTANCE:
C PRODUCES 158.213 KW LOSS FOR ENTERED TEST CONDITIONS.
$VINTAGE, 1
C <NAME><NAME><----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P
AUT508 354.92242080
$VINTAGE, 0
C TIME-STEP-DAMPING BRANCHES ACROSS MAJOR AIR GAPS FOR HARMONIC ORDER: 500.0
C <NAME><NAME><----><---->RRRRRRXXXXXXCCCCCC
AUT5A4AUT5B0 373.16 .00142
AUT5B4AUT5C0 373.16 .00142
AUT5C4AUT503 373.16 .00142
AUT501 164.84 .00322
AUT501AUT502 164.84 .00322
AUT502AUT503 164.84 .00322
C 10000-HENRY FLUX MONITORING INDUCTANCES TO INITIALIZE TACS.
C <NAME><NAME><----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX
51AUT5A5AUT5D1 3769911.185
51AUT5B5AUT5D2 3769911.185
51AUT5C5AUT5D3 3769911.185
51AUT501AUT5D4 3769911.185
51AUT502AUT5D5 3769911.185
C End of $INCLUDE. File name = aut5-brn.inc
$VINTAGE, 1
C <NAME><NAME><----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P
C CORNER GROUND DELTA
Y1 1000.
C ENERGIZE VIA 100 OHM RESISTANCES FOR DAMPING
C 0.1 ohm is now essential; 100 is too large relative to xl
SRC2A EGCT1A .1
SRC2B EGCT1B .1
SRC2C EGCT1C .1
$VINTAGE, 0
BLANK CARD ENDING ALL ATP-SORTED BRANCH CARDS (from blank.inc)
C $INCLUDE aut5-swx.inc
C WSM touches using DOS editor to make readable.
C MEAS. SWITCHES TO MONITOR CURRENT IN 10000-HENRY CORE FLUX-MONITORING INDUCTS.
C <NAME><NAME><T CLOSE ><T. OPEN ><AMP MARG> 67890123 MEASURING 56789012345678 P
AUT5D1 MEASURING
AUT5D2AUT5B0 MEASURING
AUT5D3AUT5C0 MEASURING
AUT5D4AUT5B0 MEASURING
AUT5D5AUT5C0 MEASURING
C End of $INCLUDE. File name = aut5-swx.inc
C ********************************************* STUDY-DEPENDENT CABLE SWITCHING:
C SWITCH CARD: COL 1-2 IS 0 FOR ORDINARY & GAUSS. STATISTICS SWITCHES, NOT TACS.
C CLOSING, OPENING TIMES AND STATISTICS PARAMETERS BELOW ARE IN SECONDS.
C <NAME><NAME><T CLOSE ><T. OPEN ><AMP MARG> 67890123 MEASURING 56789012345678 P
C <NAME><NAME><GAUS AVG><GAUS SIG><AMP MARG> STATISTICS P
C METERS INTERNAL CURRENT OF PHASE LEG A,B,C
AUT5AMAUT5A5 MEASURING 1
AUT5BMAUT5B5 MEASURING 1
AUT5CMAUT5C5 MEASURING 1
C ENERGIZED AT ZERO
SRCA SRC2A MEASURING 1
SRCB SRC2B MEASURING 1
SRCC SRC2C MEASURING 1
C ENERGIZED AT ZERO PLUS
C SRCA SRC2A .0001 999. 1.
C SRCB SRC2B .0001 999. 1.
C SRCC SRC2C .0001 999. 1.
BLANK CARD ENDING ALL ATP-SORTED SWITCH CARDS (from blank.inc)
$LISTON
C $INCLUDE aut5-src.inc
C WSM touches file with DOS editor to make it into a DOS file
C TACS-CONTROLLED SOURCES FOR HYSTERESIS/RESIDUAL MMF IN TRANSFORMER AUT5
C <NAME> *= VOLTAGE IF POSITIVE; CURRENT IF NEGATIVE
60AUT5A5-1
60AUT5B5-1
60AUT5C5-1
60AUT5B0-1
60AUT5C0-1
60AUT501-1
60AUT502-1
C End of $INCLUDE. File name = aut5-src.inc
C *= VOLTAGE IF POSITIVE; CURRENT IF NEGATIVE
C <NAME> <CRESTVAL><FREQ ><ANGLE > <T START ><T STOP >
C 138 KV L-L = 112676.5 V CREST or 1877.94 at 1 hertz
C 40 percent = 751.17
C 60 percent = 1126.76
C 80 percent = 1502.35
C 100 percent = 1877.94
C 120 percent = 2253.53
14SRCA 1 2253.53 1. 0. -1.0 999.
14SRCB 1 2253.53 1. -120. -1.0 999.
14SRCC 1 2253.53 1. -240. -1.0 999.
BLANK CARD ENDING ALL ATP-SORTED SOURCE CARDS (from blank.inc)
C <BUS ><VOLT><MON.><FOR ><NODE>< >< >< >< >< >< >< >< >
EGCT1A EGCT1B EGCT1C
BLANK CARD ENDING ALL ATP-SORTED OUTPUT CARDS (from blank.inc)
193.01 0.0 0.1 BRANCH
TACS ALLINATACS ALLINBTACS ALLINC
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd of 18 subcases.
C "FREQUENCY SCAN" use, with subsequent plotting of phasors vs. frequency.
C This test case uses uniform frequency spacing. It basically like DC-51,
C but involves some variations of output and batch-mode plotting.
PRINTED NUMBER WIDTH, 8, 2, { Each column of width 8 includes 2 blank bytes
FREQUENCY SCAN, 60., 10., 280., 0, { 60 < f < 280 Hz in 10-Hz increments
CHANGE PLOT FREQUENCY { Request to vary frequency of plot points and .PL4 file
5 5 10 3 20 1 { On step 5, change to 5, etc.
0.1 0.0 { dT of col 1-8 is not used; TMAX of 9-16 ==> no simulation
1 1 1 0 1 -1
5 5 10 3 20 1
C Note: 1-punch in col. 80 of the following card was replaced by LIN002 below.
C This illustrates a second way to order current output: by branch name.
C When the user does not name branches himself, NMAUTO = 1 means the
C program will provide default names. The first branch is LIN001, so
C this should be used to replace the column-80 punch. But this branch
C is in series with the 2nd branch, and has no shunt leakage, so the
C same current can be found from the second branch (name LIN002). This
C alternative is more interesting since it shows current output seems to
C be possible for the 1-phase, constant-parameter distributed model.
GEN SEND 3.0 40.
-1SEND MID .306 5.82 .012 20. { 20 miles of line from DC-37
1MID REC 34.372457.68.15781 { 10 miles of Pi from DC-3
EARTH 200. 1
C The preceding branch has a 1-punch in column 80 for current output. If a 3
C had instead been used, branch voltage would have been produced, too. Note
C it would be equal to the node voltage of node EARTH. No choice of polarity
C is possible, this way. More interesting is the "-5" request below, which
C allows polarity reversal (a minus sign will be added).
BLANK card ending branches
C Note: 1-punch in col. 80 of the following card was replaced by SWT001 below.
C This illustrates a second way to order current output: by switch name.
REC EARTH MEASURING
BLANK card ending (here, non-existent) switches
14GEN 100. 60. -1.
BLANK card ending sources
C Requests for output variables follow. About currents, there is no control
C over location in the output vector. Switch currents always precede branch
C currents, and both are in natural order (of component number). But order of
C voltages is controlled by the user. Note requests for node voltages have
C been split, and in between there is a request for a voltage difference.
C This is the order in the output vector, note. Yet, all voltage requests
C made here (after the blank card ending source cards) follow all that might
C have been made with column-80 punches on branch or switch cards.
GEN SEND MID REC
-5 EARTH { -5 ==> 2A6 name pairs for voltage differences (branch V)
-1SWT001LIN002 { -1 ==> Branch/switch current out; use A6 component names
EARTH { Node voltage. Note branch voltage 2 lines above has opposite polarity
BLANK card ending node voltage outputs
PRINTER PLOT
C Note the peculiar 60 units/inch on the following plot cards. This is not
C a mistake. With 6 lines/inch of the printer plot, & the F-scan frequency
C increment of 10 Hz, there will be exactly one line for each freq step:
14660. 60.280. GEN SEND MID REC
19660. 60.280. SEND MID EARTH REC
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 3rd of 18 subcases illustrates IF-THEN-ENDIF for conditional data
C assembly as described in the July, 1998, newsletter. Environment
C variable RESISTOR must be set prior to use, in order. to avoid
C an error stop. This can not be set internally, unfortunately.
PRINTED NUMBER WIDTH, 9, 2,
POWER FREQUENCY, 50.,
.0001 .020 50.
1 1 1 1 1 -1
5 5 20 20
C 1st of 2 identical, disconnected networks uses manually-defined branches:
C IF ( RESISTOR .EQ. 'ONEHALF' ) SET VARIABLE { Execute OS command SET RESISTOR=ONEHALF
IF ( RESISTOR .EQ. 'ONEHALF' ) THEN { If user wants to use resistive alternative:
GEN TRAN 0.5 { 1st of 3 alternatives is resistor
ELSEIF ( RESISTOR .EQ. 'ZERO' ) THEN { If user wants to use resistive alternative:
GEN TRAN 0.0 { 2nd of 3 alternatives is illegal resistor
ELSE { Alternatively (if environment variable RESISTOR is not YES):
GEN TRAN 0.5 { 3rd of 3 alternatives is inductor
ENDIF { Terminate 5-line IF statement that uses RESISTOR to choose the model
TRAN 0.5
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14GEN 1.414 50. 0.0 -1.
BLANK card terminating program source cards.
GEN TRAN
BLANK card ending program output-variable requests.
BLANK card ending all plot cards
BEGIN NEW DATA CASE --- NEWSORT ---
C 4th of 18 subcases illustrates new /-card sorting as described in a story of
C the January, 1999, newsletter. Correction of trouble sorting /INITIAL was
C involved as explained by Prof. Juan A. Martinez Velasco of the Polytechnic
C University of Catalunya in Barcelona, Spain. This is a modification of his
C FILE4.DAT DELTAT is increased and # of steps is reduced. Blank cards no
C longer are required, so have been removed. A requirement of the new logic
C is that all data after miscellaneous data cards must be sorted. This was
C not true of the old logic (see DC-8). Note that some data classes are not
C being used (no problem). Note the NEWSORT declaration on line 1, which
C was patterned after the use of NOSORT in DC-36. In effect, the NEWSORT
C request will make sure SZBED carries a minus sign. Not illustrated is
C parallel OLDSORT to make sure SZBED does _not_ carry a minus sign (this
C would force use of the old sorting logic). It is assumed STARTUP remains
C as it was (no minus sign applied to SZBED) in order to handle DC-8.
.000100 .010
1 1 -1
5 5 20 20
/OUTPUT
BUS-0 BUS-1 BUS-2
/SOURCE
14BUS-0 1.0 50.0 0.0 -1.0 1.0
/BRANCH
BUS-0 BUS-1 10.0 10.0
BUS-1 10.0 { Illustrate an in-line comment
C Comment cards are allowed anywhere, of course. Add one here after
C the first two branches just to show there is no problem.
/SWITCH
BUS-1 BUS-2 .00035 1.0 1.0 1
/BRANCH
BUS-2 10.0
/INITIAL
2BUS-2 0.5
BEGIN NEW DATA CASE
C 5th of 18 subcases is an improvement over DC4AT5.DAT of SPY @5. But that
C required SPY. Beginning 23 August 1998, parametric studies can be run
C using higher-level batch-mode commands. See October, 1998, newsletter.
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 5 1 0 { Loop five times
PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns)
.0001 .020
1 -1 1 1 1 -1 2
5 5 20 20
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be
C __MILLIHENRIES__ = 1000. - ( KNT - 1.0 ) * 100. { L =1000, 900, 800, 700, & 600 mHenry
C MILLIH = 1000. - ( KNT - 1.0 ) * 100. { L =1000, 900, 800, 700, & 600 mHenry
VALUE = 1000. - ( KNT - 1.0 ) * 100. $$ { L =1000, 900, 800, 700, & 600 mHenry
__MILLIHENRIES__ = VALUE { 1st of 2 uses of preceding intermediate variable
MILLIH = VALUE { 2nd of 2 uses of preceding intermediate variable
C Note about preceding. There are two uses of the same inductance value. The
C 2nd is 6 columns wide for use with narrow-format series R-L-C data. The 1st
C is for wide-format, which allows 16 columns. The underscore is required as
C the first character if it is not a letter. Note that the second symbol is
C contained in the first, so order is mandatory. If order were switched, an
C error message would terminate execution following a warning message:
C ++++ Notice. Symbol 1 is contained within symbol 2.
C Swap these two and try again.
C Intermediate variable VALUE was introduced 29 Dec 98. Before that, each
C of the two symbols had the right hand side of VALUE. For this trivial
C case, the difference is unimportant. But for other cases, the formula is
C much too long to be contained on a single line. So, it is broken down
C into components, with an intermediate variable used for each. Note the
C $$ on the right of the intermediate variable. By definition, such a
C variable is used within the pocket calculator only. The name is not to be
C searched for and replaced in later data cards outside of the $PARAMETER
C block (terminated by the following blank card). For industrial-strength
C use of intermediate variables, see Orlando Hevia's induction motor model
C in Rule Book documentation of $PARAMETER.
BLANK card ends $PARAMETER definitions that are processed just b4 branch cards
C Series R-L-C, narrow: RRRRRRLLLLLLCCCCCC
TRAN 10.0 3
GEN TRAN6 5.0MILLIH { $PARAMETER will define this inductance
$VINTAGE, 1, { The preceding used old, normal, narrow format. Switch to wide
C Series R-L-C, wide: RRRRRRRRRRRRRRRRLLLLLLLLLLLLLLLLCCCCCCCCCCCCCCCC
TRAN6 10.0 3
GEN TRAN 5.0__MILLIHENRIES__ { $PARAMETER will define this inductance
$VINTAGE, 0, { Return to old, normal, narrow format
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14GEN 70. 50. -1.
BLANK card terminating program source cards.
GEN TRAN TRAN6
BLANK card ending program output-variable requests.
CALCOMP PLOT { Needed to cancel PRINTER PLOT of preceding subcase
194 2. 0.0 20. -.40 .40TRAN
BLANK card ending batch-mode plot cards
-1 10.TRAN { Statistical tabulation of branch voltage (TRAN, TERRA)
BLANK card ending all statistical tabulation request cards
BEGIN NEW DATA CASE
C 6th of 18 subcases illustrates data that created the compiled TACS code
C that is built into the UTPF. As used in RUN.BAT the output will be
C CODETACS.INC and DECKTAC*.INC --- files put in the UTPF. The logic
C shown below was reproduced 9 times to produce the 60 variables mentioned
C in the January, 1998, newsletter (see TAL explanation). Below will be
C found the basic pocket calculator formulas that were used to illustrate
C how slow Dube's MODELS and TACS are for such basic mathematics.
C BENCHMARK DC-18b --- original source of formulas to be found below.
COMPILED TACS MAKE { Alternative to SET COMPTACS=MAKE used because later subcase
SINGLE STEP IF MAKE { Program will reduce TMAX to DELTAT if SET COMPTACS=MAKE
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
ABSOLUTE TACS DIMENSIONS
10 90 100 20 30 400 350 60
.00002 2.0
1 -11 0 0 1 -1
5 5 20 20 100 100 1000 1000 10000 10000
TACS STAND ALONE
99TEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99TEST2 = 1.E2 * COS ( 2.0 * 3.14159 * TEST1 / 100. )
99TEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( TEST2 ) )
99TEST5 = 2.0 * PI * TEST1 / 100.
99TEST6 = 2.0 * SIN ( TEST5 )
99TWOX = 2.E-4 * ( TEST2 * TEST2 + 2500. * TEST6 * TEST6 )
C 99TEST1 = 10. ---- This line is replacement for preceding to approximate removal of math
33TEST1 TEST2 TEST3 TEST5 TEST6 TWOX
BLANK card ending all TACS data cards
CALCOMP PLOT { Negative IPLOT = -11 will ignore the following plot cards
143 .2 0.0 2.0 TEST1 TEST2 TEST3
143 .2 0.0 2.0 TEST5 TEST6 TWOX
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 7th of 18 subcases illustrates use of pocket calculator to vary the frequency
C of a distributed-parameter transmission line within the FREQUENCY SCAN loop
C by means of TO SUPPORTING PROGRAM to provide LINE CONSTANTS rederivation.
C This use of LINE CONSTANTS for each frequency of the F-scan will not be
C seen in the .LIS file as proof of use unless /OUTPUT is added to the
C request that declares "TO SUPPORTING PROGRAM" (see illustration on comment).
C This subcase added 9 November 1999. This does work. Prior to adding the data
C here, there was an attempt to add it to end of DC-52, but that failed. It
C seems there is some conflict with LMFS (complicated). The $PUNCH statement
C was not being executed right. Ignore that problem by moving data here.
C FMIMFS--DELFFS--FMAXFS--NPD-----
C 11111111222222223333333344444444
FREQUENCY SCAN 2.0 200. 8
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETER 1 1 0
C Preceding MAXKNT is not used since FREQUENCY SCAN is involved. ATP will
C internally calculate the correct # of steps from FMINFS, FMAXFS, etc.
C Note we scan from 2 Hz to 200 Hz --- 2 decades, with 8 points/decade
C for a total of 16 plus the starting point, so 17 total. The output
C looks like exponential explosion as resonance around 300 Hz (would the
C 138 miles be 1/4 of a wavelength at 300 Hz? Probably) is approached.
C The 8 points/decade are needed to see the right edge accurately.
C 30 July 2001, preceding NOSTAT is added. This is a binary switch
C to suppress statistical tabulation ("no statistics") if value is 1.
C But for F-scan, there is no statistical tabulation, so value is
C immaterial.
$PARAMETER
C FRECU_____ = 0.01* ( EXP( LOG(10.0) /10.0 ) ) ** (KNT-1.0)
C The preceding was Tsu-huei's original, single desired line. But when first
C tried, there was an error (see April, 2000 newsletter). For about half
C a year, the following lower-order equivalent was required:
C EXPON = KNT - 1.0 $$
C BASE = EXP( LOG(10.0) / 10.0 ) $$
C FRECU_____= .01 * BASE ** EXPON
C Preceding was not nice for plotting. Switch to 8 points/decade beginning
C at 2 Hz (below which nothing much happens):
C EXPON = KNT - 1.0 $$
C The preceding line worked fine. But it was not as demanding of pocket
C calculator logic as the following. As mentioned in the newsletter, the
C leading "-" on the leading constant was mishandled prior to 13 June 2000:
C EXPON = -1.0 + KNT $$
C BASE = EXP( LOG(10.0) / 8.0 ) $$
C FRECU_____= 2.0 * BASE ** EXPON
C The preceding 3 lines were used successfully until 13 October 2007. Then,
C new logic to handle extraneous parens was added to POCKET, and it seemed
C that this also should handle the original problem. To prove it, comment
C out the preceding 3 cards and replace by the equivalent original. The .LIS
C file is esstentially unchanged. WSM, 13 Oct 2007 :
FRECU_____ = 2.0 * ( EXP( LOG(10.0) /8.0 ) ) ** (KNT-1.0)
BLANK card terminates $PARAMETER definitions
EXACT PHASOR EQUIVALENT { Switch from lumped-R to exact Pi-equiv. of distributed
C The preceding is not really necessary. It should make the representation more
C accurate, of course. At the highest frequency, if we are approaching 1/4 of
C a wavelength, there should be a difference. But that difference seems small.
C Consider the output vector at the ending frequency, with and without EPE:
C With EPE: 17 200. 303000. 0.0 371408.7751 -8.40003412 583671.5409 -30.2933187 1684.552878 49.89669008
C Without : 17 200. 303000. 0.0 371386.1316 -8.41590291 583780.5064 -30.3165151 1685.249616 49.82908871
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (if 8 output columns)
-1. -1. 60. 60. { Note dT and Tmax are ignored for FREQ SCAN
1 -1 1 0 1
SSSS1AAAAA1A 15. 1
SSSS1BBBBB1B 15.
SSSS1CCCCC1C 15.
C Add /OUTPUT anywhere on right of following request to produce output of
C the supporting program. Without this, all output will be suppressed.
C < TO SUPPORTING PROGRAM (NEXT) > /OUTPUT { Add to see LINE CONSTANTS output
< TO SUPPORTING PROGRAM (NEXT) > { Request for a jump to a supporting program
C The following data, through the blank card, is from 2nd subcase of DCNEW-3.
LINE CONSTANTS { This is the supporting program that is to be executed in-line
BRANCH AAAA1AAAAA2ABBBB1BBBBB2BCCCC1CCCCC2C
UNTRANSPOSED { Request for the line to be represented as untransposed
C LINE CONSTANTS DATA FOR JOHN DAY TO LOWER MONUMENTAL 500-KV LINE.
1.3636 .05215 4 1.602 -20.75 50. 50.
1.3636 .05215 4 1.602 -19.25 50. 50.
2.3636 .05215 4 1.602 - 0.75 77.5 77.5
2.3636 .05215 4 1.602 0.75 77.5 77.5
3.3636 .05215 4 1.602 19.25 50. 50.
3.3636 .05215 4 1.602 20.75 50. 50.
0.5 2.61 4 0.386 -12.9 98.5 98.5
0.5 2.61 4 0.386 12.9 98.5 98.5
BLANK CARD ENDING CONDUCTOR CARDS OF "LINE CONSTANTS" CASE
100. FRECU_____ 1 11 1 138. 1
BLANK CARD ENDING FREQUENCY CARDS OF "LINE CONSTANTS" CASE
$PUNCH, sup.pch ! { Output the just-created branch cards to a disk file
BLANK card ending LINE CONSTANTS data cases { Last of supporting program data
C End of data for supporting program. Branch data input now resumes:
$CLOSE, UNIT=7 STATUS=KEEP { Disconnect file containing KC Lee branch cards
$OPEN, UNIT=7 STATUS=UNKNOWN FILE=dumpch.lis FORM=FORMATTED
$INSERT, sup.pch, { Dynamic connection of disk file created by $PUNCH above
AAAA2A 1000. { Resistive load at far end of 1st of 3 phases
BBBB2B AAAA2A { 2nd of 3 phases is terminated by the same R
CCCC2C AAAA2A
BLANK card ending branch cards
BLANK card ending switch cards
POLAR OUTPUT VARIABLES { 1st of 3 alternatives for output gives mag, angle
C Preceding is one of 3 alternatives. The other two are, after commenting:
C BOTH POLAR AND RECTANGULAR { Request for (in order): mag, angle, real, imag
C RECTANGULAR OUTPUT VARIABLES { 3rd of 3 alternative outputs gives real, imag
14SSSS1A 303000.FRECU_____ 0. -1.
14SSSS1B 303000.FRECU_____ -120. -1.
14SSSS1C 303000.FRECU_____ 120. -1.
BLANK card ending source cards
SSSS1AAAAA1AAAAA2A { Names of 3 nodes (all a-phase, note) for V-node output
C Column headings for the 4 output variables follow. These are divided among the 3 possible FS variable classes as follows ....
C First 3 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 For each variable, magnitude is followed immediately by angle. Both halves of the pair are labeled identically, note.
C Step F [Hz] SSSS1A SSSS1A AAAA1A AAAA1A AAAA2A AAAA2A SSSS1A SSSS1A
C AAAA1A AAAA1A
C 1 2.0 303000. 0.0 303004.5676 -.028517799 301911.4209 -.185323276 301.7646507 1.720504686
C New parameter values follow: 1) 2.6670429
C 2 2.6670429 303000. 0.0 303008.121 -.038018478 301917.2789 -.247632726 301.7897693 2.294000532
C New parameter values follow: 1) 3.5565588
C 3 3.5565588 303000. 0.0 303014.4401 -.050681287 301931.0223 -.330715516 301.8820814 3.058524978
BLANK card ending output variable requests (here, only node voltages)
C New parameter values follow: 1) 84.339301
C 14 84.339301 303000. 0.0 311692.5772 -1.38233484 333737.5916 -8.33491659 541.8603856 48.84191334
C New parameter values follow: 1) 112.46827
C 15 112.46827 303000. 0.0 319327.0284 -2.11141826 362848.3744 -11.7652408 709.4891996 53.86022223
C New parameter values follow: 1) 149.97884
C 16 149.97884 303000. 0.0 335125.0318 -3.6799135 425532.3722 -17.6205065 1015.840142 55.61763042
C New parameter values follow: 1) 200.
C 17 200. 303000. 0.0 371408.7751 -8.40003412 583671.5409 -30.2933187 1684.552878 49.89669008
F-SCAN COMPONENTS MAG MAG { Access "mag" and "mag" next (only 1st used)
197.25 0.0 2.5 0.02000SSSS1AAAAA1A Magnitude of i-aAmperes
F-SCAN COMPONENTS ANGLE ANGLE { Access "angle" and "angle" next
147.25 0.0 2.5 -40. 0.0AAAA1AAAAA2A Angles of v-nodeDegrees
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 8th of 18 subcases illustrates two separate $PARAMETER blocks within a
C PCVP (POCKET CALCULATOR VARIES PARAMETERS) loop over time simulation.
C The data comes from Prof. Juan Martinez as should be explained in the
C April, 2001, newsletter. Prior to correction on 19 December 2000, symbol
C RESIST was not being correctly handled on the 2nd pass. It should have
C value 5 ohms on the first pass and 15 ohms on the second pass. Note that
C it is in the 1st of two blocks. Prior to correction, memory of this was
C being overlaid by memory for the 2nd block. Now, there is no overlaying.
C 25 December 2000, add variable HOLD01 to verify the new non-volatile
C storage (any variable HOLDxx that begins with the 4 letters HOLD).
C Gabor Furst asked for this extension as should be explained in the
C April, 2001 (or newer), newsletter. HOLD01 begins with value zero,
C so has value 1 on first pass and 3 on 2nd pass of loop over KNT. The
C variable is not really being used for anything in this illustration,
C but it certainly could be. Except for its non-volatile nature, it is
C a normal variable in every other respect.
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 2 0 1
.000100 .030
20 -1
$PARAMETER
RESIST = 5.0 + (KNT-1.0)*10.0
HOLD01 = HOLD01 + KNT { Non-volatile variable because 1st 4 letters "HOLD"
BLANK card ending 1st of 2 $PARAMETER blocks
$PARAMETER
AB_TON = 0.01
ABTOFF = 0.01 + ( KNT - 1.0 ) * 0.01
GND__K = KNT SERIALIZE 'GND___'
BLANK card ending 2nd of 2 $PARAMETER blocks
GND__2 100. 0
GND__1 RESIST 0
SOURCAMIDD_A RESIST 0
BLANK card ending branch cards
MIDD_AGND__K AB_TON ABTOFF
BLANK card ending switch cards
14SOURCA 0 1.5E5 50. -1. 1.
BLANK card ending source cards
MIDD_A
BLANK card ending names of nodes for voltage output
BLANK card ending batch-mode plot cards
BLANK card ending all statistical tabulation request cards
BEGIN NEW DATA CASE
C 9th of 18 subcases is copied from 5th, but MODELS data has been added
C to process and display extrema. This service begins 1 January 2001.
C Optional col. 32: MEXTR 0 or 1 ==> minimum extrema; 2 ==) all extrema
MODELS PROCESSES EXTREMA 0 { Use MODELS in simulation but only for extrema
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 5 1 1 { Loop five times
PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns)
.0001 .020
1 -1 1 1 0 -1
5 5 20 20
C Note integer miscellaneous data card has extrema flag MAXOUT = 0 as proof
C that the user does not need to remember to set it. PCVP automatically will
C set this to unity. After that, MEXTR will suppress such output for passes
C 2 onward unless MEXTR = 2 (which will show all extrema). MODELS output
C will always be seen, although note that such output will disappear if all
C write( statements are removed from the MODELS data. In the following
C data, INPUT variables are special. In place of the usual variable types,
C one of 4 new ones must be used: Xmax, Tmax, Xmin, and Tmin corresponding to
C the 4 rows of extrema output. The name within parentheses ("dummy") is not
C used, so can be any character string beginning with a letter. It can not
C be omitted, however (unfortunately). Next comes the fixed in-line comment,
C which can begin in any column. But, to the right of the colon, format is
C fixed. It is assumed I2, 1X, 2A6 information is involved. The integer is
C the same one used on a batch-mode plot card: 4 if node voltage, 8 if branch
C voltage, and 9 if branch current (or an S.M., a U.M., or a TACS variable).
C The two names that follow correspond to the printed heading. One minor
C enhancement is use of symbolic BLANK, which ATP interprets as real blanks.
C Each INPUT variable must be on its own line (see the 4 lines just below).
C 2 January 2001, WRITE( statements are left-adjusted to column 1 to take
C advantage of new logic that will hold the user's case within text (as
C delineated by apostrophes).
C 3 January 2001, add an initial $PARAMETER block that has no use. This is
C inspirated by Prof. Martinez's data TESTA1.DAT which inspired correction
C to allow such practice.
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be
UNUSED = 5.0 { Following data does not involved character string "unused".
BLANK card ends $PARAMETER definitions that are processed before MODELS data
MODELS
INPUT tran1 {Xmax(dummy)} -- Plot type and variable names: 8 TRAN6 BLANK
INPUT tran2 {Tmax(dummy)} -- Plot type and variable names: 8 TRAN6 BLANK
INPUT tran3 {Xmin(dummy)} -- Plot type and variable names: 8 TRAN6 BLANK
INPUT tran4 {Tmin(dummy)} -- Plot type and variable names: 8 TRAN6 BLANK
MODEL m1 -- Begin definition of a model named m1
VAR pass, sum -- Verify that values are kept between passes
INPUT row1 {dflt: 0} -- Connects with Xmax, the 1st of 4 extrema rows
INPUT row2 {dflt: 0} -- Connects with Tmax, the 2nd of 4 extrema rows
INPUT row3 {dflt: 0} -- Connects with Xmin, the 3rd of 4 extrema rows
INPUT row4 {dflt: 0} -- Connects with Tmin, the 4th of 4 extrema rows
INIT pass:=0 -- Begin with the pass number equal to zero
sum:= 0 ENDINIT -- Begin with the sum of Xmax equal to zero
EXEC -- For each new set of extrema, we want to do this:
pass:=pass+1 -- Count the pass number of the loop over simulations
sum := sum + row1 -- Sum all Xmax as trivial demonstration of math
write(' ************ In MODELS, pass = ', pass, 'sum =', sum )
write( ' Four rows Xmax, Tmax, Xmin, Tmin = ', row1, row2, row3, row4 )
write( ' ' ) -- Blank line separates this pass from following one
C 3 March 2001, to satisfy need of Prof. Juan Martinez in Barcelona, add the
C following copy of the preceding 3 lines. This produces a 2nd copy of the
C preceding output in disk file MODELS.1 that is connected to I/O unit 37
C as used by $CLOSE, $OS, and $OPEN commands immediately before plotting.
write1(' ************ In MODELS, pass = ', pass, 'sum =', sum )
write1( ' Four rows Xmax, Tmax, Xmin, Tmin = ', row1, row2, row3, row4 )
write1( ' ' ) -- Blank line separates this pass from following one
ENDEXEC -- End of operations for each new set of extrema
ENDMODEL -- End of the model named m1
USE m1 as m1
INPUT row1 := tran1
INPUT row2 := tran2
INPUT row3 := tran3
INPUT row4 := tran4
ENDUSE
ENDMODELS
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be
C __MILLIHENRIES__ = 1000. - ( KNT - 1.0 ) * 100. { L =1000, 900, 800, 700, & 600 mHenry
C MILLIH = 1000. - ( KNT - 1.0 ) * 100. { L =1000, 900, 800, 700, & 600 mHenry
VALUE = 1000. - ( KNT - 1.0 ) * 100. $$ { L =1000, 900, 800, 700, & 600 mHenry
__MILLIHENRIES__ = VALUE { 1st of 2 uses of preceding intermediate variable
MILLIH = VALUE { 2nd of 2 uses of preceding intermediate variable
BLANK card ends $PARAMETER definitions that are processed just b4 branch cards
C Series R-L-C, narrow: RRRRRRLLLLLLCCCCCC
TRAN 10.0 3
GEN TRAN6 5.0MILLIH { $PARAMETER will define this inductance
$VINTAGE, 1, { The preceding used old, normal, narrow format. Switch to wide
C Series R-L-C, wide: RRRRRRRRRRRRRRRRLLLLLLLLLLLLLLLLCCCCCCCCCCCCCCCC
TRAN6 10.0 3
GEN TRAN 5.0__MILLIHENRIES__ { $PARAMETER will define this inductance
$VINTAGE, 0, { Return to old, normal, narrow format
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14GEN 70. 50. -1.
BLANK card terminating program source cards.
GEN TRAN TRAN6
BLANK card ending program output-variable requests.
C 3 March 2001, to satisfy need of Prof. Juan Martinez in Barcelona, add the
C following 6 $-cards. Together with preceding WRITE1 of MODELS, this
C demonstrates the needed ability for MODELS to create a disk file that will
C be disconnected and renamed at the end of a PCVP pass. When execution is
C complete, TEMPFILE.LIS will contain the 5th and last of 5 successive
C contents, of course. Result of the TYPE command will not be detected in
C the .LIS file, unfortunately, but it can be seen on the screen during
C execution, if there is any question about content of the preceding 4. Of
C course, the WRITE1 statements are the same as WRITE, so in fact the
C output of TYPE is indirectly documented in the .LIS file.
$CLOSE, UNIT=37 STATUS=KEEP { Disconnect file connected to MODELS WRITE1(
$OS, DEL TEMPFILE.LIS { Make sure this file does not exist so RENAME is legal
$OS, RENAME MODELS.1 TEMPFILE.LIS { Save MODELS WRITE1( file to a new name
$OS, TYPE TEMPFILE.LIS { Display on screen the MODELS WRITE1( file contents
C $OS, PAUSE { For interactive use on a fast computer, allow time to inspect
C the output produced by the preceding TYPE command.
$OPEN, UNIT=37 FILE=models.1 FORM=FORMATTED STATUS=NEW ! { Connect new file for WRITE1(
BLANK card ending batch-mode plot cards
C The following request for statistical tabulation should work fine. But it
C was demonstrated in the 5th subcase, so why repeat the burden here?
C -1 10.TRAN { Statistical tabulation of branch voltage (TRAN, TERRA)
C 30 July 2001, add preceding NOSTAT = 1. Since there is not going to be
C any statistical tabulation, why burden associated storage unnecessarily?
C Because NOSTAT = 1, in fact the following card will not be read.
BLANK card ending all statistical tabulation request cards
BEGIN NEW DATA CASE
C 10th of 18 subcases is related to 8th. It illustrates the use of
C character string substitution as well as numeric substitution within
C an IF block. This need, too, comes from Prof. Martinez. Since the
C network is purely resistive, there are no real dynamics other than
C the switching. Note the first pass involves the switch to NAME2
C whereas the 2nd involves the switch to NAME1. On the first of 2
C passes, RESIST = 5 whereas on the second RESIST = 10. Increase
C step size from 100 usec to 2 msec to speed simulation. Normally,
C node names of character strings will be a full 6 bytes long. But
C just to illustrate that shorter names are handled correctly, 5-byte
C NAMEK is used in this subcase. WSM, 6 Jan 2001.
C 26 January 2001, switch the order of the two parameters within the
C IF block. Of course, the answer should not be affected. But that
C error first reported by Prof. Juan Martinez would mistreat data of
C which the first assignment was character rather than numeric. Switch
C the order to demonstrate that the problem has been cured.
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 2 0 1
.002 .030
1 -1
$PARAMETER { Begin block that defines numeric RESIST and character NAMEK
C IF( KNT .EQ. 1.0 ) THEN This standard form involves THEN on right, note.
IF( KNT .EQ. 1.0 ) { This is equivalent (THEN on right is not required)
NAMEK = 'NAME2' { Beginning 26 January 2001, character line appears first
RESIST = 5.0 + ( KNT - 1.0 ) * 10.0 { Previously, this numeric line was first
ELSE { Alternatively (if not the first pass, so for KNT = 2 onward):
NAMEK = 'NAME1'
RESIST = ( KNT - 1.0 ) * 10.0
ENDIF { Terminate 7-line block that defines symbols RESIST and NAMEK
BLANK card ending 1st of 2 $PARAMETER blocks
NAME2 100.
NAME1 RESIST
SOURC MID RESIST
BLANK card ending program branch cards.
MID NAMEK 0.009 0.020
BLANK card terminating program switch cards
14SOURC 0 100. 50. -1.
BLANK card terminating program source cards
MID
BLANK card ending program output-variable requests.
BLANK card ending batch-mode plot cards
BLANK card ending all statistical tabulation request cards
BEGIN NEW DATA CASE
C 11th of 18 subcases illustrates a PCVP loop over time simulations, with
C MODELS used as a part of the simulation. This is unlike the preceding
C 9th subcase, which involved MODELS PROCESSES EXTREMA. The 9th subcase
C uses MODELS not within the dT loop, but rather only upon completion of
C the simulation. The present use within the dT loop became operational
C 17 May 2001 when SUBR1, TACS1, and TACS2 were corrected to properly
C handle MAGVOLT.DAT from Orlando Hevia in Santa Fe, Argentina.
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 3 1 0 { 3 passes, minim print
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be
VOLTAG = 25.0 $$
REACTV = 0.5 + 0.5*KNT - 0.5 $$
ADMITC = (REACTV/VOLTAG**2)*1000000.
POWERR = REACTV
BLANK card ends $PARAMETER definitions that are processed before MODELS data
C 2.E-5 40.E-3 60.0 60.0
C 100 1 1 1 1
C Preceding miscellaneous data cards of MAGVOLT are modified to simplify and
C speed this illustration. The PCVP card involves 3 passes whereas MAGVOLT
C had 61. Reduce the number of time steps from 2000 to 50, and do not save
C the .PL4 files (ICAT = 1 in column 64 creates .001, .002, etc. files).
.0002 .010 60.0 60.0
1 1 1 1 0 -1 0
5 5
MODELS HYBRID { Begin MODELS data that is to participate in the simulation
INPUT voltage {v(bus2a)}
node0 {v(busx)}
tmax {atp(tmax)}
C Comment added 18 September 2001. The preceding line illustrates the ATP(
C function of MODELS. In this case, it is used to read the value of variable
C TMAX, which is the ending time of the study. Beware of loop counter KNT
C of PCVP use. Inside the $PARAMETER block, KNT is used, note. But if
C accessed using the ATP( function, a bad value will be obtained. This is
C because KNT is not the name of the variable used to count the PCVP loops
C within the program. Rather, MNT is. So, what the user might guess would be
C ATP(KNT) instead must be ATP(MNT). The need for this warning was inspired
C by a report of trouble by Marta Val Escudero of ESB International in Dublin,
C Ireland. As should be reported in a future newsletter, "PCVP loop index KNT
C conflicted with the energization number of STATISTICS prior to separation of
C the two as summarized in the July, 2000, issue."
MODEL maxim
INPUT voltage parame tmax
VAR maxim kvoltage
INIT
maxim:=0
ENDINIT
EXEC
kvoltage:=ABS(voltage)
IF (kvoltage > maxim) THEN maxim:=kvoltage
ENDIF
IF (t > tmax) THEN
C Only on the final time step, one line is written to the
C MODELS.1 file. This is a little inconvenient for test cases.
C We would like to redirect this output to the .LIS file by changing
C WRITE1 to WRITE. But, since we are minimizing out (see the PCVP
C declaration, output within the dT loop of the 2nd or later pass is
C being suppressed, and only a single such output (for pass 1) would
C be seen. So, we leave WRITE1. For the record, content should be:
C DCN25, 11th subcase. PARAME = 0.5 MAXIM = 1.04828954
C DCN25, 11th subcase. PARAME = 1.0 MAXIM = 1.22325171
C DCN25, 11th subcase. PARAME = 1.5 MAXIM = 1.41404906
C Note we have added labeling and preserved lower case
C by left-adjusting the WRITE1 command:
C write1(parame, ' ', maxim)
write1( ' DCN25, 11th subcase. PARAME = ', parame,
'MAXIM = ', maxim )
ENDIF
ENDEXEC
ENDMODEL
USE maxim AS sensit
INPUT voltage:=voltage parame:=node0 tmax:=tmax
ENDUSE
ENDMODELS
BLANK card ending MODELS data. Next come branch data:
BUSX POWERR
BUS0A BUS1A .4948 1.979
BUS1A BUS2A 5.562 22.25
BUS2A 800.0
BUS3A ADMITC
BLANK card ending program branch cards.
C BUS1A BUS3A 5.E-3 1.0 Switch card of MAGVOLT used multiple of dT
BUS1A BUS3A 5.001E-3 1.0 { Perturb T-close so problem with roundoff
BLANK card terminating program switch cards (none, for this case)
11BUSX -1 1.0 0.0 1.0
14BUS0A 1.0 50.0 -90.0 -1.0 1.0
BLANK card terminating program source cards.
BUS0A BUS1A BUS2A BUS3A
BLANK card ending program output-variable requests.
BLANK card ending plot cards (none for this data)
BLANK card ending statistical tabulations (none illustrated here)
BEGIN NEW DATA CASE
C 12th of 18 subcases illustrates lack of difference in the first 6 bytes
C of two variable names of $PARAMETER usage. This originally was disk
C file DISABLE0.DAT from Prof. Juan Martinez in Barcelona, Spain. It
C was sent to BPA on 14 September 2001. Because BLOCK_A1 & BLOCK_A2
C are the same through byte 6, execution fails. The first sign of trouble
C is incorrect evaluation of BLOCK_A2 (note BLOCK_A1 is correct):
C Parameter 2 defined. Value = "$ENABLE " ----- result of evaluation of the preceding IF-THEN-ELSE-ENDIF block.
C Parameter 3 defined. Value = 1.500000E+00 ----- result of evaluation of the preceding IF-THEN-ELSE-ENDIF block.
C Of course, parameter 3 should be character just as parameter 2 is. It
C as not numeric as this indicates. As a result of this error, later data
C is bad, and execution should end with "KILL = 4. The last card is a
C series R-L-C branch having zero impedance. ..."
$PARAMETER
OPTI_1=1.0
BLANK
$PARAMETER
IF(OPTI_1.EQ.1.0)
BLOCK_A1='$ENABLE '
BLOCK_A2='C '
ENDIF
BLANK
.1e-6 .2e-4
50 1
BLOCK_A1
BUS_10BUS__1 1.0 0
BLOCK_A2
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14BUS_10 1.5E5 50. -1. 1.
BLANK card terminating program source cards.
BUS__1
BLANK card ending program output-variable requests.
BLANK card ending plot cards (none for this data)
BEGIN NEW DATA CASE
C 13th of 18 subcases is a correction of the preceding subcase. The only
C change is replacement of BLOCK_A2 by A2_BLOCK (2 places). Then the
C two symbols are different through the first 6 bytes. The subcase now
C is handled properly. Note that it represents a variation of the 10th
C subcase, which had a note dated 26 January 2001. If the 10th subcase
C no longer ends its $PARAMETER block with a character variable, this
C present subcase does. Also, it illustrates the use of 2 of them within
C a single block.
$PARAMETER
OPTI_1=1.0
BLANK
$PARAMETER
IF(OPTI_1.EQ.1.0)
BLOCK_A1='$ENABLE '
A2_BLOCK='C '
ENDIF
BLANK
.1e-6 .2e-4
50 1
BLOCK_A1
BUS_10BUS__1 1.0 0
A2_BLOCK
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14BUS_10 1.5E5 50. -1. 1.
BLANK card terminating program source cards.
BUS__1
BLANK card ending program output-variable requests.
BLANK card ending plot cards (none for this data)
BEGIN NEW DATA CASE
C 14th of 18 subcases illustrates dynamic as opposed to static $DISABLE
C use. This originally was disk file DISABLE3.DAT from Prof. Juan
C Martinez in Barcelona, Spain. It was sent to BPA on 14 September 2001.
C Prior to a change to CIMAGE on 17 Sept 01, execution failed.
$PARAMETER
OPTI_1=1.0
BLANK
$PARAMETER
C IF(OPTI_1.EQ.1.0)
BLOCK_A1='$ENABLE '
BLOCK_A2='C '
C ENDIF
BLANK card ending $PARAMETER block
$PARAMETER
C IF(OPTI_1.EQ.1.0)
BLOCK_B1='$DISABLE'
BLOCK_B2='$ENABLE '
C ENDIF
BLANK card ending $PARAMETER block
.1e-6 .2e-4
50 1
BLOCK_A1
BUS_10BUS__1 1.0 0
BLOCK_A2
BLOCK_B1
BUS_10BUS__1 1.0 0
BLOCK_B2
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14BUS_10 1.5E5 50. -1. 1.
BLANK card terminating program source cards.
BUS__1
BLANK card ending program output-variable requests.
BLANK card ending plot cards (none for this data)
BEGIN NEW DATA CASE
C 15th of 18 subcases shows that a FREQUENCY SCAN with multi-part output
C is possible with no output voltage. This addition is made 16 June 2002
C to demonstrate that a problem, first observed by Orlando Hevia, has been
C corrected. The change to HEADL4 was simple: (N2 .EQ. NV) was replaced
C by (N3 .EQ. NC). Before this correction, execution died in plotting as
C ATP attempted to read too much of the .PL4 file. Now, everything is OK,
C as the following shows. The plot illustrates a mixture of magnitude and
C real part on the same plot. At the peak, which occurs at 160 Hz, both
C have the same value, of course. 160 Hz is the resonant frequency of the
C series R-L-C branch. The variation is nice and smooth and simple using
C 10 Hz steps. Here, using 20 Hz, the curves are slightly lumpy.
PRINTED NUMBER WIDTH, 10, 2, { Each column of width 10 includes 2 blank bytes
FREQUENCY SCAN 20. 20. 200. { F in Hz = 20, 40, ... 200
1.E-6 .001
1 -1 1 1
XX 1. 1.0 1.E3 1
BLANK card ending branch cards (just this one)
BLANK card ending switch cards (none)
BOTH POLAR AND RECTANGULAR { Any of the 3 possible multi-part outputs works
14XX 100. 60. -1.
BLANK card ends source cards
C XX { Prior to correction of HEADL4 on 16 June 2002, node voltage was needed
BLANK card ending names of nodes for node-voltage output
F-SCAN COMPONENTS MAG REAL { Access magnitude and real parts next
19620. 0. 200. -40.120.XX XX { Plot the magnitude and the real part
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 16th of 18 subcases is added 18 February 2003 to demonstrate that the
C complaint by Prof. Juan Martinez has been corrected. The preceding day,
C he had sent REPLACE.DAT to demonstrate that a $PARAMETER block within
C a dynamic $DISABLE block was being mishandled. The 3rd pass, execution
C was being stopped by an error message within TACS (which is not even
C being used).
C Following request carries params MAXKNT IOPCVP NOSTAT
POCKET CALCULATOR VARIES PARAMETERS 3 0
$PARAMETER
C IF(MNT.GT.1.0) -- Original line until 16 April 2009 involved erroneous MNT
IF( KNT .GT. 1.0 ) { Correct MNT --> KNT to agree with all other DC*.DAT
KNT1____='$DISABLE'
ELSE
KNT1____='$ENABLE '
ENDIF
BLANK card ends 1st $PARAM block
.0002 .02
50 -1
$PARAMETER { SIMPLE PARAMETER DEFINITION, JUST TO ILLUSTRATE THE CASE
PARAM2=1.0
BLANK card ends 2nd $PARAM block
KNT1____ { This will be replaced by $ENABLE on 1st pass, & $DISABLE thereafter
$PARAMETER { THIS BLOCK IS CAUSING CERTAIN DAMAGE WITHIN PCVP ROUTINE
PARAM3=100.0
BLANK card ends 3rd $PARAM block
VOLTSA PARAM3
$ENABLE { End of either $DISABLE or $ENABLE block that begins with KNT1____
$PARAMETER { THIS COULD BE ANY MATHEMATICAL COMPUTATION TO OBTAIN THE CONDI
TIPO_F=MNT
BLANK card ends 4th $PARAM block
$PARAMETER { $PARAMETER BLOCK TO COMMENT OR UNCOMMENT NEXT BRANCH CARD
IF(TIPO_F.LT.1.0)
OPTION=' VOLT'
ELSE
OPTION='C VOLT'
ENDIF
BLANK card ends 5th $PARAM block
OPTIONSA PARAM2 { BRANCH CARD RELATED TO PREVIOUS PARAMETER
VOLTSA 1.0
BLANK card ends branch cards
VOLT_AVOLTSA MEASURING 1
BLANK card ends switch cards
14VOLT_A 0 1. 50. -1. 1.
BLANK card ends source cards
VOLTSA
BLANK card ends names for node voltage output
BLANK card ends batch-mode plot cards
BLANK card ends requests for statistical tabulation
BEGIN NEW DATA CASE
C 17th of 18 subcases is added 11 May 2003. This data comes from Ricardo
C Tenorio of ABB in Vasteras, Sweden. The 2 symbol names are not distinct
C through 6 bytes. Unlike the 12th subcase, ATP did not automatically
C shut itself off with an unrelated error message. Mr. Tenorio found that
C execution continued, and results were wrong. In E-mail dated May 5th,
C he explained: I've tried to simplify the case as far as I could. The
C case involves an EDM module designed to apply generic faults, i.e 1-ph,
C 2-ph & 3-ph, to a named bus." Data interpretation provided the first
C sign of trouble: "Only the first variable BLOCk_11 is correctly
C interpreted as a string! The others are interpreted as number ...
C Line 22 ends this IF-THEN-ELSE-ENDIF block. |ENDIF
C Parameter 1 defined. Value = "$DISABLE" ----- result of evaluation ...
C Parameter 2 defined. Value = 0.000000E+00 ----- result of evaluation ...
C Yes, clearly wrong. Most serious was continued execution. So, a special
C error trap was designed for MATDAT. Beginning May 11th, it halts this
C case as well as the 12th subcase (which now is stopped sooner, and with
C a relevant message for the user). What is needed (in addition to the
C ambiguous names)? An IF( block and text (as opposed to just math).
C This uses the pocket calculator, which in turn requires distinct names.
C The first 6 bytes of each is BLOCK_ so data must be rejected.
C EATS testing of RUNEATS.BAT required a modification. Immediately
C below the $PARAMETER block will be found the two uses. Note that
C these two cards have been commented out. The problem of EATS is
C simple: switch cards must be read for counting, yet $PARAMETER is
C skipped during counting. So, if not commented out, execution would
C die in OVER4 with an illegal decode (BLOCK... is not a switch card).
C So, we modify to avoid death at the hands of EATS. For normal (non-
C EATS) use, the change makes no difference. Execution will be halted
C before any data below the $PARAMETER loop is read. It makes no
C difference what data is down there.
50.E-6 10.E-3 50. 50.
1 -1 1 1 1 -1
5 5 20 20
51SOUR1AEQUIVA 1.003 9.5900
52SOUR1BEQUIVB 0.319 6.9810
53SOUR1CEQUIVC
BLANK card ends branch data
$PARAMETER
IF( 3.0 .EQ. 1.0 )
BLOCK_11='C '
BLOCK_12='C '
ELSEIF( 3.0 .EQ. 3.0 )
BLOCK_11='$DISABLE'
BLOCK_12='C '
ENDIF
BLANK card ends $PARAMETER block
C BLOCK_11 { This symbol is to be replaced by either $DISABLE or a comment card
C Note: the preceding card was commented to prevent death in OVER4 using EATS
EQUIVA 0.0100 9999.
C BLOCK_12 { This symbol is to be replaced by a comment card (either condition)
C Note: the preceding card was commented to prevent death in OVER4 using EATS
EQUIVAEQUIVB 0.0100 9999.
EQUIVA 0.0100 9999.
EQUIVB 0.0100 9999.
EQUIVC 0.0100 9999.
$ENABLE { End possible preceding $DISABLE after last switch card
BLANK card ends switch data
14SOUR1A 428660.705 60. 0.00 -1. 1.
14SOUR1B 428660.705 60. -120.00 -1. 1.
14SOUR1C 428660.705 60. 120.00 -1. 1.
BLANK card ending source cards
C Since there will be an error message, remaining data is immaterial; it
C will not be used, anyway. So let's avoid waste by eliminating remaining
C data through the end of the subcase.
BEGIN NEW DATA CASE
C 18th of 18 subcases is added 8 October 2006. This data involves a PCVP
C loop that produces 5 time-simulations with appended power and energy output
C due to column-80 punches > 4. Original complaint was from Orlando Hevia of
C Santa Fe, Argentina on 18 Aug 2006. Prior to correction in October of 2006,
C column-80 punches above "4" were unusable with "FIND" code of statistical
C tabulations. Punches 5-16 (use hex if above 9) serve to append power P or
C energy "E" or both. The printed dT-loop output was believed to be correct.
C But the subsequent "FIND" command offered no class for power & no class for
C energy. Now it does. As for UTPF changes, look for WSM06AUG since the
C work began in August. For uniformity, no SEP or OCT were used. The
C vast majority of changes are confined to DICTAB. WSM, 6 October 2006
C Controls of following request card: MAXKNT IOPCVP NOSTAT { Loop 5 times with
POCKET CALCULATOR VARIES PARAMETERS 5 1 0 { reduced printout
$PARAMETER { Each loop KNT=1, 2, ... MAXKNT will involve different TIMECLOSEX
TIMECLOSEX = (KNT)*0.0004
C 13 October 2007, the preceding line is modified by the addition of
C extraneous parentheses around the KNT. Orlando Hevia had sent data to
C show that RTR6__= 0.684+0.054*(KNT)**1.75 is wrong (Salford ATP dies)
C whereas RTR6__= 0.684+0.054*KNT**1.75 works fine. This resulted
C in a correction to POCKET. To build verification into a test cases, WSM
C now simply adds unused parens () to this 18th subcase. Answer is unchanged.
C However, this did _not_ stress the new logic as intended. Well, keep the
C preceding addition but add a second, unused variable by copying Orlando :
RTR6__= 0.684+0.054*(KNT)**1.75 { Orlando's formula that was fatal b4 Oct. 2007
BLANK card ending parameter block
C DELTAT TMAX
1.E-4 .200 { Orlando Hevia had had 20 usec and .25 sec here
C IOUT IPLOT IDOUBL KSSOUT MAXOUT
200 11 1
C Replace Orlando's J. Marti line by 42 constant-parameter miles from DC-38:
-1FASEA1FASEA2 .305515.8187.01210 42. 0 { 42-mile, constant-
-2FASEB1FASEB2 .031991.5559.01937 42. 0 { parameter, 3-phase
-3FASEC1FASEC2 { transmission line.
C For list of all column-80 punches for power or energy, see DC-37, subcase 6.
C In the following, "4" is old whereas "G" and "D" date to 1999. Outputs
C will be as shown farther below:
FASEA2 1.0E01 G
FASEB2 1.0E01 D
FASEC2 1.0E01 4
GENEA OPENXX RTR6__ { Do-nothing branch will confirm value of RTR6
BLANK card ending branches
GENEA FASEA1TIMECLOSEX 0.05 3
GENEB FASEB1TIMECLOSEX 0.05 0
GENEC FASEC1TIMECLOSEX 0.05 0
BLANK card ending switches
14GENEA 100000. 50. 0.0 -1.0
14GENEB 100000. 50. -120. -1.0
14GENEC 100000. 50. 120. -1.0
BLANK card ending sources
FASEA1
C Column headings for the 9 EMTP output variables follow. These are divided among the 5 possible classes as follows ....
C First 4 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 Next 2 output variables are either power or energy or both as a pair (column 80 punches > 4).
C Remember that 1 of the voltages really are powers, and currents are energies; the 4-punches will overlay.
C Step Time FASEB2 FASEC2 GENEA FASEA1 GENEA FASEB2 FASEC2 FASEA2 FASEB2
C TERRA TERRA FASEA1 FASEA1 TERRA TERRA TERRA TERRA
C End actual dT-loop output. For clarity, let's label each of the preceding
C output variables, showing the variable class and how it was requested:
C Variable class bran V power bran V node V bran I bran I energy energy energy
C Produced by: "D" "4" "3" A6 name "3" "D" "4" "G" "D"
C Summary class <----- 4 voltages appear first -----> <--- 3 currents next ---> <-- appended P,E
C Note that power and energy of the 4-punch are counted among voltages and
C currents whereas any due to punches 4-16 are appended on the right in a
C separate, final variable class that here has 2 entries. 9 = 4 + 3 + 2
BLANK card ending node voltage outputs
C The user who does not want to see any batch-mode plotting can preserve the
C following 3 lines (no need to remove) while setting misc. data IPLOT < 0.
PRINTER PLOT { Character mode always works. CALCOMP PLOT might not (unknown)
193.01 0.0 .06 BRANCH { Show the 3 energy signals on a printer plot
FASEA2 FASEB2 FASEC2
BLANK blank card ending plot cards
C BASE-------BUS1--BUS2-- ...
-4 1.E6 FASEA2 { Compartments 93-95 will be used for base energy = 1.E6
-4 1.E6 FASEB2 { Compartments 103-105 will be used for ...
-4 2.E6 FASEC2 { Compartments 44-49 will be used for base energy = 2.E6
-3 2.E8 FASEC2 { Compartments 18-24 will be used for base power = 2.E8
C Note about preceding. Three of the 4 tabulations lost a compartment when
C the number of shots was reduced from 10 to 5. If the user wants to see a
C smoother tabulation, he can increase shots to 10 or even 20. Of course,
C TIMECLOSEX = KNT*0.0004 then would have its final factor correspondingly
C reduced to .0002 or .0001, respectively. Orlando used 20.
-1 1.E4 FASEB2 { Compartments 101-104 will be used for base voltage = 1.E4
-2 1.E3 GENEA FASEA1 { Boxes 90-97 will be used for base current = 1.E3
FIND { Search all extrema vectors to find the worst of some set of variables
ALL { Rather than list variables separately, search a set of "all something"
6 { In response to "all what?" this answers: "all energy" outputs
EXCLUDE { Repeat the search after excluding the worst shot, which is shot # 4
DISK /LIST { Create deterministic data file for last excluded shot, which is # 4
C Note about preceding. If no EXCLUDE, shot #4 still would be punched. If
C EXCLUDE use, it is data for the worst excluded shot that will be created.
BLANK card ending statistical tabulations
BEGIN NEW DATA CASE
BLANK
|