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
|
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
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
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=9 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:
C $OPEN, UNIT=22 STATUS=SCRATCH FORM=UNFORMATTED { Disconnect *.BIN on LUNIT2
$CLOSE, UNIT=22 STATUS=UNKNOWN { Disconnect to prevent damage by stray WRITE
$OPEN, UNIT=22 STATUS=SCRATCH FILE=dumxx.bin ! { Reconnect dummy .BIN on LUNIT2
-1 100. GENA A1
$CLOSE, UNIT=9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
$OPEN, UNIT=9 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
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 4.37603226E+02 This occurred during
C energization 3 for the variable having names "B1 " and "BSW1 ".
C This was variable 4 of 4 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 3 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=8 FILE=dc24lun8.dat FORM=FORMATTED STATUS=OLD RECL=80 ! { Case
C $OPEN, UNIT=22 FORM=UNFORMATTED STATUS=SCRATCH RECL=16000
$CLOSE, UNIT=22 STATUS=UNKNOWN { Disconnect to prevent damage by stray WRITE
$OPEN, UNIT=22 STATUS=SCRATCH FILE=dumxx.bin ! { Reconnect dummy .BIN on LUNIT2
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=4 STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP"
$OPEN, UNIT=4 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
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=4 STATUS=KEEP
$OPEN, UNIT=4 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 3 subcases repeats the first subcase exactly, only instead of
C dumping full tables, MEMSAV = 2 (columns 59-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
1 1 0 0 1 0 { STATIST. misc. data card } 1
2 2 10 10 { Printout frequency change for all soutions
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=9 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=9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
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
BLANK
|
