diff options
| author | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
|---|---|---|
| committer | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
| commit | b18347ffc9db9641e215995edea1c04c363b2bdf (patch) | |
| tree | f3908dc911399f1a21e17d950355ee56dc0919ee /benchmarks/dc16.dat | |
Initial commit.
Diffstat (limited to 'benchmarks/dc16.dat')
| -rw-r--r-- | benchmarks/dc16.dat | 230 |
1 files changed, 230 insertions, 0 deletions
diff --git a/benchmarks/dc16.dat b/benchmarks/dc16.dat new file mode 100644 index 0000000..3650a6a --- /dev/null +++ b/benchmarks/dc16.dat @@ -0,0 +1,230 @@ +BEGIN NEW DATA CASE
+C BENCHMARK DC-16
+C Test of "SYSTEMATIC" option which allows the user to have switch closing
+C times varied regularly. There are 3 independent, systematically-varied
+C switches taking 2, 4, and 3 steps, respectively. In addition to the base
+C case, then, there are to be 2 * 4 * 3 = -NENERG = 24 energizations.
+PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
+99.03E-6 .025 60.
+ 1 1 1 2 1 -1 -24
+C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
+ 1 1 { ISW=1 ==> Printed switching times; ITEST=1 ==> no bias
+ 5 5 20 20 { Printout frequency changes for base case
+ 0GENA A1 7. 4
+ 0GENB B1 7. 3
+ 0GENC C1 7.
+ 0ENDA A10 7.
+ 0ENDB B10 7.
+ 0ENDC C10 7.
+-1ASW1 A5 .3 2.1146 0.645 50. 0
+-2BSW1 B5 .0268 .5397 0.021 50. 0
+-3CSW1 C5
+ 0A5 A5F 1.
+ 0B5 B5F 1.
+ 0C5 C5F 1.
+-1A5F ASW10 ASW1 A5
+-2B5F BSW10
+-3C5F CSW10
+BLANK card ending branch cards
+ B10 BSW10 10.E-3 1.0
+ B1 BSW1 4.E-3 .2E-3 2 SYSTEMATIC
+ C1 CSW1 6.E-3 .3E-3 4 SYSTEMATICTARGET
+ A10 ASW10 2.E-3 .4E-3 3 SYSTEMATIC
+ A1 ASW1 1.E-3 SYSTEMATICB1 BSW1
+ C10 CSW10 12.E-3 1.0
+BLANK card ending switch cards
+14GENA 303. 60. 0.0 -1.
+14GENB 303. 60. -120.0 -1.
+14GENC 303. 60. 120.0 -1.
+14ENDA 303. 60. - 10.0 -1.
+14ENDB 303. 60. -130.0 -1.
+14ENDC 303. 60. 110.0 -1.
+BLANK card ending source cards
+ ASW1 BSW1
+BLANK card ending output variable requests (node voltages only, here)
+C *** Close switch "B1 " to "BSW1 " after 6.04083000E-03 sec.
+C *** Close switch "C1 " to "CSW1 " after 6.04083000E-03 sec.
+C *** Close switch "A10 " to "ASW10 " after 6.04083000E-03 sec.
+C *** Close switch "A1 " to "ASW1 " after 7.03113000E-03 sec.
+C 80 .0079224 3.639449214 12.5053685 -283.627168 177.6693051 .5124656E-3
+C 100 .009903 1.299612647 -58.9506992 -271.930029 38.31674092 .6510477E-4
+C *** Close switch "B10 " to "BSW10 " after 1.00020300E-02 sec.
+C 120 .0118836 -13.4793168 7.386747874 -54.0094052 -227.85024 .0068001249
+C *** Close switch "C10 " to "CSW10 " after 1.20816600E-02 sec.
+C
+C Last step: 252 .02495556 -2.83666996 -5.98940516 -327.625802 161.8641485
+C Last step continued ....... .1629402E-3 .3686585491
+ PRINTER PLOT
+ 144 4. 0.0 50. ASW1 { Axis limits: (-3.276, 3.176)
+C GENA GENB ASW1 BSW1 GENA
+C A1 B1 A1
+C First shot: Random switching times for energization number 1 :
+C 2 5.9000000E-03 3 6.0000000E-03 4 5.6000000E-03 5 6.9000000E-03
+C 58.93128751 192.0984488 -326.075521 336.1486991 .0315970504
+C Times of maxima : .0138642 .00604083 .02495556 .00722919 .01554771
+C
+C 2nd shot: Random switching times for energization number 2 :
+C 2 6.1000000E-03 3 6.0000000E-03 4 5.6000000E-03 5 7.1000000E-03
+C 56.04192203 165.7654174 -322.282369 -326.291206 .0329276386
+C Times of maxima : .01267584 .00623889 .02495556 .01346808 .01574577
+BLANK card ending plot cards (for base case only)
+C 23rd shot: Random switching times for energization number 23 :
+C 23rd shot: 2 6.8000000E-3 3 6.9000000E-3 4 7.3000000E-3 5 7.8000000E-3
+C 23rd shot: -39.0656922 169.6381869 -379.444302 -311.343589 .0253668269
+C Times of maxima: .01742928 .0069321 .00812046 .01346808 .01673607
+C
+C Last shot: Random switching times for energization number 24 :
+C Last shot: 2 7.0000000E-3 3 6.9000000E-3 4 7.3000000E-3 5 8.0000000E-3
+C Last shot: 49.82358949 165.4706377 -382.803147 310.7780459 .0299093498
+C Last shot: .01356711 .00713016 .00812046 .02188563 .01673607
+C
+C 1 ) --------------------------------------------------------------------------
+C Statistical distribution of peak voltage for branch "GENB " to "B1 ".
+C Interval voltage voltage in Frequency Cumulative
+C number in per unit physical units (density) frequency
+C 57 2.8500000 0.14250000E+03 0 0
+C 58 2.9000000 0.14500000E+03 1 1
+C < < Etc. (rows 59 through 75 omitted to save space) > >
+C 76 3.8000000 0.19000000E+03 0 21
+C 77 3.8500000 0.19250000E+03 3 24
+C Summary for following request: Mean = 3.48750000E+00 3.49384237E+00
+C Variance = 5.61413043E-02 5.91771084E-02
+C Standard deviation = 2.36941563E-01 2.43263455E-01
+-1 50. GENB B1
+C 2 ) --------------------------------------------------------------------------
+C Statistical distribution of peak voltage at node "ASW1 ". The base voltage
+C Interval voltage voltage in Frequency Cumulative
+C number in per unit physical units (density) frequency
+C 21 1.0500000 0.31815000E+03 0 0
+C 22 1.1000000 0.33330000E+03 5 5
+C 23 1.1500000 0.34845000E+03 5 10
+C 24 1.2000000 0.36360000E+03 3 13
+C 25 1.2500000 0.37875000E+03 5 18
+C 26 1.3000000 0.39390000E+03 3 21
+C 27 1.3500000 0.40905000E+03 3 24
+C Distribution parameters for .. data follow. Grouped data Ungrouped data
+C Mean = 1.18541667E+00 1.18028303E+00
+C Variance = 7.38677536E-03 6.74976204E-03
+C Standard deviation = 8.59463517E-02 8.21569354E-02
+ 0 ASW1 BSW1 { See preceding comments for 1st of 3 tables produced
+C MODTAB AINCR XMAXMX
+STATISTICS DATA 0 -15. 0.0 ---- Fixed 15 boxes used now
+C 5 ) --------------------------------------------------------------------------
+C Statistical distribution of peak current for branch "GENB " to "B1 ".
+C Interval current current in Frequency Cumulative
+C number in per unit physical units (density) frequency
+C 19 0.9976152 0.95970790E+00 0 0
+C 20 1.0501213 0.10102188E+01 2 2
+C 21 1.1026274 0.10607298E+01 0 2
+C < < Etc. (rows 22 through 32 omitted to save space) > >
+C 33 1.7327002 0.16668611E+01 0 22
+C 34 1.7852062 0.17173720E+01 2 24
+C Distribution parameters for .. data follow. Grouped data Ungrouped data
+C Mean = 1.40672501E+00 1.39937500E+00
+C Variance = 3.96702046E-02 3.93384525E-02
+C Standard deviation = 1.99173805E-01 1.98339236E-01
+-2 GENB B1 { See preceding comment cards for resulting single table
+-3 GENA A1 { Ungrouped: Mean= 1.63806617E0, Variance = 1.01307955E-1
+-4 GENA A1 { Ungrouped: Mean= 1.30997561E0, Variance = 1.73383337E-2
+C MODTAB AINCR XMAXMX 15 Aug 03, add Sturges use:
+STATISTICS DATA 0 Sturges 0.0 ---- Sturges gives # of boxes
+-4 GENA A1 { Ungrouped:
+C 8 ) -------------------------------------------------------------------------------------------------------------------------------
+C Statistical distribution of peak energy for branch "GENA " to "A1 ". Base energy for per unit output = 2.33417187E-02
+C Interval energy energy in Frequency Cumulative Per cent
+C number in per unit physical units (density) frequency .GE. current value
+C 9 .9050250 2.11248398E-02 0 0 100.000000
+C 10 1.0055834 2.34720442E-02 1 1 95.833333
+C 11 1.1061417 2.58192486E-02 2 3 87.500000
+C 12 1.2067000 2.81664531E-02 2 5 79.166667
+C 13 1.3072584 3.05136575E-02 5 10 58.333333
+C 14 1.4078167 3.28608619E-02 9 19 20.833333
+C 15 1.5083751 3.52080663E-02 5 24 .000000
+C Summary of preceding table follows: Grouped data Ungrouped data
+C Mean = 1.29887852E+00 1.30997561E+00
+C Variance = 1.92713806E-02 1.73383337E-02
+C Standard deviation = 1.38821398E-01 1.31675107E-01
+BLANK card ending "SYSTEMATIC" output-variable requests
+BEGIN NEW DATA CASE
+C 2nd of 3 subcases is related to 1st in that it demonstrates SYSTEMATIC use.
+C But network is single phase and much simpler. Basically, a 50-Hz generator
+C will charge a capacitor upon switch closure. Illustrate different options:
+C a) Misc. data KSSOUT of column 32 is 3 whereas in DC-16 it was 2. The
+C difference is a phasor branch flow (now seen) for the base case.
+C b) Misc. data MAXOUT of column 40 is zero whereas in DC-16 it was 1. But
+C this makes no difference. The base case show extrema either way.
+C c) Misc. data MEMSAV of column 56 is 2 meaning that program tables will
+C be partially dumped, as 3rd subcase of DC-24 (STATISTICS). For ease of
+C documentation, this will be FORMATTED rather than UNFORMATTED (note
+C FORM= of the $OPEN statement).
+C d) IMAX, read from cols. 25-32 of the STATISTICS/SYSTEMATIC misc. data
+C card, has value 1 to request the addition of normal, non-statistic
+C extrema.
+C e) There is a single SYSTEMATIC switch. NENERG of the misc. data card
+C can be found without the minus sign in column 44 of the switch card.
+C f) Illustrate PARTIAL TABLE DUMPING for faster table dumping. This
+C illustration is added 5 October 2001.
+PARTIAL TABLE DUMPING { Faster table dumping (dump only what is used, not limit)
+PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
+ .000500 .010 { Take 20 steps over half a cycle of the 50-Hz power frequency
+C IOUT IPLOT IDOUBL KSSOUT MAXOUT NENERG
+ 1 1 1 3 0 -1 -4
+C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
+ 1 1 0 1 1 1
+ 5 1 { Printout frequency remains unchanged; print each time step
+ GEN SWIT 1.0 1
+ SWIT 1.0
+ CAP 0.0 0.0 1.0E3 3
+BLANK card ending branch cards
+ SWIT CAP 5.E-3 1.E-3 4 SYSTEMATIC
+BLANK card ending switch cards
+14GEN 400. 50. 0.0 -1.
+BLANK card ending source cards
+$OPEN, UNIT=LUNIT9 FILE=systemat.lu9 STATUS=UNKNOWN FORM=FORMATTED !
+ GEN SWIT CAP
+BLANK card ending output variable requests (node voltages only, here)
+ PRINTER PLOT
+ 184 2. 4.0 10. CAP GEN SWIT
+BLANK card ending plot cards (for base case only)
+$CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
+$OPEN, UNIT=LUNIT9 FILE=dum.lu9 STATUS=UNKNOWN FORM=FORMATTED !
+-2 500. CAP { Tabulate branch current of (CAP, TERRA) using I-base = 500.
+BLANK card ending statistical tabulate (none possible, since no output variable)
+BEGIN NEW DATA CASE
+C 3rd of 3 subcases is same as 2nd except for the following changes:
+C a) Add OMIT BASE CASE to drop the base case (no normal dT-loop output
+C preceding the 4 energizations indicated by NENERG.
+C b) IMAX, read from cols. 25-32 of the STATISTICS/SYSTEMATIC misc. data
+C card, has value 2 to request the replacement of statistical extrema
+C by normal, non-statistic extrema.
+C c) Misc. data MAXOUT of column 40 is 1. But this makes no difference.
+C d) Misc. data KSSOUT of column 32 is 2 as in DC-16. This omits the
+C phasor branch flow of the preceding subcase.
+OMIT BASE CASE { Energizations will not be preceded by normal dT-loop output
+ .000500 .010
+C IOUT IPLOT IDOUBL KSSOUT MAXOUT NENERG
+ 1 1 1 3 1 -1 -4
+C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED
+ 1 1 0 2 1 1
+ 5 5 { Less dT-loop output after time step number 5
+ GEN SWIT 1.0 1
+ SWIT 1.0
+ CAP 0.0 0.0 1.0E3 3
+BLANK card ending branch cards
+ SWIT CAP 5.E-3 1.E-3 4 SYSTEMATIC
+BLANK card ending switch cards
+14GEN 400. 50. 0.0 -1.
+BLANK card ending source cards
+C $OPEN, UNIT=LUNIT9 FILE=systemat.lu9 STATUS=UNKNOWN FORM=FORMATTED !
+ GEN SWIT CAP
+BLANK card ending output variable requests (node voltages only, here)
+$DISABLE { Comment out this card, if a base case solution is to be added
+ PRINTER PLOT
+BLANK card ending plot cards (for base case only)
+$ENABLE { Comment out this card, if a base case solution is to be added
+C $CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE
+-2 500. CAP { Tabulate branch current of (CAP, TERRA) using I-base = 500.
+BLANK card ending statistical tabulate (none possible, since no output variable)
+BEGIN NEW DATA CASE
+BLANK
+EOF
|