path: root/benchmarks/dc3.dat

BEGIN NEW DATA CASE
$BEGIN PL4 COMMENTS
C     BENCHMARK DC-3                                                            
C              Energization of 180-mile 3-phase line represented by 18 identical
C              Pi-sections.  Transposed at 60 and 120 miles, note.  XOPT = 3000.
c              Note nice  DELTAT.   For exact "M39." agreement,  use   0.4630E-4
c              It is test case DC-54 that uses the plot file for "REPLOT" usage.
c              Note this and preceding 2 comment cards illustrate lower-case "$$
C c" (the third letter of the alphabet).   End comment text.
$END PL4 COMMENTS
CHANGE PLOT FREQUENCY
C  Beginning  21 March 2010,  the ten integers that define either plotting or
C  printing frequency may be keyed as a mixture of I8  and  E8.0  FORMAT.  There
C  exist several special cases, however,  for large values.  Interpretation will
C  continue to use 6I6 integers if all can be encoded this way. Examples follow:
       5       5      10       1     { The original alternative has all integers
C Interpretation of this: Plot out :      5     5    10     1     0     0   | ...
C      5    5.E0      10       1    1.E4    9.E5   { Near the limit of interpret
C Interpretation of this: Plot out :      5     5    10     1 10000900000   | ...
C      5    5.E0      10       1    1.E4    1.E9 { Too big a value for interpret
C Interpretation of this: 5 pairs.  But a value is too big for I6 display.  | ...
C      5    5.E0      10       1    1.E4    2.E9 { Verify largest 32-bit integer
C Interpretation of this: 5 pairs.  But a value is too big for I6 display.  | ...
C      5    5.E0      10       1    1.E4   2.1E9   { Too big a value for 32 bits
C rejection of this:  ===  Error defining one or more integers of printing or
C    plotting frequency.  Too big for 32-bit integer storage.  Halt in SUBR1.
CUSTOM PLOT FILE     { Request for REAL*8 .PL4 file (see July, 1995, newsletter)
$MONITOR      { Request time (no other special effect, since not LUNIT6 to disk)
C DISK PLOT DATA   { Toggle the Apollo default of LUNIT4 = -4  to  +4 (use disk)
C $CLOSE,  UNIT=4  STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP"
C $OPEN,  UNIT=4  FILE=dc3to54.pl4  FORM=UNFORMATTED STATUS=UNKNOWN RECL=8000  !
C$ text   Beginning 31 August 1994, "C$" in columns 1-2 is taken to be a comment
PRINTED NUMBER WIDTH, 13, 2,
 .000050    .010   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
C      1       1       1       1       1      -1       0       2                
    1.E0    1.E0       1       1       1      -1       0       2                
C       The preceding integer misc. data card illustrates that either the number
C       of steps for printing (cols. 1-8) or the number of steps for plotting
C       (cols. 9-16) may be E-field.  This is to accomodate values > 1.E8.
C   See previous explanation of  CHANGE PLOT FREQUENCY.  The same procedure is
C   applicable to the following card,  which modifies printout frequency:
C      5       5      10      10      20      20 { Escalating printout frequency
       5    5.E0    1.E1      10    2.E1      20   { Illustrate some E-field use
$COMMENT                { Begin discarding (making invisible) all comment cards
C     Ok, immediately we have a comment card, which should not be seen on LUNIT6
 1GEN-A 1-A               34.372457.68.15781                                    
 2GEN-B 1-B               35.735164.43-.031538.002451.79.16587                  
 3GEN-C 1-C               35.735164.43-.031537.455151.72-.021938.002451.79.16587
C                     2nd comment card.  This, too, should not be seen on LUNIT6
$COMMENT   { End discarding of comment cards.  Once again such cards are visible
C                     Is this comment card visible?  It should be seen on LUNIT6
 11-A   2-A   GEN-A 1-A          { Sections 2 through 18 are copies of the first
 21-B   2-B                      { which has just been inputted.                
 31-C   2-C                                                                     
 12-A   3-A   GEN-A 1-A                                                         
 22-B   3-B                                                                     
 32-C   3-C                                                                     
 13-A   4-A   GEN-A 1-A                                                         
 23-B   4-B                                                                     
 33-C   4-C                                                                     
 14-A   5-A   GEN-A 1-A                                                         
 24-B   5-B                                                                     
 34-C   5-C                                                                     
 15-A   6-A   GEN-A 1-A                                                         
 25-B   6-B                                                                     
 35-C   6-C                                                                     
 16-C   7-C   GEN-A 1-A      { Note transposition:  /C/A/B/  rather than  /A/B/C
 26-A   7-A                                                                     
 36-B   7-B                                                                     
 17-C   8-C   GEN-A 1-A                                                         
 27-A   8-A                                                                     
 37-B   8-B                                                                     
 18-C   9-C   GEN-A 1-A                                                         
 28-A   9-A                                                                     
 38-B   9-B                                                                     
 19-C   10-C  GEN-A 1-A                                                         
 29-A   10-A                                                                    
 39-B   10-B                                                                    
 110-C  11-C  GEN-A 1-A                                                         
 210-A  11-A                                                                    
 310-B  11-B                                                                    
 111-C  12-C  GEN-A 1-A                                                         
 211-A  12-A                                                                    
 311-B  12-B                                                                    
 112-B  13-B  GEN-A 1-A  { Note 2nd transposition:  /B/C/A/  rather than  /C/A/B
 212-C  13-C                                                                    
 312-A  13-A                                                                    
 113-B  14-B  GEN-A 1-A                                                         
 213-C  14-C                                                                    
 313-A  14-A                                                                    
 114-B  15-B  GEN-A 1-A                                                         
 214-C  15-C                                                                    
 314-A  15-A                                                                    
 115-B  16-B  GEN-A 1-A                                                         
 215-C  16-C                                                                    
 315-A  16-A                                                                    
 116-B  17-B  GEN-A 1-A                                                         
 216-C  17-C                                                                    
 316-A  17-A                                                                    
 117-B  18-B  GEN-A 1-A                                                         
 217-C  18-C                                                                    
 317-A  18-A                                                                    
 0M-A   GEN-A              400.0   { 400 Ohm closing resistors, to be shorted by
 0M-B   GEN-B              400.0   { breaker poles at times  9.98,  14,  and  14
 0M-C   GEN-C              400.0   { msec, respectively.}                      1
 0POLE-AM-A                15.0                                                 
 0POLE-BM-B                15.0                                                 
 0POLE-CM-C                15.0                                                 
BLANK card ending branch cards
  E-A   POLE-A    0.      20.0                                                 1
  E-B   POLE-B  0.00398   20.0  { Closing will be at 4.0 msec, all computer }  3
  E-C   POLE-C  0.00398   20.0  { This backoff from 4.0 was needed by PRIME }  1
  M-A   GEN-A   0.00998   20.0            { Breaker poles across 400 Ohm closing
  M-B   GEN-B   0.013998  20.0            { resistors.   Note artificial opening
  M-C   GEN-C   0.013998  20.0            { time (in fact, there is no opening).
BLANK card ending switches
14E-A        -1.0      60.0    -90.0                                            
14E-B        -1.0      60.0   -210.0                                            
14E-C        -1.0      60.0    30.0                                             
BLANK card ending sources
  18-C  18-B  18-A                                                              
$MONITOR      { Request time (no other special effect, since not LUNIT6 to disk)
C  Step      Time        E-B          18-C         18-B         18-A         E-A
C                     POLE-B                                              POLE-A
C ***        Switch  "E-A   "  to  "POLE-A"  closed  after  0.00000000E+00  sec.
C    0     0.0           0.0          0.0          0.0          0.0          0.0
C    1   .5E-4   .8771390038  .735053E-17  .754138E-17  -.12623E-16  -.277151E-4
C    2   .1E-3   .8898949351  .408898E-15  .419476E-15  -.70276E-15   -.42636E-4
C    3  .15E-3   .9023577997  .111958E-13  .114844E-13  -.19259E-13  -.657706E-4
C    4   .2E-3   .9141912694  .201088E-12  .206252E-12  -.34623E-12  -.902274E-4
C    5  .25E-3   .9254563554  .266406E-11  .273221E-11  -.45918E-11   -.11094E-3
C   10   .5E-3   .9821092355  .5616974E-7  .5757735E-7  -.975638E-7  -.211746E-3
C   20    .001   1.079796132  .0030080496  .0030756136  -.005532677  -.393566E-3
BLANK card ending output variables requests (node voltages, here)
$MONITOR      { Request time (no other special effect, since not LUNIT6 to disk)
C Final step:      200     .01     0.0   .7226176276   -.887792385   .0782571591
C Final step continued:       .0011948668  -.449544E-3  -.687051E-3  -.687051E-3
C Variable max : 1.188718074  .8796065581  .1512528782  .0782571591  .0012022221
C Times of max :        .002        .0085       .00525          .01       .00995
C Variable min :         0.0  -.133205909  -.887792385  -1.27643205   -.79278E-3
C Times of min :         0.0        .0048          .01       .00525        .0041
  PRINTER PLOT                                                                  
 144 1. 0.0 10.         18-A  18-B  18-C         { Axis limits: (-1.276,  0.880)
 194 1. 0.0 10.         GEN-C M-C   E-A   POLE-A { Axis limits: (-1.105,  1.202)
BLANK card ending plot cards
BEGIN NEW DATA CASE                                                             
C      2nd of 6 subcases illustrates  $PI  which becomes available 14 July 2001.
C      This is for scaling Pi-circuit data  --  in effect, varying the length of
C      the Pi-circuits.  Originator of the idea was Dr. Sven Demmig of Bewag
C      (power utility of city of Berlin, Germany) as explained in the October,
C      2001 newsletter.  Although a cascade connection of sections is typical,
C      in fact there is no mandatory association of $PI and cascading.  In this
C      illustration,  consider 4 different and disconnected cascade connections
C      of 2 Pi-circuits.  The first is unscaled whereas the following 3 involve
C      different combinations of scaling.  All 4 disconnected subnetworks give
C      exactly the same answer because data is either normal or multiplied by
C      2;  and if the later, it subsequently is scaled by 1/2, thereby restoring
C      the original values.
POWER FREQUENCY                    3000.   { Avoid warning about unusual XOPT
PRINTED NUMBER WIDTH, 11, 2,
 .000050    .005   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
       1       1       1       1       1      -1                
       5       5      10      10      20      20 { Escalating printout frequency
C   Begin cascading of 3-phase Pi-circuits from DC-3.  To establish the right 
C   answer,  we define data params using the 1st section, and then use reference
C   branch for the second.  So, a total of 6 cells are occupied in List 3:
 1GENA  MIDA              34.372457.68.15781                                    
 2GENB  MIDB              35.735164.43-.031538.002451.79.16587                  
 3GENC  MIDC              35.735164.43-.031537.455151.72-.021938.002451.79.16587
 1MIDA  ENDA  GENA  MIDA 
 2MIDB  ENDB                             
 3MIDC  ENDC                                                                     
C   Preceding establishes right answer without scaling.  Next, scale.  There are  
C   3 illustrations.  Each is a disconnected circuit that should give the same 
C   answer as the first.  Each of the copies uses scaling in a different way.
C   Begin with Pi-circuit data that is double what it should be, so will require
C   a scaling factor of a half. Like preceding, this occupies 6 cells of List 3:
$PI, 0.5,            { Multiply R, L, and C values of Pi-circuits by this factor
 1GENA  MIDD              68.744915.36.31562                                    
 2GENB  MIDE              71.470328.86-.063076.004903.58.33174                  
 3GENC  MIDF              71.470328.86-.063074.910303.44-.043876.004903.58.33174
$PI, 1.0,            { Cancel preceding scaling;  return to normal, unscaled use
 1MIDD  ENDD  GENA  MIDD 
 2MIDE  ENDE                             
 3MIDF  ENDF                                                                     
C   2nd of 3 alternatives to the original circuit will define but not use the
C   doubled data.  Nothing else will be connected to far end (MIDG, H, I).
C   Instead,  the alternative will consist of 2 scaled copies of this data.
C   That makes a total of 3 copies of data, so 18 cells of List 3 are occupied:
 1GENA  MIDG              68.744915.36.31562
 2GENB  MIDH              71.470328.86-.063076.004903.58.33174                  
 3GENC  MIDI              71.470328.86-.063074.910303.44-.043876.004903.58.33174
$PI, 0.5,            { Multiply R, L, and C values of Pi-circuits by this factor
 1GENA  MIDJ  GENA  MIDG 
 2GENB  MIDK                             
 3GENC  MIDL                                                                     
 1MIDJ  ENDJ  GENA  MIDG 
 2MIDK  ENDK                             
 3MIDL  ENDL                                                                     
C   3rd of 3 alternatives begins (1st half) the same as the preceding second,
C   but it ends with the 2nd section a copy of this 1st section.  This adds 6  
C   more cells for List 3 (for the 1st section only).  So, the total burden on
C   List 3 is the sum of the 4 parts: 6 + 6 + 18 + 6 = 36 cells.
 1GENA  MIDM  GENA  MIDG 
 2GENB  MIDN                             
 3GENC  MIDO                                                                     
$PI, 1.0,            { Cancel preceding scaling;  return to normal, unscaled use
 1MIDM  ENDM  GENA  MIDM 
 2MIDN  ENDN                             
 3MIDO  ENDO                                                                     
BLANK card ending branch cards
BLANK card ending switches
C 14GENA        -1.0      60.0     -90.0                                            
C   Note we omit phase "a" (preceding source) to produce great imbalance.  Also,
C   we will add a steady-state solution, for more generality.  As a result, the
C   simulation should be pursely sinusoidal.  If  PRINTER PLOT  below is changed
C   to  CALCOMP PLOT  &  HPI = 1.  msec/inch  is changed to  .5  (columns 6-7),
C   this will be obvious visually.
14GENB          -1.0      60.0    -210.0                           -1.
14GENC          -1.0      60.0      30.0                           -1.                  
$WIDTH, 79,   { Request narrow, 80-column LUNIT6 output to minimize phasor print
BLANK card ending sources
C     Total network loss  P-loss  by summing injections =   5.813026183430E-07
C GENB                  1.0   .6605121612E-3  -.4572787068E-4   .3302560806E-3
C               -30.0000000       67.9588504  -.3270749771E-3       -0.1384619
$WIDTH, 132, { Done with 80 columns, so return to wide, 132-column LUNIT6 output
  ENDA  ENDD  ENDJ  ENDM  ENDB  ENDE  ENDK  ENDN  ENDC                                                            
C  Step     Time      ENDA       ENDD       ENDJ       ENDM       ENDB       ENDE       ENDK       ENDN       ENDC
C    0       0.0   -.0023158  -.0023158  -.0023158  -.0023158  .86301625  .86301625  .86301625  .86301625  -.8705073
C    1     .5E-4   .00145168  .00145168  .00145168  .00145168  .87229895  .87229895  .87229895  .87229895  -.8609192
C    2     .1E-3   .00521861  .00521861  .00521861  .00521861  .88127179  .88127179  .88127179  .88127179  -.8510252
C    3    .15E-3   .00898371  .00898371  .00898371  .00898371  .88993164  .88993164  .88993164  .88993164  -.8408287
BLANK card ending output variables requests (node voltages, here)
C  100      .005   .19079849  .19079849  .19079849  .19079849  .20945009  .20945009  .20945009  .20945009  .74501057
C Variable maxima: .19987826  .19987826  .19987826  .19987826  .99771368  .99771368  .99771368  .99771368  .74501057
C Times of maxima:     .0042      .0042      .0042      .0042      .0014      .0014      .0014      .0014       .005
C Variable minima: -.0023158  -.0023158  -.0023158  -.0023158  .20945009  .20945009  .20945009  .20945009  -.8705073
C Times of minima:       0.0        0.0        0.0        0.0       .005       .005       .005       .005        0.0
  PRINTER PLOT    { Conversion to CALCOMP PLOT would show that curves are smooth
 144 1. 0.0 5.0 -1. 1.0 ENDA  ENDB  ENDC
BLANK card ending plot cards
BEGIN NEW DATA CASE                                                             
C      3rd of 6 subcases illustrates  HIGH ORDER PI CIRCUIT  as first mentioned
C      in the April, 1998, newsletter.   Prior to addition on  22 November 2001,
C      this data was stored at home as  f:\data\highpi.dat    Data of this
C      subcase is identical to the first;  the solution should be identical.
C      But instead of data for the 1st Pi-circuit being read from in-line card
C      images, it is to be read from a separate disk file.  Of course, for 3
C      phases,  this is artificially complicated.  But for really high order
C      (e.g., 400), there should be a big saving when the C-like alternative is
C      used.  But here, initially, we just illustrate the universal, FORMATTED
C      alternative,  which stores Pi-circuit data in separate  DC3HIGH.DAT  The
C      use is artificial.  Practical use almost always will be C-like.
PRINTED NUMBER WIDTH, 13, 2,
 .000050    .010   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
       1       1       0       0       0      -1                
       5       5      10      10      20      20 { Escalating printout frequency
C 1GEN-A 1-A               34.372457.68.15781
C 2GEN-B 1-B               35.735164.43-.031538.002451.79.16587                  
C 3GEN-C 1-C               35.735164.43-.031537.455151.72-.021938.002451.79.16587
C   The following is FORMATTED.  For the C-like alternative,  see the 4th subcase.
C HIGH ORDER PI CIRCUIT        3  File=[]dc3high.dat  FORM=FORMATTED C-LIKE OUTPUT
C     Preceding single card is replaced by following 2-card alternative on
C     9 September 2003.  Not only was  C-LIKE OUTPUT  being chopped on the
C     right,  there was no space for additional attributes such as  REAL*4
C     for single-precision data  and  $UNITS  to apply XOPT and COPT.  So now
C     2 or more cards are allowed.  Only the last is to carry the file name
C     after the  File=  tag.  Tags on the optional preface card are free-
C     format but they must be to the right of column 32 (for # of phases):
  HIGH ORDER PI CIRCUIT           FORM=FORMATTED   C-LIKE OUTPUT { Optional card
  HIGH ORDER PI CIRCUIT        3  File=[]dc3high.dat 
C     Note that preceding creates REAL*8 output file  dc3high.clk  by default.
C     The alternative single-precision file  dc3high.444  which will be used in
C     the following subcase # 5 was created by adding  REAL*4  anywhere to the
C     right of the input data file name (separate by 1 or more blanks).  The
C     resulting file then will be 108 bytes in size rather than 180.  The
C     difference is  3 arrays * 6 elements/array * 4 bytes/element = 72 bytes.
C     But 5th subcase also used  $UNITS  and this requires a change in the code
C     see comments near the bottom of  SUBR3.SPL
 11-A   2-A   GEN-A 1-A          { Sections 2 through 18 are copies of the first
 21-B   2-B                      { which has just been inputted.                
 31-C   2-C                                                                     
 12-A   3-A   GEN-A 1-A                                                         
 22-B   3-B                                                                     
 32-C   3-C                                                                     
 13-A   4-A   GEN-A 1-A                                                         
 23-B   4-B                                                                     
 33-C   4-C                                                                     
 14-A   5-A   GEN-A 1-A                                                         
 24-B   5-B                                                                     
 34-C   5-C                                                                     
 15-A   6-A   GEN-A 1-A                                                         
 25-B   6-B                                                                     
 35-C   6-C                                                                     
 16-C   7-C   GEN-A 1-A      { Note transposition:  /C/A/B/  rather than  /A/B/C
 26-A   7-A                                                                     
 36-B   7-B                                                                     
 17-C   8-C   GEN-A 1-A                                                         
 27-A   8-A                                                                     
 37-B   8-B                                                                     
 18-C   9-C   GEN-A 1-A                                                         
 28-A   9-A                                                                     
 38-B   9-B                                                                     
 19-C   10-C  GEN-A 1-A                                                         
 29-A   10-A                                                                    
 39-B   10-B                                                                    
 110-C  11-C  GEN-A 1-A                                                         
 210-A  11-A                                                                    
 310-B  11-B                                                                    
 111-C  12-C  GEN-A 1-A                                                         
 211-A  12-A                                                                    
 311-B  12-B                                                                    
 112-B  13-B  GEN-A 1-A  { Note 2nd transposition:  /B/C/A/  rather than  /C/A/B
 212-C  13-C                                                                    
 312-A  13-A                                                                    
 113-B  14-B  GEN-A 1-A                                                         
 213-C  14-C                                                                    
 313-A  14-A                                                                    
 114-B  15-B  GEN-A 1-A                                                         
 214-C  15-C                                                                    
 314-A  15-A                                                                    
 115-B  16-B  GEN-A 1-A                                                         
 215-C  16-C                                                                    
 315-A  16-A                                                                    
 116-B  17-B  GEN-A 1-A                                                         
 216-C  17-C                                                                    
 316-A  17-A                                                                    
 117-B  18-B  GEN-A 1-A                                                         
 217-C  18-C                                                                    
 317-A  18-A                                                                    
 0M-A   GEN-A              400.0   { 400 Ohm closing resistors, to be shorted by
 0M-B   GEN-B              400.0   { breaker poles at times  9.98,  14,  and  14
 0M-C   GEN-C              400.0   { msec, respectively.}                      1
 0POLE-AM-A                15.0                                                 
 0POLE-BM-B                15.0                                                 
 0POLE-CM-C                15.0                                                 
BLANK card ending branch cards
  E-A   POLE-A    0.      20.0                                                 1
  E-B   POLE-B  0.00398   20.0  { Closing will be at 4.0 msec, all computer }  3
  E-C   POLE-C  0.00398   20.0  { This backoff from 4.0 was needed by PRIME }  1
  M-A   GEN-A   0.00998   20.0            { Breaker poles across 400 Ohm closing
  M-B   GEN-B   0.013998  20.0            { resistors.   Note artificial opening
  M-C   GEN-C   0.013998  20.0            { time (in fact, there is no opening).
BLANK card ending switches
14E-A        -1.0      60.0    -90.0                                            
14E-B        -1.0      60.0   -210.0                                            
14E-C        -1.0      60.0    30.0                                             
BLANK card ending sources
  18-C  18-B  18-A                                                              
BLANK card ending output variables requests (node voltages, here)
C  100      .005       0.0  -.084110964  -.076310793   -.92825821  -.621001E-3  .1845237E-3  .9370288E-3  .9370288E-3
C  120      .006       0.0  .4539217657  .0096691741  -1.01404035   -.10259E-3   -.52007E-3  .0010499755  .0010499755
C  140      .007       0.0  .6915082663   -.20985407  -.843083014  .6173076E-3  -.951807E-3  .5938305E-3  .5938305E-3
C  160      .008       0.0  .8512280695  -.465621295  -.605281204  .9350157E-3   -.96575E-3  .1846072E-3  .1846072E-3
C  180      .009       0.0  .8583902573   -.70840222  -.288190483  .0011711805  -.794874E-3  -.285103E-3  -.285103E-3
C ***        Close switch  "M-A   "  to  "GEN-A "  after  1.00000000E-02  sec.
C  200       .01       0.0  .7226176276  -.887792385  .0782571591  .0011948668  -.449544E-3  -.687051E-3  -.687051E-3
  PRINTER PLOT                                                                  
BLANK card ending plot cards
BEGIN NEW DATA CASE                                                             
C      4th of 6 subcases is the same as the preceding 3rd subcase except for
C      two changes.  First and most importantly,  the C-like alternative is used 
C      for the data rather than the FORMATTED alternative.  Second, the ending
C      time  T-max  has been halved.  Solution through  T = T-max  is identical.
PRINTED NUMBER WIDTH, 13, 2,
 .000050    .005   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
       1       1       0       0       0      -1                
       5       5      10      10      20      20 { Escalating printout frequency
C 1GEN-A 1-A               34.372457.68.15781
C 2GEN-B 1-B               35.735164.43-.031538.002451.79.16587                  
C 3GEN-C 1-C               35.735164.43-.031537.455151.72-.021938.002451.79.16587
  HIGH ORDER PI CIRCUIT        3  File=[]dc3high.clk  { Use default type: C-like 
 11-A   2-A   GEN-A 1-A          { Sections 2 through 18 are copies of the first
 21-B   2-B                      { which has just been inputted.                
 31-C   2-C                                                                     
 12-A   3-A   GEN-A 1-A                                                         
 22-B   3-B                                                                     
 32-C   3-C                                                                     
 13-A   4-A   GEN-A 1-A                                                         
 23-B   4-B                                                                     
 33-C   4-C                                                                     
 14-A   5-A   GEN-A 1-A                                                         
 24-B   5-B                                                                     
 34-C   5-C                                                                     
 15-A   6-A   GEN-A 1-A                                                         
 25-B   6-B                                                                     
 35-C   6-C                                                                     
 16-C   7-C   GEN-A 1-A      { Note transposition:  /C/A/B/  rather than  /A/B/C
 26-A   7-A                                                                     
 36-B   7-B                                                                     
 17-C   8-C   GEN-A 1-A                                                         
 27-A   8-A                                                                     
 37-B   8-B                                                                     
 18-C   9-C   GEN-A 1-A                                                         
 28-A   9-A                                                                     
 38-B   9-B                                                                     
 19-C   10-C  GEN-A 1-A                                                         
 29-A   10-A                                                                    
 39-B   10-B                                                                    
 110-C  11-C  GEN-A 1-A                                                         
 210-A  11-A                                                                    
 310-B  11-B                                                                    
 111-C  12-C  GEN-A 1-A                                                         
 211-A  12-A                                                                    
 311-B  12-B                                                                    
 112-B  13-B  GEN-A 1-A  { Note 2nd transposition:  /B/C/A/  rather than  /C/A/B
 212-C  13-C                                                                    
 312-A  13-A                                                                    
 113-B  14-B  GEN-A 1-A                                                         
 213-C  14-C                                                                    
 313-A  14-A                                                                    
 114-B  15-B  GEN-A 1-A                                                         
 214-C  15-C                                                                    
 314-A  15-A                                                                    
 115-B  16-B  GEN-A 1-A                                                         
 215-C  16-C                                                                    
 315-A  16-A                                                                    
 116-B  17-B  GEN-A 1-A                                                         
 216-C  17-C                                                                    
 316-A  17-A                                                                    
 117-B  18-B  GEN-A 1-A                                                         
 217-C  18-C                                                                    
 317-A  18-A                                                                    
 0M-A   GEN-A              400.0   { 400 Ohm closing resistors, to be shorted by
 0M-B   GEN-B              400.0   { breaker poles at times  9.98,  14,  and  14
 0M-C   GEN-C              400.0   { msec, respectively.}                      1
 0POLE-AM-A                15.0                                                 
 0POLE-BM-B                15.0                                                 
 0POLE-CM-C                15.0                                                 
BLANK card ending branch cards
  E-A   POLE-A    0.      20.0                                                 1
  E-B   POLE-B  0.00398   20.0  { Closing will be at 4.0 msec, all computer }  3
  E-C   POLE-C  0.00398   20.0  { This backoff from 4.0 was needed by PRIME }  1
  M-A   GEN-A   0.00998   20.0            { Breaker poles across 400 Ohm closing
  M-B   GEN-B   0.013998  20.0            { resistors.   Note artificial opening
  M-C   GEN-C   0.013998  20.0            { time (in fact, there is no opening).
BLANK card ending switches
14E-A        -1.0      60.0    -90.0                                            
14E-B        -1.0      60.0   -210.0                                            
14E-C        -1.0      60.0    30.0                                             
BLANK card ending sources
  18-C  18-B  18-A                                                              
BLANK card ending output variables requests (node voltages, here)
C   80      .004   .8173610312  -.086429015  -.083753604   -.67093214  -.487946E-3          0.0          0.0  .555112E-18
C  100      .005       0.0  -.084110964  -.076310793   -.92825821  -.621001E-3  .1845237E-3  .9370288E-3  .9370288E-3
BLANK card ending plot cards
BEGIN NEW DATA CASE
C      5th of 6 subcases is the same as the preceding 4th subcase except for
C      two changes.  First, single precision is used for the C-like Pi-circuit
C      data rather than the default REAL*8 alternative.  Second, XOPT and COPT
C      are to be applied to control units just as for any other branch.  Answers
C      are very close (REAL*4 input data change 7th or 8 digit of some numbers).
PRINTED NUMBER WIDTH, 13, 2,
 .000050    .005   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
       1       1       0       0       0      -1                
       5       5      10      10      20      20 { Escalating printout frequency
C 1GEN-A 1-A               34.372457.68.15781
C 2GEN-B 1-B               35.735164.43-.031538.002451.79.16587                  
C 3GEN-C 1-C               35.735164.43-.031537.455151.72-.021938.002451.79.16587
  HIGH ORDER PI CIRCUIT           REAL*4   $UNITS  { Optional card of attributes
C HIGH ORDER PI CIRCUIT           REAL*4           { Optional card of attributes
C HIGH ORDER PI CIRCUIT        3  File=[]dc3high.444   
  HIGH ORDER PI CIRCUIT        3  File=[]dc3high.opt   
 11-A   2-A   GEN-A 1-A          { Sections 2 through 18 are copies of the first
 21-B   2-B                      { which has just been inputted.                
 31-C   2-C                                                                     
 12-A   3-A   GEN-A 1-A                                                         
 22-B   3-B                                                                     
 32-C   3-C                                                                     
 13-A   4-A   GEN-A 1-A                                                         
 23-B   4-B                                                                     
 33-C   4-C                                                                     
 14-A   5-A   GEN-A 1-A                                                         
 24-B   5-B                                                                     
 34-C   5-C                                                                     
 15-A   6-A   GEN-A 1-A                                                         
 25-B   6-B                                                                     
 35-C   6-C                                                                     
 16-C   7-C   GEN-A 1-A      { Note transposition:  /C/A/B/  rather than  /A/B/C
 26-A   7-A                                                                     
 36-B   7-B                                                                     
 17-C   8-C   GEN-A 1-A                                                         
 27-A   8-A                                                                     
 37-B   8-B                                                                     
 18-C   9-C   GEN-A 1-A                                                         
 28-A   9-A                                                                     
 38-B   9-B                                                                     
 19-C   10-C  GEN-A 1-A                                                         
 29-A   10-A                                                                    
 39-B   10-B                                                                    
 110-C  11-C  GEN-A 1-A                                                         
 210-A  11-A                                                                    
 310-B  11-B                                                                    
 111-C  12-C  GEN-A 1-A                                                         
 211-A  12-A                                                                    
 311-B  12-B                                                                    
 112-B  13-B  GEN-A 1-A  { Note 2nd transposition:  /B/C/A/  rather than  /C/A/B
 212-C  13-C                                                                    
 312-A  13-A                                                                    
 113-B  14-B  GEN-A 1-A                                                         
 213-C  14-C                                                                    
 313-A  14-A                                                                    
 114-B  15-B  GEN-A 1-A                                                         
 214-C  15-C                                                                    
 314-A  15-A                                                                    
 115-B  16-B  GEN-A 1-A                                                         
 215-C  16-C                                                                    
 315-A  16-A                                                                    
 116-B  17-B  GEN-A 1-A                                                         
 216-C  17-C                                                                    
 316-A  17-A                                                                    
 117-B  18-B  GEN-A 1-A                                                         
 217-C  18-C                                                                    
 317-A  18-A                                                                    
 0M-A   GEN-A              400.0   { 400 Ohm closing resistors, to be shorted by
 0M-B   GEN-B              400.0   { breaker poles at times  9.98,  14,  and  14
 0M-C   GEN-C              400.0   { msec, respectively.}                      1
 0POLE-AM-A                15.0                                                 
 0POLE-BM-B                15.0                                                 
 0POLE-CM-C                15.0                                                 
BLANK card ending branch cards
  E-A   POLE-A    0.      20.0                                                 1
  E-B   POLE-B  0.00398   20.0  { Closing will be at 4.0 msec, all computer }  3
  E-C   POLE-C  0.00398   20.0  { This backoff from 4.0 was needed by PRIME }  1
  M-A   GEN-A   0.00998   20.0            { Breaker poles across 400 Ohm closing
  M-B   GEN-B   0.013998  20.0            { resistors.   Note artificial opening
  M-C   GEN-C   0.013998  20.0            { time (in fact, there is no opening).
BLANK card ending switches
14E-A        -1.0      60.0    -90.0                                            
14E-B        -1.0      60.0   -210.0                                            
14E-C        -1.0      60.0    30.0                                             
BLANK card ending sources
  18-C  18-B  18-A                                                              
BLANK card ending output variables requests (node voltages, here)
C   80   .004  .8173610397   -.08642902  -.083753607  -.670932147  -.487946E-3          0.0          0.0  .555112E-18
C  100   .005          0.0  -.084110972  -.076310803  -.928258207  -.621001E-3  .1845237E-3  .9370288E-3  .9370288E-3
BLANK card ending plot cards
BEGIN NEW DATA CASE
C      6th of 6 subcases is the same as the preceding except that the master
C      Pi-circuit is defined in user-supplied  SUBROUTINE HOPCOD  rather than
C      in the HOPC disk file.  The file name must be  "File=HOPCOD."  where the
C      following file type (001) thus far is being ignored.  Later, if 2 or more
C      different definitions might be provided in  SUBROUTINE HOPCOD,  the file
C      type will allow selection among 999 of them.   Since 5th subcase uses
C      REAL*4  data whereas here we use REAL*8 data,  expect about 7 digits of
C      agreement in the answers.
PRINTED NUMBER WIDTH, 13, 2,
 .000050    .005   3000.    { XOPT = 3 KHz means reactance in ohms at this freq.
       1       1       0       0       0      -1                
       5       5      10      10      20      20 { Escalating printout frequency
C 1GEN-A 1-A               34.372457.68.15781
C 2GEN-B 1-B               35.735164.43-.031538.002451.79.16587                  
C 3GEN-C 1-C               35.735164.43-.031537.455151.72-.021938.002451.79.16587
  HIGH ORDER PI CIRCUIT        3  File=CODE.001  { Define data in code (not disk)  
 11-A   2-A   GEN-A 1-A          { Sections 2 through 18 are copies of the first
 21-B   2-B                      { which has just been inputted.                
 31-C   2-C                                                                     
 12-A   3-A   GEN-A 1-A                                                         
 22-B   3-B                                                                     
 32-C   3-C                                                                     
 13-A   4-A   GEN-A 1-A                                                         
 23-B   4-B                                                                     
 33-C   4-C                                                                     
 14-A   5-A   GEN-A 1-A                                                         
 24-B   5-B                                                                     
 34-C   5-C                                                                     
 15-A   6-A   GEN-A 1-A                                                         
 25-B   6-B                                                                     
 35-C   6-C                                                                     
 16-C   7-C   GEN-A 1-A      { Note transposition:  /C/A/B/  rather than  /A/B/C
 26-A   7-A                                                                     
 36-B   7-B                                                                     
 17-C   8-C   GEN-A 1-A                                                         
 27-A   8-A                                                                     
 37-B   8-B                                                                     
 18-C   9-C   GEN-A 1-A                                                         
 28-A   9-A                                                                     
 38-B   9-B                                                                     
 19-C   10-C  GEN-A 1-A                                                         
 29-A   10-A                                                                    
 39-B   10-B                                                                    
 110-C  11-C  GEN-A 1-A                                                         
 210-A  11-A                                                                    
 310-B  11-B                                                                    
 111-C  12-C  GEN-A 1-A                                                         
 211-A  12-A                                                                    
 311-B  12-B                                                                    
 112-B  13-B  GEN-A 1-A  { Note 2nd transposition:  /B/C/A/  rather than  /C/A/B
 212-C  13-C                                                                    
 312-A  13-A                                                                    
 113-B  14-B  GEN-A 1-A                                                         
 213-C  14-C                                                                    
 313-A  14-A                                                                    
 114-B  15-B  GEN-A 1-A                                                         
 214-C  15-C                                                                    
 314-A  15-A                                                                    
 115-B  16-B  GEN-A 1-A                                                         
 215-C  16-C                                                                    
 315-A  16-A                                                                    
 116-B  17-B  GEN-A 1-A                                                         
 216-C  17-C                                                                    
 316-A  17-A                                                                    
 117-B  18-B  GEN-A 1-A                                                         
 217-C  18-C                                                                    
 317-A  18-A                                                                    
 0M-A   GEN-A              400.0   { 400 Ohm closing resistors, to be shorted by
 0M-B   GEN-B              400.0   { breaker poles at times  9.98,  14,  and  14
 0M-C   GEN-C              400.0   { msec, respectively.}                      1
 0POLE-AM-A                15.0                                                 
 0POLE-BM-B                15.0                                                 
 0POLE-CM-C                15.0                                                 
BLANK card ending branch cards
  E-A   POLE-A    0.      20.0                                                 1
  E-B   POLE-B  0.00398   20.0  { Closing will be at 4.0 msec, all computer }  3
  E-C   POLE-C  0.00398   20.0  { This backoff from 4.0 was needed by PRIME }  1
  M-A   GEN-A   0.00998   20.0            { Breaker poles across 400 Ohm closing
  M-B   GEN-B   0.013998  20.0            { resistors.   Note artificial opening
  M-C   GEN-C   0.013998  20.0            { time (in fact, there is no opening).
BLANK card ending switches
14E-A        -1.0      60.0    -90.0                                            
14E-B        -1.0      60.0   -210.0                                            
14E-C        -1.0      60.0    30.0                                             
BLANK card ending sources
  18-C  18-B  18-A                                                              
BLANK card ending output variables requests (node voltages, here)
C   80   .004   .8173610312  -.086429015  -.083753604   -.67093214  -.487946E-3          0.0          0.0  .555112E-18
C  100   .005           0.0  -.084110964  -.076310793   -.92825821  -.621001E-3  .1845237E-3  .9370288E-3  .9370288E-3
  PRINTER PLOT                                                                  
BLANK ---     extraneous; added here to prove that one or 2 extras are tolerated
BLANK ---     extraneous; added here to prove that one or 2 extras are tolerated
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK