path: root/benchmarks/dc33.dat

BEGIN NEW DATA CASE
C     BENCHMARK DC-33
C      Type-96 hysteretic inductor with TACS logic to give flux by integrating v
C      Batch-mode EMTP CalComp plotting, including  "X-Y PLOT",  is illustrated.
C      For split version with TACS postponed until postprocessing, see DC-45 and
C      DC-46.   As  for data ordering,  this is intentionally shuffled from what
C      would be natural in order to illustrate the use of "/" cards for sorting.
C      Initial condition & load flow classes were added during  November of 1991 
C      after  Prof. Juan Martinez  complained from  Barcelona, Spain,  that such
C      features were not working.   He was right.  One year later, he complained
C      again (it still did not work right.   Somehow, work was lost;  DC-33  had 
C      not been changed;  the optimistic report of the January, 1992, newsletter
C      was incorrect).  So, sorting by class was repaired again on 28 Novemb 92.
C      1st of 9 stacked subcases total.
 .000050    .020
       1       1       1       1       1      -1
       5       5      20      20
TACS HYBRID
/SOURCE
14GEN           377.       60.                                 -1.
/TACS
 1FLUX    +GEN
       1.0
       0.0       1.0
90GEN
91GROUND
99CURR          -1.0*     GROUND
33FLUX  GEN   GROUNDCURR
/OUTPUT
  GEN
