From b18347ffc9db9641e215995edea1c04c363b2bdf Mon Sep 17 00:00:00 2001 From: Angelo Rossi Date: Wed, 21 Jun 2023 12:04:16 +0000 Subject: Initial commit. --- benchmarks/dc24.dat | 878 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 878 insertions(+) create mode 100644 benchmarks/dc24.dat (limited to 'benchmarks/dc24.dat') diff --git a/benchmarks/dc24.dat b/benchmarks/dc24.dat new file mode 100644 index 0000000..bd2c786 --- /dev/null +++ b/benchmarks/dc24.dat @@ -0,0 +1,878 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-24 +C Monte Carlo simulation that is run in two parts. The original six +C energizations of "M39." stand-alone case have been reduced to 3, & +C the remaining 3 are performed in companion case DC-40, which is a +C "START AGAIN" continuation of this MEMSAV = 1 case. The use of +C standard random numbers and fixed seed gives results that should be +C repeatable & same for all computers of interest. But M39 differs! +C Note that this data case has no base case. However, if the user +C wants to add this, he need merely comment out the "OMIT BASE CASE" +C request card, and do likewise for the $DISABLE/ENABLE cards that +C precede statistics tabulation requests. Save the resulting binary +C disk file TPTABLES.BIN under the name dc24at40.bin for use by +C DC-40. For a "STATISTICS" random-opening data case, see DC-48. +C Beware of PARALLEL.LIS entry for same USERID that might modify NENERG! +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) +C TRULY RANDOM NUMBERS { Same as NSEED=0 (this overrides stat misc. data card) +REPEATABLE RANDOM NUMBERS { Same as NSEED=1 (this overrides zero value below) +C CENTRAL STATISTICS FILE { Request word that precedes the disk file name +C %/wsm/parallel.atp ! { Leading percent sign avoids confusion of "/" in column 1 +C USER IDENTIFICATION Warsaw ! { Overrides A6 USERID of STARTUP for plots +C DIAGNOSTIC 0 0 0 0 0 0 0 0 0 0 0 9 9 9 9 +OMIT BASE CASE { Comment out this card, if a base case solution is to be added +C Request special output for MS Excel. EOFE service begins 20 July 2007. The +C output file name will be parallel to the input data, with distinctive type : +C Type --- Ruler for file type of Excel output file: +EXTREMA OUTPUT FOR EXCEL .exc { Columns 29-32 carry optional file type +C EXTREMA OUTPUT FOR EXCEL { Blank 29-32 will use default file type .xcl + 100.E-6 20.E-3 60. + 1 1 1 2 1 -1 1 0 3 +C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED + 1 1 0 0 1 { KSTOUT is blank, not 0! } 0 + 2 2 10 10 { Printout frequency change only if base case + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +C Next come 3 switches whose closing times are to be varied by rolling dice. +C The 1st of the 3 STATISTICS switches has type code "76" which means that +C the distribution is uniform with mean of 2 msec and standard deviation of +C 0.1 msec (columns 15 through 34). The 2nd and 3rd switches have type code +C zero, so use Gaussian distribution. The means are 4 and 6 msec, and the +C standard deviations are 1/2 and 1 msec, respectively. +76A1 ASW1 2.E-3 .1E-3 STATISTICS + B1 BSW1 4.E-3 0.5E-3 STATISTICS 4 + C1 CSW1 6.E-3 1.E-3 STATISTICS + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 +C The following is DIAGNOSTIC from "RANDNZ" usage of overlay 12. Not a +C normal case output, this is useful in case switching times are incorrect: +C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow .... +C 0.439822971502571000E+02 0.296996800000000000E+07 +C 0.296996800000000000E+07 0.000691499561071 +C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.205132719793000000E+12 +C 0.326925688100000000E+10 0.761183183873072 +C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.225804303513790000E+15 +C 0.692893886000000000E+09 0.161326929461211 +C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.478574878121350000E+14 +C 0.296220010300000000E+10 0.689690956613049 +C "RFUNL1/RANDNZ". KNT = 1. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.204596198914108000E+15 +C 0.113680185200000000E+10 0.264682306908071 +C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14 +C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565 +C The following unit-9 connection could come earlier on this first subcase +C because of FORM=UNFORMATTED, which is compatible with VECRSV/VECISV of +C MS-DOS (which uses unit 9 for scratch storage). But the present location +C is always safe, so its use is good practice. It is mandatory for 3rd case. +$OPEN, UNIT=LUNIT9 FILE=dc24at40.ext STATUS=UNKNOWN FORM=UNFORMATTED RECL=1000 ! + ASW10 BSW10 CSW10 { Request for these node voltage outputs +C GENA GENC ENDA B1 ASW10 +C Reference angle A1 C1 A10 BSW1 +C +C GENB ENDA +C B1 A10 +C Random switching times for energization number 1 : +C 1 1.8826802E-03 2 4.2472809E-03 3 5.3713577E-03 +C 0.0 141.9164096 -72.1873524 -35.0598126 372.1266197 -429.52043 +C 1.257471517 -1.85970093 +C Times of maxima : .002 .0141 .0127 .0051 .0075 +C .02 .017 +C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow .... +C 0.439822971502571000E+02 0.785177671157890000E+14 +C 0.146997761300000000E+10 0.342255833791569 +C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.101529883752298000E+15 +C 0.115184215400000000E+10 0.268184150103480 +C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.795565857346270000E+14 +C 0.906510819000000000E+09 0.211063497466967 +C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.626117957575120000E+14 +C 0.405748372000000000E+10 0.944706546142697 +C "RFUNL1/RANDNZ". KNT = 2. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.280246343056681000E+15 +C 0.402195997700000000E+10 0.936435530195013 +C Random switching times for energization number 2 : +C 1 1.8999095E-03 2 4.7979502E-03 3 7.5258324E-03 +C 0.0 141.9164096 61.55446116 -41.9189588 371.6253603 334.5616569 +C 1.207725131 -1.83695142 +C Times of maxima : .002 .0154 .0139 .0059 .0028 +C .0058 .017 +BLANK card ending the specification of program outputs (node voltages, here) +C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow .... +C 0.439822971502571000E+02 0.277792753651414000E+15 +C 0.285888072600000000E+10 0.665635039564222 +C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.197460032864095000E+15 +C 0.320639779100000000E+10 0.746547661488876 +C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.221462689026580000E+15 +C 0.129034293200000000E+10 0.300431375391781 +C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.891226959703090000E+14 +C 0.212457830900000000E+10 0.494666935177520 +C "RFUNL1/RANDNZ". KNT = 3. SY, D6, SEEDRN, D18 (result) follow .... +C 0.000000000000000000E+00 0.146742499224322000E+15 +C 0.646589186000000000E+09 0.150545776356012 +C Random switching times for energization number 3 : +C 1 1.9308674E-03 2 3.9933351E-03 3 4.9659084E-03 +C 0.0 136.7926646 -89.6601448 -53.6359842 437.6032257 -510.753735 +C 1.532437698 -1.74295212 +C Times of maxima : .0021 .0051 .0083 .0045 .0074 +C .0044 .0169 +C Cards associated with the base-case solution must not be present: +$DISABLE { Comment out this card, if a base case solution is to be added + PRINTER PLOT + 144 3. 0.0 20. ASW10 +BLANK card ending base-case plot cards +$ENABLE { Comment out this card, if a base case solution is to be added +C 1 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 27 1.3500000 0.13500000E+03 0 0 +C 28 1.4000000 0.14000000E+03 1 1 +C 29 1.4500000 0.14500000E+03 2 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.40833333E+00 1.40208495E+00 +C Variance = 8.33333333E-04 8.75092111E-04 +C Standard deviation = 2.88675135E-02 2.95819558E-02 +C The following disconnects DC24.BIN to protect; automatic result of reopening: +$OPEN, UNIT=LUNIT2 STATUS=SCRATCH FORM=UNFORMATTED { Disconnect *.BIN on LUNIT2 +-1 100. GENA A1 +$CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE +$OPEN, UNIT=LUNIT9 FILE=dc24dum.dum STATUS=SCRATCH ! +FIND { Enter interactive search for which shot and which variable gave extremum +HELP { Show the user that there is a summary explanation of "FIND" on-line +C Peak extremum of subset has value 5.62169583E+02 This occurred during +C energization 2 for the variable having names "CSW10 " and " ". +C This was variable 3 of 3 in the subset; 3 shots were searched. + ASW10 BSW10 CSW10 { Search these 3 node voltages for extremum +EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot +C Peak extremum of subset has value 5.59424445E+02 This occurred during +C energization 1 for the variable having names "CSW10 " and " ". +C This was variable 3 of 3 in the subset; 2 shots were searched. +C Remember excluded shots are numbered 2 +C Add the following DISK command on 29 September 2002. Look in the January, +C 2003 or later newsletter for a summary. Only the framework of the feature +C now exists in DICTAB, but it is being illustrated here in order that it not +C be forgotten. The idea is this. Each EXCLUDE command excludes another shot +C (number 2 for the preceding). Deterministic simulation of this often is of +C interest, and DISK is designed to create such a data file. The name will +C be SHOT0002.DAT and this file will be located along with the .LIS file, +C assuming parallel naming. +C DISK { Create a disk file of data for deterministic simulation of EXCLUDEd shot +C 6 October 2002, the RESIMULATE command is working at home as illustrated +C by f:\data\DUM24.DAT First, the preceding DISK is enhanced by /LIST: +C 8 October 2002, /FULL is created to use full 10 bytes of T-close or open +C as opposed to the default 9, with one blank on the left as separator. +DISK /LIST /FULL { Create disk file of data for deterministic simulation of EXCLUDEd shot +C RESIMULATE { Abort this data case. Begin with SHOT0002.DAT just created by DISK +C If the preceding line is activated, all following data will be ignored as +C the present data is aborted in favor of the DISK-created SHOT0002.DAT But +C we do not want to lose all the following output. So, comment this line +C here and note that a new 4th subcase of DC-48 illustrates RESIMULATE use. +EXCLUDE { Excluding the shot of previous output means only one is left for use +C Peak extremum of subset has value 5.53096221E+02 This occurred during +C energization 3 for the variable having names "BSW10 " and " ". +C This was variable 2 of 3 in the subset; 1 shots were searched. +C Remember excluded shots are numbered 1 2 +RESET { Reinitialize by erasing history of previous selections +ALL { Next, consider searching all output variables of some class +8 { As with batch-mode plotting, "8" indicates branch voltages (search ALL) +C Peak extremum of subset has value 1.41916410E+02 This occurred during +C energization 1 for the variable having names "GENA " and "A1 ". +C This was variable 1 of 3 in the subset; 3 shots were searched. +EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot +EXCLUDE { Excluding the shot of previous output means only one is left for use +EXCLUDE { Final attempt to exclude should lead to error message: nothing is left +ALL { Next, consider searching all output variables of some class +9 { As with batch-mode plotting, "9" indicates branch currents (search ALL) +EXCLUDE { Repeat the preceding search after 1st excluding the just-found shot +C $KEY { Un-comment to allow keyboard input (until "END") at this point +QUIT { Done locating extrema, so exit the "FIND" alternative +C &&&&&&&&&&&&& End request and documentation. Begin 3 coupled cards ...... +C 2 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 9 0.4500000 0.13635000E+03 0 0 +C 10 0.5000000 0.15150000E+03 3 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 4.75000000E-01 4.62734306E-01 +C Variance = 0.00000000E+00 9.53165933E-05 +C Standard deviation = 0.00000000E+00 9.76302173E-03 +-1 GENA A1 CONT. +C 3 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENC " to "C1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 4 0.2000000 0.60600000E+02 0 0 +C 5 0.2500000 0.75750000E+02 2 2 +C 6 0.3000000 0.90900000E+02 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 2.41666667E-01 2.45766731E-01 +C Variance = 8.33333333E-04 2.19347831E-03 +C Standard deviation = 2.88675135E-02 4.68345845E-02 +-1 GENC C1 CONT. +C 4 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "ENDA " to "A10 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 2 0.1000000 0.30300000E+02 0 0 +C 3 0.1500000 0.45450000E+02 2 2 +C 4 0.2000000 0.60600000E+02 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.41666667E-01 1.43690600E-01 +C Variance = 8.33333333E-04 9.61072666E-04 +C Standard deviation = 2.88675135E-02 3.10011720E-02 +-1 ENDA A10 +C 5 ) -------------------------------------------------------------------------- +C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR +C 5 ) -------------------------------------------------------------------------- +C A distribution of peak values among all output branches of the last data card +C statistical distribution is for the maximum of the peaks at all of these outpu +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 9 0.4500000 0.13635000E+03 0 0 +C 10 0.5000000 0.15150000E+03 3 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 4.75000000E-01 4.62734306E-01 +C Variance = 0.00000000E+00 9.53165933E-05 +C Standard deviation = 0.00000000E+00 9.76302173E-03 +C 3456789012345678901234567890123456789012345678 +C MODTAB AINCR XMAXMX +STATISTICS DATA 1 0.0 0.0 ---- Only individual +-1 GENC C1 { =.245766731, Variance=2.19347831E-3} CONT. +-1 ENDA A10 { 1.43690600E-01, 9.61072666E-04 +C 8 ) -------------------------------------------------------------------------- +C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR +C 8 ) -------------------------------------------------------------------------- +C A distribution of peak values among all output branches of the last data card +C statistical distribution is for the maximum of the peaks at all of these outpu +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 4 0.2000000 0.60600000E+02 0 0 +C 5 0.2500000 0.75750000E+02 2 2 +C 6 0.3000000 0.90900000E+02 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 2.41666667E-01 2.45766731E-01 +C Variance = 8.33333333E-04 2.19347831E-03 +C Standard deviation = 2.88675135E-02 4.68345845E-02 +STATISTICS DATA 2 0.0 0.0 ------ Only union of 2 curves +-1 GENC C1 CONT. +-1 ENDA A10 { See comments for solution +STATISTICS DATA 3 0.0 0.0 ---- Both individual & union +-1 GENC C1 { =.245766731, Variance=2.19347831E-3} CONT. +-1 ENDA A10 { 1.43690600E-01, 9.61072666E-04 +C 12 ) ------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENC " to "C1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 8 0.9131973 0.56211367E+02 0 0 +C 9 1.0273470 0.63237788E+02 1 1 +C 10 1.1414966 0.70264209E+02 0 1 +C 11 1.2556463 0.77290630E+02 1 2 +C 12 1.3697959 0.84317051E+02 0 2 +C 13 1.4839456 0.91343472E+02 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.19857145E+00 1.20977941E+00 +C Variance = 5.21205808E-02 5.31495428E-02 +C Standard deviation = 2.28299323E-01 2.30541846E-01 +STATISTICS DATA 3 -5.0 0.0 ---- Fixed number of 5 boxes +-1 GENC C1 { See preceding comments for 1st of 3} CONT. +-1 ENDA A10 +STATISTICS DATA 3 .05 0.0 ---- Restore STARTUP values +-1 GENC C1 ENDA A10 { =2.45766731E-1, Variance = 2.19347831E-3 + 0 ASW10 BSW10 CSW10 { Mean=1.40245965E+0, Variance = 8.47042357E-2 +C 1st of 4 tables is documented on preceding request; the 4th of 4 follows: +C 21 ) ------------------------------------------------------------------------- +C SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMARY SUMMAR +C 21 ) ------------------------------------------------------------------------- +C The following is a distribution of peak overvoltages among all output nodes of +C This distribution is for the maximum of the peaks at all output nodes with V- +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 36 1.8000000 0.54540000E+03 0 0 +C 37 1.8500000 0.56055000E+03 2 2 +C 38 1.9000000 0.57570000E+03 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.84166667E+00 1.84234351E+00 +C Variance = 8.33333333E-04 2.35830395E-04 +C Standard deviation = 2.88675135E-02 1.53567703E-02 +C &&&&&&&&& End documentation of node voltage tabulations (4th of 4 tables) +STATISTICS DATA 3 -4.0 0.0 ---- Fixed number of 4 boxes +C 22 ) ------------------------------------------------------------------------- +C Statistical distribution of peak power for branch "B1 " to "BSW1 ". +C Interval power power in Frequency Cumulative +C number in per unit physical units (density) frequency +C 16 0.9468728 0.35188195E+03 0 0 +C 17 1.0060524 0.37387457E+03 2 2 +C 18 1.0652319 0.39586719E+03 0 2 +C 19 1.1244115 0.41785981E+03 0 2 +C 20 1.1835910 0.43985244E+03 1 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.03564214E+00 1.05962916E+00 +C Variance = 1.05066576E-02 1.04274409E-02 +C Standard deviation = 1.02501988E-01 1.02114842E-01 +-3 B1 BSW1 { See preceding comments for this branch power table +-2 0. GENB B1 ENDA A10 { =1.10335104E+0, Variance = 2.09698061E-2 +-4 0. B1 BSW1 { Ungroup: Mean=1.16595031E+0, Variance = 3.80542557E-2 +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +C 2nd of 6 subcases repeats the preceding problem exactly, only instead +C of rolling dice to find the switch closing times, such times will be +C read from a user-supplied disk file dc24lun8.dat (see $OPEN use). +C The solution of this 2nd subcase should exactly agree with the first. +C However , for variety , a base case solution (missing from the first +C subcase) and a character plot have been added to this 2nd subcase. +C Also, the ".PL4" plot file (actually, just the header) is saved +C permanently for later re-connection using the 4th subcase of DC-40. +PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP +C $STARTUP, dc37star.dat { Use disk file for re-initialization at this point +C This is a Monte Carlo case, for which TENERG is critical. But this +C parameter is initialized only in the STARTUP file. For this 2nd or +C later data case, we reinitialize via the preceding $STARTUP request +USER SUPPLIED SWITCH TIMES { Bypass the dice, reading TCLOSE & TOPEN from LUNIT8 +C The following connection of unit 2 is required only because this is the +C 2nd subcase. It follows a case in which MEMSAV=1 dumps tables via I/O +C unit LUNIT2=2, and then this file is disconnected. Without the manual +C connection that follows, nothing would be connected, and execution would +C die for Apollo in "SUBR5" where REWIND LUNIT2 is found. So connect it: +$OPEN, UNIT=LUNIT8 FILE=dc24lun8.dat FORM=FORMATTED STATUS=OLD RECL=80 ! { Case +$OPEN, UNIT=LUNIT2 FORM=UNFORMATTED STATUS=SCRATCH RECL=16000 +C $STARTUP, dc37star.dat { Re-initialize is required to set TENERG for table dump +C { This is a peculiarity of having "STATISTICS" as the +C { 2nd or later data subcase within a single disk file. +C DISK PLOT DATA { Toggle the Apollo default of LUNIT4 = -4 to +4 (use disk) +$CLOSE, UNIT=LUNIT4 STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP" +$OPEN, UNIT=LUNIT4 FILE=dc24b40d.pl4 ! { All we need for C-like case +C Before continuing with case, let's document the contents of that unit-8 file: +C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 1: +C 1.882680E-03 4.247281E-03 5.371358E-03 7.95E-3 9.95E-3 +C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line: +C .1E-3 0.5E-3 1.E-3 1.0 1.0 +C 1.0 +C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 2: +C 1.899909E-03 4.797950E-03 7.525832E-03 7.95E-3 9.95E-3 +C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line: +C .1E-3 0.5E-3 1.E-3 1.0 1.0 +C 1.0 +C Marker preceding TCLOSE(1:6) and TOPEN(1:6) for energization 3: +C 1.930867E-03 3.993335E-03 4.965908E-03 7.95E-3 9.95E-3 +C 11.95E-3 { End of TCLOSE; TOPEN begin on the next line: +C .1E-3 0.5E-3 1.E-3 1.0 1.0 +C 1.0 +C 345678901234567890123456789012345678901234567890123456789012345678901234567890 +C Cancel special output for MS Excel of 1st subcase. Addition on 20 July 2007: +C Type --- Ruler for file type of Excel output file: +EXTREMA OUTPUT FOR EXCEL .off { Columns 29-32 carry optional file type + 100.E-6 20.E-3 60. + 1 1 1 2 1 -1 0 2 3 +C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED + 1 1 0 0 1 { KSTOUT is blank, not 0! } 1 + 2 2 10 10 { Printout frequency change only if base case + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +76A1 ASW1 1.95E-3 .1E-3 STATISTICS + B1 BSW1 3.95E-3 0.5E-3 STATISTICS 4 + C1 CSW1 5.95E-3 1.E-3 STATISTICS + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 +C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14 +C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565 + ASW10 BSW10 CSW10 +C Step Time GENA GENC ENDA B1 ASW10 BSW10 +C A1 C1 A10 BSW1 +C 0 0.0 -.5684E-13 .28422E-13 -.5684E-13 -151.5 0.0 0.0 +C 1 .1E-3 .56843E-13 .29559E-11 -.2842E-12 0.0 0.0 0.0 +C 2 .2E-3 0.0 -.2927E-11 .28422E-12 0.0 0.0 0.0 +BLANK card ending the specification of program outputs (node voltages, here) +C Last step: 200 .02 1.60684633 5.46144739 15.5569682 209.921255 126.692915 +C Variable max:136.792665 67.5308064 48.9224823 437.603226 320.835627 553.096221 +C Times of max: .0021 .0061 .0093 .0045 .0029 .0049 +C Variable min:-100.37615 -85.515782 -58.824356 -292.03358 -548.79093 -353.02344 +C Times of min: .0033 .0127 .0084 .0056 .0074 .0141 + PRINTER PLOT + 144 1. 0.0 10. ASW10 { Plot limits: (-5.488, 3.208) +$CLOSE, UNIT=LUNIT4 STATUS=KEEP +$OPEN, UNIT=LUNIT4 FILE=dc24bdum.pl4 ! { All we need for C-like case +BLANK card ending base-case plot cards +C GENA GENC ENDA B1 ASW10 +C Reference angle A1 C1 A10 BSW1 +C Random switching times for energization number 1 : +C 1 1.8826800E-03 2 4.2472810E-03 3 5.3713580E-03 +C 0.0 141.91641 -72.187352 -35.059813 372.12662 -429.52043 422.745884 +C Times of maxima : .002 .0141 .0127 .0051 .0075 .0051 +C +C Random switching times for energization number 2 : +C 1 1.8999090E-03 2 4.7979500E-03 3 7.5258320E-03 +C 0.0 141.91641 61.5544612 -41.918959 371.62536 334.561657 497.525193 +C Times of maxima : .002 .0154 .0139 .0059 .0028 .0065 +C +C Random switching times for energization number 3 : +C 1 1.9308670E-03 2 3.9933350E-03 3 4.9659080E-03 +C 0.0 136.792665 -89.660145 -53.635984 437.603226 -510.75374 +C cont. 553.096221 513.648909 .96346791 1.5324377 -1.7429521 +C .0021 .0051 .0083 .0045 .0074 +C cont. .0049 .011 .02 .0044 .0169 +C Cards associated with the base-case solution must not be present: +C 1 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 27 1.3500000 0.13500000E+03 0 0 +C 28 1.4000000 0.14000000E+03 1 1 +C 29 1.4500000 0.14500000E+03 2 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.40833333E+00 1.40208495E+00 +C Variance = 8.33333333E-04 8.75092111E-04 +C Standard deviation = 2.88675135E-02 2.95819558E-02 +-1 100. GENA A1 { See preceding comments that document this tabulation +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +C 3rd of 6 subcases repeats the first subcase exactly, only instead of +C dumping full tables, MEMSAV = 2 (columns 49-56) means that only +C the extrema and switching times are to be dumped, and +C this will be FORMATTED rather than UNFORMATTED (for +C this change, only $OPEN is changed). The resulting card-image file +C named dc24ct40.lis is used by the 2nd and later subcases of DC-40. +C Also, as an illustration, time-step loop output of all energizations +C will be shown (because column-48 variable KSTOUT is zero not blank), +C and the miniature character plots of switch times will be suppressed. +PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP +C $STARTUP, dc37star.dat { Use disk file for re-initialization at this point +C This is a Monte Carlo case, for which TENERG is critical. But this +C parameter is initialized only in the STARTUP file. For this 2nd or +C later data case, we reinitialize via the preceding $STARTUP request +OMIT BASE CASE { Comment out this card, if a base case solution is to be added + 100.E-6 20.E-3 60. + 1 1 1 2 1 -1 2 0 3 +C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED +C 1 1 0 0 1 0 { STATIST. misc. data card } 1 + 1 1 0 0 2309 0 { STATIST. misc. data card } 1 + 2 2 10 10 { Printout frequency change for all soutions +C Preceding STATISTICS miscellaneous data card is handled differently beginning +C 23 September 2005 as new reserved value IDICE = 2309 is introduced as a +C request for alternative, computer-dependent dice. This is to satisfy Apple +C Macintosh G5 as used by Stu Cook of JUST Services in Ontario, Canada. As of +C this date, Mac G5 ATP is the only version to honor the 2309 alternative. +C For all versions, any other IDICE value (including the value unity shown) +C is treated the same as zero or blank. I.e., any other value does nothing. +C Henceforth, IDICE = 2309 is a request for alternative, computer-dependent +C random numbers (dice) of the operating system's library. So, we illustrate +C use here to document the two associated warning lines immediately before the +C input of branch data begins. If Stu Cook wants to compare his results with +C everybody else's, he should restore the value IDICE = 1 (alternatively, +C zero or blank would be equivalent). Salford ATP output follows: +C Statistics data. 1 1 0 0 2309 1 | 1 1 ... +C Printout : 2 2 10 10 0 0 | 2 2 ... +C Computer-dependent (not universal) dice have been requested for this +C Monte Carlo study by reserved value IDICE = 2309. +C Sorry, no alternative computer-dependent dice are available. ATP will +C ignore this request. Universal dice will be used. +C Series R-L-C. 0.000E+00 1.857E-02 0.000E+00 | 0GENA A1 ... + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +76A1 ASW1 2.E-3 .1E-3 STATISTICS + B1 BSW1 4.E-3 0.5E-3 STATISTICS 4 + C1 CSW1 6.E-3 1.E-3 STATISTICS + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 +C Last gen: ENDC -103.6321034277 303. -.710542736E-14 .7324106878E-14 +C Last gen: 284.72686409813 110.000 -.177635684E-14 -165.9637565 +C The following unit-9 connection must not be made earlier, since the same +C I/O channel is used as UNFORMATTED by VECRSV/VECISV of MS-DOS. This is +C different than the first subcase (read those comments), note: +$OPEN, UNIT=LUNIT9 FILE=dc24ct40.lis STATUS=UNKNOWN FORM=FORMATTED ! + ASW10 BSW10 CSW10 { Request for these node voltage outputs +C GENA GENC ENDA B1 ASW10 BSW10 +C Reference angle A1 C1 A10 BSW1 +C Random switching times for energization number 1 : +C 1 1.8826802E-03 2 4.2472809E-03 3 5.3713577E-03 +C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 +C 1 .1E-3 0.0 0.0 0.0 0.0 0.0 0.0 +C 2 .2E-3 0.0 0.0 0.0 0.0 0.0 0.0 +C 4 .4E-3 0.0 0.0 0.0 0.0 0.0 0.0 +C 6 .6E-3 0.0 0.0 0.0 0.0 0.0 0.0 +C 8 .8E-3 0.0 0.0 0.0 0.0 0.0 0.0 +C 10 .001 0.0 0.0 0.0 0.0 0.0 0.0 +C ==== Table dumping for all subsequent restorations. Time [sec] = 1.8000E-03 +C *** Close switch "A1 " to "ASW1 " after 1.90000000E-03 sec. +C 20 .002 141.91641 0.0 0.0 0.0 0.0 0.0 +C 30 .003 -34.73764 0.0 0.0 0.0 318.450266 -159.22513 +BLANK card ending the specification of program outputs (node voltages, here) +C 190 .019 -8.1503339 19.8420817 .549065534 24.9565797 230.197042 85.0661304 +C 200 .02 .352725683 14.1822201 13.4881293 207.346558 128.761754 181.144509 +C 201 .0201 -.23343325 23.5035485 8.73828958 199.126365 123.327149 204.30575 +C Extrema of output variables follow. Order and column positioning are the sam +C Variable max:136.792665 56.3388361 40.5581876 437.603226 320.835627 553.096221 +C Times of max: .0021 .0053 .0092 .0045 .0029 .0049 +C Variable min:-100.37615 -89.660145 -53.635984 -272.03803 -510.75374 -342.9791 +C Times of min: .0033 .0051 .0083 .0055 .0074 .014 +C Warning! No statistical tabulations are allowed here because of our use +C ------- of output file LUNIT9 = 9 being "FORMATTED". Any attempted +C tabulation will die in "MCBANK" called by "DICTAB". 7 Dec 88 +C -1 100. GENA A1 +$CLOSE, UNIT=LUNIT9 STATUS=KEEP { Disconnect to prevent damage by stray WRITE +$OPEN, UNIT=LUNIT9 FILE=dc24dum.dum STATUS=UNKNOWN FORM=FORMATTED RECL=1000 ! +C The preceding $OPEN is added 11 October 2010 as required by new 4th subcase +C 1 ) -------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". +C Interval voltage voltage in Frequency Cumulative +C number in per unit physical units (density) frequency +C 27 1.3500000 0.13500000E+03 0 0 +C 28 1.4000000 0.14000000E+03 1 1 +C 29 1.4500000 0.14500000E+03 2 3 +C Distribution parameters for .. data follow. Grouped data Ungrouped data +C Mean = 1.40833333E+00 1.40208495E+00 +C Variance = 8.33333333E-04 8.75092111E-04 +C Standard deviation = 2.88675135E-02 2.95819558E-02 +C Column positioning of "NO SWITCH PLOTS" on blank terminator is arbitrary: +BLANK card ending statistical tabulation requests --- NO SWITCH PLOTS --- +BEGIN NEW DATA CASE +C 4th of 6 subcases produces simulations that are identical to the 2nd. +C However, output looks substantially different because STATISTICS +C switches are not used. Instead of USER SUPPLIED SWITCH TIMES to +C define the closing times of switches, the SELECT function of the +C pocket calculator will be used. Structural differences include lack +C of a base case (the 2nd subcase had 1) preceding the 3 energizations +C Output of connectivity and the phasor solution will be suppressed so +C as not to be repeated 3 times (boring). Batch-mode plotting will be +C omitted for the same reason. Finally, output of the dT-loop will be +C be made less frequent. Data is added 10 October 2010. WSM. +PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP +C Controls of following request card: MAXKNT IOPCVP NOSTAT { Loop 3 times with +POCKET CALCULATOR VARIES PARAMETERS 3 0 0 { full printout + 100.E-6 20.E-3 60. + 20 1 0 0 1 + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +TCLOSE_AAA = SELECT ( KNT ) { Closing time of phase "a" is taken from following: + .00188268.001899909.001930867 { Data is 3E10.0 for 3 energizations of phase a +C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase a +TCLOSE_BBB = SELECT ( KNT ) { Closing time of phase "b" is taken from following: +.004247281 .00479795.003993335 { Data is 3E10.0 for 3 energizations of phase b +C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase b +TCLOSE_CCC = SELECT ( KNT ) { Closing time of phase "c" is taken from following: +.005371358.007525832.004965908 { Data is 3E10.0 for 3 energizations of phase c +C Pass 1>< Pass 2 >< Pass 3 > Ruler for preceding 3 T-close values of phase c +BLANK card ends $PARAMETER definitions that are processed just b4 branch card +C Switch cards begin with the 3 STATISTICS switches of the 2nd subcase, which +C are being converted to deterministic switches with T-close of columns 15- +C 24 defined by the preceding $PARAMETER block. T-open has been set to 99 +C seconds (never reached) to ensure that the switches never open: +C 76A1 ASW1 1.95E-3 .1E-3 STATISTICS + A1 ASW1 TCLOSE_AAA 99.0 +C B1 BSW1 3.95E-3 0.5E-3 STATISTICS 4 + B1 BSW1 TCLOSE_BBB 99.0 4 +C C1 CSW1 5.95E-3 1.E-3 STATISTICS + C1 CSW1 TCLOSE_CCC 99.0 + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 +C The following start of the dT-loop output is for the 1st of the 3 passes: +C Column headings for the 10 EMTP output variables follow. These are divided among the 5 possible classes as follows .... +C First 7 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Branch power consumption (power flow, if a switch) is treated like a branch voltage for this grouping (1 variables); +C Branch energy consumption (energy flow, if a switch) is treated like a branch current for this grouping (1 variables). +C Step Time GENA GENC ENDA B1 ASW10 BSW10 CSW10 B1 GENB ENDA +C A1 C1 A10 BSW1 BSW1 B1 A10 +C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 +C *** Close switch "A1 " to "ASW1 " after 1.90000000E-03 sec. +C 20 .002 141.91641 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 +C 40 .004 2.37028418 0.0 0.0 0.0 -46.220806 23.1104031 23.1104031 0.0 0.0 0.0 +C *** Close switch "B1 " to "BSW1 " after 4.30000000E-03 sec. +C *** Close switch "C1 " to "CSW1 " after 5.40000000E-03 sec. +C 60 .006 -15.795936 10.618273 0.0 -123.04522 -168.96416 422.15175 -253.18759 .357709401 -.40883323 0.0 + ASW10 BSW10 CSW10 +BLANK card ending the specification of program outputs (node voltages, here) +C 180 .018 7.07846743 32.9539555 -12.904568 -.92594731 299.739293 -60.943514 -240.37687 .62556771 .018267223 -1.498213 +C 200 .02 5.74620332 -11.67979 16.6256782 248.357719 125.624205 168.195467 -283.39188 .769701668 1.25747152 -1.2455527 +C Variable maxima : 141.91641 56.029476 33.4968354 372.12662 334.561657 422.745884 559.424445 .769701668 1.25747152 .049072408 +C Times of maxima : .002 .0122 .0092 .0051 .0028 .0051 .0114 .02 .02 .0095 +C Variable minima : -103.91401 -72.187352 -35.059813 -191.64938 -429.52043 -351.293 -316.83918 0.0 -1.0806309 -1.8597009 +C Times of minima : .0032 .0141 .0127 .0159 .0075 .0142 .0191 0.0 .0109 .017 +BLANK card ending batch-mode plot cards +C Following statistical tabulation is identical to that of 2nd subcase of DC-24: +C 1 ) ------------------------------------------------------------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". Base voltage for per unit output = 1.00000000E+02 +C Interval voltage voltage in Frequency Cumulative Per cent +C number in per unit physical units (density) frequency .GE. current value +C 27 1.3500000 1.35000000E+02 0 0 100.000000 +C 28 1.4000000 1.40000000E+02 1 1 66.666667 +C 29 1.4500000 1.45000000E+02 2 3 .000000 +C Summary of preceding table follows: Grouped data Ungrouped data +C Mean = 1.40833333E+00 1.40208495E+00 +C Variance = 8.33333333E-04 8.75092111E-04 +C Standard deviation = 2.88675135E-02 2.95819558E-02 +-1 100. GENA A1 { See preceding comments that document this tabulation +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +C 5th of 6 subcases produces simulations that are identical to the 4th. +C However, instead of the SELECT function (which is a generalization +C of USER SUPPLIED SWITCH TIMES), dice will be rolled just as was the +C case for the 1st subcase. Yes, the same dice, with the RAN function +C and the GAUSS function of the pocket calculator replacing internal +C rolling that ATP performs in support of STATISTICS switches. Within +C the PCVP loop, the user must roll his own dice after first ensuring +C that seeding of the random number generator is identical to what ATP +C does for the 1st subcase. Data is added 16 October 2010. WSM. +PRINTED NUMBER WIDTH, 11, 1, { Set dT-loop column width as done b4 by STARTUP +REPEATABLE RANDOM NUMBERS { Needed since no NSEED of STATISTICS misc. data card +C Note about preceding line. Without it, dice no longer will be deterministic +C and the resulting .LIS file should change every time data is simulated. +C That assumes KOMPAR of STARTUP is small enough (e.g., value unity) so that +C the program date and time are real. But for KOMPAR = 4, all simulations +C will be identical because seeding of the number generator then will depend +C on a fixed time (the WW I Armistice time) and zero elapsed CP time. I.e., +C program output will be deterministic although different from the output for +C the data as shown. That is because the fixed seed of RRN differs from the +C fixed seed that results from the WW I time. Complicated, eh? +C Controls of following request card: MAXKNT IOPCVP NOSTAT { Loop 3 times with +POCKET CALCULATOR VARIES PARAMETERS 3 0 0 { full printout + 100.E-6 20.E-3 60. + 20 1 0 0 1 + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +C Diagnostic printout of the 1st subcase reveals that dice are rolled 5 times +C prior to each energization. There are 3 STATISTICS switches, so 3 of the 5 +C are used for them. But prior to these, there are 2 more. First, there is +C a wasted roll that is used for nothing. This has been known for many years +C but was not removed for reasons of historical continuity. Second, there is +C a roll that is used for the uniform bias that might, if the user wants it, +C be added to the closing times of all switches. So, here are the 5 rolls : +UNUSE1of2 = GAUSS ( 0.0 ) { OVER12 wastes 1 use for no good reason if STATISTIC +UNUSE2of2 = GAUSS ( 0.0 ) { OVER12 rolls a 2nd time for bias if STATISTICS case +C Now we are ready for the 3 switches. The 1st has type code 76, which means +C a uniform distribution. This is a little involved because it is low level. +C The 1st subcase uses a mean of 2 msec and a standard deviation of 0.1 msec. +C But the RAN function produces a uniform distribution over (0, 1). To map +C the RAN output RNG into the desired closing time, OVER12 uses this code: +C RNG = TIMEV + BIAS + SQRT3 * ABSZ ( TOPEN(NDXI) ) * ( 2.0 * RNG - 1.0 ) +C For this case, there is no bias BIAS, TIMEV is the mean equal to 2 msec, +C SQRT3 is the square root of 3, and TOPEN is the standard deviation. So: +TCLOSE_AAA = 2.0E-3 + 0.17320508E-3 * ( 2.0 * RAN ( 0.0 ) - 1.0 ) { Uniform closing time of phase "a" +C This is slightly more efficient, using the constant 0.17320508E-3 rather +C than scaling an ATP variable. But the answers should be identical. ATP is +C smart enough to recognize SQRT3 as a program variable, and use the full, +C internal precision associated with it. Program output will be unchanged if +C the preceding were replaced by: +C TCLOSE_AAA = 2.0E-3 + SQRT3 * 1.E-4 * ( 2.0 * RAN ( 0.0 ) - 1.0 ) { Uniform closing time of phase "a" +C Switches 2 and 3 (phases b and c) are easier since they are Gaussian: +TCLOSE_BBB = 4.0E-3 + 0.5E-3 * GAUSS ( 0.0 ) { Gauss closing time of phase "b" +TCLOSE_CCC = 6.0E-3 + 1.0E-3 * GAUSS ( 0.0 ) { Gauss closing time of phase "c" +BLANK card ends $PARAMETER definitions that are processed just b4 branch card +C 76A1 ASW1 2.E-3 .1E-3 STATISTICS + A1 ASW1 TCLOSE_AAA 99.0 +C B1 BSW1 4.E-3 0.5E-3 STATISTICS 4 + B1 BSW1 TCLOSE_BBB 99.0 4 +C C1 CSW1 6.E-3 1.E-3 STATISTICS + C1 CSW1 TCLOSE_CCC 99.0 + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 + ASW10 BSW10 CSW10 +BLANK card ending the specification of program outputs (node voltages, here) +BLANK card ending batch-mode plot cards +C Following statistical tabulation is identical to that of 4th subcase of DC-24: +C 1 ) ------------------------------------------------------------------------------------------------------------------------------- +C Statistical distribution of peak voltage for branch "GENA " to "A1 ". Base voltage for per unit output = 1.00000000E+02 +C Interval voltage voltage in Frequency Cumulative Per cent +C number in per unit physical units (density) frequency .GE. current value +C 27 1.3500000 1.35000000E+02 0 0 100.000000 +C 28 1.4000000 1.40000000E+02 1 1 66.666667 +C 29 1.4500000 1.45000000E+02 2 3 .000000 +C Summary of preceding table follows: Grouped data Ungrouped data +C Mean = 1.40833333E+00 1.40208495E+00 +C Variance = 8.33333333E-04 8.75092111E-04 +C Standard deviation = 2.88675135E-02 2.95819558E-02 +-1 100. GENA A1 { See preceding comments that document this tabulation +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +C 6th of 6 subcases is added 24 July 2011 to illustrate the new +C protection against overflow of List 23. Data is the same as for +C subcase 1 although the base case has been activated to prove that +C omission is not a requirement. Note the NENERG = 90K shots of +C the integer misc. data card. For other translations, this is a +C case of overkill. But for Salford ATP, ordered COMMON blocks +C result in an enormous List 23 in excess of a million words. For +C Salford ATP, such huge NENERG is necessary. WSM. + 100.E-6 20.E-3 60. + 50000 20001 0 2 1 -1 1 0 90000 +C Interpretation of the preceding integer misc. data card is improved on +C 7 August 2011 when alternative variable-width encoding and the use of "K" +C allow the handling of much larger integers. The old format was (2I5, 8I3) +C which resulted in overflow for NENERG > 999. The new is not troubled : +C Misc. data. 1.000E-04 2.000E-02 6.000E+01 | 100.E-6 20.E-3 60. +C Misc. data. 50K 20K 0 2 1 -1 1 0 90K 0 | 50000 20001 ... +C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED + 1 1 0 0 1 { KSTOUT is blank, not 0! } 0 + 2 2 10 10 { Printout frequency change only if base case + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 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 +76A1 ASW1 2.E-3 .1E-3 STATISTICS + B1 BSW1 4.E-3 0.5E-3 STATISTICS 4 + C1 CSW1 6.E-3 1.E-3 STATISTICS + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-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 + ASW10 BSW10 CSW10 { Request for these node voltage outputs +BLANK card ending the specification of program outputs (node voltages, here) +C Overflow of List 23. General information follows. Specifically, here in SUBR15 prior to entry into the time-step loop, the demand +C for storage that is required for a Monte Carlo (or a SYSTEMATIC) study has been computed, and has been found to be inadequate. +C Total working space is L23TOT = 1502629 cells whereas each shot requires NDX1 = 24 cells. For all NENERG = 90000 shots +C requested, at least 2160000 will be needed for List 23. But recall that a Monte Carlo study can be simulated in pieces using two +C or more separate ATP executions. It is safer not to place all eggs in one basket. Also, simulation should be faster if two or +C more central processors are used simultaneously. CENTRAL STATISTICS FILE is the request word for such operation as illustrated by +C BENCHMARK DC-66. +C ------------------------------------------------------------------------------------------------------------------------------------ +C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ +C Etc. (the familiar KILL = 1 error termination). +BLANK card ending base-case plot cards +-1 GENA A1 +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +BLANK + -- cgit v1.2.3