BEGIN NEW DATA CASE C BENCHMARK DCNEW-27 C 1st of 10 subcases that verify correct use of manually-specified initial C conditions for problems involving various distributed line models. Prior C to beginning on 22 November 1998, there was no such illustration. Manual C specification of initial conditions can, however, be found in many test C cases (e.g., DC-34 that involves Prof. Hian Lauw's U.M.). But none used C a distributed line. Begin with 1-phase, constant-parameter distributed. PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) .000050 .002 1 1 1 1 C Create two parallel, independent, identical circuits of line plus resistor. C The (SEND, REC) alternative will involve no phasor solution whereas the C alternative (GEN, END) will have a phasor solution. Answers agree: -1SEND REC .306 5.82 .012 100. { Half the length of original -1GEN END .306 5.82 .012 100. { Half the length of original C 19 July 2009, we add branch name "RESIST" for use with initial conditions. C 34567890123 Name RESIST --- Equivalent, alternative request is left-adjusted REC Name:RESIST 100. { Name branch using right-adjusted " Name:" } 1 END 100. 1 BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. { Excitation of (SEND, REC) is transient only 14GEN 100. 60. { Note phasor presence } -1. BLANK card ending source data C (SEND, REC) had no phasor solution, initialize it using values of (GEN, END): C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 60. 2REC 19.956813332748-34.69967901024 C Distributed: I-k real I-k imag I-m real I-m imag 3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239 C 19 July 2009, convert the following linear branch current from the original C use of terminal node names (REC, TERRA) to the use of branch name RESIST. C Just today, such alternative use of branch names is being added to ATP. WSM C 3REC .19956813332748 { Unnecessary, since purely resistive branch 3Name: RESIST.19956813332748 { Unnecessary, since purely resistive branch SEND GEN REC END { List of nodes for node-voltage output C Travel time is less than 17 steps, so document the solution this far: C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Next 2 output variables are branch currents (flowing from the upper node to the lower node); C Step Time SEND GEN REC END REC END C TERRA TERRA C 0 0.0 100. 100. 19.95681 19.95681 .1995681 .1995681 C 1 .5E-4 99.98224 99.98224 20.60655 20.60655 .2060655 .2060655 C 2 .1E-3 99.92895 99.92895 21.2497 21.2497 .212497 .212497 C 3 .15E-3 99.84016 99.84016 21.88529 21.88529 .2188529 .2188529 C 4 .2E-3 99.71589 99.71589 22.51311 22.51311 .2251311 .2251311 C 5 .25E-3 99.5562 99.5562 23.13293 23.13293 .2313293 .2313293 C 6 .3E-3 99.36113 99.36113 23.74453 23.74453 .2374453 .2374453 C 7 .35E-3 99.13076 99.13076 24.3477 24.3477 .243477 .243477 C 8 .4E-3 98.86517 98.86517 24.94221 24.94221 .2494221 .2494221 C 9 .45E-3 98.56446 98.56446 25.52786 25.52786 .2552786 .2552786 C 10 .5E-3 98.22873 98.22873 26.10444 26.10444 .2610444 .2610444 C 11 .55E-3 97.85809 97.85809 26.67175 26.67175 .2667175 .2667175 C 12 .6E-3 97.45269 97.45269 27.22958 27.22958 .2722958 .2722958 C 13 .65E-3 97.01266 97.01266 27.77774 27.77774 .2777774 .2777774 C 14 .7E-3 96.53816 96.53816 28.31603 28.31603 .2831603 .2831603 C 15 .75E-3 96.02937 96.02937 28.84425 28.84425 .2884425 .2884425 C 16 .8E-3 95.48645 95.48645 29.36223 29.36223 .2936223 .2936223 C 17 .85E-3 94.90961 94.90961 29.86972 29.86972 .2986972 .2986972 BLANK card ending node voltage outputs BLANK card ending plot cards BEGIN NEW DATA CASE C 2nd of 10 subcases is like 1st, but more general. The line is 3-phase, and C excitation and connection are unbalanced in order that all line modes have C strong phasor history. The line is half the length originaly used in DC-38 C and nonlinear elements at the far end have been replaced by linear, 200-ohm C resistors. Phase c of excitation is omitted, with the voltage source is C replaced by a 4th 200-ohm resistor to ground. All voltages and currents C are substantial. PRINTED NUMBER WIDTH, 11, 2, { Reduced precision POWER FREQUENCY, 60., { Definition here allows omision on initial condition card .000050 .00100 { TMAX = 1 msec slightly exceeds the travel time of the line 1 -1 1 1 -1SENDA RECA .305515.8187.01210 100. 0 { 100-mile, constant- -2SENDB RECB .031991.5559.01937 100. 0 { parameter, 3-phase -3SENDC RECC { transmission line. RECA 200. RECB 200. RECC 200. SENDC 200. { Replace phase c source by this resistor BLANK card follows the last branch card BLANK line terminates the last (here, nonexistent) switch 14SENDA 100. 60. 0.0 14SENDB 100. 60. -120. BLANK card ending source data C The following initial-condition cards were constructed from branch flow C output with sources present for a phasor solution. Full precision of C the ATP output is being used, note. Purely resistive branches are not C supposed to require initialization, so these have been omitted. C 2SENDA 100. 0.0 60. C Note omission of frequency on right of preceding card. Since equal to C the power frequency, its appearance on the 1st node voltage initial C condition card is not required. The power frequency will be used. C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SENDA 100.0 0.0 2RECA 77.881089324669-16.76951438169 2SENDB -50.-86.60254037844 2RECB -84.53015167415-56.11022779303 2SENDC 8.9524884789137-4.001192729873 2RECC -7.4834609042614.3043134435609 C Distributed: I-k real I-k imag I-m real I-m imag 3SENDA RECA .38833508724138-.020723270343 -.3894054466233.08384757190847 3SENDB RECB -.3776059443784-.3319155592773.42265075837076.28055113896513 3SENDC RECC -.0447624423946.02000596364937 .0374173045213-.0215215672178 C Step Time RECC RECB RECA SENDC SENDA SENDB C 0 0.0 -7.483461 -84.53015 77.881089 8.9524885 100. -50. C 1 .5E-4 -7.563122 -83.45714 78.183167 9.0261691 99.982235 -48.35879 C 2 .1E-3 -7.640246 -82.35489 78.457621 9.0967927 99.928947 -46.70041 C 3 .15E-3 -7.714656 -81.22339 78.704198 9.1641842 99.840155 -45.02543 C 4 .2E-3 -7.786324 -80.06303 78.922812 9.2283197 99.71589 -43.33445 C 5 .25E-3 -7.855227 -78.87422 79.113386 9.2891764 99.556196 -41.62808 C 6 .3E-3 -7.921338 -77.65739 79.27585 9.3467328 99.361131 -39.90692 C 7 .35E-3 -7.984635 -76.41297 79.410149 9.4009683 99.130763 -38.17157 C 8 .4E-3 -8.045095 -75.14139 79.516233 9.4518636 98.865174 -36.42267 C 9 .45E-3 -8.102697 -73.84312 79.594066 9.4994008 98.56446 -34.66082 C 10 .5E-3 -8.15742 -72.51862 79.643619 9.5435629 98.228725 -32.88666 C 11 .55E-3 -8.209246 -71.16835 79.664875 9.5843341 97.85809 -31.10082 C 12 .6E-3 -8.258233 -69.7927 79.657818 9.6216977 97.452687 -29.30393 C 13 .65E-3 -8.304204 -68.39235 79.622466 9.6556448 97.01266 -27.49662 C 14 .7E-3 -8.347225 -66.9677 79.558825 9.6861614 96.538164 -25.67954 C 15 .75E-3 -8.387279 -65.51925 79.466917 9.7132365 96.029369 -23.85335 C 16 .8E-3 -8.424354 -64.04753 79.346775 9.7368605 95.486454 -22.01867 C 17 .85E-3 -8.458403 -62.55302 79.198474 9.7570132 94.909614 -20.17617 C 18 .9E-3 -8.489362 -61.03619 79.022118 9.7736699 94.299054 -18.32651 C 19 .95E-3 -8.517428 -59.49781 78.817564 9.786901 93.654989 -16.47033 C 20 .1E-2 -8.542467 -57.93828 78.585006 9.7966549 92.977649 -14.6083 C Note about preceding solution. This not only is smooth, it was verified C by adding the phasor solution and omitting initial conditions. Agreement C is exact. 1 { Request node voltage output of all nodes BLANK card ending plot cards BEGIN NEW DATA CASE C 3rd of 10 subcases is like 1st, but without the built-in phasor solution, C and with the constant-parameter line replaced by JMarti frequency-dependent C line modeling. Line length has changed, unfortunately. Rather than build C data to derive Marti branch cards, we borrowed those already punched by C the 2nd subcase of DCNEW-5. Line length is 108.7 km (7 time steps) C Reduce 2 in col. 31 to 0 so all Marti data card input is seen. PRINTED NUMBER WIDTH, 10, 2, .000050 .001 1 1 1 1 C <++++++> Cards punched by support routine on 11-Nov-18 11.00.00 <++++++> C ***** TRANSPOSED JMARTI line segment ****** C JMARTI SETUP C BRANCH JDA LMA JDB LMB JDC LMC C LINE CONSTANTS C METRIC C 1.3161 .1151 4 2.355 -6.5 22. 8.0 C 0.5 4.189 4 0.914 4.5 27.3 15.7 C BLANK card ending conductor cards of imbedded "LINE CONSTANTS" data C 1500. 60. 1 108.7 C 1500. .01 1 108.7 9 10 C BLANK card ending frequency cards of inbedded "LINE CONSTANTS" data C BLANK card ending "LINE CONSTANTS" cases (back to "JMARTI SETUP") C 1 { Note use of default conductance G to ground } C .30 30 0 1 3 0 0 -1SEND REC 0. 0.00 -2 1 25 4.6006674814379022100E+02 -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 3.72990420164290220E+07 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 2.89354778735701405E+07 14 3.8038812947515656600E-04 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 -7.93497908986478425E+10 3.97029750605793992E+10 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 1.82589462762646523E+05 1.82772052225409076E+05 0.58951119 0.00000000 REC 100. 1 BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. BLANK card ending source data C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678 C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865 C C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945 C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15 C The following initial condition cards are built from preceding phasor flow. C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 60. 2REC 14.563049373597-25.89164864838 C Distributed: I-k real I-k imag I-m real I-m imag 3SEND REC .14695409358625-.2529237293961-.145630493736 .25891648648378 C First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Next 1 output variables are branch currents (flowing from the upper node to the lower node); C Step Time REC SEND REC C TERRA C 0 0.0 14.56305 100. .1456305 C 1 .5E-4 15.04761 99.98224 .1504761 C 2 .1E-3 15.52757 99.92895 .1552757 C 3 .15E-3 16.00208 99.84016 .1600208 C 4 .2E-3 16.47093 99.71589 .1647093 C 5 .25E-3 16.93392 99.5562 .1693392 C 6 .3E-3 17.3909 99.36113 .173909 C 7 .35E-3 17.8417 99.13076 .178417 C 8 .4E-3 18.28609 98.86517 .1828609 C 9 .45E-3 18.72385 98.56446 .1872385 C 10 .5E-3 19.155 98.22873 .19155 C 11 .55E-3 19.57944 97.85809 .1957944 C 12 .6E-3 19.99694 97.45269 .1999694 C 13 .65E-3 20.40735 97.01266 .2040735 C 14 .7E-3 20.8105 96.53816 .208105 C 15 .75E-3 21.20625 96.02937 .2120625 C 16 .8E-3 21.59436 95.48645 .2159436 C Showing 16 steps is more than enough. Without manual initial conditions, the C surge arrives at the receiving end on step 8. Above exceeds 2 travel times. C Furthermore, above agrees exactly to output that includes a phasor solution. 1 { Request output of all node voltages BLANK card ending plot cards BEGIN NEW DATA CASE C 4th of 10 subcases is like 2nd, but with the constant-parameter line C replaced by JMarti frequency-dependent line modeling. Line length has C changed, unfortunately. Rather than build data to derive Marti branch C cards, we borrowed those used by DCNEW-4. Line length is 138 miles. C Reduce 2 in col. 31 to 1 so some Marti data card input is seen. PRINTED NUMBER WIDTH, 10, 2, { Reduced precision .000050 .001 { Note TMAX exceeds travel time of .74 msec by about 40% 1 1 1 1 { Request phasor output to prove there is none C PUNCHED CARD OUTPUT OF "JMARTI SETUP" WHICH BEGAN AT 14.32.51 09/14/82 C 1.3636 .05215 4 1.602 -20.75 50. 50. C 1.3636 .05215 4 1.602 -19.25 50. 50. C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5 C 2.3636 .05215 4 1.602 0.75 77.5 77.5 C 3.3636 .05215 4 1.602 19.25 50. 50. C 3.3636 .05215 4 1.602 20.75 50. 50. C 0.5 2.61 4 0.386 -12.9 98.5 98.5 C 0.5 2.61 4 0.386 12.9 98.5 98.5 C C 100. 5000. 1 138. 1 3-2 C 100. 60.00 1 138. 1 3 C 100. .01 1 138. 1 9 10 3 C -1SENDA RECA 1. -2 3 1 18 0.47487144790587633736E+03 -0.117967342612603676E+01 -0.316671345021416945E+01 0.255662317832426673E+02 -0.601654631032706213E+02 -0.118825306534453652E+03 0.128169142098635430E+04 0.589610442538186533E+04 0.292891488173549769E+05 0.115665757459954655E+06 0.456394278488149546E+06 0.153681652821077613E+07 0.681848054466802243E+07 0.155477562610202006E+08 0.231789897204889809E+08 0.143827766961221036E+08 0.241524708800130989E+08 0.245049034784953832E+08 0.476808364145524129E+08 0.321645453133055891E+00 0.897616891499021768E+00 0.165632367032913183E+01 0.166669419493737986E+01 0.253073682990780913E+01 0.200014655383098283E+02 0.124257951002701473E+03 0.657345842950598183E+03 0.278623017775457413E+04 0.117145310573639467E+05 0.427410884649070877E+05 0.202553264108715757E+06 0.968789502295299084E+06 0.296908522006270138E+07 0.737968316276072920E+07 0.121009167337626640E+08 0.132891445947112725E+08 0.252097199908134253E+08 14 0.86849922355034181931E-03 0.580543547692016874E-01 0.269744993210885436E+00 0.358646519088149721E+00 0.775351125410302933E+00 0.512138885239557462E+01 0.838265420474258915E+01 0.700229782114550385E+02 0.393084113996385128E+03 0.138208543918962840E+04 0.706160326497018531E+03 0.527246586822120662E+04 0.528343251309079615E+05 0.178307801042203880E+08 -0.178914532139174500E+08 0.228360173022402533E+02 0.102750699821169636E+03 0.139222638455845377E+03 0.269608322018776320E+03 0.389220463867018722E+03 0.427793721751318500E+03 0.132148883486994549E+04 0.277178340638843952E+04 0.449341296469570784E+04 0.596603161966094024E+04 0.114714152641659587E+05 0.268314807408425590E+05 0.204027813729699328E+05 0.204231841543429041E+05 -2SENDB RECB 1. -2 3 1 13 0.28580875773089406522E+03 0.381761604345214016E+04 -0.269529222495194318E+04 0.445018224900628539E+03 0.214658130219351907E+03 0.897325461203551775E+02 0.130480129629523390E+03 0.683037485206990933E+02 0.746626715384045969E+02 0.133938729387751110E+04 0.249954495902361043E+04 0.187764116219525090E+05 0.391031620802855759E+06 0.562668431190421223E+07 0.358668999787749571E+01 0.376679186833120605E+01 0.706439034439901437E+01 0.109721634064781350E+02 0.140506334140286961E+02 0.234292596911152251E+02 0.407686998675834982E+02 0.797598814775078289E+02 0.129106739969348078E+04 0.242138030939135859E+04 0.182362646041542116E+05 0.379814243183306855E+06 0.548808385606378759E+07 13 0.74123751236187229426E-03 0.192189595814863363E+02 0.210970481418782463E+01 0.383406546226866283E+02 0.307999912886148377E+02 0.634013410317139090E+03 -0.253243195247046863E+03 0.210540381301789816E+05 0.755533947306502050E+05 0.229576557459444489E+07 0.607777584135737231E+13 -0.156620677810598129E+13 0.153143851607186420E+13 -0.604300997216751013E+13 0.325736126030671755E+04 0.366694913707514260E+03 0.675727697588602189E+04 0.516024976370228876E+04 0.249456113770459378E+05 0.115067772934717399E+06 0.105205466546493019E+06 0.240930799961836921E+06 0.633276538663731131E+06 0.975820481990592802E+06 0.976796302472583498E+06 0.974176721795814519E+06 0.975150898517610403E+06 -3SENDC RECC 1. -2 3 1 13 0.27248404992352924125E+03 0.928429573805075023E+03 0.205756155237142632E+03 0.538106630554037608E+03 0.105199818376817531E+03 0.821382476778671116E+02 0.133356875032521778E+03 0.654774478488469374E+02 0.612546121562179922E+02 0.905544916173387175E+02 0.973244316091135204E+03 0.267481871266201051E+04 0.138566411938941123E+05 0.646075558955914967E+06 0.297277477955690506E+01 0.413177069940672514E+01 0.739410817322770197E+01 0.103072536754267949E+02 0.135874091964585515E+02 0.234792992956042887E+02 0.384075708508238121E+02 0.692539345057236435E+02 0.997101922255402791E+02 0.941460937860312015E+03 0.259295974388396348E+04 0.134538093446956691E+05 0.628668738677182439E+06 20 0.74871165651714274120E-03 0.249456855794677017E-01 0.516637614002139312E+01 0.435227595525510835E+01 0.696293545843379624E+01 0.113102519337405665E+02 0.107077756312168120E+02 0.136526319661978881E+02 0.103057494642840792E+03 0.647491597081196815E+03 0.221213037293261823E+04 0.164010494681972204E+05 0.517776934139508076E+05 0.393233870054795443E+06 -0.101893253559274247E+05 -0.173413454816459129E+06 0.524218185885897590E+06 0.711044562955687729E+12 -0.559640645251973480E+12 0.582892127161358521E+12 -0.734296849907947113E+12 0.887741958857318814E+01 0.183177145564733675E+04 0.151655674045953921E+04 0.245910936284398730E+04 0.383273362703089043E+04 0.399285585783856465E+04 0.496339281485561617E+04 0.253464455063045348E+05 0.272156609769823021E+05 0.475215596432249440E+05 0.149791320991727040E+06 0.333884409922466803E+06 0.541583000670874084E+06 0.670984155731562481E+06 0.569769858667429973E+06 0.132622344670840449E+07 0.440534234233555838E+07 0.440974768467789434E+07 0.436787151891064062E+07 0.437223939042955148E+07 0.57154434 0.70710678 -0.41761362 0.00000000 0.00000000 0.00000000 0.58879039 0.00000000 0.80696823 0.00000000 0.00000000 0.00000000 0.57154434 -0.70710678 -0.41761362 0.00000000 0.00000000 0.00000000 RECA 200. RECB 200. RECC 200. SENDC 200. { Replace phase c source by this resistor BLANK card follows the last branch card BLANK line terminates the last (here, nonexistent) switch 14SENDA 100. 60. 0.0 14SENDB 100. 60. -120. BLANK card ending source data C The following phasor solution was obtained by setting T-start = -1.0 for C each of the two preceding sources. This provides basis for initialization. C SENDA 100. 100. .37357344098553 .37380400389422 18.678672049276 3.7819577671008 C 0.0 0.0 -.0131269767114 -2.0124852 .65634883557087 0.6563488 C C RECA 74.582077276654 77.192523685071 -.3729103863833 .38596261842535 -14.89671428218 C -19.90475978176 -14.9430397 .09952379890879 165.0569603 -.451244944E-14 C C SENDB -50. 100. -.3833667890536 .47733933630699 21.899004304761 -1.819030186807 C -86.60254037844 -120.0000000 -.2843989223551 -143.4302949 9.4902958555069 9.4902959 C C RECB -87.9019678687 97.402329523617 .43950983934349 .48701164761808 -23.71803449157 C -41.95781025551 -154.4836884 .20978905127755 25.5163116 -.213197515E-14 C C SENDC 8.9664534487851 11.539874097595 -.0448322672439 .05769937048798 -.3329217354709 -.5808639378985 C -7.264393074381 -39.0135136 .03632196537191 140.9864864 .2136962482E-15 0.0000000 C C RECC -7.027068223793 9.9587590075801 .03513534111896 .0497937950379 -.2479422024276 C 7.0567126304834 134.8794005 -.0352835631524 -45.1205995 .3134699531E-16 C The following initial-condition cards were constructed from branch flow C output with sources present for a phasor solution. Full precision of C the ATP output is being used, note. Purely resistive branches are not C supposed to require initialization, so these were omitted. 2SENDA 100. 0.0 60. 2RECA 74.582077276654-19.90475978176 2SENDB -50.-86.60254037844 2RECB -87.9019678687-41.95781025551 2SENDC 8.9664534487851-7.264393074381 2RECC -7.0270682237937.0567126304834 C Distributed: I-k real I-k imag I-m real I-m imag 3SENDA RECA .37357344098553-.0131269767114-.3729103863833.09952379890879 3SENDB RECB -.3833667890536-.2843989223551.43950983934349.20978905127755 3SENDC RECC -.0448322672439.03632196537191.03513534111896-.0352835631524 1 { Request node voltage outputs for each network node C First 6 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 Step Time RECC RECB RECA SENDC SENDA SENDB SENDA SENDB SENDC C RECA RECB RECC C 0 0.0 -7.02707 -87.902 74.58208 8.966453 100. -50. .3735734 -.383367 -.044832 C 1 .5E-4 -7.15772 -87.0944 74.94189 9.101568 99.98224 -48.3588 .3737522 -.377936 -.045508 C 2 .1E-3 -7.28683 -86.2569 75.27724 9.233558 99.92895 -46.7004 .3738012 -.372372 -.046168 C 3 .15E-3 -7.41337 -85.3888 75.58585 9.362285 99.84016 -45.0254 .3737174 -.366676 -.046811 C 4 .2E-3 -7.5373 -84.4904 75.86758 9.487714 99.71589 -43.3345 .3735005 -.360851 -.047439 C 5 .25E-3 -7.65858 -83.5619 76.12234 9.609793 99.5562 -41.6281 .3731507 -.354897 -.048049 C 6 .3E-3 -7.77715 -82.6038 76.35003 9.728472 99.36113 -39.9069 .3726683 -.348817 -.048642 C 7 .35E-3 -7.89297 -81.6164 76.55059 9.843704 99.13076 -38.1716 .3720534 -.342613 -.049219 C 8 .4E-3 -8.00599 -80.6 76.72395 9.955443 98.86517 -36.4227 .3713062 -.336287 -.049777 C 9 .45E-3 -8.11617 -79.5549 76.87004 10.06365 98.56446 -34.6608 .3704271 -.329843 -.050318 C 10 .5E-3 -8.22347 -78.4816 76.98882 10.16828 98.22873 -32.8867 .3694163 -.32328 -.050841 C 11 .55E-3 -8.32786 -77.3804 77.08024 10.2693 97.85809 -31.1008 .3682743 -.316604 -.051347 C 12 .6E-3 -8.42928 -76.2517 77.14427 10.36668 97.45269 -29.3039 .3670014 -.309814 -.051833 C 13 .65E-3 -8.52771 -75.0959 77.1809 10.46037 97.01266 -27.4966 .3655981 -.302915 -.052302 C 14 .7E-3 -8.62311 -73.9135 77.1901 10.55035 96.53816 -25.6795 .364065 -.295908 -.052752 C 15 .75E-3 -8.71543 -72.7047 77.17185 10.63653 96.02937 -23.8533 .362402 -.288795 -.053183 C 16 .8E-3 -8.80466 -71.47 77.12599 10.71883 95.48645 -22.0187 .3606083 -.281579 -.053594 C 17 .85E-3 -8.89082 -70.21 77.05293 10.79745 94.90961 -20.1762 .3586886 -.274264 -.053987 C 18 .9E-3 -8.97382 -68.925 76.95252 10.87225 94.29905 -18.3265 .3566416 -.266852 -.054361 C 19 .95E-3 -9.05361 -67.6156 76.82477 10.94318 93.65499 -16.4703 .354468 -.259345 -.054716 C 20 .1E-2 -9.13015 -66.2822 76.66974 11.01022 92.97765 -14.6083 .3521686 -.251745 -.055051 BLANK card ending plot cards BEGIN NEW DATA CASE C BENCHMARK DC-9 C 5th of 10 subcases is like preceding, but with JMarti frequency C dependence replaced by Semlyen branch cards from DC-31. This example C comes from BPA's Dr. Tsu-huei Liu, who first made it work. PRINTED NUMBER WIDTH, 11, 2, { Reduced precision .000050 .00100 1 1 1 1 1 C L= 138.0MILES, RHO= 27.0, SS FREQ= 60.00, NSS=0, KFIT= 0, KPS=2, KYC=10 C 1.3636 .05215 4 1.602 -20.75 50. 50. C 1.3636 .05215 4 1.602 -19.25 50. 50. C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5 C 2.3636 .05215 4 1.602 0.75 77.5 77.5 C 3.3636 .05215 4 1.602 19.25 50. 50. C 3.3636 .05215 4 1.602 20.75 50. 50. C 0.5 2.61 4 0.386 -12.9 98.5 98.5 C 0.5 2.61 4 0.386 12.9 98.5 98.5 C C 27. 5000. 1 138. C 27. 60.00 1 138. C 27. 6.00 1 138. 6 20 C -1SENDA RECA 0.55456E-02 0.77998E-03 1 1 2 2 3 0.14938602E+02 0.90516313E+02-0.16261865E-03 0.14500324E-02 0.60000000E+02 0.00000E+00 0.15161E+05 0.75119E+00 0.00000E+00 0.17105E+04 0.24881E+00 0.00000E+00 0.59584E+03-0.11954E-02 0.00000E+00 0.39933E+05-0.74162E-03 -1SENDB RECB 0.74392E-02 0.74149E-03 2 2 2 2 3 0.83801231E+00 0.38634735E+02-0.44536501E-04 0.20637999E-02 0.60000000E+02 0.00000E+00 0.26608E+06 0.83767E+00 0.00000E+00 0.35627E+05 0.16233E+00 0.00000E+00 0.13630E+03-0.49171E-03 0.00000E+00 0.32940E+05-0.75349E-04 -1SENDC RECC 0.42823E-02 0.74017E-03 3 3 2 2 3 0.26535168E+01 0.67894100E+02-0.45225808E-04 0.11682248E-02 0.60000000E+02 0.00000E+00 0.61698E+06 0.90903E+00 0.00000E+00 0.15239E+05 0.90969E-01 0.00000E+00 0.25667E+03-0.34465E-03 0.00000E+00 0.11831E+05-0.28233E-04 0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 0.71114E+00 0.00000E+00-0.41983E-15 0.00000E+00 0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 0.42027E+00 0.00000E+00 0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 0.22471E+00 0.00000E+00-0.27373E-15 0.00000E+00 0.84054E+00 0.00000E+00 0.42027E+00 0.00000E+00-0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 BLANK card follows the last branch card BLANK line terminates the last (here, nonexistent) switch 14SENDA 408000. 60. 0.0 14SENDB 408000. 60. -120. 14SENDC 408000. 60. 120. BLANK card ending source data C The following initial-condition cards were constructed from branch flow C output with sources present for a phasor solution. Full precision of C the ATP output is being used, note. Purely resistive branches are not C supposed to require initialization, so these were omitted. C C SENDA 408000. 408000. -14.35539062039 427.20215266352 -.29284996866E7 -.29284996866E7 C 0.0 0.0 426.96089048118 91.9256883 -.87100021658E8 87100021.6581596 C C RECA 426275.19929468 426279.30459602 .1509903313E-13 .9664560746E-12 .23142460224E-8 C 1870.8268620134 0.2514568 -.966338121E-12 -89.1048263 .20597711132E-6 C C SENDB -204000. 408000. 390.21858493854 442.08847622165 -.30943350074E7 -.30943350074E7 C -353338.3647441 -120.0000000 -207.7779506506 -28.0337631 -.90132949314E8 90132949.3138173 C C RECB -211353.2339786 423078.09601739 .6536993169E-12 .7963152328E-12 .14252537328E-7 C -366503.5959121 -119.9709473 -.454747351E-12 -34.8244892 -.1678477368E-6 C C SENDC -204000. 408000. -374.9796048273 439.89286922318 -.2384468398E7 -.2384468398E7 C 353338.36474405 120.0000000 -229.9913745272 -148.4773916 -.89706460393E8 89706460.3928106 C C RECC -210971.8157184 425738.55951138 -.170530257E-12 .190658623E-12 .33753610037E-7 C 369789.41849017 119.7055709 .8526512829E-13 153.4349488 -.2253587274E-7 2SENDA 408000. 0.0 60. 2RECA 426275.199295 1870.82686201 2SENDB -204000. -353338.36474 2RECB -211353.23398 -366503.59591 2SENDC -204000. 353338.36474 2RECC -210971.81572 369789.41849 3SENDA RECA -14.35539062039426.96089048118.1509903313E-13-.966338121E-12 3SENDB RECB 390.21858493854-207.7779506506.6536993169E-12-.454747351E-12 3SENDC RECC -374.9796048273-229.9913745272-.170530257E-12.8526512829E-13 1 { Request for the output of each node voltage C First 6 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Step Time RECC RECB RECA SENDA SENDB SENDC C 0 0.0 -210971.8 -211353.2 426275.2 408000. -204000. -204000. C 1 .5E-4 -217903.1 -204391.8 426140.86 407927.52 -197303.9 -210623.6 C 2 .1E-3 -224764. -197371.2 425871.55 407710.1 -190537.7 -217172.4 C 3 .15E-3 -231544. -190283.3 425452.69 407347.83 -183703.7 -223644.1 C 4 .2E-3 -238240.9 -183128.9 424883.53 406840.83 -176804.6 -230036.3 C 5 .25E-3 -244852.7 -175910. 424163.92 406189.28 -169842.6 -236346.7 C 6 .3E-3 -251377.2 -168628.7 423293.95 405393.41 -162820.2 -242573.2 C 7 .35E-3 -257812.1 -161287.4 422273.86 404453.51 -155740. -248713.5 C 8 .4E-3 -264155.2 -153888.8 421103.95 403369.91 -148604.5 -254765.4 C 9 .45E-3 -270404.2 -146435.4 419784.61 402142.99 -141416.2 -260726.8 C 10 .5E-3 -276557. -138929.9 418316.3 400773.2 -134177.6 -266595.6 C 11 .55E-3 -282611.4 -131374.9 416699.51 399261.01 -126891.3 -272369.7 C 12 .6E-3 -288565.2 -123773.2 414934.81 397606.96 -119560. -278046.9 C 13 .65E-3 -294416.4 -116127.4 413022.82 395811.65 -112186.2 -283625.4 C 14 .7E-3 -300162.8 -108440.2 410964.19 393875.71 -104772.5 -289103.2 C 15 .75E-3 -305803.1 -100717.9 408762.02 391799.82 -97321.65 -294478.2 C 16 .8E-3 -311336.8 -92970.64 406424.38 389584.73 -89836.18 -299748.6 C 17 .85E-3 -316755.1 -85178.69 403935.32 387231.23 -82318.79 -304912.4 C 18 .9E-3 -322059.1 -77353.3 401302.31 384740.14 -74772.16 -309968. C 19 .95E-3 -327248.4 -69499.21 398526.63 382112.35 -67198.95 -314913.4 C 20 .1E-2 -332321.6 -61619.94 395609.19 379348.81 -59601.88 -319746.9 BLANK card ending plot cards BEGIN NEW DATA CASE PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) .000050 .001 1 1 1 1 C 6th of 10 subcases illustrates multiple frequencies for initial condit. C Create 3 parallel, independent, identical circuits of line plus resistor. C The (SEND, REC) alternative will be initialized for 60 Hz whereas the C alternative (GEN, END) and (GEN, MARTI) will use 50 Hz --- with the C former being constant-parameter and the latter JMARTI. It works. There is C perfect agreement with related solution that begins with phasors. The C Marti line is identical to what already was used (see 3rd subcase), so C documentation will be minimized. Comment cards have been dropped, and the C 1 of col. 31 is restored to 2. -1SEND REC .306 5.82 .012 100. { Half the length of original -1GEN END .306 5.82 .012 100. { Half the length of original -1GEN MARTI 2. 0.00 -2 1 25 4.6006674814379022100E+02 -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 3.72990420164290220E+07 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 2.89354778735701405E+07 14 3.8038812947515656600E-04 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 -7.93497908986478425E+10 3.97029750605793992E+10 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 1.82589462762646523E+05 1.82772052225409076E+05 0.58951119 0.00000000 REC 100. END 100. MARTI 100. BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. 14GEN 100. 50. C 14GEN2 100. 50. -1. BLANK card ending source data C GEN 100. 100. .26567063730001 .43724454992683 13.283531865001 3.0530027479459 C 0.0 0.0 -.3472778554952 -52.5836910 17.363892774758 17.3638928 C C END 25.859996851553 45.233901262338 -.2585999685155 .45233901262338 -10.23052911705 C -37.11288706432 -55.1314894 .37112887064316 124.8685106 -.737257477E-17 C C GEN 100. 100. .18468803864924 .33562381569663 9.2344019324622 3.4793989521433 C 0.0 0.0 -.2802386019853 -56.6135784 14.011930099264 14.0119301 C C MARTI 18.346330639813 33.926399692036 -.1834633063981 .33926399692036 -5.755002980319 C -28.53791772569 -57.2640259 .28537917725693 122.7359741 .7014788056E-15 C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 2REC 19.956813332748-34.69967901024 2GEN 100. 0.0 2END 25.859996851553-37.11288706432 2MARTI 18.346330639813-28.53791772569 C Distributed: I-k real I-k imag I-m real I-m imag [Hz] 3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239 60. 3GEN END .26567063730001-.3472778554952-.2585999685155.37112887064316 50. 3GEN MARTI .18468803864924-.2802386019853-.1834633063981.28537917725693 C Note that preceding branch card has blank cols. 75-80 for the frequency. By C definition, frequency will remain unchanged at 50 Hz (2nd Type-3 card). SEND GEN REC END MARTI { Names of nodes for voltage output BLANK card ending node voltage outputs C First 5 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Step Time SEND GEN REC END MARTI C 0 0.0 100. 100. 19.956813 25.859997 18.346331 C 1 .5E-4 99.982235 99.987663 20.606552 26.439082 18.791563 C 2 .1E-3 99.928947 99.950656 21.249697 27.012298 19.232816 C 3 .15E-3 99.840155 99.888987 21.885291 27.578849 19.669386 C 4 .2E-3 99.71589 99.802673 22.51311 28.138595 20.101117 C 5 .25E-3 99.556196 99.691733 23.13293 28.691398 20.527892 C 6 .3E-3 99.361131 99.556196 23.744531 29.237122 20.949603 C 7 .35E-3 99.130763 99.396096 24.347696 29.775632 21.366146 C 8 .4E-3 98.865174 99.21147 24.94221 30.306796 21.777336 C 9 .45E-3 98.56446 99.002366 25.527863 30.830482 22.183024 C 10 .5E-3 98.228725 98.768834 26.104445 31.346561 22.583287 C 11 .55E-3 97.85809 98.510933 26.671752 31.854905 22.97808 C 12 .6E-3 97.452687 98.228725 27.229584 32.35539 23.367227 C 13 .65E-3 97.01266 97.922281 27.77774 32.847892 23.750616 C 14 .7E-3 96.538164 97.591676 28.316028 33.332289 24.128147 C 15 .75E-3 96.029369 97.236992 28.844255 33.808461 24.499706 C 16 .8E-3 95.486454 96.858316 29.362233 34.276292 24.865139 C 17 .85E-3 94.909614 96.455742 29.869718 34.735588 25.224341 C 18 .9E-3 94.299054 96.029369 30.366492 35.186194 25.577322 C 19 .95E-3 93.654989 95.579301 30.852618 35.628304 25.924072 C 20 .1E-2 92.977649 95.105652 31.327783 36.061624 26.264482 BLANK card ending plot cards BEGIN NEW DATA CASE C BENCHMARK DC-9 C 7th of 10 subcases extends the preceding. It has two 3-phase Semlyen C lines with different frequencies (50 and 60 Hz). Semlyen branch cards C are copied from the 5th subcase (from DC-31), so omit the comments. PRINTED NUMBER WIDTH, 11, 2, { Allow output line nicely within 80 columns .000050 .001 { Just 20 steps. Extend to 40 msec if plot card is enabled 1 -1 1 1 { Note IPLOT of cols. 9-16 disables plotting -1SENDA RECA 0.55456E-02 0.77998E-03 1 1 2 2 3 0.14938602E+02 0.90516313E+02-0.16261865E-03 0.14500324E-02 0.60000000E+02 0.00000E+00 0.15161E+05 0.75119E+00 0.00000E+00 0.17105E+04 0.24881E+00 0.00000E+00 0.59584E+03-0.11954E-02 0.00000E+00 0.39933E+05-0.74162E-03 -1SENDB RECB 0.74392E-02 0.74149E-03 2 2 2 2 3 0.83801231E+00 0.38634735E+02-0.44536501E-04 0.20637999E-02 0.60000000E+02 0.00000E+00 0.26608E+06 0.83767E+00 0.00000E+00 0.35627E+05 0.16233E+00 0.00000E+00 0.13630E+03-0.49171E-03 0.00000E+00 0.32940E+05-0.75349E-04 -1SENDC RECC 0.42823E-02 0.74017E-03 3 3 2 2 3 0.26535168E+01 0.67894100E+02-0.45225808E-04 0.11682248E-02 0.60000000E+02 0.00000E+00 0.61698E+06 0.90903E+00 0.00000E+00 0.15239E+05 0.90969E-01 0.00000E+00 0.25667E+03-0.34465E-03 0.00000E+00 0.11831E+05-0.28233E-04 0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 0.71114E+00 0.00000E+00-0.41983E-15 0.00000E+00 0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 0.42027E+00 0.00000E+00 0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 0.22471E+00 0.00000E+00-0.27373E-15 0.00000E+00 0.84054E+00 0.00000E+00 0.42027E+00 0.00000E+00-0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 RECA 200. RECB 200. RECC 200. SENDC 200. { Replace phase c source by this resistor C Illustrate reference branch copying for Semlyen model. 2nd circuit follows. C But because the 2nd will be excited at 50 Hz whereas the 1st has its phasor C parameters derived at 60 Hz, there will be warning messages about 2nd. Note C that the phasor solution of the second is substantially in error (50 Hz is C not close to 60 Hz). Yet, the result is amazingly smooth. To see this, C extend TMAX to a couple of cycles (40 msec), and enable plotting (now, C IPLOT = -1 will ignore the plot card, which involves RECA, RECB, RECC). -1GENA ENDA SENDA RECA -1GENB ENDB SENDB RECB -1GENC ENDC SENDC RECC ENDA 200. ENDB 200. ENDC 200. GENC 200. { Replace phase c source by this resistor BLANK card follows the last branch card BLANK line terminates the last (here, nonexistent) switch 14SENDA 100. 60. 0.0 14SENDB 100. 60. -120. 14GENA 100. 50. 0.0 14GENB 100. 50. -120. BLANK card ending source data C The following initial-condition cards were constructed from branch flow C output with sources present for a phasor solution. Full precision of C the ATP output is being used, note. Purely resistive branches are not C supposed to require initialization, so these were omitted. C The following 3 branch flows are for 60 Hz: C SENDA 100. 100. .36998975926495 .37180820218231 18.499487963247 2.9848546654506 C 0.0 0.0 -.0367276088128 -5.6689829 1.8363804406416 1.8363804 C C RECA 74.690024760283 78.777238585258 -.3734501238014 .39388619292629 -15.5146332978 C -25.04502985478 -18.5373124 .12522514927392 161.4626876 -.253833413E-14 C C SENDB -50. 100. -.3685220242714 .49798629886772 23.71618850994 -1.806979904073 C -86.60254037844 -120.0000000 -.3349356228992 -137.7334962 7.5840811711722 7.5840812 C C RECB -86.37515570642 101.04091926346 .43187577853211 .5052045963173 -25.52316841401 C -52.4270907289 -148.7435235 .26213545364451 31.2564765 -.116668827E-13 C C SENDC 9.9402547548508 10.505728992727 -.0497012737743 .05252864496364 -.2759258541716 -.4312032208709 C -3.40024662007 -18.8841944 .01700123310035 161.1158056 -.204060401E-15 0.0000000 C C RECC -7.010739885979 7.8810498462914 .0350536994299 .03940524923146 -.1552773666993 C 3.6000656842441 152.8191734 -.0180003284212 -27.1808266 -.154803741E-15 C The following 3 branch flows are for 50 Hz: C GENA 100. 100. .36998975926495 .37180820218231 18.499487963247 2.9848546654506 C 0.0 0.0 -.0367276088128 -5.6689829 1.8363804406416 1.8363804 C C ENDA 74.690024760283 78.777238585258 -.3734501238014 .39388619292629 -15.5146332978 C -25.04502985478 -18.5373124 .12522514927392 161.4626876 -.865843855E-14 C C GENB -50. 100. -.3685220242714 .49798629886772 23.71618850994 -1.806979904073 C -86.60254037844 -120.0000000 -.3349356228992 -137.7334962 7.5840811711722 7.5840812 C C ENDB -86.37515570642 101.04091926346 .43187577853211 .5052045963173 -25.52316841401 C -52.4270907289 -148.7435235 .26213545364451 31.2564765 .1027303242E-13 C C GENC 9.9402547548508 10.505728992728 -.0497012737743 .05252864496364 -.2759258541716 -.4312032208709 C -3.40024662007 -18.8841944 .01700123310035 161.1158056 -.319839641E-16 0.0000000 C C ENDC -7.010739885979 7.8810498462914 .0350536994299 .03940524923146 -.1552773666993 C 3.6000656842441 152.8191734 -.0180003284212 -27.1808266 -.56249764E-16 C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SENDA 100.0 0.0 2RECA 74.690024760283-25.04502985478 2SENDB -50.-86.60254037844 2RECB -86.37515570642 -52.4270907289 2SENDC 9.9402547548508 -3.40024662007 2RECC -7.0107398859793.6000656842441 2GENA 100.0 0.0 2ENDA 74.690024760283-25.04502985478 2GENB -50.-86.60254037844 2ENDB -86.37515570642 -52.4270907289 2GENC 9.9402547548508 -3.40024662007 2ENDC -7.0107398859793.6000656842441 C Distributed: I-k real I-k imag I-m real I-m imag [Hz] 3SENDA RECA .36998975926495-.0367276088128-.3734501238014.12522514927392 60. 3SENDB RECB -.3685220242714-.3349356228992.43187577853211.26213545364451 3SENDC RECC -.0497012737743.01700123310035 .0350536994299-.0180003284212 3GENA ENDA .36998975926495-.0367276088128-.3734501238014.12522514927392 50. 3GENB ENDB -.3685220242714-.3349356228992.43187577853211.26213545364451 3GENC ENDC -.0497012737743.01700123310035 .0350536994299-.0180003284212 C About phasor frequency on the right, note that only the first of each group C of 3 coupled conductors requires a nonzero frequency. A blank for following C phases means that the declared value will remain in effect. As for the 1st C frequency (60 Hz), this could be omitted provided it is equal to the power C frequency. But different places in the world have different power frequency C so it is prudent practice to declare this frequency explicitely. RECA ENDA RECB ENDB RECC ENDC { Names of nodes for voltage output BLANK card ending node names for voltage output C Step Time RECA ENDA RECB ENDB RECC ENDC C 0 0.0 74.690025 74.690025 -86.37516 -86.37516 -7.01074 -7.01074 C 1 .5E-4 75.163363 75.132945 -85.37169 -84.34315 -7.059456 -7.969316 C 2 .1E-3 75.595736 75.435454 -84.33743 -82.81562 -7.123114 -8.190566 C 3 .15E-3 76.001094 75.753382 -83.27324 -81.54915 -7.18451 -8.234526 C 4 .2E-3 76.379406 76.079746 -82.17945 -80.40443 -7.243416 -8.226221 C 5 .25E-3 76.730556 76.398672 -81.05645 -79.31142 -7.299777 -8.200803 C 6 .3E-3 77.054425 76.700684 -79.90465 -78.23596 -7.353565 -8.170835 C 7 .35E-3 77.350898 76.981153 -78.72446 -77.16148 -7.404759 -8.141389 C 8 .4E-3 77.619873 77.238057 -77.51631 -76.07989 -7.453337 -8.114506 C 9 .45E-3 77.861255 77.470657 -76.28063 -74.98724 -7.499282 -8.090883 C 10 .5E-3 78.07496 77.678819 -75.01785 -73.88158 -7.542576 -8.070602 C 11 .55E-3 78.260914 77.862692 -73.72842 -72.76194 -7.583202 -8.05347 C 12 .6E-3 78.419051 78.02254 -72.41281 -71.62787 -7.621145 -8.039181 C 13 .65E-3 78.549315 78.158673 -71.07147 -70.47916 -7.65639 -8.027389 C 14 .7E-3 78.651663 78.271409 -69.70489 -69.31574 -7.688924 -8.017745 C 15 .75E-3 78.726915 78.367856 -68.31324 -68.14096 -7.71979 -8.014505 C 16 .8E-3 78.777534 78.482038 -66.89676 -66.94954 -7.75238 -8.021203 C 17 .85E-3 78.796385 78.578276 -65.45863 -65.69399 -7.778268 -7.991112 C 18 .9E-3 78.78656 78.66135 -63.99679 -64.41226 -7.800598 -7.945947 C 19 .95E-3 78.749002 78.71769 -62.51169 -63.1208 -7.819641 -7.907277 C 20 .1E-2 78.683581 78.746939 -61.00415 -61.8202 -7.835692 -7.875092 C ----------------- Preceding was before Tsu-huei's OVER8 correction. C Step Time RECA ENDA RECB ENDB RECC ENDC C 0 0.0 74.690025 74.690025 -86.37516 -86.37516 -7.01074 -7.01074 C 1 .5E-4 75.163363 74.744884 -85.37169 -84.65358 -7.059456 -8.300563 C 2 .1E-3 75.595736 75.003289 -84.33743 -83.20507 -7.123114 -8.557454 C 3 .15E-3 76.001094 75.325776 -83.27324 -81.97004 -7.18451 -8.597739 C 4 .2E-3 76.379406 75.66166 -82.17945 -80.84254 -7.243416 -8.582616 C 5 .25E-3 76.730556 75.989351 -81.05645 -79.76068 -7.299777 -8.551501 C 6 .3E-3 77.054425 76.299083 -79.90465 -78.69329 -7.353565 -8.516962 C 7 .35E-3 77.350898 76.58653 -78.72446 -77.62521 -7.404759 -8.483715 C 8 .4E-3 77.619873 76.849919 -77.51631 -76.54911 -7.453337 -8.453539 C 9 .45E-3 77.861255 77.088671 -76.28063 -75.46146 -7.499282 -8.42697 C 10 .5E-3 78.07496 77.302754 -75.01785 -74.36051 -7.542576 -8.403987 C 11 .55E-3 78.260914 77.492377 -73.72842 -73.24544 -7.583202 -8.384334 C 12 .6E-3 78.419051 77.657849 -72.41281 -72.11583 -7.621145 -8.367661 C 13 .65E-3 78.549315 77.799509 -71.07147 -70.97152 -7.65639 -8.353591 C 14 .7E-3 78.651663 77.917697 -69.70489 -69.81243 -7.688924 -8.341754 C 15 .75E-3 78.726915 77.996637 -68.31324 -68.63631 -7.71979 -8.317196 C 16 .8E-3 78.777534 78.003092 -66.89676 -67.41624 -7.75238 -8.217882 C 17 .85E-3 78.796385 78.078973 -65.45863 -66.12507 -7.778268 -8.151926 C 18 .9E-3 78.78656 78.16892 -63.99679 -64.82237 -7.800598 -8.093319 C 19 .95E-3 78.749002 78.231147 -62.51169 -63.5241 -7.819641 -8.050547 C 20 .1E-2 78.683581 78.263977 -61.00415 -62.22435 -7.835692 -8.017481 144 4. 0.0 40. ENDA ENDB ENDC BLANK card ending plot cards BEGIN NEW DATA CASE C 8th of 10 subcases demonstrates manual initialization of a JMarti line. C Use two different frequencies. Single phase is used for simplicity. C One frequency will be very low (.01 Hz) to prove no need for true dc. C This subcase will use phasor solution for .01 Hz. The following subcase C will replace the phasor solution here by initial condition cards. PRINTED NUMBER WIDTH, 13, 2, .000050 .001 1 1 1 1 -1SEND REC 1. 0.00 -2 1 25 4.6006674814379022100E+02 -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 3.72990420164290220E+07 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 2.89354778735701405E+07 14 3.8038812947515656600E-04 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 -7.93497908986478425E+10 3.97029750605793992E+10 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 1.82589462762646523E+05 1.82772052225409076E+05 0.58951119 0.00000000 REC 100. -1GEN END SEND REC { 2nd JMarti line is disconnected from the first END 100. BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. { No phasor presence; I.C. cards will be used 14GEN 100. .010 { Low-freq phasor solution } -1. BLANK card ending source data C The following initial-condition cards were constructed from branch flow C output with the source present for a phasor solution. Full precision of C the ATP output is being used, note. C The following branch flow is for 60 Hz: C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678 C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865 C C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945 C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15 C The following initial condition cards are built from preceding phasor flow. C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 2REC 14.563049373597-25.89164864838 C Distributed: I-k real I-k imag I-m real I-m imag F [Hz] 3SEND REC .14695409358625-.2529237293961-.145630493736 .25891648648378 60. C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Step Time SEND REC GEN END C 0 0.0 100. 14.56304937 100. 74.15408515 C 1 .5E-4 99.98223524 15.04761164 100. 74.15408201 C 2 .1E-3 99.92894726 15.52757104 100. 74.15407898 C 3 .15E-3 99.84015501 16.00208444 100. 74.15407596 C 4 .2E-3 99.71589003 16.47092745 99.99999999 74.15407293 C 5 .25E-3 99.55619646 16.93392289 99.99999999 74.1540699 C 6 .3E-3 99.36113105 17.39090307 99.99999998 74.15406686 C 7 .35E-3 99.13076311 17.84170463 99.99999998 74.15406383 C 8 .4E-3 98.86517447 18.28608879 99.99999997 74.15406078 C 9 .45E-3 98.56445951 18.72384749 99.99999996 74.1540577 C 10 .5E-3 98.22872507 19.15499653 99.99999995 74.15405464 C 11 .55E-3 97.85809043 19.57943512 99.99999994 74.1540516 C 12 .6E-3 97.45268728 19.99693921 99.99999993 74.15404857 C 13 .65E-3 97.01265965 20.40734552 99.99999992 74.15404553 C 14 .7E-3 96.53816388 20.81050376 99.9999999 74.15404249 C 15 .75E-3 96.02936857 21.20624627 99.99999989 74.15403945 C 16 .8E-3 95.48645447 21.59436163 99.99999987 74.1540364 C Showing 16 steps is more than enough. Without manual initial conditions, the C surge arrives at the receiving end on step 8. Above exceeds 2 travel times. C Furthermore, above agrees exactly to output that includes a phasor solution. SEND REC GEN END { Node names for node voltage output BLANK card ending output requests BLANK card ending plot cards BEGIN NEW DATA CASE C 9th of 10 subcases has the same two JMarti lines as the preceding. C It replaces the low-frequency phasor solution of the preceding by C by initial condition cards. This is manual, multi-frequency JMarti. PRINTED NUMBER WIDTH, 13, 2, .000050 .001 1 1 1 1 -1SEND REC 1. 0.00 -2 1 25 4.6006674814379022100E+02 -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 3.72990420164290220E+07 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 2.89354778735701405E+07 14 3.8038812947515656600E-04 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 -7.93497908986478425E+10 3.97029750605793992E+10 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 1.82589462762646523E+05 1.82772052225409076E+05 0.58951119 0.00000000 REC 100. -1GEN END SEND REC C GEN END 100. END 100. BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. C 14GEN 100. 60. -1. 14GEN 100. .010 BLANK card ending source data C The following initial-condition cards were constructed from branch flow C output with the source present for a phasor solution. Full precision of C the ATP output is being used, note. C The following branch flow is for 60 Hz: C SEND 100. 100. .14695409358625 .29251652690644 7.3477046793123 2.9354052953678 C 0.0 0.0 -.2529237293961 -59.8424965 12.646186469807 12.6461865 C C REC 14.563049373597 29.706226229343 -.145630493736 .29706226229343 -4.412299383945 C -25.89164864838 -60.6439325 .25891648648378 119.3560675 -.152655666E-15 C The following branch flow is for .01 Hz: C GEN 100. 100. .74154511249386 .74160775929634 37.077255624693 9.5784693158927 C 0.0 0.0 .00963923155659 0.7447371 -.4819615778294 -0.4819616 C C END 74.154085146507 74.16034831202 -.7415408514651 .7416034831202 -27.4987863088 C .96380384134885 0.7446492 -.0096380384135 -179.2553508 .5416403003E-15 C C The following initial condition cards are built from preceding phasor flow. C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 2GEN 100. 0.0 2REC 14.563049373597-25.89164864838 2END 74.154085146507.96380384134885 C Distributed: I-k real I-k imag I-m real I-m imag F [Hz] 3SEND REC .14695409358625-.2529237293961 -.145630493736.25891648648378 60. 3GEN END .74154511249386.00963923155659-.7415408514651-.0096380384135 .01 C Showing 16 steps is more than enough. Without manual initial conditions, the C surge arrives at the receiving end on step 8. Above exceeds 2 travel times. C Furthermore, above agrees exactly to output that includes a phasor solution. SEND REC GEN END { Node names for node voltage output C Step Time SEND REC GEN END C 0 0.0 100. 14.56304937 100. 74.15408515 C 1 .5E-4 99.98223524 15.04761164 100. 74.15408201 C 2 .1E-3 99.92894726 15.52757104 100. 74.15407898 C 3 .15E-3 99.84015501 16.00208444 100. 74.15407596 C 4 .2E-3 99.71589003 16.47092745 99.99999999 74.15407293 C 5 .25E-3 99.55619646 16.93392289 99.99999999 74.1540699 C 6 .3E-3 99.36113105 17.39090307 99.99999998 74.15406686 C 7 .35E-3 99.13076311 17.84170463 99.99999998 74.15406383 C 8 .4E-3 98.86517447 18.28608879 99.99999997 74.15406078 C 9 .45E-3 98.56445951 18.72384749 99.99999996 74.1540577 C 10 .5E-3 98.22872507 19.15499653 99.99999995 74.15405464 C 11 .55E-3 97.85809043 19.57943512 99.99999994 74.1540516 C 12 .6E-3 97.45268728 19.99693921 99.99999993 74.15404857 C 13 .65E-3 97.01265965 20.40734552 99.99999992 74.15404553 C 14 .7E-3 96.53816388 20.81050376 99.9999999 74.15404249 C 15 .75E-3 96.02936857 21.20624627 99.99999989 74.15403945 C 16 .8E-3 95.48645447 21.59436163 99.99999987 74.1540364 BLANK card ending output requests BLANK card ending plot cards BEGIN NEW DATA CASE C 10th of 10 subcases is like the 1st except that a repeated initial C condition has been added to illustrate the new trap for these. Use C 1st subcase data prior to 19 July 2009 (see comments in 1st subcase). C WSM appends this subcase 1 August 2009. .000050 .002 1 1 0 0 { No need to see connectivity or phasor output -1SEND REC .306 5.82 .012 100. { Half the length of original -1GEN END .306 5.82 .012 100. { Half the length of original REC 100. 1 END 100. 1 BLANK card terminating branch data BLANK card terminating all (in this case, nonexistent) switches 14SEND 100. 60. { Excitation of (SEND, REC) is transient only 14GEN 100. 60. { Note phasor presence } -1. BLANK card ending source data C (SEND, REC) had no phasor solution, initialize it using values of (GEN, END): C Name K V-k real V-k imag ----- Node V initial condition is type 2 2SEND 100. 0.0 60. 2REC 19.956813332748-34.69967901024 C Distributed: I-k real I-k imag I-m real I-m imag 3SEND REC .20753030414706-.3196504046191-.1995681333275.34699679010239 3REC .19956813332748 { Unnecessary, since purely resistive branch 3REC .19956813332748 { This is the duplicate initial condition SEND GEN REC END { List of nodes for node-voltage output BLANK card ending node voltage outputs BLANK card ending plot cards BEGIN NEW DATA CASE BLANK