summaryrefslogtreecommitdiffstats
path: root/benchmarks/dcn23.dat
diff options
context:
space:
mode:
authorAngelo Rossi <angelo.rossi.homelab@gmail.com>2023-06-21 12:04:16 +0000
committerAngelo Rossi <angelo.rossi.homelab@gmail.com>2023-06-21 12:04:16 +0000
commitb18347ffc9db9641e215995edea1c04c363b2bdf (patch)
treef3908dc911399f1a21e17d950355ee56dc0919ee /benchmarks/dcn23.dat
Initial commit.
Diffstat (limited to 'benchmarks/dcn23.dat')
-rw-r--r--benchmarks/dcn23.dat484
1 files changed, 484 insertions, 0 deletions
diff --git a/benchmarks/dcn23.dat b/benchmarks/dcn23.dat
new file mode 100644
index 0000000..756fd44
--- /dev/null
+++ b/benchmarks/dcn23.dat
@@ -0,0 +1,484 @@
+BEGIN NEW DATA CASE
+C BENCHMARK DCNEW-23
+C Illustrate modeling of Type-51,52 elements to model unsymmetric series
+C connection of [R] and [L] as described in the October, 1997, newsletter
+C Begin with uncoupled branches of 1 ohm resistance and just 2 phases (no
+C Type-53). The answer is obvious by inspection: the sinusoidal source
+C voltage is exactly split in half. Very low frequency and small time
+C step means that this varies very slowly --- in the 7th decimal place.
+ .00001 .00005 { Note just 5 very small time steps T-max = 5 * DELTAT
+ 1 -1 1 0
+ SENDA RECA 1.0 1
+ SENDB RECB SENDA RECA 1
+51RECA MODEL [R][L] { Note special request in cols. 15-26
+52RECB { Cols. 27-44 of preceding card is blank
+ 1.0 0.0
+ 0.0 0.0
+ 0.0 1.0
+ 0.0 0.0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 2 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA SENDB RECB
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 2nd of 15. Identical to preceding except that more compact alternative
+C for 2-phase (Type-51,52 with no 53) is illustrated. Nonblank col 27-44
+C distinguishes this case, which has no following, separate cards for the
+C [R] and [L] matrices.
+ .00001 .00005
+ 1 -1 1 0
+ SENDA RECA 1.0 1
+ SENDB RECB SENDA RECA 1
+C 789012345678901234 890123456789012
+51RECA MODEL [R][L] 1.0 0.0 0.0 0.0
+52RECB 0.0 0.0 1.0 0.0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA SENDB RECB
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 3rd of 15. Identical to 1st except that here the Type-51,52 branch is
+C inputted before rather than after the two series R-L-C branches. Also,
+C current output using column-80 punches is illustrated. This replaces
+C the current output of the Type-0 branches in series
+ .00001 .00005
+ 1 -1 1 0
+51RECA MODEL [R][L] 1
+52RECB 1
+ 1.0 0.0
+ 0.0 0.0
+ 0.0 1.0
+ 0.0 0.0
+ SENDA RECA 1.0 0
+ SENDB RECB SENDA RECA 0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 2 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA SENDB RECB
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 4th of 15 subcases progresses from preceding 2-phase to 3-phase. Begin
+C with special Z0Z1Z2 model that assumes decoupled symmetrical component
+C impedances Zo for the zero sequence, Z1 for the positive sequence,
+C and finally Z2 for the negative sequence. Otherwise, the problem is
+C the same: 5 very small steps of uncoupled, resistive voltage dividers.
+C About current outputs, this subcases introduces selective use of col.
+C 80 punches. Note phase C output of the series R-L-C is missing, and
+C that this is supplied by the punch on the Type-53 card. Switching the
+C zero and the 1 of these two phase-C cards does not change the numbers
+C (it only changes the heading names). That was for the 3rd of 3. One
+C can change the 2nd (phase B) in comparable fashion without changing
+C numbers of the time-step loop output. This illustrates selective
+C branch current output.
+C diagnostic 0 0 9 9 9 9
+ .00001 .00005
+ 1 -1 1 0
+ SENDA ENDA 1.0 1
+ SENDB ENDB SENDA ENDA 1
+ SENDC ENDC SENDA ENDA 0
+51ENDA MODEL Z0Z1Z2 1.0 { Ro, Lo right-adjusted in cols. 32 & 44
+52ENDB 1.0 { R1 and L1 .... }
+53ENDC 1.0 { R2 and L2 .... } 1
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 0.001 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA ENDA SENDB ENDB SENDC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 5th of 15. Identical to 1st except that here the Type-51,52,53 branch
+C is inputted before rather than after the three series R-L-C branches.
+ .00001 .00005
+ 1 -1 1 0
+51ENDA MODEL Z0Z1Z2 1.0
+52ENDB 1.0
+53ENDC 1.0
+ SENDA ENDA 1.0 0
+ SENDB ENDB SENDA ENDA 0
+ SENDC ENDC SENDA ENDA 0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 0.001 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 6th of 15. Identical to case before preceding subcase except that here
+C MODEL [R][L] is used in place of MODEL Z0Z1Z2 for the Type-51,52,53.
+C Also, column-80 requests for branch current output have been added.
+ .00001 .00005
+ 1 -1 1 0
+ SENDA RECA 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+51RECA MODEL [R][L] { Cols. 27-44 of this card is blank } 1
+52RECB { This is signal for separate, later } 1
+53RECC { matrices [R] and [L]. } 1
+ 1.0 0.0 0.0 { R(1,1), R(1,2), R(1,3)
+ 0.0 0.0 0.0 { L(1,1), L(1,2), L(1,3)
+ 0.0 1.0 0.0 { Etc. row 2
+ 0.0 0.0 0.0
+ 0.0 0.0 1.0 { Etc. row 3
+ 0.0 0.0 0.0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 0.001 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 7th of 15. Identical to 1st except that here the Type-51,52,53 branch
+C is inputted before rather than after the three series R-L-C branches.
+ .00001 .00005
+ 1 -1 1 0
+51RECA MODEL [R][L]
+52RECB
+53RECC
+ 1.0 0.0 0.0
+ 0.0 0.0 0.0
+ 0.0 1.0 0.0
+ 0.0 0.0 0.0
+ 0.0 0.0 1.0
+ 0.0 0.0 0.0
+ SENDA RECA 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 0.001 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 8th of 15. Combination of preceding subcases involves Type-51,52,53
+C branch specified with both MODEL [R][L] and MODEL Z0Z1Z2. This is
+C a double-size problem; each half has the same solution as previously.
+ .00001 .00005
+ 1 -1 1 0
+ SENDA ENDA 1.0 0
+ SENDB ENDB SENDA ENDA 0
+ SENDC ENDC SENDA ENDA 0
+51ENDA MODEL Z0Z1Z2 1.0
+52ENDB 1.0
+53ENDC 1.0
+ SENDA RECA 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+51RECA MODEL [R][L]
+52RECB
+53RECC
+ 1.0 0.0 0.0
+ 0.0 0.0 0.0
+ 0.0 1.0 0.0
+ 0.0 0.0 0.0
+ 0.0 0.0 1.0
+ 0.0 0.0 0.0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.001 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 0.001 -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 0.001 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 9th of 15. Mix new Type-51,52,53 with the old. With diagonal matrices,
+C the data can be represented either way, of course (it is symmetric).
+C Z1 = Zo (actually, Z1 = Z0) so the answer can be shown to be
+C identical to that using Type-51,52,53 modeling with sequence impedances.
+C In fact, there are 2 identical, uncoupled networks driven from the same
+C balanced 3-phase sources at SENDA, SENDB, and SENDC. Another first:
+C switch from resistance to inductance.
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+51ENDA MODEL Z0Z1Z2 0.0 1.0 { Sequence Ro, Lo in [ohms, mHenry]
+52ENDB 0.0 1.0 { Sequence R1, L1 in [ohms, mHenry]
+53ENDC 0.0 1.0 { Note Z2 = Z1 so [Z] is symmetric
+ SENDA ENDA 0.0 1.0 0
+ SENDB ENDB SENDA ENDA 0
+ SENDC ENDC SENDA ENDA 0
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+51RECA 0.0 1.0 { Ro, Lo in [ohms, mHenry]
+52RECB 0.0 1.0 { Ro, Lo in [ohms, mHenry]
+53RECC { Blank 27-44 on this 3rd card ==> sequence data
+ SENDA RECA 0.0 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 50. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 50. -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 50. 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C diagnostic 0 0 9 9 9 9 0 9 0 0 0 9 9 9 9 9 9 9 9
+C 10th of 15 subcases. Generalize preceding by changing from R only to
+C both R and L. Also, the data order has been shuffled, to illustrate
+C that the answer is unaffected.
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+C 22 Feb 2000, switch the following 3 data cards from narrow to wide format to
+C illustrate that this works following correction a few days earlier:
+C 51RECA 0.3 1.0 { Ro, Lo in [ohms, mHenry]
+C 52RECB 0.3 1.0 { R1, L1 in [ohms, mHenry]
+C 53RECC
+$VINTAGE, 1,
+C 34567890123456789012345678901234567890123456789012345678
+51RECA 0.3 1.0 { Ro, Lo in [ohms, mHenry]
+52RECB 0.3 1.0 { R1, L1 in [ohms, mHenry]
+ 0.0 0.0 { Unused cell 2,2
+53RECC 0.0 0.0 { Unused cell 3,1
+ 0.0 0.0 { Unused cell 3,2
+ 0.0 0.0 { Unused cell 3,3
+$VINTAGE, 0, { Done with wide format for 51, 52, 53, so return to narrow
+51ENDA MODEL Z0Z1Z2 0.3 1.0 { Sequence Ro, Lo in [ohms, mHenry]
+52ENDB 0.3 1.0 { Sequence R1, L1 in [ohms, mHenry]
+53ENDC 0.3 1.0 { Note Z2 = Z1 so [Z] is symmetric
+ SENDA ENDA 0.3 1.0 0
+ SENDB ENDB SENDA ENDA 0
+ SENDC ENDC SENDA ENDA 0
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+ SENDA RECA 0.3 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 50. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 50. -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 50. 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 11th of 15 subcases. Generalize preceding by adding inter-phase coupling
+C for the first time. The negative-sequence impedance still is equal to
+C the positive, so both old and new Type-51,52,53 branches can be used.
+C But Zo not equal to Z1 means that phase-domain matrices no longer
+C are diagonal.
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+ SENDA ENDA 0.3 1.0 0
+ SENDB ENDB SENDA ENDA 0
+ SENDC ENDC SENDA ENDA 0
+51ENDA MODEL Z0Z1Z2 0.3 1.0 { Sequence Ro, Lo in [ohms, mHenry]
+52ENDB 0.1 0.5 { Sequence R1, L1 in [ohms, mHenry]
+53ENDC 0.1 0.5 { Note Z2 = Z1 so [Z] is symmetric
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+ SENDA RECA 0.3 1.0 0
+ SENDB RECB SENDA RECA 0
+ SENDC RECC SENDA RECA 0
+51RECA 0.3 1.0 { Ro, Lo in [ohms, mHenry]
+52RECB 0.1 0.5 { Ro, Lo in [ohms, mHenry]
+53RECC
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 50. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 50. -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 50. 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 12th of 15 subcases. Generalize preceding by making the negative-sequenc
+C impedance different from the positive, so the old Type-51,52,53 branch
+C no longer can be used. Instead, both MODEL [R][L] and MODEL Z0Z1Z2
+C are illustrated. There are two identical halves with two identical
+C solutions. As for the phase-domain matrix data, it was copied from
+C near the end of DCN22.DBG with enough precision to make the effect of
+C truncation unobservable in the printed output.
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+ SENDA ENDA 0.3 1.0 1
+ SENDB ENDB SENDA ENDA 1
+ SENDC ENDC SENDA ENDA 1
+ ENDA ENDB 1.E7
+ ENDB ENDC ENDA ENDB
+ ENDC ENDA ENDA ENDB
+51ENDA MODEL Z0Z1Z2 0.3 1.0 { Sequence Ro, Lo in [ohms, mHenry]
+52ENDB 0.1 0.5 { Sequence R1, L1 in [ohms, mHenry]
+53ENDC .101 0.8
+ SENDA RECA 0.3 1.0 1
+ SENDB RECB SENDA RECA 1
+ SENDC RECC SENDA RECA 1
+ RECA RECB 1.E7
+ RECB RECC RECA RECB
+ RECC RECA RECA RECB
+51RECA MODEL [R][L]
+52RECB MODEL [R][L]
+53RECC MODEL [R][L]
+$DISABLE
+ .1670000000 .0665866025 .0664133975 { R(1,1), R(1,2), R(1,3)
+ .7666666667-.1720084679 .4053418013 { L(1,1), L(1,2), L(1,3)
+ .0664133975 .1670000000 .0665866025 { Row 2 of [R]
+ .4053418013 .7666666667-.1720084679 { Row 2 of [L]
+ .0665866025 .0664133975 .1670000000 { Row 3 of [R]
+-.1720084679 .4053418013 .7666666667 { Row 3 of [L]
+$ENABLE
+C Preceding disabled data was from years past. Correct this 2 April 2002.
+C The following is copied from DCN22c (the compensation alternative):
+$UNITS, 60.0, 60.0,
+C USERNL begins with Lo, L1, L2 [H] = 1.00000E-03 5.00000E-04 8.00000E-04
+C Converted to Xo, X1, X2 = 3.76991E-01 1.88496E-01 3.01593E-01
+C 3x3 phase-domain impedance matrix in ohms follow. For each row I, X(I,J) is below R(I,J). w = 3.769911E+02 rad/sec.
+ .1670000000 .0991483886 .0338516114 { R(1,1), R(1,2), R(1,3) in ohms
+ .2890265241 .0436936220 .0442709723 { X(1,1), X(1,2), X(1,3) in ohms
+ .0338516114 .1670000000 .0991483886 { Row 2 of [R]
+ .0442709723 .2890265241 .0436936220 { Row 2 of [X]
+ .0991483886 .0338516114 .1670000000 { Row 3 of [R]
+ .0436936220 .0442709723 .2890265241 { Row 3 of [X]
+$UNITS, -1.0, -1.0, { Done with ohms and micromhos at 60 Hz, so restore original
+C Whereas the preceding could not be used in DCN22c, it can be used here
+C because the conversion is not delayed until the dT loop. For Type-51,
+C 52, 53 branches, the conversion is perfomed at data-input time.
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 50. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 2.0 50. -120. { three-phase, sinusoidal excitation
+14SENDC 2.0 50. 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA RECA ENDA { SENDB RECB ENDB SENDC RECC ENDC
+BLANK card ending node voltage outputs
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C diagnostic 0 0 9 9 9 9 0 9 9 9
+C 13th of 15. 1st use of USE AR for [A] and [R] input rather than
+C the original [R] and [L]. The 2 phases in fact are uncoupled for
+C simplicity. Series L of each phase divide the voltage. Series R-L-C
+C is 1 mHenry, and A(1,1) = A(2,2) correspond to the inverse of this.
+C So, L(1,1) = L(2,2) = 1.E-3 ===> A(1,1) = A(2,2) = 1000. This is
+C entered instead of resistance. Instead of inductance, enter R = 0
+C (the [R] matrix is zero; no resistance in this problem)
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+ USE AR
+51RECA MODEL [R][L]
+52RECB
+ 1000. 0.0
+ 0.0 0.0
+ 0.0 1000.
+ 0.0 0.0
+ USE RL
+ SENDA RECA 0.0 1.0 0
+ SENDB RECB SENDA RECA 1
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+ USE AR
+51ENDA 1000. 0.0
+52ENDB 0.0 0.0 1000. 0.0
+ USE RL
+ SENDA ENDA 0.0 1.0 0
+ SENDB ENDB SENDA ENDA 1
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.1 0.0 { 1st of 2 sources. Note balanced,
+14SENDB 2.0 0.1 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 14th of 15. Like preceding, except resistance of 1 ohm is added to
+C each series, uncoupled inductance. There still is a voltage divider.
+C Another change is use of col-80 output for 2nd phase (Type-52).
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+ USE AR
+51RECA MODEL [R][L]
+52RECB 1
+ 1000. 0.0
+ 1.0 0.0
+ 0.0 1000.
+ 0.0 1.0
+ USE RL
+ SENDA RECA 1.0 1.0 0
+ SENDB RECB SENDA RECA 1
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+ USE AR
+51ENDA 1000. 1.0
+52ENDB 0.0 0.0 1000. 1.0
+ USE RL
+ SENDA ENDA 1.0 1.0 0
+ SENDB ENDB SENDA ENDA 1
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.1 0.0 { 1st of 2 sources. Note balanced,
+14SENDB 2.0 0.1 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 15th of 15. Like preceding, except coupling is added to [R] and [L]
+C matrices. Both halves of problem still should have the same answer
+C (note matrices remain symmetric), although voltage division now
+C longer will be exactly half and half.
+ .0001 .0005
+ 1 1 1 0 1 -1
+ 5 5 20 20 100 100
+ USE AR
+51RECA MODEL [R][L]
+52RECB
+ 1000. 200.
+ 1.0 0.2
+ 200. 1000.
+ 0.2 1.0
+ USE RL
+ SENDA RECA 1.0 1.0 0
+ SENDB RECB SENDA RECA 1
+C Next, build a copy of this, but using the old (Type-51,52,53) modeling:
+ USE AR
+51ENDA 1000. 1.0
+52ENDB 200. 0.2 1000. 1.0
+ USE RL
+ SENDA ENDA 1.0 1.0 0
+ SENDB ENDB SENDA ENDA 1
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 2.0 0.1 0.0 { 1st of 2 sources. Note balanced,
+14SENDB 2.0 0.1 -120. { three-phase, sinusoidal excitation
+BLANK card follows the last source card
+ SENDA RECA ENDA SENDB RECB ENDB
+BLANK card ending node voltage outputs
+ CALCOMP PLOT
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+BLANK