/PLOT
C Last step:  400      .02   116.4994069   2.765761424  -2.76576141  .9416257618
C Last step continued .....:               116.4994069  2.765761424  -2.76576142
C Final 5 max. : 3.499909256  .9905492104  376.9925583  3.499613141  3.499909256
C Associated times :   .0125       .00415       .01665       .00415        .0125
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
 194 2. 0.0 20.         TACS  CURR              ABCDEFGHIJKLMNOP1234567890123456
  X-Y PLOT    Horizontal Axis label123
     10.    -1.1     1.1
      8.    -4.0     4.0
 194 4. 0.0 20.         TACS  FLUX  TACS  CURR  Graph heading---Vertical axis la
  X-Y PLOT    Horizontal Axis label123
   9999.    -1.1     1.1     "9999." returns to conventional (vs. time) plotting
      8.    -4.0     4.0
  PRINTER PLOT
 194 4. 0.0 20.         TACS  CURR                { Axis limits: (-3.500, 3.500)
C    Now that all plotting is done with, let's illustrate the three declarations
C    that will choose the destination  of  any  subsequent  "CALCOMP PLOT"  use:
  SCREEN PLOT  { If a vector plot were to follow, it would go only to the screen
  PEN PLOT     { If a vector plot were to follow, ... go only to CalComp plotter
  SCREEN PEN   { If a vector plot were to follow, .. go to both screen & plotter
C         The preceding 3 declarations really could only be tested by a user who
C         had both a vector-graphic screen and a CalComp plotter.  Hence no use.
/LOAD FLOW
C   Comment of /LOAD FLOW.  Since  DC-33  does not really involve the load flow, 
C   Comment of /LOAD FLOW.  we can not use real load flow data cards.   In their 
C   Comment of /LOAD FLOW.  place,  we use these three comment cards.   End test 
/BRANCH
96GROUNDGEN                8888.       1.E-9                                   1
             1.0            -0.7
             2.0             0.9
             3.5             1.0
           9999.
/REQUEST
PRINTED NUMBER WIDTH, 13, 2,  { Request maximum precision (for 8 output columns)
/INITIAL
 2GROUND             0. { Node voltage initial condition in fact changes nothing
C  Initial conditions really are not a part of this problem, although we do want
C  to illustrate that they,  too,  can be sorted,  and inserted after  /LOADFLOW 
C  The preceding redefinition of node voltage  at  GROUND  changes zero to zero.
/SWITCH
  GROUND                                              MEASURING                1
BLANK card terminates all TACS data
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
C       Total network loss  P-loss  by summing injections =   7.106450000000E-06
C Inject:  GEN       377.             377.             .377E-7   2.6786345332877
C Inject:            0.0              0.0      -2.678634533288       -89.9999992
C           ---- Initial flux of coil  "GROUND"  to  "GEN   "  =  1.00000000E-09
C
C  Step     Time        GEN          GROUND       GROUND       TACS         TACS
C                                     TERRA        GEN          FLUX         GEN
C ***      Phasor I(0) =  0.3770000E-07   Switch  "GROUND"  to  "      "  closed
C   0      0.0          377.      .377E-7          0.0          0.0          0.0
C   1    .5E-4   376.9330268  .5235646391  -.523564601  .0094233257  376.9330268
C   2    .1E-3   376.7321312   1.04694322  -1.04694318  .0282649546  376.7321312
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C     BENCHMARK DC-33
C      2nd  of  9  subcases  is the same as the first except that control system
C      modeling is to be done using  MODELS  rather than the older  TACS.  Jerry
C      Almos of BPA converted the data on 18 Nov 1991 after WSM changed the code
C      to allow such sorting using the same  /TACS  declaration. 
 .000050    .020
       1       1       1       1       1      -1
       5       5      20      20
MODELS              { Note the change;  the 1st subcase had  "TACS HYBRID"  here
/SOURCE
14GEN           377.       60.                                 -1.
/MODELS        { 9 July 1995,  the former  /TACS  was changed to this new class
INPUT GEN {V(GEN)}, GROUND {I(GROUND)}
MODEL DC33
INPUT gen, ground
VAR   flux, curr
HISTORY gen {DFLT: 0}, flux {DFLT: 0}, ground {DFLT: 0} 
INIT
  curr:=0
ENDINIT
EXEC
  COMBINE AS INTEGRATOR
    LAPLACE(flux/gen):=(1.0|S0)/(1.0|S1)
  ENDCOMBINE
  curr:=-ground
ENDEXEC
ENDMODEL
USE DC33 AS DC33
  INPUT GEN:=GEN, GROUND:=GROUND
ENDUSE
RECORD
  DC33.FLUX AS FLUX
  DC33.GEN AS GEN
  DC33.GROUND AS GROUND
  DC33.CURR AS CURR
ENDRECORD
ENDMODELS 
/OUTPUT
  GEN
/PLOT
C Last step:  400      .02   116.4994069   2.765761424  -2.76576141  .9604757618  
C Variable maxima :                 377.   3.499613141  3.499909256   1.00939921         
C Times of maxima :                  0.0        .00415        .0125       .00415          
C Variable minima :          -376.992558   -3.49990926  -3.49961314  -.990568949  
C Times of minima :               .00835         .0125       .00415        .0125       
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
  CALCOMP PLOT
  SCREEN PLOT
 194 2. 0.0 20.         MODELSCURR              ABCDEFGHIJKLMNOP1234567890123456
  X-Y PLOT    Horizontal Axis label123
     10.    -1.1     1.1
      8.    -4.0     4.0
 194 4. 0.0 20.         MODELSFLUX  MODELSCURR  Graph heading---Vertical axis la
  X-Y PLOT    Horizontal Axis label123
   9999.    -1.1     1.1     "9999." returns to conventional (vs. time) plotting
      8.    -4.0     4.0
/BRANCH
96GROUNDGEN                8888.       1.E-9                                   1
             1.0            -0.7
             2.0             0.9
             3.5             1.0
           9999.
/REQUEST
PRINTED NUMBER WIDTH, 13, 2,  { Request maximum precision (for 8 output columns)
/SWITCH
  GROUND                                              MEASURING                1
BLANK card terminates all MODELS data  { 9 July 1995,  MODELS  replaces old
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
C       Total network loss  P-loss  by summing injections =   7.106450000000E-06
C Inject:  GEN       377.             377.             .377E-7   2.6786345332877
C Inject:            0.0              0.0      -2.678634533288       -89.9999992
C           ---- Initial flux of coil  "GROUND"  to  "GEN   "  =  1.00000000E-09
C
C  Step     Time        GEN          GROUND       GROUND       TACS         TACS
C                                     TERRA        GEN          FLUX         GEN
C ***      Phasor I(0) =  0.3770000E-07   Switch  "GROUND"  to  "      "  closed
C   0      0.0          377.      .377E-7  -.25696E-16      .009425         377.      .377E-7     -.377E-7
C   1    .5E-4   376.9330268  -1.42571976  1.425719797  .0282733257  376.9330268  -1.42571976  1.425719759
C   2    .1E-3   376.7321312  -1.41394374  1.413943779  .0471149546  376.7321312  -1.41394374  1.413943741
C   3   .15E-3   376.3973844  -1.40217609   1.40217613  .0659431925  376.3973844  -1.40217609  1.402176093
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C      3rd  of  9  subcases ends the original series.  4-9 added 5 March 2007
C     Illustration of user-supplied residual flux of Type-96 hysteretic
C     inductor was supplied by Prof. Mustafa Kizilcay of FH Osnabrueck in
C     Germany (see E-mail dated 23 March 2001).  This data is added to show
C     that non-zero (but very small) flux of the phasor solution now is
C     tolerated.  Prior to the change several days ago,  it was not.  Prior
C     to the change,  residual flux needed to be exactly zero for the user-
C     supplied flux to be used.  No longer.  Phasor flux < EPSILN  is all that
C     now is required.
POWER FREQUENCY                      50.
PRINTED NUMBER WIDTH, 13, 2,  { dT loop columns have width 13 including 2 blanks
ZERO FLUX TOLERANCE                1.E-8 { Columns 33-40 carry Type-96 tolerance
C     Note about preceding declaration, which was added  29 March 2001.  After
C     reading later discussion about trouble with the  IDEAL TRANSFORMER  data,
C     Prof. Kizilcay agreed that definition of the flux tolerance independent of
C     EPSILN  would be a good idea.  The user can define whatever value he wants
C     although here nothing will change (1.E-8 = value of EPSILN in STARTUP).
 .000200    .200     50.
      40       1       1       1       1
C        Following TACS data is not required.  It is used only to integrate the
C        inductor voltage, thereby indicating the change in inductor flux  ---
C        change with respect to the user-supplied initial value of 0.2 volt-sec.
TACS HYBRID
 1FLUSSA                          +MA                 1.
        1.                   { Numerator   is :    1                                                   
                  1.         { Denominator is :    0 + 1 * s                                                   
90MA          { Type-96 voltage of electric network defines this TACS source                                                                1.
33FLUSSA      { Output of integrator,  the change of flux,  is only TACS output
BLANK card ending TACS data
  QA    P1A                .0172 .1722                                         0
  P3A                      .0172 .1722                                         0
  MA                         1.8                                               0
96MA                        8888       0.20                                    1
         -0.1199         -0.3754  { Preceding 0.2 is residual flux in [volt-sec]
        -0.06264         -0.3643  { 2nd of many current-flux pairs
         -0.0353         -0.3536  { 3rd of many.  These define hysteresis loop.
        -0.02025         -0.3428
        -0.00732         -0.3214
       -0.002869            -0.3
       -0.000538         -0.2786
         0.01345             0.0
         0.02744          0.2786
         0.02977             0.3
         0.03422          0.3214
         0.04715          0.3428
          0.0622          0.3536
         0.08954          0.3643
          0.1468          0.3754
          0.2254          0.3826
          0.5494          0.3968
            9999    { Terminate last point of Type-96 hysteresis element
  P2A                      1.E10
BLANK card ending all branch cards
  P2A   P3A         .020       1.0     { Switch closes at 20 msec }            1
BLANK card ending switches
14QA           8165.       50.      -90.                           -1.
11P1A          1E-10   { 1st half of ungrounded voltage source and ideal Xformer 
18P2A          .0025MA { 2nd half of ungrounded voltage source and ideal Xformer        
C   Note about preceding two source cards, which represent an ideal transformer.  
C   Yes, there is the ungrounded voltage source,  but for engineering purposes,
C   this has value zero (the  1.E-10  in columns 11-20 of the 1st source card).
C   The second source card,  a Type-18 source, defines turns ratio and 4 names.
C   In theory,  these two source cards could be replaced by the following three  
C   branch cards,  which are more readable:
C  IDEAL TRANSFORMER
C 1P1A   P2A                            400.0    { Rated voltage of high side
C 2MA                                     1.0    { Rated voltage of low  side 
C   The resulting phasor solution would be nearly the same, too.  Nearly,  
C   but not quite.  Prof. Kizilcay's Type-11 source value is small, note   
C   ---  small enough to make the phasor flux at time zero < EPSILN.  If
C   instead the  IDEAL TRANSFORMER  is used,  it would appear that the EPSILN
C   threshold is exceeded.  In place of the user-desired residual flux of 0.2
C   volt-seconds,  the .LIS file will show the noise from the phasor solution:
C        ---- Initial flux of coil  "MA    "  to  "      "  =  -1.87122597E-06
C   This then is carried into the dT loop.  Without meaningful residual flux, 
C   there is no significant in-rush current.  Instead of Prof. Kizilcay's 
C   30 amps, the  PRINTER PLOT  shows only about 1.243E-2 amps.  For engineering
C   purposes,  the simulation will be completely different (and wrong).
BLANK card ending sources
C  MA                   .1797146122E-11  .36742498645E-8      -.998414547E-12  .9984146668E-12      .8600609455E-40  .8600609455E-40
C                      -.3674249425E-8      -89.9719755      -.488343774E-15     -179.9719755      .1834212477E-20        0.0000000
C       Total network loss  P-loss  by summing injections =   3.333364056743E-03
C       ---- Initial flux of coil  "MA    "  to  "      "  =  2.00000000E-01
C  The preceding result of phasor solution confirms correct operation.  The
C  voltage across the Type-96 branch is not zero,  but it is small enough
C  compared with floating-point miscellaneous data parameter EPSILN so no
C  special precautions must be taken.  Understand what is happening.  The
C  initial flux is equal to the imaginary part of the phasor voltage divided
C  by the frequency,  and this is small:  .367E-8 / OMEGA < EPSILN 
  QA    MA   { Names of nodes for node-voltage output
C  First  2     output variables are electric-network voltage differences (upper voltage minus lower voltage);
C  Next   2     output variables are branch currents (flowing from the upper node to the lower node);
C  Next   1     output variables belong to TACS (with "TACS" an internally-added upper name of pair).
C  Step    Time        QA           MA           P2A          MA           TACS
C                                                P3A          TERRA        FLUSSA
C     0     0.0   .499946E-12  .179715E-11          0.0  -.99841E-12          0.0
C    40    .008   4799.266585  .2154661E-8          0.0  .278303E-11  .211415E-10
C    80    .016   -7765.37646  -.349521E-8          0.0  .438944E-12  .803576E-11
C   ***        Close switch  "P2A   "  to  "P3A   "  after  2.02000000E-02  sec.
C   120    .024   7765.376456  70.13133651  15591.82659  .0177128714  .0990291974
C   160    .032   -4799.26658  110.5331176  33643.29393  22.70094726  1.167639417
BLANK card ending node voltage outputs
C  1000     0.2   .9430083E-8  -105.223163  -23383.3158   -.97653E-3  .0484350031
C Variable maxima :     8165.   159.787186  40889.62376  30.82522183  1.523703303
C Times of maxima :      .005        .0288        .0294         .054         .054
C Variable minima :    -8165.  -118.702745  -23490.9996   -.16538292  -.198576963
C Times of minima :      .015        .0792        .1996        .1794        .1846
  CALCOMP PLOT
 193.02 0.0 0.2 -10. 30.MA      { Plot Type-96 current over entire time span
  PRINTER PLOT
 194 4. 20. 60.         MA                         { Axis limits: (0.000  3.083)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C      4th  of  9  subcases begins a series using  /UM  on  DC-35 U.M. data
C     1st of 6 subcases that have same solution but slightly different data.
C     Begin with  /SOURCE  (not /UM for U.M.) and full Class-1 U.M. data.
C     Output has been reduced from 2 lines per dT to one which consists of
C     just the U.M. variables.  For the 10 dT that are taken,  signals agree
C     exactly with  DC37.LIS  and are documented on comment cards of  /OUTPUT
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
/OUTPUT
C  The following shows dT-loop output for each of the subcases of the /UM set:
C Step   Time     UM-1        UM-1        UM-1        UM-1        UM-1        UM-1        UM-1        UM-1        UM-1
C                 TQGEN       OMEGM       THETAM      IPA         IPB         IPC         IE1         IE2         IE3
C   0    0.0        0.0         0.0         0.0         0.0         0.0         0.0         0.0         0.0         0.0
C   5   .002    -.6898285         0.0         0.0  -150.98163  16.0278591  134.953771  15.5970274    131.9385  -147.53553
C  10   .004   -11.779008  -.21882419  -.42361E-3  -203.70113  -91.734138   295.43527  -90.126277  287.718333  -197.59206
/SOURCE
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
/BRANCH
 0      BUS-A0            1.0
 0      BUS-B0            1.0
 0      BUS-C0            1.0
 0BUS-A1BUS-A0            0.005 1.0
 0BUS-B1BUS-B0BUS-A1BUS-A0
 0BUS-C1BUS-C0BUS-A1BUS-A0
 0BUSRA1                   1.E-5
 0BUSRB1                   1.E-5
 0BUSRC1                   1.E-5
 0BUS-OM                              6.00E4
 0      BUS-OM            33.333
/REQUEST
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 8 output columns)
/SOURCE
19     { Data for U.M. begins with source type 19 keyed in columns 1-2
 0   { The zero does nothing.  In fact, a blank card would suffice for DC-35
BLANK card ending U.M. header (3rd of 3 non-comment cards of Class-1 U.M. data)
 4 1 1111BUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     BUS-A1            1 0.0           { 1st winding card
0.063         0.0003925     BUS-B1            1 0.0
0.063         0.0003925     BUS-C1            1 0.0
0.083         0.0003925     BUSRB1            1 0.0
0.083         0.0003925     BUSRC1            1 0.0
0.083         0.0003925     BUSRA1            1 0.0 { 6th and final winding card
/PLOT
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C      5th  of  9  subcases continues series using  /UM  on  DC-35 U.M. data
C     2nd of 6 examples that have the same solution but slightly different data
C     Replace  /SOURCE  by  /UM  and omit all 3 non-comment header cards.  But
C     then  U.M. HEADER CARD (UMHC)  must be added to define Class-1 data.
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
/OUTPUT
C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM   { Shorten output by eliminating node voltages
/SOURCE
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
/BRANCH
 0      BUS-A0            1.0
 0      BUS-B0            1.0
 0      BUS-C0            1.0
 0BUS-A1BUS-A0            0.005 1.0
 0BUS-B1BUS-B0BUS-A1BUS-A0
 0BUS-C1BUS-C0BUS-A1BUS-A0
 0BUSRA1                   1.E-5
 0BUSRB1                   1.E-5
 0BUSRC1                   1.E-5
 0BUS-OM                              6.00E4
 0      BUS-OM            33.333
/REQUEST
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 8 output columns)
/UM
C   Note that there is no Class-1 U.M. data;  ATP will supply the Type-19 card.
 4 1 1111BUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     BUS-A1            1 0.0           { 1st winding card
0.063         0.0003925     BUS-B1            1 0.0
0.063         0.0003925     BUS-C1            1 0.0
0.083         0.0003925     BUSRB1            1 0.0
0.083         0.0003925     BUSRC1            1 0.0
0.083         0.0003925     BUSRA1            1 0.0 { 6th and final winding card
C 14BUSRA1  180.0      60.0     0.0
/PLOT
  PRINTER PLOT  { Batch-mode plotting is for the line printer not vector-graphic
/REQUEST
U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1
 0                             { In fact, a blank card would suffice for DC-35
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE        ---  NEWSORT  ---
C      6th  of  9  subcases continues series using  /UM  on  DC-35 U.M. data
C     3rd of 6 examples that have the same solution but slightly different data
C     Return to use of  /SOURCE  and full  Class-1  U.M. data as in 1st example.
C     But add  U.M. HEADER CARD (UMHC)  which overrides the normal Class-1 data.
C       See explanation of the preceding  NEWSORT  within  DCNEW-25.  There 
C       is old (normal) data sorting logic and there is new data sorting 
C       logic as summarized in the January, 1999, Can/Am newsletter.  Putting
C       NEWSORT  on the case-separator card as shown has the effect of placing
C       a minus sign on  SZBED  of  STARTUP,  and this selects the newer sorting
C       alternative.  Either should handle  /UM  sorting.  The  NEWSORT  could
C       be added to each subcase without significant change to the .LIS file.
C       Text of  BLANK  cards will vary a little,  but that should be all.
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
/OUTPUT
C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM   { Shorten output by eliminating node voltages
/SOURCE
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
/BRANCH
 0      BUS-A0            1.0
 0      BUS-B0            1.0
 0      BUS-C0            1.0
 0BUS-A1BUS-A0            0.005 1.0
 0BUS-B1BUS-B0BUS-A1BUS-A0
 0BUS-C1BUS-C0BUS-A1BUS-A0
 0BUSRA1                   1.E-5
 0BUSRB1                   1.E-5
 0BUSRC1                   1.E-5
 0BUS-OM                              6.00E4
 0      BUS-OM            33.333
/REQUEST
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 8 output columns)
U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1
 0                             { In fact, a blank card would suffice for DC-35
/SOURCE
C    Any number of comments should be legal anywhere within Class-1 U.M. data
C    as now begins.  Because of UMHC (U.M. HEADER CARD) use,  only the first of
C    the 3 non-comment lines of Class-1 U.M. data ever will be seen. The crucial
C    2nd will be replaced by the card following the UMHC declaration,  and the
C    blank 3rd card is read and discarded,  so will not been seen in the .LIS
C    file.  Yet, it is required.  Even though not-used,  the 3rd is checked
C    for being blank,  and execution will be halted if it is not so.  Add "XXX"
C    at the start of the 2nd to emphasize that content is arbitrary: 
19     { Data for U.M. begins with source type 19 keyed in columns 1-2
XXX 0  { 2nd of 3 non-comment header lines for U.M. will be replaced by UMHC crd
BLANK card ending U.M. header (3rd of 3 never will been seen in .LIS file)
 4 1 1111BUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     BUS-A1            1 0.0           { 1st winding card
0.063         0.0003925     BUS-B1            1 0.0
0.063         0.0003925     BUS-C1            1 0.0
0.083         0.0003925     BUSRB1            1 0.0
0.083         0.0003925     BUSRC1            1 0.0
0.083         0.0003925     BUSRA1            1 0.0 { 6th and final winding card
/PLOT
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C      7th  of  9  subcases continues series using  /UM  on  DC-35 U.M. data
C     4th of 6 examples that have the same solution but slightly different data
C     Data of the preceding 3rd example has been massaged.  First, "/" cards
C     for sorting have been removed to demonstrate that UMHC use does not
C     require sorting.  The full 3-card Class-1 U.M. data has been replaced
C     by the special Type-19 summary line which will refer to UMHC for the
C     missing 2nd nonblank card of Class-1 data.  Finally,  since no sorting
C     to insert the blank card ending U.M. data,  this must be added manually.
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 8 output columns)
U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1
 0                             { In fact, a blank card would suffice for DC-35
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
 0      BUS-A0            1.0
 0      BUS-B0            1.0
 0      BUS-C0            1.0
 0BUS-A1BUS-A0            0.005 1.0
 0BUS-B1BUS-B0BUS-A1BUS-A0
 0BUS-C1BUS-C0BUS-A1BUS-A0
 0BUSRA1                   1.E-5
 0BUSRB1                   1.E-5
 0BUSRC1                   1.E-5
 0BUS-OM                              6.00E4
 0      BUS-OM            33.333
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
19 { Class-1 U.M. data has been omitted WSM
 4 1 1111BUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     BUS-A1            1 0.0           { 1st winding card
0.063         0.0003925     BUS-B1            1 0.0
0.063         0.0003925     BUS-C1            1 0.0
0.083         0.0003925     BUSRB1            1 0.0
0.083         0.0003925     BUSRC1            1 0.0
0.083         0.0003925     BUSRA1            1 0.0 { 6th and final winding card
BLANK card ends U.M. data
BLANK terminates the last SOURCE card
C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM   { Shorten output by eliminating node voltages
BLANK card ends OUTPUT variable requests
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C      8th  of  9  subcases continues series using  /UM  on  DC-35 U.M. data
C     5th of 6 subcases that have the same solution but slightly different data
C     Data is like the 2nd example except that a second and third disconnected
C     machine has been added to illustrate multi-machine use for 3 machines.
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
/OUTPUT
C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM   { Shorten output by eliminating node voltages
/SOURCE
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
/BRANCH
 0      BUS-A0            1.0
 0      BUS-B0            1.0
 0      BUS-C0            1.0
 0BUS-A1BUS-A0            0.005 1.0
 0BUS-B1BUS-B0BUS-A1BUS-A0
 0BUS-C1BUS-C0BUS-A1BUS-A0
 0BUSRA1                   1.E-5
 0BUSRB1                   1.E-5
 0BUSRC1                   1.E-5
 0BUS-OM                              6.00E4
 0      BUS-OM            33.333
/REQUEST
U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1
 0                             { In fact, a blank card would suffice for DC-35
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 8 output columns)
/SOURCE
19 { Class-1 U.M. data has been omitted (this card is needed because of /SOURCE)
 4 1 1111BUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     BUS-A1            1 0.0           { 1st winding card
0.063         0.0003925     BUS-B1            1 0.0
0.063         0.0003925     BUS-C1            1 0.0
0.083         0.0003925     BUSRB1            1 0.0
0.083         0.0003925     BUSRC1            1 0.0
0.083         0.0003925     BUSRA1            1 0.0 { 6th and final winding card
C 14BUSRA1  180.0      60.0     0.0
/PLOT
  PRINTER PLOT  { Batch-mode plotting is for the line printer not vector-graphic
C  ----------------  Done with data of 1st (real) U.M.  Begin data of 2nd:
C  2nd, unused U.M. has the same data as the 1st,  but it is disconnected,
C  so should not affect the solution.  Outputs have been suppressed.  For
C  variety,  this U.M. still uses  /UM  whereas the first has been converted
C  to  /SOURCE  followed by the special Type-19 header card.  Note that the
C  Type-19 card is mandatory because any source might follow  /SOURCE  whereas
C  only U.M. data can follow  /UM  (i.e., the "19" is needed for recognition
C  as U.M. data).  For the 2nd, disconnected machine,  note that the original
C  A3  text  BUS  has become  GUS  everywhere.  To restore the 9 outputs of
C  the 2nd machine,  add  "111"  to cols. 7-9 of the card immediately following
C  the  /UM  declaaration,  and change  "0"  to  "1"  in column 47 of each of
C  the 6 winding cards that end machine data.  If this is done,  output consists
C  of 2 lines each time step,  and the two will be identical.
/UM
 4 1 1   GUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     GUS-A1            0 0.0           { 1st winding card
0.063         0.0003925     GUS-B1            0 0.0
0.063         0.0003925     GUS-C1            0 0.0
0.083         0.0003925     GUSRB1            0 0.0
0.083         0.0003925     GUSRC1            0 0.0
0.083         0.0003925     GUSRA1            0 0.0 { 6th and final winding card
/BRANCH
 0      GUS-A0            1.0
 0      GUS-B0            1.0
 0      GUS-C0            1.0
 0GUS-A1GUS-A0            0.005 1.0
 0GUS-B1GUS-B0GUS-A1GUS-A0
 0GUS-C1GUS-C0GUS-A1GUS-A0
 0GUSRA1                   1.E-5
 0GUSRB1                   1.E-5
 0GUSRC1                   1.E-5
 0GUS-OM                              6.00E4
 0      GUS-OM            33.333
/SOURCE
14GUS-A0  180.0      60.0     0.0
14GUS-B0  180.0      60.0     -120.0
14GUS-C0  180.0      60.0     +120.0
11GUS-OM-1 -6.00
C  ----------------  Done with data of 2nd U.M.  Begin data of 3rd:
C  3rd, unused U.M. has the same data as the 1st,  but it is disconnected,
C  so should not affect the solution.  Outputs have been suppressed.  Etc.
C  (like the 2nd except that the original  BUS  --->  MUS  here).
/UM
 4 1 1   MUS-OM       2 0.06                       0.18850   { 1st machine table
0.0           0.020291      0
0.0           0.020291      0
0.063         0.0001911     MUS-A1            0 0.0           { 1st winding card
0.063         0.0003925     MUS-B1            0 0.0
0.063         0.0003925     MUS-C1            0 0.0
0.083         0.0003925     MUSRB1            0 0.0
0.083         0.0003925     MUSRC1            0 0.0
0.083         0.0003925     MUSRA1            0 0.0 { 6th and final winding card
/BRANCH
 0      MUS-A0            1.0
 0      MUS-B0            1.0
 0      MUS-C0            1.0
 0MUS-A1MUS-A0            0.005 1.0
 0MUS-B1MUS-B0MUS-A1MUS-A0
 0MUS-C1MUS-C0MUS-A1MUS-A0
 0MUSRA1                   1.E-5
 0MUSRB1                   1.E-5
 0MUSRC1                   1.E-5
 0MUS-OM                              6.00E4
 0      MUS-OM            33.333
/SOURCE
14MUS-A0  180.0      60.0     0.0
14MUS-B0  180.0      60.0     -120.0
14MUS-C0  180.0      60.0     +120.0
11MUS-OM-1 -6.00
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C      9th  of  9  subcases ends the series using  /UM  on  DC-35 U.M. data
C     6th of 6 examples that have the same solution but slightly different data
C     Data is the same as the 2nd example except that  /UM  data has been moved
C     to a  $INCLUDE  file.  Use or non-use of  $INCLUDE  should have no effect,
C     and this demonstrates none for this simplest, first level of  $INCLUDE.
C       Enhancement by Gabor Furst.  WSM's single  $INCLUDE  was replaced by
C       two for more generality.  A second induction motor was added to show
C       that this is possible.  Mr. Furst observed the following on  9 April
C       2007:  "I thought we could have this sub case with two rather than
C       one  $INCLUDE file.  I attach my version of this, which demonstrates
C       how adding another machine is simply adding another $INCLUDE file."
C       Mr. Furst simulated to  Tmax = 1.044 seconds  but program developers
C       truncate this drastically (to 10 steps) to speed execution.  This is
C       just enough to start the rotors and show approximate balance.  To see
C       the entire settling into the steady state,  restore the longer T-max,
C       remove the minus sign from  IPLOT  immediately below T-max,  and
C       activate the batch-mode plot card (now omitted by $DISABLE).  Finally,
C       control printout (every 5th step for IOUT is a lot, for realistic use).
$PREFIX, []    { $INCLUDE files are located in same place as this main data file
 .000400    .004   { Take only 10 time steps -- enough to begin motor rotation
       5      -1   { No plotting (which has nothing special to do with /UM use)
/OUTPUT
C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM   { Shorten output by eliminating node voltages
/SOURCE
14BUS-A0  180.0      60.0     0.0
14BUS-B0  180.0      60.0     -120.0
14BUS-C0  180.0      60.0     +120.0
11BUS-OM-1 -6.00
14BUS-A3  180.0      60.0     0.0
14BUS-B3  180.0      60.0     -120.0
14BUS-C3  180.0      60.0     +120.0
11BUS-3M-1 -6.00
/BRANCH
C the following branch records are transferred to $INCLUDE  dc33inc1.dat
C  0      BUS-A0            1.0    
C  0      BUS-B0            1.0
C 0      BUS-C0            1.0
C 0BUS-A1BUS-A0            0.005 1.0C
C 0BUS-B1BUS-B0BUS-A1BUS-A0
C 0BUS-C1BUS-C0BUS-A1BUS-A0
C 0BUSRA1                   1.E-5
C 0BUSRB1                   1.E-5
C 0BUSRC1                   1.E-5
C 0BUS-OM                              6.00E4
C 0      BUS-OM            33.33
C
 0      BUS-A3            1.0
 0      BUS-B3            1.0
 0      BUS-C3            1.0
 0BUS-A4BUS-A3            0.005 1.0
 0BUS-B4BUS-B3BUS-A4BUS-A3
 0BUS-C4BUS-C0BUS-A4BUS-A3
 0BUSRA4                   1.E-5
 0BUSRB4                   1.E-5
 0BUSRC4                   1.E-5
 0BUS-3M                              6.00E4
 0      BUS-3M            33.333
$INCLUDE,  dc33inc1.dat, { Separate file for /UM followed by all data for U.M. 1
$INCLUDE,  dc33inc2.dat, { Separate file for /UM followed by all data for U.M. 2
C   About the preceding data modularization,  there is lack of balance.  This
C   is intentional.  Whereas  DC33INC1  has  /BRANCH  data,  DC33INC2  has none.
C   Each also involves  /REQUEST  which precedes a  UMHC  request.  Since only
C   one has meaning, expect rejection accompanied by a warning message for the
C   2nd of the two.
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
  CALCOMP PLOT
$DISABLE    { The following plot card assumes use of Gabor Furst's longer T-max
 193 .1 0.0 1.0         UM-1  OMEGM UM-2  OMEGM Induction motor w in [rad/sec]
$ENABLE
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
BLANK