BEGIN NEW DATA CASE C BENCHMARK DC-13 C Illustration of SUBROUTINE SATURA, the magnetic saturation routine C that converts an RMS saturation curve of voltage vs. current into the C corresponding locus of instantaneous flux vs. current. Also shown C are B-H hysteresis generation and ZnO conversion from old formats. C 12 data subcases are present, so many other features are tested, too SATURATION $ERASE -1. 10. .001 1 0.0 5.0 2.0 5.0 3.0 3.5 4.0 2.0 5.0 1.0 10.0 1.0 9999 C Final row of resulting printout: 6 100.0000000000 0.2350000000 60. .664 .01 1 1.0 100. 1.6 105. 2.75 110. 4.0 113. 6.0 116. 10. 119. 44. 131. 9999 C Final row of resulting printout: 8 1601.6892070492 326.3046435192 $PUNCH, dc13a.pch ! { Exclamation holds lower case -1. 10. .001 1 2.0 5.0 9999 60. 303.11 300. 1 1 .000306 0.7 .000415 0.8 .00076 0.9 .00198 1.0 .0108 1.1 9999 C Final row of resulting printout: 6 26.6610642630 1250.7681935392 $PUNCH 60. 34.5 300.0 1 .000306 0.7 .000415 0.8 .00076 0.9 .00198 1.0 .0108 1.1 9999 C Final row of resulting printout: 6 234.2387011236 142.3625174923 $PUNCH, dc13b.pch ! { Exclamation holds lower case C 3.76303783E+00 9.05943293E+01 C 6.10062069E+00 1.03536376E+02 C 1.33222222E+01 1.16478423E+02 C 3.87652515E+01 1.29420470E+02 C 2.34238701E+02 1.42362517E+02 BLANK card ending all "SATURATION" data cases BEGIN NEW DATA CASE C 2nd of 12 subcases will illustrate use of another supporting program: OLD TO NEW ZNO { Special request card to convert old ZnO data to new format .000050 .020000 { Note: this data case once was numbered DC-38B. 1 1 1 0 1 -1 5 5 20 1 30 5 50 50 -1SENDA RECA .305515.8187.01210 200. 0 -2SENDB RECB .031991.5559.01937 200. 0 -3SENDC RECC C 34567890123456789012345678901234567890 92RECA 5555. -1. -1. { Old ZnO data formats involved a dummy 1. 1. { characteristic (these two data cards) 9999. 92RECB RECA 5555. { 2nd arrester of 3-phase group is copy of 1st 92RECC RECA 5555. { 3rd arrester of 3-phase group is copy of 1st BLANK card ending branch data BLANK card terminating nonexistent switches 14SENDA 408000. 60. 0.0 14SENDB 408000. 60. -120. 14SENDC 408000. 60. 120. BLANK card ending source cards C 345678901234567890123456789012345678901234567890123456789012345678901234567890 C The following 3 cards are ZnO characteristics, using the old formats: -1 2500. 26. 0.5 778000. -1 2500. 26. 0.5 778000. -1 2500. 26. 0.5 778000. 1 { 1st of remaining cards to be discarded searching for next case beginning PRINTER PLOT 144 3. 0.0 20. RECA BLANK card ending plot cards { Final card to be discarded b4 new case starts BEGIN NEW DATA CASE $PUNCH C 92RECA 5555. C 0.778000000000000E+06 -0.100000000000000E+03 0.000000000000000E+00 C 0.250000000000000E+04 0.260000000000000E+02 0.500000000000000E+00 C 9999 C 92RECB 5555. C 0.778000000000000E+06 -0.100000000000000E+03 0.000000000000000E+00 C 0.250000000000000E+04 0.260000000000000E+02 0.500000000000000E+00 C 9999 C 92RECC 5555. C 0.778000000000000E+06 -0.100000000000000E+03 0.000000000000000E+00 C 0.250000000000000E+04 0.260000000000000E+02 0.500000000000000E+00 C 9999 C 3rd of 12 subcases will illustrate use of another supporting program. C This came from Prof. Ned Mohan and graduate student Jim Frame at the C University of Minnesota in Minneapolis around 1981. HYSTERESIS C ITYPE LEVEL { Request Armco M4 oriented silicon steel -- only 1 available 1 2 { That was ITYPE=1. As for LEVEL=2, moderate accuracy output 500. 1.0 { Current and flux coordinates of positive saturation point $PUNCH C -1.25000000E+02 -9.70588235E-01 C -1.56250000E+01 -9.19411765E-01 C 9.37500000E+00 -8.47058824E-01 C 2.18750000E+01 -6.58823529E-01 C 4.21875000E+01 5.70588235E-01 C 6.56250000E+01 7.52941176E-01 C 1.12500000E+02 8.64705882E-01 C 2.07812500E+02 9.41176471E-01 C 5.00000000E+02 1.00000000E+00 C 6.87500000E+02 1.00588235E+00 C 9999. C A quarter of a century late, Orlando Hevia in Santa Fe, Argentina, enhances C both the code and the illustrations. He adds to the original ITYPE = 1 a C new ITYPE = 2 alternative, which is for Armco M4 oriented silicon steel. C Next, illustrate the remaining 3 levels for ITYPE = 1 (these are 1, 3, and C 4) followed by the 4 levels for ITYPE = 2. Orlando had $PUNCH requests C after each, although WSM removes them to minimize such redundant output. C In E-mail to BPA's Dr. Tsu-huei Liu on February 21, 2006, Orlando wrote: C "I received a curve of the B-H characteristics of ARMCO M-6 steel. It is C for 0.014-inch thickness (0.35 mm), with a bit more losses than ARMCO M-4." C ITYPE LEVEL { Request Armco M4 oriented silicon steel 1 1 { That was ITYPE=1 and LEVEL=2. 1st of 4 levels is LEVEL = 1 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M4 oriented silicon steel 1 3 { That was ITYPE=1 and LEVEL=2. 3rd of 4 levels is LEVEL = 3 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M4 oriented silicon steel 1 4 { That was ITYPE=1 and LEVEL=2. 4th of 4 levels is LEVEL = 4 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M6 oriented silicon steel 2 1 { 1st of 4 illustrations of M6 steel is for LEVEL = 1 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M6 oriented silicon steel 2 2 { 2nd of 4 illustrations of M6 steel is for LEVEL = 2 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M6 oriented silicon steel 2 3 { 3rd of 4 illustrations of M6 steel is for LEVEL = 3 500. 1.0 { Current and flux coordinates of positive saturation point C ITYPE LEVEL { Request Armco M6 oriented silicon steel 2 4 { 4th of 4 illustrations of M6 steel is for LEVEL = 4 500. 1.0 { Current and flux coordinates of positive saturation point BLANK card ending stacked "HYSTERESIS" data subcases BEGIN NEW DATA CASE C 4th of 12 subcases will illustrate the printing of error messages. C Were there not at least five data cards remaining at this point, the C data subcase would not even be read. The program would instead C discard remaining cards and jump out to issue the initial prompt C once again. But the present set of comment cards prevent this. KILL CODES, 13, 14, { Just look at two very short error messages: KILL = 13, 14 C The following 3 lines will never be executed because of the nature of the C preceding KILL CODES. However, they are legal non-comment cards. Once C we switch to the destruction of comment cards using NOCOMM= 1 in STARTUP, C the 5th subcase was missing because it was not separated from the 4th by at C least 5 data cards (variable KASEND = 5 in STARTUP). So, 14 Dec 94, we C add 3 noncomments: PRINTED NUMBER WIDTH, 15, 2, { 1st of 3 never-processed noncomment lines PRINTED NUMBER WIDTH, 15, 2, { 2nd of 3 never-processed noncomment lines PRINTED NUMBER WIDTH, 15, 2, { 3rd of 3 never-processed noncomment lines BEGIN NEW DATA CASE C 5th of 12 subcases is like the 4th except that output will be C truncated to 125 columns --- short of the 132 columns needed. TRUNCATE OUTPUT LINES, 125, { 132-column lines are formed, but then chopped off KILL CODES, 13, 14, { The same two error messages as preceding subcase, note C The following 3 lines will never be executed because of the nature of the C preceding KILL CODES. However, they are legal non-comment cards. Once C we switch to the destruction of comment cards using NOCOMM= 1 in STARTUP, C this 5th subcase was missing because it did not have 5 or more data cards C (variable KASEND = 5 in STARTUP). So, 14 Dec 94, we add 2 noncomments: PRINTED NUMBER WIDTH, 15, 2, { 1st of 3 never-processed noncomment lines PRINTED NUMBER WIDTH, 15, 2, { 2nd of 3 never-processed noncomment lines PRINTED NUMBER WIDTH, 15, 2, { 3rd of 3 never-processed noncomment lines BEGIN NEW DATA CASE C 6th of 12 subcases illustrates batch-mode usage of Orlando C Hevia's fitter to approximate piecewise-line (I, Psi) points C of an ordinary, old, Type-93 nonlinear reactor by a smooth, C hyperbolic tangent. See story in January, 1998, newsletter. TRUNCATE OUTPUT LINES, 132, { Cancel preceding truncation; return to normal 132 $ERASE C 345678901234567890123456 C 345678901234567890123456789012 ITMAX KPL DXL2 SMOOTH SATURATION USING TANH 600 40 .08 C Preceding parameters to the right of col. 32 are optional changes to defaults: C ITMAX = the maximum number of iterations of parameter fitter (default = 3000) C KPL = the number of uniform current steps for curve plotting (default = 50) C DXL2 = radius in inches of marking circle to show data points (default = .05) Arbitrary line will not be seen (illustrate how program ignores all but Type 93) Arbitrary line will not be seen (illustrate how program ignores all but Type 93) C AAAAAABBBBBBCCCCCCDDDDDDCURR FLUX I 93BUS1 BUS2 0.005 30.0 1 -.005 30.0 { Illegal point illustrates rejection 0.0 0.0 { 2nd illegal point (special case of origin) .005 30.0 .010 40.0 .020 45.0 .030 47.0 9999. Arbitrary line will not be seen (illustrate how program ignores all but Type 93) Arbitrary line will not be seen (illustrate how program ignores all but Type 93) 93BUS1 BUS2 0.005 30.0 .002 11.0 .005 30.0 .010 40.0 .020 43.0 .030 45.0 .040 46.0 9999. Arbitrary line will not be seen (illustrate how program ignores all but Type 93) Arbitrary line will not be seen (illustrate how program ignores all but Type 93) 93BUS1 BUS2 0.005 30.0 .002 11.0 .005 30.0 9999. Arbitrary line will not be seen (illustrate how program ignores all but Type 93) Arbitrary line will not be seen (illustrate how program ignores all but Type 93) $PUNCH { Create .PCH file as well as show contents in the printed output BLANK card ending last SSUT data case BEGIN NEW DATA CASE C 7th of 12 subcases illustrates batch-mode usage of Orlando C Hevia's supporting program HYSTER. This serves to convert C a Type-98 element into a Type-96 element. For background, C see the story in the July, 2000, newsletter. C 345678901234567890123456789012345678901234567890123456 C NOZOOM FSCALE DXL2 HYSTERESIS HEVIA 1 1.1 .04 C The 3 parameters on the preceding request card govern the screen plot C that automatically will document hysteresis of the result: C NOZOOM is a binary switch that indicates whether the entire C curve or just the hysteresis loop itself is to be plotted. C Value 1 will plot everything (no zoom) whereas value zero C will drop the final point, which is outside the loop proper. C FSCALE is a scaling factor to waste space above and to the right C of the 1st and 2nd-quadrant plot. For example, value 1.1 C will waste 10%, providing this much margin. C DXL2 is the "radius" in screen inches of the squares that mark data C points. Value zero means that there will be no such marking. C So much for the special-request card. On to reactor. Everything about C the following is normal Type 98 except for columns 27-74, from which 3 C new floating-point parameters are read. Of course, if currents are in C amperes and voltage is in volts, power will be in watts and residual C flux will be in volt-seconds: C <--Losses in W-><--Freq in Hz--><-Residual flux> I 98NODE1 NODE2 40774.59 50.0 400.0 3 0.122419130E-01 0.563980835E+03 0.628200099E-01 0.782295959E+03 0.215062916E+00 0.100061115E+04 0.576934814E+00 0.121892639E+04 0.791255713E+00 0.129949500E+04 0.131488657E+01 0.143724146E+04 0.266661191E+01 0.165555664E+04 0.386263156E+01 0.178290710E+04 0.550489426E+01 0.192845056E+04 9999 $PUNCH, dc13c.pch ! { Exclamation holds lower case C NOZOOM FSCALE DXL2 REPLOT 0 1.0 .07 C The following subcase illustrates generalization of the preceding hysteresis C derivation by Orlando Hevia. WSM adds 26 May 2001. C The following case produces Type-96 hysteretic branch cards that will result C in automatic initialization. There will be 8888. for I-steady and blank C fields for flux and remanent flux. Note negative losses on the following C Type-98 card (this is the way the user requests this option): C <--Losses in W-><--Freq in Hz--><-Residual flux> I 98NODE1 NODE2 -40774.59 50.0 400.0 3 0.122419130E-01 0.563980835E+03 0.628200099E-01 0.782295959E+03 0.215062916E+00 0.100061115E+04 0.576934814E+00 0.121892639E+04 0.791255713E+00 0.129949500E+04 0.131488657E+01 0.143724146E+04 0.266661191E+01 0.165555664E+04 0.386263156E+01 0.178290710E+04 0.550489426E+01 0.192845056E+04 9999 $PUNCH C NOZOOM FSCALE DXL2 REPLOT 0 1.0 .07 BLANK card terminates stacked Type-98 elements within HYSTERESIS HEVIA BEGIN NEW DATA CASE C 8th of 12 subcases illustrates in-line usage of Orlando Hevia's C hysteresis-creating program. The effect should be identical to C replacement of the Type-98 branch by the Type-96 branch that was C punched by the preceding subcase. The basic network was created C by taking one phase of DC-7, and replacing the reactor using data C of the preceding subcase. Simulation continues only long enough C for the first pulse of current inrush of the reactor. .000025 .001 60. 60. 1 1 1 0 0 -1 5 5 30 1 GENA SWA 89352. -1SENDA RECA .137681.07755.6806124.27 RECA TRANA GENA SWA 1 TRANA GNDA 7.5398 1 GNDA .0001 C Note special request in 15:26 of the following branch card. This is C what provides the in-line connection to Orlando Hevia's fitter, C replacing the Type-98 branch by an equivalent (better) Type-96 branch: C <--Losses in W-><--Freq in Hz--><-Residual flux> I 98TRANA HEVIA HYSTER 40774.59 50.0 400.0 3 0.122419130E-01 0.563980835E+03 0.628200099E-01 0.782295959E+03 0.215062916E+00 0.100061115E+04 0.576934814E+00 0.121892639E+04 0.791255713E+00 0.129949500E+04 0.131488657E+01 0.143724146E+04 0.266661191E+01 0.165555664E+04 0.386263156E+01 0.178290710E+04 0.550489426E+01 0.192845056E+04 9999 BLANK card ends all branch cards SWA SENDA 1.0 BLANK card ends all switch cards 14GENA 188000. 60. 0.0 -1. BLANK card ending source cards GENA TRANA RECA C Step Time TRANA GENA TRANA RECA TRANA RECA TRANA C TERRA TERRA TRANA GNDA C *** Switch "SWA " to "SENDA " closed after 0.00000000E+00 sec. C 0 0.0 0.0 188000. 0.0 0.0 0.0 0.0 0.0 C 1 .25E-4 .52901E-14 187991.65 .52901E-14 .52894E-14 -.2049E-16 -.1203E-16 .84641E-17 C 2 .5E-4 .96607E-14 187966.602 .96607E-14 .96583E-14 -.2049E-16 -.2196E-16 -.1471E-17 BLANK card ending output variable requests (just node voltages, here) C 32 .8E-3 .90572E-14 179514.534 .90572E-14 .89899E-14 -.2044E-16 -.2044E-16 .11666E-20 C 33 .825E-3 81947.1269 178980.254 81947.1269 81954.0422 .007179619 131.122181 131.114993 C 34 .85E-3 284030.522 178430.075 284030.522 284054.493 .039243952 192.256743 192.217471 C 35 .875E-3 383163.494 177864.047 383163.494 383195.837 .097698701 -33.507413 -33.60515 C 36 .9E-3 365690.044 177282.221 365690.044 365720.926 .16330787 5.81105035 5.64770591 C 37 .925E-3 368481.981 176684.647 368481.981 368513.116 .153316524 -1.0272643 -1.1806177 C 38 .95E-3 367719.904 176071.379 367719.904 367750.991 .18045259 .141786623 -.03870274 C 39 .975E-3 367528.209 175442.471 367528.209 367559.301 .2075535 -.06041831 -.26800857 C 40 1.E-3 367198.987 174797.979 367198.987 367230.074 .234635209 -.02407335 -.25874528 BLANK card terminating plot cards BEGIN NEW DATA CASE C 9th of 12 subcases illustrates batch-mode use of Orlando Hevia's C supporting program to find parameters A, B, and C of the Type-15 C Standler surge function. See the January, 2000, newsletter. The C alternative of in-line fitting rather than batch-mode fitting can C be found in the 3rd subcase of DC-19. STANDLER SURGE FUNCTION 10.0 50.0 0.0 1000. S10500 -1 $PUNCH, dc13d.pch ! { Exclamation holds lower case C C STANDLER SURGE FUNCTION C C 10.0 50.0 0.0 1000. S10500 -1 C C Surge peak : 1.0000E+03 C C Time to peak : 1.0000E+01 us C C Time to half value : 5.0000E+01 us C 15S10500-1 1913.48 .34489E-4 .28995155 Standler 10.0 50.0 10.0 1000. S10501 0 $PUNCH C C 10.0 50.0 10.0 1000. S10501 0 C C Surge peak : 1.0000E+03 C C Virtual front time : 1.0000E+01 us C C Virtual half time : 5.0000E+01 us C C Virtual front defined by 10.0 and 90.0 % C 15S10501 2543.7879 .23168E-4 .65859353 Standler 10.0 50.0 30.0 1000. S10503 1 10.0 50.0 50.0 1000. S10505 Voltage 8.0 20.0 0.0 1000. S82000 Current 8.0 20.0 10.0 1000. S82010 voltage 8.0 20.0 30.0 1000. S82030 current 8.0 20.0 50.0 1000. S82050 VOLTAGE $PUNCH C C 10.0 50.0 30.0 1000. S10503 1 C C Surge peak : 1.0000E+03 C C Virtual front time : 1.0000E+01 us C C Virtual half time : 5.0000E+01 us C C Virtual front defined by 30.0 and 90.0 % C 15S10503 2372.9809 .25096E-4 .52633653 Standler C This is 1st of 6. Etc. for 2nd through 6th illustrations. BLANK card ending Standler surges BEGIN NEW DATA CASE C 10th of 12 subcases illustrates use of LOSSY SATURATION from C Orlando Hevia. This is added 11 May 2001 as feature becomes C available. In general terms, this supporting program is similar C to SATURATION (see 1st subcase), but losses are taken into C account. LOSSY SATURATION $ERASE C C FREQ: AS SATURATION C C VBASE: AS SATURATION C C PBASE: AS SATURATION C C KCON : CONNECTION C -1, BLANK OR 0: SINGLE PHASE (THE SAME AS SATURATION) C VOLTAGE, FLUX AND CURRENT WITHOUT SCALING. C C 1 STAR (Y, WYE) CONNECTION NEUTRAL GROUNDED C VOLTAGE AND FLUX SCALED DIVIDING BY SQRT(3) C C 3 DELTA (TRIANGLE) CONNECTION C CURRENT SCALED DIVIDING BY SQRT(3) C TRIPLEN HARMONICS REMOVED C C KTHIRD : AS SATURATION C C IOUT : TO BE PASSED TO .PCH FILE IN COLUMN 80 C C BUS1, BUS2 : TO BE PASSED TO .PCH FILE FOR NON-LINEAR RESISTOR C BUS3, BUS4 : TO BE PASSED TO .PCH FILE FOR NON-LINEAR INDUCTOR C IF BUS3=BUS4=' ', BUS1 AND BUS2 WILL BE USED C C FREQ VBASE PBASE IPUNCH KTHIRD KCON IOUT BUS1 BUS2 BUS3 BUS4 C ------________--------________--------________--------______------______------ C -1. 10. .001 0 0 C I V C --------------________________ 0.0 5.0 2.0 5.0 3.0 3.5 4.0 2.0 5.0 1.0 10.0 1.0 9999 -1. 10. .001 0 1 0.0 5.0 2.0 5.0 3.0 3.5 4.0 2.0 5.0 1.0 10.0 1.0 9999 -1. 10. .001 1 0 0.0 5.0 2.0 5.0 3.0 3.5 4.0 2.0 5.0 1.0 10.0 1.0 9999 -1. 10. .001 1 1 0.0 5.0 2.0 5.0 3.0 3.5 4.0 2.0 5.0 1.0 10.0 1.0 9999 C 60. .664 .01 0 0 1.0 100. 1.6 105. 2.75 110. 4.0 113. 6.0 116. 10. 119. 44. 131. 9999 C $PUNCH, dc13e.pch ! { Exclamation holds lower case -1. 10. .001 1 2.0 5.0 9999 60. 303.11 300. 1 1 .000306 0.7 .000415 0.8 .00076 0.9 .00198 1.0 .0108 1.1 9999 C $PUNCH 60. 34.5 300.0 1 .000306 0.7 .000415 0.8 .00076 0.9 .00198 1.0 .0108 1.1 9999 C $PUNCH, dc13f.pch ! { Exclamation holds lower case C INI EPSIR CONTROLS 1024 1.E-8 { # of points per quarter cycle; convergence tolerance C C THE FOLLOWING IS A REAL LIFE CASE C C C FREQ VBASE PBASE IPUNCH KTHIRD KCON IOUT BUS1 BUS2 BUS3 BUS4 C ------________--------________--------________--------______------______------ 50. 0.001 1.0E-6 0 3 1DOMAR DOMBR DOMAL DOMBL C C IRMS----------VRMS------------P--------------- 3.01 29107.00 50640.0 4.39 31000.00 60640.0 7.51 33200.00 75760.0 14.45 35387.0 94960.0 26.09 37280.00 114000.0 53.73 39213.00 135360.0 9999. $PUNCH C C SAMPLES OF EACH CONNECTION C CURRENT, VOLTAGE AND POWER TO OBTAIN RESISTANCE=REACTANCE C C CURRENT, VOLTAGE AND LOSSES SCALED ACCORDINGLY C C SINGLE PHASE C C FREQ VBASE PBASE IPUNCH KTHIRD KCON IOUT BUS1 BUS2 BUS3 BUS4 C ------________--------________--------________--------______------______------ 50. 0.001 1.0E-6 0 -1 1 C C IRMS----------VRMS------------P--------------- 1.414213562 10000.0 10000.000 3.39411255 12000.0 28800.000 7.071067812 14000.0 70000.000 9999. C C Y NEUTRAL CONNECTED TO GROUND C C FREQ VBASE PBASE IPUNCH KTHIRD KCON IOUT BUS1 BUS2 BUS3 BUS4 C ------________--------________--------_______------______------______------ 50. 0.001 1.0E-6 0 0 1 1 C C IRMS----------VRMS------------P--------------- 1.414213562 17320.50808 30000.000 3.39411255 20784.60969 86400.000 7.071067812 24248.71131 210000.000 9999. C C DELTA CONNECTED C C FREQ VBASE PBASE IPUNCH KTHIRD KCON IOUT BUS1 BUS2 BUS3 BUS4 C ------________--------________--------________--------______------______------ 50. 0.001 1.0E-6 0 0 3 C C IRMS----------VRMS------------P--------------- 2.449489743 10000.0 30000.000 5.878775383 12000.0 86400.000 12.24744871 14000.0 210000.000 9999. $PUNCH BLANK card ends stacked subcases of LOSSY SATURATION data BEGIN NEW DATA CASE C 11th of 12 subcases illustrates batch-mode use of Orlando Hevia's C supporting program to find parameters A, B, and C of the Type-15 C Heidler surge function. See the January, 2001, newsletter. The C alternative of in-line fitting rather than batch-mode fitting can C be found in the ?? subcase of DC-??. Add this subcase 16 Nov 02 HEIDLER SURGE FUNCTION C FRONT TAIL PERC AMPLIT BUS IV N C ------========--------========XX------========-- 1.0E-6 50.0E-6 10.0 100.0 BUS 1 5 1.0E-6 50.0E-6 0.0 100.0 BUS -1 5 1.0E-6 50.0E-6 30.0 100.0 BUS Voltage 5 1.0E-6 50.0E-6 50.0 100.0 BUS current 5 C The preceding data was within HEIDLER.DAT as supplied by Orlando Hevia $PUNCH C Essence of Orlando Hevia's HEIDSLOP.DAT is inserted here 17 August 2003. C This verifies a new fitting option that allows the user to specify his C desired peak value and slope. C NEW OPTION FOR HEIDLER SURGE SOURCE C IF PERC IS .LT. 0.0, THE SOURCE WILL HAVE THE C SLOPE AS REQUIRED IN FIRST FIELD (NOT T1) C SLOPE> IV N 1.0E+10 50.0E-6 -1.0 1.0E4 OFFLIN -1 4 $PUNCH BLANK card ending input for Heidler fitter BEGIN NEW DATA CASE C 12th of 12 subcases illustrates batch-mode use of Orlando Hevia's C supporting program TWOEXP to find parameters of Dommel's Type-15 C exponential surge funct. See the January, 2001, newsletter (end C of the Standler story mentions "other alternatives"). This does C not imply more complicated (even Dommel's original 2-exponential C surge is benefitted). Final note: the alternative of in-line C fitting rather than batch-mode fitting can C be found in the ?? subcase of DC-??. Add this subcase 16 Nov 02 TWO EXP SURGE FUNCTION C FRONT TAIL PERC AMPLIT BUS IV N C ------========--------========XX------========-- 1.0E-6 50.0E-6 10.0 100.0 BUS 1 5 1.0E-6 50.0E-6 0.0 100.0 BUS -1 5 1.0E-6 50.0E-6 30.0 100.0 BUS Voltage 5 1.0E-6 50.0E-6 50.0 100.0 BUS current 5 $PUNCH BLANK card ending input for 2-exponential fitter BEGIN NEW DATA CASE BLANK