summaryrefslogtreecommitdiffstats
path: root/benchmarks/dc40.dat
blob: dea4726fe239232310cfc1bb53302f930d326aba (plain)
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
BEGIN NEW DATA CASE
C     BENCHMARK DC-40
C     This present data case restarts the  "STATISTICS" (Monte Carlo) simulation
C     of  DC-24,  which ran only  3 energizations.  When this 2nd half finishes,
C     6 will have been completed, to be compared with  "M39."  answers of DC-24.
C     Yet answers differ due to different random number generator and use of the
C     seed (now, initialization only occurs once).  As for the disk file that is
C     connected to  UNIT=2  immediately below,  note  that  it is disconnected 3
C     lines later,  and  DUM.BIN  is reconnected.   This is to protect  DC24.BIN
C     from a stray write.  But beware!  The replacement file  DUM.BIN  must be a
C     perfect copy of the original DC24TO40.BIN to ensure correct answers, since
C     it  is  DUM.BIN  from which restoration will actually occur for the second
C     and later simulations of this restarted case.  Finally, be careful of OPMC
C     Beware of  PARALLEL.LIS  entry for same  USERID  that might modify NENERG!
$OPEN,  UNIT=LUNIT2  FILE=dc24.bin  STATUS=OLD  FORM=UNFORMATTED ! { Hold case
START AGAIN                      { Request for awakening of hibernating solution
    9999               { Terminate modifications (none, here) to switching times
C Cancel special output for MS Excel of DC-24 tables.  Addition on 21 July 2007:
C                           Type  ---  Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL    .off        { Columns 29-32 carry optional file type
$OPEN, UNIT=LUNIT2  FILE=dc24at40.ext  STATUS=OLD  FORM=UNFORMATTED  RECL=16000  !
LOAD MORE SHOTS    { Load the separate extrema produced by 3 shots made in DC-24
MISCELLANEOUS DATA CARDS
 100.E-6  20.E-3     60.
       1       1       1       1       1               0       0       6
C    ISW   ITEST   IDIST    IMAX   IDICE  KSTOUT  KNTRPT                   NSEED
       1       1       0       0       1         { KSTOUT is blank, not 0! }   1
C                     GENA         GENC          ENDA         B1           ASW10
C Reference angle       A1           C1           A10        BSW1
C
C                     GENB         ENDA
C                      B1           A10
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 =   4.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.439822971502571000E+02   0.446592684878350000E+14
C                                   0.198544027000000000E+09   0.046227133600041
C      "RFUNL1/RANDNZ".   KNT =   4.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.137132374008640000E+14
C                                   0.370179203200000000E+10   0.861890621483326
C      "RFUNL1/RANDNZ".   KNT =   4.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.255679073858209000E+15
C                                   0.396569462500000000E+10   0.923335232073441
C      "RFUNL1/RANDNZ".   KNT =   4.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.273906562054126000E+15
C                                   0.361268631800000000E+10   0.841144080739468
C      "RFUNL1/RANDNZ".   KNT =   4.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.249524631297943000E+15
C                                   0.421126952700000000E+10   0.980512594571337
C                           Random switching times for energization number  4  :
C                       1  2.1466476E-03    2  4.4995694E-03    3  8.0648931E-03
C          0.0   125.9695623   -58.920794  52.58141134  370.8688937  -424.516718
C                1.249136867  -1.87200012
C Times of maxima :    .0023        .0114        .0094        .0049        .0078
C                      .0049        .0171
TIME STEP LOOP   { Transfer control to overlay 16 for a resumption of simulation
C      "RFUNL1/RANDNZ".   KNT =   5.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.439822971502571000E+02   0.290868174960364000E+15
C                                   0.104773356000000000E+09   0.024394447915256
C      "RFUNL1/RANDNZ".   KNT =   5.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.723659092556500000E+13
C                                   0.386599910100000000E+10   0.900123059051111
C      "RFUNL1/RANDNZ".   KNT =   5.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.267020691906970000E+15
C                                   0.257511465000000000E+10   0.599565601442009
C      "RFUNL1/RANDNZ".   KNT =   5.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.177860593760851000E+15
C                                   0.170306619500000000E+10   0.396525998366997
C      "RFUNL1/RANDNZ".   KNT =   5.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.117629079022456000E+15
C                                   0.280968690400000000E+10   0.654181210324168
C                           Random switching times for energization number  5  :
C                       1  2.0344905E-03    2  3.8690347E-03    3  6.3961916E-03
C          0.0   131.4745297  -76.1213478  -90.4743742  472.0143285  -556.117151
C                1.594127662  -1.83850686
C Times of maxima :    .0022        .0127        .0083        .0139        .0074
C                      .0044        .0172
C      "RFUNL1/RANDNZ".   KNT =   6.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.439822971502571000E+02   0.194062264772377000E+15
C                                   0.275743720900000000E+10   0.642015880206600
C      "RFUNL1/RANDNZ".   KNT =   6.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.190453430588422000E+15
C                                   0.169578189400000000E+10   0.394829989876598
C      "RFUNL1/RANDNZ".   KNT =   6.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.117125959636687000E+15
C                                   0.220147476700000000E+10   0.512570786988363
C      "RFUNL1/RANDNZ".   KNT =   6.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.152053660681924000E+15
C                                   0.322846893200000000E+10   0.751686499454081
C      "RFUNL1/RANDNZ".   KNT =   6.   SY, D6, SEEDRN, D18 (result)  follow ....
C                            0.000000000000000000E+00   0.222987120664309000E+15
C                                   0.100859058100000000E+10   0.234830794157460
C                           Random switching times for energization number  6  :
C                       1  2.0043546E-03    2  4.3397550E-03    3  5.2772330E-03
C          0.0   131.4745297  -102.270651   45.1085527  370.2365789  -461.601966
C                 1.25947653  -1.88310377
C Times of maxima :    .0022        .0054        .0094        .0052        .0076
C                        .02         .017
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         25           1.2500000      0.12500000E+03             0             0
C         26           1.3000000      0.13000000E+03             1             1
C         27           1.3500000      0.13500000E+03             2             3
C         28           1.4000000      0.14000000E+03             1             4
C         29           1.4500000      0.14500000E+03             2             6
C Summary for following request:        Mean = 1.35833333E+00     1.34924018E+00
C                                   Variance = 3.66666667E-03     4.10518277E-03
C                         Standard deviation = 6.05530071E-02     6.40717002E-02
-1  100.      GENA  A1
-1            GENC  C1    ENDA  A10   { =2.53418455E-1, Variance = 3.02359400E-3
              ASW10 BSW10 CSW10   { Mean=1.49453887E+0, Variance = 6.41172862E-2
C                         MODTAB   AINCR  XMAXMX
STATISTICS DATA                1    -10.     0.0   ---- Fixed number of 10 boxes
-3            B1    BSW1 { Ungroup: Mean=1.07790313E+0, Variance = 1.44614228E-2
STATISTICS DATA                1     .05     0.0  --- Return to p.u. box of 1/20
-2    0.      GENB  B1    ENDA  A10   { =1.11785582E+0, Variance = 1.92038029E-2
STATISTICS DATA                1    -5.0     0.0   ---- Fixed number of  5 boxes
-4    0.      B1    BSW1 { Ungroup: Mean=1.27016791E+0, Variance = 8.50686724E-2
C 12 ) -------------------------------------------------------------------------
C Statistical distribution of peak energy   for branch  "B1    "  to   "BSW1  ".
C   Interval           energy          energy   in      Frequency     Cumulative
C     number      in per unit      physical units       (density)      frequency
C          5        0.9692791      0.66169358E+00               0              0
C          6        1.1631349      0.79403229E+00               2              2
C          7        1.3569907      0.92637101E+00               2              4
C          8        1.5508465      0.10587097E+01               1              5
C          9        1.7447023      0.11910484E+01               0              5
C         10        1.9385582      0.13233872E+01               1              6
C Distribution parameters for .. data follow.    Grouped data     Ungrouped data
C                                       Mean = 1.32468140E+00     1.27016791E+00
C                                   Variance = 8.51815081E-02     8.50686724E-02
C                         Standard deviation = 2.91858713E-01     2.91665343E-01
C   Switch  "A1    "  to  "ASW1  "
C     8.33333377E-03   switch closings per column
C              Columns
C          5   10   15   20   25   30
C     +----+----+----+----+----+----+            (Note: all "A" below did line
C     A                             *                   up before left shift,
C     A                             *                   in about column 127).
C     A                             *
C     A                             *
C     |                             *             Pull in to col 80  ---->     A
C     |                             *             Pull in to col 80  ---->     A
C     A                             *
C     |                             *             Pull in to col 80  ---->     A
C     A                             *
C     A                             *
C     A                             *
C     A                             *
C     |  2.0000E-03                 *             Pull in to col 80  ---->     A
C     A                             *
C     |                             *             Pull in to col 80  ---->     A
C     A                             *
C     A                             *
C     A                             *
C     A                             *
C     A                             *
C     A                             *
C     A                             *
C     |                             *             Pull in to col 80  ---->     A
C     A                             *
C     A                             *
C   Time (scale = 1.44337567E-05  Sec/line)
BLANK card ending statistical output variable requests
BEGIN NEW DATA CASE
C         ======================================================================
C         2nd of 4 subcases will follow an explanation of how to produce  vector
C         plots of statistical distributions in Apollo window using input of the
C         preceding solution.  See   "STATPLOT"  as described in Section I-F-1-g
C         of the rule book.   Execute the following in a 2/3-height window  with
C         font  F7x13.B  and nearly full width:
C                $ STATPLOT.EXE
C                Y 5  .9  1.9
C                HELP
C                GO
C                  <<<  Pasted data beginning with 1st row of the
C                  <<<  statistical tabulation of peak energy for
C                  <<<  branch  "B1    "  to  "BSW1  ".   This begins
C                  <<<  in row 10 and extends through row 20,  followed
C                  <<<  4 lines of summary statistics.   See Rule Book.
C                LINE  1.4  0.1
C                LINE  1.8  0.1
C                LINE  1.0  1.0
C         ======================================================================
C         On  to  the 2nd subcase,  which will illustrate statistical tabulation
C         of previously-generated energizations without the addition of any  new
C         solutions.   It also illustrates the combination of 2 sets of extrema,
C         with the 1st coming from the 1st subcase  of  DC-24  (just as with the
C         preceding data subcase)  and  the  2nd  coming from the 3rd subcase of
C         DC-24.   This illustrates the use of  "LOAD NEXT PART"  that announces
C         or requests concatenation of the 2nd set of extrema.   Finally,  there
C         is illustration of the erasure of energizations from each set.   Since
C         both sets of data are really identical,  killing one extrema in 1 file
C         and the remaining 2 in the other,  in any combination,  will result in
C         3 remaining extrema that should tabulate identically to the 3 of DC-24
$OPEN,  UNIT=LUNIT2  FILE=dc24.bin  STATUS=OLD  FORM=UNFORMATTED ! { Hold case
START AGAIN            { Load hibernating solution (except for separate extrema)
    9999               { Terminate modifications (none, here) to switching times
C Cancel special output for MS Excel of DC-24 tables.  Addition on 21 July 2007:
C                           Type  ---  Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL    .off        { Columns 29-32 carry optional file type
$OPEN, UNIT=LUNIT2  FILE=dc24at40.ext  STATUS=OLD  FORM=UNFORMATTED  RECL=16000  !
C 345678901234567890 KIL KIL KIL KIL  ---  Ruler for  (20X, 15I4)  of next card:
LOAD MORE SHOTS        2             { Load extrema of DC-24 minus shot number 2
TOGGLE EXTREMA MODE        { Change from default UNFORMATTED mode to card images
$OPEN, UNIT=LUNIT2  FILE=dc24ct40.lis  STATUS=OLD  FORM=FORMATTED !
LOAD MORE SHOTS        1   3        { Load additional history, minus shots 1 & 3
C         Following card means next case will have  L2FORM = 0  within  "REQUES"
C         whether code is overlaid or not.   Otherwise,  virtual computers would
C         have  kept  L2FORM = 1  whereas overlaid would have reset  L2FORM = 0.
TOGGLE EXTREMA MODE  { Restore the default UNFORMATTED mode of "LOAD MORE SHOTS"
MISCELLANEOUS DATA CARDS
 100.E-6  20.E-3     60.
       1       1       1       1       1               0       0       3
C    ISW   ITEST   IDIST    IMAX   IDICE  KSTOUT  KNTRPT                   NSEED
       1       1       0       0       1         { KSTOUT is blank, not 0! }   1
TIME STEP LOOP   { Transfer control to overlay 16 for a resumption of simulation
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
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}
C    -----------  Following is 2nd mini-switch plot, which is the 1st Gaussian :
C   Switch  "B1    "  to  "BSW1  "
C     1.32420007E-02   switch closings per column
C              Columns
C          5   10   15   20   25   30
C     +----+----+----+----+----+----+
C     X
C     X
C     X
C     X
C     A*
C     A *
C     A   *
C     A       *
C     A           *
C     A                 *
C     A                       *
C     A                           *
C     |  4.0000E-03                 *              Pulled 3 columns left --->  A
C     |                           *                Pulled 3 columns left --->  A
C     A                       *
C     A                 *
C     A           *
C     |       *                                    Pulled 3 columns left --->  A
C     A   *
BLANK card ending statistical tabulation requests
BEGIN NEW DATA CASE
C         3rd  of 4 total subcases, which will illustrate statistical tabulation
C         of previously-generated energizations without the addition of any  new
C         solutions.   It is like 2nd, only with a different erasure of 3 shots.
C         Tabulations are the same.  Finally, it illustrates the omission of the
C         miniature character plots of switch closing times via NO SWITCH PLOTS.
$OPEN,  UNIT=LUNIT2  FILE=dc24.bin  STATUS=OLD  FORM=UNFORMATTED ! { Hold case
START AGAIN            { Load hibernating solution (except for separate extrema)
    9999               { Terminate modifications (none, here) to switching times
C Cancel special output for MS Excel of DC-24 tables.  Addition on 21 July 2007:
C                           Type  ---  Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL    .off        { Columns 29-32 carry optional file type
$OPEN, UNIT=LUNIT2  FILE=dc24at40.ext  STATUS=OLD  FORM=UNFORMATTED  RECL=16000  !
C 345678901234567890 KIL KIL KIL KIL  ---  Ruler for  (20X, 15I4)  of next card:
LOAD MORE SHOTS        3             { Load extrema of DC-24 minus shot number 3
TOGGLE EXTREMA MODE        { Change from default UNFORMATTED mode to card images
$OPEN, UNIT=LUNIT2  FILE=dc24ct40.lis  STATUS=OLD  FORM=FORMATTED  !
LOAD MORE SHOTS        1   2        { Load additional history, minus shots 1 & 2
C         Following card means next case will have  L2FORM = 0  within  "REQUES"
C         whether code is overlaid or not.   Otherwise,  virtual computers would
C         have  kept  L2FORM = 1  whereas overlaid would have reset  L2FORM = 0.
TOGGLE EXTREMA MODE  { Restore the default UNFORMATTED mode of "LOAD MORE SHOTS"
MISCELLANEOUS DATA CARDS
 100.E-6  20.E-3     60.
       1       1       1       1       1               0       0       3
C    ISW   ITEST   IDIST    IMAX   IDICE  KSTOUT  KNTRPT                   NSEED
       1       1       0       0       1         { KSTOUT is blank, not 0! }   1
TIME STEP LOOP   { Transfer control to overlay 16 for a resumption of simulation
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
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 4 total subcases will illustrate the repeated simulation of the
C         already-solved 2nd energization.   By definition,  only  a single shot
C         will occur, and by means of  KSTOUT = 0  we can look at time-step loop
C         output just as though there were no  "STATISTICS".   Such  output  now
C         can be followed by plotting (activated on  30 March 1989).   Note that
C         the plot file header was created by the 2nd subcase of  DC-24.   It is
C         necessary to connect this since only time-step-loop output is produced
C         by  "START AGAIN"  usage such as here.   Miscellaneous  data parameter
C         ICAT = 2  will save the ".PL4" plot file,  so subsequent,  independent
C         plotting is possible  (and is illustrated by the 2nd plot of the first
C         subcase of DC-54).  There is no statistical tabulation for the present
C         subcase, of course.  Variable  KNTRPT = 2  repeats the 2nd shot of the
C         previously-solved 1st subcase of DC-24.  The zero (not blank!)  KSTOUT
C         provides desired time-step loop output.
$OPEN,  UNIT=LUNIT2  FILE=dc24.bin  STATUS=OLD  FORM=UNFORMATTED ! { Hold case
START AGAIN                      { Request for awakening of hibernating solution
    9999               { Terminate modifications (none, here) to switching times
C Cancel special output for MS Excel of DC-24 tables.  Addition on 21 July 2007:
C                           Type  ---  Ruler for file type of Excel output file:
EXTREMA OUTPUT FOR EXCEL    .off        { Columns 29-32 carry optional file type
$OPEN, UNIT=LUNIT2  FILE=dc24at40.ext  STATUS=OLD  FORM=UNFORMATTED  RECL=16000  !
C 345678901234567890 KIL KIL KIL KIL  ---  Ruler for  (20X, 15I4)  of next card:
LOAD MORE SHOTS                  { Load extrema produced by 1st subcase of DC-24
MISCELLANEOUS DATA CARDS
 100.E-6  20.E-3     60.
       1       1       1       1       1               0       2       3
C               Preceding KNTRPT = 3 was changed from former 6 on 23 Sept 1998.
C    ISW   ITEST   IDIST    IMAX   IDICE  KSTOUT  KNTRPT                   NSEED
       1       1       0       0       1       0       2                       1
CHANGE PRINTOUT FREQUENCY
       1       1      20      10
$CLOSE,  UNIT=LUNIT4  STATUS=DELETE   { Destroy empty date/time plot file of "SYSDEP"
$OPEN,  UNIT=LUNIT4  FILE=dc24b40d.pl4  STATUS=OVERLAY  ! { C-like uses only header
C  READ PL4 HEADER     { Messes up C-like usage of PL4
C Transfer control to the time-step loop.           |TIME STEP LOOP   { Transfer
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  2  :
C                       1  1.8999095E-03    2  4.7979502E-03    3  7.5258324E-03
C   18   .0018   .284217E-13  -.26148E-11  .284217E-12          0.0          0.0
C                .192639E-17  .966554E-16
C  ***         Close switch  "A1    "  to  "ASW1  "  after  1.90000000E-03  sec.
C   19   .0019   -.28422E-13   .26148E-11  -.28422E-12          0.0          0.0
C                .210598E-16  .966554E-16
TIME STEP LOOP   { Transfer control to overlay 16 for a resumption of simulation
C  200     .02   18.40117142   -5.2655535  13.97952075  79.19748583  128.2703628
C                .4457755752  -1.28521406
C Variable max : 141.9164096  61.55446116  25.22399908  371.6253603  334.5616569
C                1.207725131  .5125980517
C Times of max :        .002        .0154        .0132        .0059        .0028
C                      .0058        .0105
C Variable min : -103.914014  -46.8051457  -41.9189588  -288.448852  -318.274018
C                -1.09393135  -1.83695142
C Times of min :       .0032        .0144        .0139        .0071        .0088
C                      .0071         .017
  PRINTER PLOT
 144 1. 0.0 10.         ASW10                    { Plot limits: (-3.183,  3.346)
BLANK card ending batch-mode plot cards
BEGIN NEW DATA CASE
BLANK