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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-24
C Monte Carlo simulation that is run in two parts. The original six
C energizations of "M39." stand-alone case have been reduced to 3, &
C the remaining 3 are performed in companion case DC-40, which is a
C "START AGAIN" continuation of this MEMSAV = 1 case. The use of
C standard random numbers and fixed seed gives results that should be
C repeatable & same for all computers of interest. But M39 differs!
C Note that this data case has no base case. However, if the user
C wants to add this, he need merely comment out the "OMIT BASE CASE"
C request card, and do likewise for the $DISABLE/ENABLE cards that
C precede statistics tabulation requests. Save the resulting binary
C disk file TPTABLES.BIN under the name dc24at40.bin for use by
C DC-40. For a "STATISTICS" random-opening data case, see DC-48.
C Beware of PARALLEL.LIS entry for same USERID that might modify NENERG!
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
C TRULY RANDOM NUMBERS { Same as NSEED=0 (this overrides stat misc. data card)
REPEATABLE RANDOM NUMBERS { Same as NSEED=1 (this overrides zero value below)
C CENTRAL STATISTICS FILE { Request word that precedes the disk file name
C %/wsm/parallel.atp ! { Leading percent sign avoids confusion of "/" in column 1
C USER IDENTIFICATION Warsaw ! { Overrides A6 USERID of STARTUP for plots
C DIAGNOSTIC 0 0 0 0 0 0 0 0 0 0 0 9 9 9 9
OMIT BASE CASE { Comment out this card, if a base case solution is to be added
C Request special output for MS Excel. EOFE service begins 20 July 2007. The
C output file name will be parallel to the input data, with distinctive type :
C Type --- Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL .exc { Columns 29-32 carry optional file type
C EXTREMA OUTPUT FOR EXCEL { Blank 29-32 will use default file type .xcl
100.E-6 20.E-3 60.
1 1 1 2 1 -1 1 0 3
C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
1 1 0 0 1 { KSTOUT is blank, not 0! } 0
2 2 10 10 { Printout frequency change only if base case
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
C Next come 3 switches whose closing times are to be varied by rolling dice.
C The 1st of the 3 STATISTICS switches has type code "76" which means that
C the distribution is uniform with mean of 2 msec and standard deviation of
C 0.1 msec (columns 15 through 34). The 2nd and 3rd switches have type code
C zero, so use Gaussian distribution. The means are 4 and 6 msec, and the
C standard deviations are 1/2 and 1 msec, respectively.
76A1 ASW1 2.E-3 .1E-3 STATISTICS
B1 BSW1 4.E-3 0.5E-3 STATISTICS 4
C1 CSW1 6.E-3 1.E-3 STATISTICS
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
C The following is DIAGNOSTIC from "RANDNZ" usage of overlay 12. Not a
C normal case output, this is useful in case switching times are incorrect:
C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow ....
C 0.439822971502571000E+02 0.296996800000000000E+07
C 0.296996800000000000E+07 0.000691499561071
C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.205132719793000000E+12
C 0.326925688100000000E+10 0.761183183873072
C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.225804303513790000E+15
C 0.692893886000000000E+09 0.161326929461211
C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.478574878121350000E+14
C 0.296220010300000000E+10 0.689690956613049
C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.204596198914108000E+15
C 0.113680185200000000E+10 0.264682306908071
C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14
C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565
C The following unit-9 connection could come earlier on this first subcase
C because of FORM=UNFORMATTED, which is compatible with VECRSV/VECISV of
C MS-DOS (which uses unit 9 for scratch storage). But the present location
C is always safe, so its use is good practice. It is mandatory for 3rd case.
$OPEN, UNIT=LUNIT9 FILE=dc24at40.ext STATUS=UNKNOWN FORM=UNFORMATTED RECL=1000 !
ASW10 BSW10 CSW10 { Request for these node voltage outputs
C GENA GENC ENDA B1 ASW10
C Reference angle A1 C1 A10 BSW1
C
C GENB ENDA
C B1 A10
C Random switching times for energization number 1 :
C 1 1.8826802E-03 2 4.2472809E-03 3 5.3713577E-03
C 0.0 141.9164096 -72.1873524 -35.0598126 372.1266197 -429.52043
C 1.257471517 -1.85970093
C Times of maxima : .002 .0141 .0127 .0051 .0075
C .02 .017
C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow ....
C 0.439822971502571000E+02 0.785177671157890000E+14
C 0.146997761300000000E+10 0.342255833791569
C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.101529883752298000E+15
C 0.115184215400000000E+10 0.268184150103480
C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.795565857346270000E+14
C 0.906510819000000000E+09 0.211063497466967
C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.626117957575120000E+14
C 0.405748372000000000E+10 0.944706546142697
C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.280246343056681000E+15
C 0.402195997700000000E+10 0.936435530195013
C Random switching times for energization number 2 :
C 1 1.8999095E-03 2 4.7979502E-03 3 7.5258324E-03
C 0.0 141.9164096 61.55446116 -41.9189588 371.6253603 334.5616569
C 1.207725131 -1.83695142
C Times of maxima : .002 .0154 .0139 .0059 .0028
C .0058 .017
BLANK card ending the specification of program outputs (node voltages, here)
C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow ....
C 0.439822971502571000E+02 0.277792753651414000E+15
C 0.285888072600000000E+10 0.665635039564222
C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.197460032864095000E+15
C 0.320639779100000000E+10 0.746547661488876
C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.221462689026580000E+15
C 0.129034293200000000E+10 0.300431375391781
C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.891226959703090000E+14
C 0.212457830900000000E+10 0.494666935177520
C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow ....
C 0.000000000000000000E+00 0.146742499224322000E+15
C 0.646589186000000000E+09 0.150545776356012
C Random switching times for energization number 3 :
C 1 1.9308674E-03 2 3.9933351E-03 3 4.9659084E-03
C 0.0 136.7926646 -89.6601448 -53.6359842 437.6032257 -510.753735
C 1.532437698 -1.74295212
C Times of maxima : .0021 .0051 .0083 .0045 .0074
C .0044 .0169
C Cards associated with the base-case solution must not be present:
$DISABLE { Comment out this card, if a base case solution is to be added
PRINTER PLOT
144 3. 0.0 20. ASW10
BLANK card ending base-case plot cards
$ENABLE { Comment out this card, if a base case solution is to be added
C 1 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 27 1.3500000 0.13500000E+03 0 0
C 28 1.4000000 0.14000000E+03 1 1
C 29 1.4500000 0.14500000E+03 2 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.40833333E+00 1.40208495E+00
C Variance = 8.33333333E-04 8.75092111E-04
C Standard deviation = 2.88675135E-02 2.95819558E-02
C The following disconnects DC24.BIN to protect; automatic result of reopening:
$OPEN, UNIT=LUNIT2 STATUS=SCRATCH FORM=UNFORMATTED { Disconnect *.BIN on LUNIT2
-1 100. GENA A1
$CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
$OPEN, UNIT=LUNIT9 FILE=dc24dum.dum STATUS=SCRATCH !
FIND { Enter interactive search for which shot and which variable gave extremum
HELP { Show the user that there is a summary explanation of "FIND" on-line
C Peak extremum of subset has value 5.62169583E+02 This occurred during
C energization 2 for the variable having names "CSW10 " and " ".
C This was variable 3 of 3 in the subset; 3 shots were searched.
ASW10 BSW10 CSW10 { Search these 3 node voltages for extremum
EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot
C Peak extremum of subset has value 5.59424445E+02 This occurred during
C energization 1 for the variable having names "CSW10 " and " ".
C This was variable 3 of 3 in the subset; 2 shots were searched.
C Remember excluded shots are numbered 2
C Add the following DISK command on 29 September 2002. Look in the January,
C 2003 or later newsletter for a summary. Only the framework of the feature
C now exists in DICTAB, but it is being illustrated here in order that it not
C be forgotten. The idea is this. Each EXCLUDE command excludes another shot
C (number 2 for the preceding). Deterministic simulation of this often is of
C interest, and DISK is designed to create such a data file. The name will
C be SHOT0002.DAT and this file will be located along with the .LIS file,
C assuming parallel naming.
C DISK { Create a disk file of data for deterministic simulation of EXCLUDEd shot
C 6 October 2002, the RESIMULATE command is working at home as illustrated
C by f:\data\DUM24.DAT First, the preceding DISK is enhanced by /LIST:
C 8 October 2002, /FULL is created to use full 10 bytes of T-close or open
C as opposed to the default 9, with one blank on the left as separator.
DISK /LIST /FULL { Create disk file of data for deterministic simulation of EXCLUDEd shot
C RESIMULATE { Abort this data case. Begin with SHOT0002.DAT just created by DISK
C If the preceding line is activated, all following data will be ignored as
C the present data is aborted in favor of the DISK-created SHOT0002.DAT But
C we do not want to lose all the following output. So, comment this line
C here and note that a new 4th subcase of DC-48 illustrates RESIMULATE use.
EXCLUDE { Excluding the shot of previous output means only one is left for use
C Peak extremum of subset has value 5.53096221E+02 This occurred during
C energization 3 for the variable having names "BSW10 " and " ".
C This was variable 2 of 3 in the subset; 1 shots were searched.
C Remember excluded shots are numbered 1 2
RESET { Reinitialize by erasing history of previous selections
ALL { Next, consider searching all output variables of some class
8 { As with batch-mode plotting, "8" indicates branch voltages (search ALL)
C Peak extremum of subset has value 1.41916410E+02 This occurred during
C energization 1 for the variable having names "GENA " and "A1 ".
C This was variable 1 of 3 in the subset; 3 shots were searched.
EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot
EXCLUDE { Excluding the shot of previous output means only one is left for use
EXCLUDE { Final attempt to exclude should lead to error message: nothing is left
ALL { Next, consider searching all output variables of some class
9 { As with batch-mode plotting, "9" indicates branch currents (search ALL)
EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot
C $KEY { Un-comment to allow keyboard input (until "END") at this point
QUIT { Done locating extrema, so exit the "FIND" alternative
C &&&&&&&&&&&&& End request and documentation. Begin 3 coupled cards ......
C 2 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 9 0.4500000 0.13635000E+03 0 0
C 10 0.5000000 0.15150000E+03 3 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 4.75000000E-01 4.62734306E-01
C Variance = 0.00000000E+00 9.53165933E-05
C Standard deviation = 0.00000000E+00 9.76302173E-03
-1 GENA A1 CONT.
C 3 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENC " to "C1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 4 0.2000000 0.60600000E+02 0 0
C 5 0.2500000 0.75750000E+02 2 2
C 6 0.3000000 0.90900000E+02 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 2.41666667E-01 2.45766731E-01
C Variance = 8.33333333E-04 2.19347831E-03
C Standard deviation = 2.88675135E-02 4.68345845E-02
-1 GENC C1 CONT.
C 4 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "ENDA " to "A10 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 2 0.1000000 0.30300000E+02 0 0
C 3 0.1500000 0.45450000E+02 2 2
C 4 0.2000000 0.60600000E+02 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.41666667E-01 1.43690600E-01
C Variance = 8.33333333E-04 9.61072666E-04
C Standard deviation = 2.88675135E-02 3.10011720E-02
-1 ENDA A10
C 5 ) --------------------------------------------------------------------------
C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR
C 5 ) --------------------------------------------------------------------------
C A distribution of peak values among all output branches of the last data card
C statistical distribution is for the maximum of the peaks at all of these outpu
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 9 0.4500000 0.13635000E+03 0 0
C 10 0.5000000 0.15150000E+03 3 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 4.75000000E-01 4.62734306E-01
C Variance = 0.00000000E+00 9.53165933E-05
C Standard deviation = 0.00000000E+00 9.76302173E-03
C 3456789012345678901234567890123456789012345678
C MODTAB AINCR XMAXMX
STATISTICS DATA 1 0.0 0.0 ---- Only individual
-1 GENC C1 { =.245766731, Variance=2.19347831E-3} CONT.
-1 ENDA A10 { 1.43690600E-01, 9.61072666E-04
C 8 ) --------------------------------------------------------------------------
C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR
C 8 ) --------------------------------------------------------------------------
C A distribution of peak values among all output branches of the last data card
C statistical distribution is for the maximum of the peaks at all of these outpu
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 4 0.2000000 0.60600000E+02 0 0
C 5 0.2500000 0.75750000E+02 2 2
C 6 0.3000000 0.90900000E+02 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 2.41666667E-01 2.45766731E-01
C Variance = 8.33333333E-04 2.19347831E-03
C Standard deviation = 2.88675135E-02 4.68345845E-02
STATISTICS DATA 2 0.0 0.0 ------ Only union of 2 curves
-1 GENC C1 CONT.
-1 ENDA A10 { See comments for solution
STATISTICS DATA 3 0.0 0.0 ---- Both individual & union
-1 GENC C1 { =.245766731, Variance=2.19347831E-3} CONT.
-1 ENDA A10 { 1.43690600E-01, 9.61072666E-04
C 12 ) -------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENC " to "C1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 8 0.9131973 0.56211367E+02 0 0
C 9 1.0273470 0.63237788E+02 1 1
C 10 1.1414966 0.70264209E+02 0 1
C 11 1.2556463 0.77290630E+02 1 2
C 12 1.3697959 0.84317051E+02 0 2
C 13 1.4839456 0.91343472E+02 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.19857145E+00 1.20977941E+00
C Variance = 5.21205808E-02 5.31495428E-02
C Standard deviation = 2.28299323E-01 2.30541846E-01
STATISTICS DATA 3 -5.0 0.0 ---- Fixed number of 5 boxes
-1 GENC C1 { See preceding comments for 1st of 3} CONT.
-1 ENDA A10
STATISTICS DATA 3 .05 0.0 ---- Restore STARTUP values
-1 GENC C1 ENDA A10 { =2.45766731E-1, Variance = 2.19347831E-3
0 ASW10 BSW10 CSW10 { Mean=1.40245965E+0, Variance = 8.47042357E-2
C 1st of 4 tables is documented on preceding request; the 4th of 4 follows:
C 21 ) -------------------------------------------------------------------------
C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR
C 21 ) -------------------------------------------------------------------------
C The following is a distribution of peak overvoltages among all output nodes of
C This distribution is for the maximum of the peaks at all output nodes with V-
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 36 1.8000000 0.54540000E+03 0 0
C 37 1.8500000 0.56055000E+03 2 2
C 38 1.9000000 0.57570000E+03 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.84166667E+00 1.84234351E+00
C Variance = 8.33333333E-04 2.35830395E-04
C Standard deviation = 2.88675135E-02 1.53567703E-02
C &&&&&&&&& End documentation of node voltage tabulations (4th of 4 tables)
STATISTICS DATA 3 -4.0 0.0 ---- Fixed number of 4 boxes
C 22 ) -------------------------------------------------------------------------
C Statistical distribution of peak power for branch "B1 " to "BSW1 ".
C Interval power power in Frequency Cumulative
C number in per unit physical units (density) frequency
C 16 0.9468728 0.35188195E+03 0 0
C 17 1.0060524 0.37387457E+03 2 2
C 18 1.0652319 0.39586719E+03 0 2
C 19 1.1244115 0.41785981E+03 0 2
C 20 1.1835910 0.43985244E+03 1 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.03564214E+00 1.05962916E+00
C Variance = 1.05066576E-02 1.04274409E-02
C Standard deviation = 1.02501988E-01 1.02114842E-01
-3 B1 BSW1 { See preceding comments for this branch power table
-2 0. GENB B1 ENDA A10 { =1.10335104E+0, Variance = 2.09698061E-2
-4 0. B1 BSW1 { Ungroup: Mean=1.16595031E+0, Variance = 3.80542557E-2
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
C 2nd of 6 subcases repeats the preceding problem exactly, only instead
C of rolling dice to find the switch closing times, such times will be
C read from a user-supplied disk file dc24lun8.dat (see $OPEN use).
C The solution of this 2nd subcase should exactly agree with the first.
C However , for variety , a base case solution (missing from the first
C subcase) and a character plot have been added to this 2nd subcase.
C Also, the ".PL4" plot file (actually, just the header) is saved
C permanently for later re-connection using the 4th subcase of DC-40.
PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP
C $STARTUP, dc37star.dat { Use disk file for re-initialization at this point
C This is a Monte Carlo case, for which TENERG is critical. But this
C parameter is initialized only in the STARTUP file. For this 2nd or
C later data case, we reinitialize via the preceding $STARTUP request
USER SUPPLIED SWITCH TIMES { Bypass the dice, reading TCLOSE & TOPEN from LUNIT8
C The following connection of unit 2 is required only because this is the
C 2nd subcase. It follows a case in which MEMSAV=1 dumps tables via I/O
C unit LUNIT2=2, and then this file is disconnected. Without the manual
C connection that follows, nothing would be connected, and execution would
C die for Apollo in "SUBR5" where REWIND LUNIT2 is found. So connect it:
$OPEN, UNIT=LUNIT8 FILE=dc24lun8.dat FORM=FORMATTED STATUS=OLD RECL=80 ! { Case
$OPEN, UNIT=LUNIT2 FORM=UNFORMATTED STATUS=SCRATCH RECL=16000
C $STARTUP, dc37star.dat { Re-initialize is required to set TENERG for table dump
C { This is a peculiarity of having "STATISTICS" as the
C { 2nd or later data subcase within a single disk file.
C DISK PLOT DATA { Toggle the Apollo default of LUNIT4 = -4 to +4 (use disk)
$CLOSE, UNIT=LUNIT4 STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP"
$OPEN, UNIT=LUNIT4 FILE=dc24b40d.pl4 ! { All we need for C-like case
C Before continuing with case, let's document the contents of that unit-8 file:
C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 1:
C 1.882680E-03 4.247281E-03 5.371358E-03 7.95E-3 9.95E-3
C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line:
C .1E-3 0.5E-3 1.E-3 1.0 1.0
C 1.0
C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 2:
C 1.899909E-03 4.797950E-03 7.525832E-03 7.95E-3 9.95E-3
C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line:
C .1E-3 0.5E-3 1.E-3 1.0 1.0
C 1.0
C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 3:
C 1.930867E-03 3.993335E-03 4.965908E-03 7.95E-3 9.95E-3
C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line:
C .1E-3 0.5E-3 1.E-3 1.0 1.0
C 1.0
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
C Cancel special output for MS Excel of 1st subcase. Addition on 20 July 2007:
C Type --- Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL .off { Columns 29-32 carry optional file type
100.E-6 20.E-3 60.
1 1 1 2 1 -1 0 2 3
C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
1 1 0 0 1 { KSTOUT is blank, not 0! } 1
2 2 10 10 { Printout frequency change only if base case
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
76A1 ASW1 1.95E-3 .1E-3 STATISTICS
B1 BSW1 3.95E-3 0.5E-3 STATISTICS 4
C1 CSW1 5.95E-3 1.E-3 STATISTICS
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14
C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565
ASW10 BSW10 CSW10
C Step Time GENA GENC ENDA B1 ASW10 BSW10
C A1 C1 A10 BSW1
C 0 0.0 -.5684E-13 .28422E-13 -.5684E-13 -151.5 0.0 0.0
C 1 .1E-3 .56843E-13 .29559E-11 -.2842E-12 0.0 0.0 0.0
C 2 .2E-3 0.0 -.2927E-11 .28422E-12 0.0 0.0 0.0
BLANK card ending the specification of program outputs (node voltages, here)
C Last step: 200 .02 1.60684633 5.46144739 15.5569682 209.921255 126.692915
C Variable max:136.792665 67.5308064 48.9224823 437.603226 320.835627 553.096221
C Times of max: .0021 .0061 .0093 .0045 .0029 .0049
C Variable min:-100.37615 -85.515782 -58.824356 -292.03358 -548.79093 -353.02344
C Times of min: .0033 .0127 .0084 .0056 .0074 .0141
PRINTER PLOT
144 1. 0.0 10. ASW10 { Plot limits: (-5.488, 3.208)
$CLOSE, UNIT=LUNIT4 STATUS=KEEP
$OPEN, UNIT=LUNIT4 FILE=dc24bdum.pl4 ! { All we need for C-like case
BLANK card ending base-case plot cards
C GENA GENC ENDA B1 ASW10
C Reference angle A1 C1 A10 BSW1
C Random switching times for energization number 1 :
C 1 1.8826800E-03 2 4.2472810E-03 3 5.3713580E-03
C 0.0 141.91641 -72.187352 -35.059813 372.12662 -429.52043 422.745884
C Times of maxima : .002 .0141 .0127 .0051 .0075 .0051
C
C Random switching times for energization number 2 :
C 1 1.8999090E-03 2 4.7979500E-03 3 7.5258320E-03
C 0.0 141.91641 61.5544612 -41.918959 371.62536 334.561657 497.525193
C Times of maxima : .002 .0154 .0139 .0059 .0028 .0065
C
C Random switching times for energization number 3 :
C 1 1.9308670E-03 2 3.9933350E-03 3 4.9659080E-03
C 0.0 136.792665 -89.660145 -53.635984 437.603226 -510.75374
C cont. 553.096221 513.648909 .96346791 1.5324377 -1.7429521
C .0021 .0051 .0083 .0045 .0074
C cont. .0049 .011 .02 .0044 .0169
C Cards associated with the base-case solution must not be present:
C 1 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 27 1.3500000 0.13500000E+03 0 0
C 28 1.4000000 0.14000000E+03 1 1
C 29 1.4500000 0.14500000E+03 2 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.40833333E+00 1.40208495E+00
C Variance = 8.33333333E-04 8.75092111E-04
C Standard deviation = 2.88675135E-02 2.95819558E-02
-1 100. GENA A1 { See preceding comments that document this tabulation
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
C 3rd of 6 subcases repeats the first subcase exactly, only instead of
C dumping full tables, MEMSAV = 2 (columns 49-56) means that only
C the extrema and switching times are to be dumped, and
C this will be FORMATTED rather than UNFORMATTED (for
C this change, only $OPEN is changed). The resulting card-image file
C named dc24ct40.lis is used by the 2nd and later subcases of DC-40.
C Also, as an illustration, time-step loop output of all energizations
C will be shown (because column-48 variable KSTOUT is zero not blank),
C and the miniature character plots of switch times will be suppressed.
PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP
C $STARTUP, dc37star.dat { Use disk file for re-initialization at this point
C This is a Monte Carlo case, for which TENERG is critical. But this
C parameter is initialized only in the STARTUP file. For this 2nd or
C later data case, we reinitialize via the preceding $STARTUP request
OMIT BASE CASE { Comment out this card, if a base case solution is to be added
100.E-6 20.E-3 60.
1 1 1 2 1 -1 2 0 3
C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
C 1 1 0 0 1 0 { STATIST. misc. data card } 1
1 1 0 0 2309 0 { STATIST. misc. data card } 1
2 2 10 10 { Printout frequency change for all soutions
C Preceding STATISTICS miscellaneous data card is handled differently beginning
C 23 September 2005 as new reserved value IDICE = 2309 is introduced as a
C request for alternative, computer-dependent dice. This is to satisfy Apple
C Macintosh G5 as used by Stu Cook of JUST Services in Ontario, Canada. As of
C this date, Mac G5 ATP is the only version to honor the 2309 alternative.
C For all versions, any other IDICE value (including the value unity shown)
C is treated the same as zero or blank. I.e., any other value does nothing.
C Henceforth, IDICE = 2309 is a request for alternative, computer-dependent
C random numbers (dice) of the operating system's library. So, we illustrate
C use here to document the two associated warning lines immediately before the
C input of branch data begins. If Stu Cook wants to compare his results with
C everybody else's, he should restore the value IDICE = 1 (alternatively,
C zero or blank would be equivalent). Salford ATP output follows:
C Statistics data. 1 1 0 0 2309 1 | 1 1 ...
C Printout : 2 2 10 10 0 0 | 2 2 ...
C Computer-dependent (not universal) dice have been requested for this
C Monte Carlo study by reserved value IDICE = 2309.
C Sorry, no alternative computer-dependent dice are available. ATP will
C ignore this request. Universal dice will be used.
C Series R-L-C. 0.000E+00 1.857E-02 0.000E+00 | 0GENA A1 ...
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
76A1 ASW1 2.E-3 .1E-3 STATISTICS
B1 BSW1 4.E-3 0.5E-3 STATISTICS 4
C1 CSW1 6.E-3 1.E-3 STATISTICS
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14
C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565
C The following unit-9 connection must not be made earlier, since the same
C I/O channel is used as UNFORMATTED by VECRSV/VECISV of MS-DOS. This is
C different than the first subcase (read those comments), note:
$OPEN, UNIT=LUNIT9 FILE=dc24ct40.lis STATUS=UNKNOWN FORM=FORMATTED !
ASW10 BSW10 CSW10 { Request for these node voltage outputs
C GENA GENC ENDA B1 ASW10 BSW10
C Reference angle A1 C1 A10 BSW1
C Random switching times for energization number 1 :
C 1 1.8826802E-03 2 4.2472809E-03 3 5.3713577E-03
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 1 .1E-3 0.0 0.0 0.0 0.0 0.0 0.0
C 2 .2E-3 0.0 0.0 0.0 0.0 0.0 0.0
C 4 .4E-3 0.0 0.0 0.0 0.0 0.0 0.0
C 6 .6E-3 0.0 0.0 0.0 0.0 0.0 0.0
C 8 .8E-3 0.0 0.0 0.0 0.0 0.0 0.0
C 10 .001 0.0 0.0 0.0 0.0 0.0 0.0
C ==== Table dumping for all subsequent restorations. Time [sec] = 1.8000E-03
C *** Close switch "A1 " to "ASW1 " after 1.90000000E-03 sec.
C 20 .002 141.91641 0.0 0.0 0.0 0.0 0.0
C 30 .003 -34.73764 0.0 0.0 0.0 318.450266 -159.22513
BLANK card ending the specification of program outputs (node voltages, here)
C 190 .019 -8.1503339 19.8420817 .549065534 24.9565797 230.197042 85.0661304
C 200 .02 .352725683 14.1822201 13.4881293 207.346558 128.761754 181.144509
C 201 .0201 -.23343325 23.5035485 8.73828958 199.126365 123.327149 204.30575
C Extrema of output variables follow. Order and column positioning are the sam
C Variable max:136.792665 56.3388361 40.5581876 437.603226 320.835627 553.096221
C Times of max: .0021 .0053 .0092 .0045 .0029 .0049
C Variable min:-100.37615 -89.660145 -53.635984 -272.03803 -510.75374 -342.9791
C Times of min: .0033 .0051 .0083 .0055 .0074 .014
C Warning! No statistical tabulations are allowed here because of our use
C ------- of output file LUNIT9 = 9 being "FORMATTED". Any attempted
C tabulation will die in "MCBANK" called by "DICTAB". 7 Dec 88
C -1 100. GENA A1
$CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
$OPEN, UNIT=LUNIT9 FILE=dc24dum.dum STATUS=UNKNOWN FORM=FORMATTED RECL=1000 !
C The preceding $OPEN is added 11 October 2010 as required by new 4th subcase
C 1 ) --------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ".
C Interval voltage voltage in Frequency Cumulative
C number in per unit physical units (density) frequency
C 27 1.3500000 0.13500000E+03 0 0
C 28 1.4000000 0.14000000E+03 1 1
C 29 1.4500000 0.14500000E+03 2 3
C Distribution parameters for .. data follow. Grouped data Ungrouped data
C Mean = 1.40833333E+00 1.40208495E+00
C Variance = 8.33333333E-04 8.75092111E-04
C Standard deviation = 2.88675135E-02 2.95819558E-02
C Column positioning of "NO SWITCH PLOTS" on blank terminator is arbitrary:
BLANK card ending statistical tabulation requests --- NO SWITCH PLOTS ---
BEGIN NEW DATA CASE
C 4th of 6 subcases produces simulations that are identical to the 2nd.
C However, output looks substantially different because STATISTICS
C switches are not used. Instead of USER SUPPLIED SWITCH TIMES to
C define the closing times of switches, the SELECT function of the
C pocket calculator will be used. Structural differences include lack
C of a base case (the 2nd subcase had 1) preceding the 3 energizations
C Output of connectivity and the phasor solution will be suppressed so
C as not to be repeated 3 times (boring). Batch-mode plotting will be
C omitted for the same reason. Finally, output of the dT-loop will be
C be made less frequent. Data is added 10 October 2010. WSM.
PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP
C Controls of following request card: MAXKNT IOPCVP NOSTAT { Loop 3 times with
POCKET CALCULATOR VARIES PARAMETERS 3 0 0 { full printout
100.E-6 20.E-3 60.
20 1 0 0 1
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b
TCLOSE_AAA = SELECT ( KNT ) { Closing time of phase "a" is taken from following:
.00188268.001899909.001930867 { Data is 3E10.0 for 3 energizations of phase a
C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase a
TCLOSE_BBB = SELECT ( KNT ) { Closing time of phase "b" is taken from following:
.004247281 .00479795.003993335 { Data is 3E10.0 for 3 energizations of phase b
C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase b
TCLOSE_CCC = SELECT ( KNT ) { Closing time of phase "c" is taken from following:
.005371358.007525832.004965908 { Data is 3E10.0 for 3 energizations of phase c
C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase c
BLANK card ends $PARAMETER definitions that are processed just b4 branch card
C Switch cards begin with the 3 STATISTICS switches of the 2nd subcase, which
C are being converted to deterministic switches with T-close of columns 15-
C 24 defined by the preceding $PARAMETER block. T-open has been set to 99
C seconds (never reached) to ensure that the switches never open:
C 76A1 ASW1 1.95E-3 .1E-3 STATISTICS
A1 ASW1 TCLOSE_AAA 99.0
C B1 BSW1 3.95E-3 0.5E-3 STATISTICS 4
B1 BSW1 TCLOSE_BBB 99.0 4
C C1 CSW1 5.95E-3 1.E-3 STATISTICS
C1 CSW1 TCLOSE_CCC 99.0
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
C The following start of the dT-loop output is for the 1st of the 3 passes:
C Column headings for the 10 EMTP output variables follow. These are divided among the 5 possible classes as follows ....
C First 7 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 Branch power consumption (power flow, if a switch) is treated like a branch voltage for this grouping (1 variables);
C Branch energy consumption (energy flow, if a switch) is treated like a branch current for this grouping (1 variables).
C Step Time GENA GENC ENDA B1 ASW10 BSW10 CSW10 B1 GENB ENDA
C A1 C1 A10 BSW1 BSW1 B1 A10
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C *** Close switch "A1 " to "ASW1 " after 1.90000000E-03 sec.
C 20 .002 141.91641 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 40 .004 2.37028418 0.0 0.0 0.0 -46.220806 23.1104031 23.1104031 0.0 0.0 0.0
C *** Close switch "B1 " to "BSW1 " after 4.30000000E-03 sec.
C *** Close switch "C1 " to "CSW1 " after 5.40000000E-03 sec.
C 60 .006 -15.795936 10.618273 0.0 -123.04522 -168.96416 422.15175 -253.18759 .357709401 -.40883323 0.0
ASW10 BSW10 CSW10
BLANK card ending the specification of program outputs (node voltages, here)
C 180 .018 7.07846743 32.9539555 -12.904568 -.92594731 299.739293 -60.943514 -240.37687 .62556771 .018267223 -1.498213
C 200 .02 5.74620332 -11.67979 16.6256782 248.357719 125.624205 168.195467 -283.39188 .769701668 1.25747152 -1.2455527
C Variable maxima : 141.91641 56.029476 33.4968354 372.12662 334.561657 422.745884 559.424445 .769701668 1.25747152 .049072408
C Times of maxima : .002 .0122 .0092 .0051 .0028 .0051 .0114 .02 .02 .0095
C Variable minima : -103.91401 -72.187352 -35.059813 -191.64938 -429.52043 -351.293 -316.83918 0.0 -1.0806309 -1.8597009
C Times of minima : .0032 .0141 .0127 .0159 .0075 .0142 .0191 0.0 .0109 .017
BLANK card ending batch-mode plot cards
C Following statistical tabulation is identical to that of 2nd subcase of DC-24:
C 1 ) -------------------------------------------------------------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ". Base voltage for per unit output = 1.00000000E+02
C Interval voltage voltage in Frequency Cumulative Per cent
C number in per unit physical units (density) frequency .GE. current value
C 27 1.3500000 1.35000000E+02 0 0 100.000000
C 28 1.4000000 1.40000000E+02 1 1 66.666667
C 29 1.4500000 1.45000000E+02 2 3 .000000
C Summary of preceding table follows: Grouped data Ungrouped data
C Mean = 1.40833333E+00 1.40208495E+00
C Variance = 8.33333333E-04 8.75092111E-04
C Standard deviation = 2.88675135E-02 2.95819558E-02
-1 100. GENA A1 { See preceding comments that document this tabulation
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
C 5th of 6 subcases produces simulations that are identical to the 4th.
C However, instead of the SELECT function (which is a generalization
C of USER SUPPLIED SWITCH TIMES), dice will be rolled just as was the
C case for the 1st subcase. Yes, the same dice, with the RAN function
C and the GAUSS function of the pocket calculator replacing internal
C rolling that ATP performs in support of STATISTICS switches. Within
C the PCVP loop, the user must roll his own dice after first ensuring
C that seeding of the random number generator is identical to what ATP
C does for the 1st subcase. Data is added 16 October 2010. WSM.
PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP
REPEATABLE RANDOM NUMBERS { Needed since no NSEED of STATISTICS misc. data card
C Note about preceding line. Without it, dice no longer will be deterministic
C and the resulting .LIS file should change every time data is simulated.
C That assumes KOMPAR of STARTUP is small enough (e.g., value unity) so that
C the program date and time are real. But for KOMPAR = 4, all simulations
C will be identical because seeding of the number generator then will depend
C on a fixed time (the WW I Armistice time) and zero elapsed CP time. I.e.,
C program output will be deterministic although different from the output for
C the data as shown. That is because the fixed seed of RRN differs from the
C fixed seed that results from the WW I time. Complicated, eh?
C Controls of following request card: MAXKNT IOPCVP NOSTAT { Loop 3 times with
POCKET CALCULATOR VARIES PARAMETERS 3 0 0 { full printout
100.E-6 20.E-3 60.
20 1 0 0 1
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b
C Diagnostic printout of the 1st subcase reveals that dice are rolled 5 times
C prior to each energization. There are 3 STATISTICS switches, so 3 of the 5
C are used for them. But prior to these, there are 2 more. First, there is
C a wasted roll that is used for nothing. This has been known for many years
C but was not removed for reasons of historical continuity. Second, there is
C a roll that is used for the uniform bias that might, if the user wants it,
C be added to the closing times of all switches. So, here are the 5 rolls :
UNUSE1of2 = GAUSS ( 0.0 ) { OVER12 wastes 1 use for no good reason if STATISTIC
UNUSE2of2 = GAUSS ( 0.0 ) { OVER12 rolls a 2nd time for bias if STATISTICS case
C Now we are ready for the 3 switches. The 1st has type code 76, which means
C a uniform distribution. This is a little involved because it is low level.
C The 1st subcase uses a mean of 2 msec and a standard deviation of 0.1 msec.
C But the RAN function produces a uniform distribution over (0, 1). To map
C the RAN output RNG into the desired closing time, OVER12 uses this code:
C RNG = TIMEV + BIAS + SQRT3 * ABSZ ( TOPEN(NDXI) ) * ( 2.0 * RNG - 1.0 )
C For this case, there is no bias BIAS, TIMEV is the mean equal to 2 msec,
C SQRT3 is the square root of 3, and TOPEN is the standard deviation. So:
TCLOSE_AAA = 2.0E-3 + 0.17320508E-3 * ( 2.0 * RAN ( 0.0 ) - 1.0 ) { Uniform closing time of phase "a"
C This is slightly more efficient, using the constant 0.17320508E-3 rather
C than scaling an ATP variable. But the answers should be identical. ATP is
C smart enough to recognize SQRT3 as a program variable, and use the full,
C internal precision associated with it. Program output will be unchanged if
C the preceding were replaced by:
C TCLOSE_AAA = 2.0E-3 + SQRT3 * 1.E-4 * ( 2.0 * RAN ( 0.0 ) - 1.0 ) { Uniform closing time of phase "a"
C Switches 2 and 3 (phases b and c) are easier since they are Gaussian:
TCLOSE_BBB = 4.0E-3 + 0.5E-3 * GAUSS ( 0.0 ) { Gauss closing time of phase "b"
TCLOSE_CCC = 6.0E-3 + 1.0E-3 * GAUSS ( 0.0 ) { Gauss closing time of phase "c"
BLANK card ends $PARAMETER definitions that are processed just b4 branch card
C 76A1 ASW1 2.E-3 .1E-3 STATISTICS
A1 ASW1 TCLOSE_AAA 99.0
C B1 BSW1 4.E-3 0.5E-3 STATISTICS 4
B1 BSW1 TCLOSE_BBB 99.0 4
C C1 CSW1 6.E-3 1.E-3 STATISTICS
C1 CSW1 TCLOSE_CCC 99.0
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
ASW10 BSW10 CSW10
BLANK card ending the specification of program outputs (node voltages, here)
BLANK card ending batch-mode plot cards
C Following statistical tabulation is identical to that of 4th subcase of DC-24:
C 1 ) -------------------------------------------------------------------------------------------------------------------------------
C Statistical distribution of peak voltage for branch "GENA " to "A1 ". Base voltage for per unit output = 1.00000000E+02
C Interval voltage voltage in Frequency Cumulative Per cent
C number in per unit physical units (density) frequency .GE. current value
C 27 1.3500000 1.35000000E+02 0 0 100.000000
C 28 1.4000000 1.40000000E+02 1 1 66.666667
C 29 1.4500000 1.45000000E+02 2 3 .000000
C Summary of preceding table follows: Grouped data Ungrouped data
C Mean = 1.40833333E+00 1.40208495E+00
C Variance = 8.33333333E-04 8.75092111E-04
C Standard deviation = 2.88675135E-02 2.95819558E-02
-1 100. GENA A1 { See preceding comments that document this tabulation
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
C 6th of 6 subcases is added 24 July 2011 to illustrate the new
C protection against overflow of List 23. Data is the same as for
C subcase 1 although the base case has been activated to prove that
C omission is not a requirement. Note the NENERG = 90K shots of
C the integer misc. data card. For other translations, this is a
C case of overkill. But for Salford ATP, ordered COMMON blocks
C result in an enormous List 23 in excess of a million words. For
C Salford ATP, such huge NENERG is necessary. WSM.
100.E-6 20.E-3 60.
50000 20001 0 2 1 -1 1 0 90000
C Interpretation of the preceding integer misc. data card is improved on
C 7 August 2011 when alternative variable-width encoding and the use of "K"
C allow the handling of much larger integers. The old format was (2I5, 8I3)
C which resulted in overflow for NENERG > 999. The new is not troubled :
C Misc. data. 1.000E-04 2.000E-02 6.000E+01 | 100.E-6 20.E-3 60.
C Misc. data. 50K 20K 0 2 1 -1 1 0 90K 0 | 50000 20001 ...
C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
1 1 0 0 1 { KSTOUT is blank, not 0! } 0
2 2 10 10 { Printout frequency change only if base case
0GENA A1 7. 2
0GENB B1 7. 1
0GENC C1 7. 2
0ENDA A10 7. 3
0ENDB B10 7.
0ENDC C10 7.
-1ASW1 A5 .3 2.1146 0.645 50. 0
-2BSW1 B5 .0268 .5397 0.021 50. 0
-3CSW1 C5
0A5 A5F 1.
0B5 B5F 1.
0C5 C5F 1.
-1A5F ASW10 ASW1 A5
-2B5F BSW10
-3C5F CSW10
BLANK card ending branch cards
76A1 ASW1 2.E-3 .1E-3 STATISTICS
B1 BSW1 4.E-3 0.5E-3 STATISTICS 4
C1 CSW1 6.E-3 1.E-3 STATISTICS
A10 ASW10 7.95E-3 1.0
B10 BSW10 9.95E-3 1.0
C10 CSW10 11.95E-3 1.0
BLANK card ending switch cards
14GENA 303. 60. 0.0 -1.
14GENB 303. 60. -120.0 -1.
14GENC 303. 60. 120.0 -1.
14ENDA 303. 60. - 10.0 -1.
14ENDB 303. 60. -130.0 -1.
14ENDC 303. 60. 110.0 -1.
BLANK card ending source cards
ASW10 BSW10 CSW10 { Request for these node voltage outputs
BLANK card ending the specification of program outputs (node voltages, here)
C Overflow of List 23. General information follows. Specifically, here in SUBR15 prior to entry into the time-step loop, the demand
C for storage that is required for a Monte Carlo (or a SYSTEMATIC) study has been computed, and has been found to be inadequate.
C Total working space is L23TOT = 1502629 cells whereas each shot requires NDX1 = 24 cells. For all NENERG = 90000 shots
C requested, at least 2160000 will be needed for List 23. But recall that a Monte Carlo study can be simulated in pieces using two
C or more separate ATP executions. It is safer not to place all eggs in one basket. Also, simulation should be faster if two or
C more central processors are used simultaneously. CENTRAL STATISTICS FILE is the request word for such operation as illustrated by
C BENCHMARK DC-66.
C ------------------------------------------------------------------------------------------------------------------------------------
C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/
C Etc. (the familiar KILL = 1 error termination).
BLANK card ending base-case plot cards
-1 GENA A1
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
BLANK
|
