BEGIN NEW DATA CASE C BENCHMARK DCNEW-6 C Derive Marti representation of 400-kV overhead Portuguese line as used C in subsequent simulation of DCNEW-7. From J. Allen Lima of EDP, 1981. C Discontinuity of April, 1992: the default conductance to ground was made C smaller, which changes the answers. To reproduce the previous answers C (we do not want to change these, for historical reasons), we have replaced C the previously-blank columns 33-40 of the Marti fitting card by 0.48D-7. C If erased, the smaller current default value of 0.32D-8 would be used. C Generally, users are advised to use the new default value. The old one C was found to be poorer (less stable) at preserving trapped charge, which C had a tendency to either decay or grow, or both. See July, 92, newsletter C 1st of 4 data cases. This overhead line case is followed by 3 for cables. JMARTI SETUP $ERASE C LINHA CEDILLO-RIO MAIOR (EDP2.DAT) C EXEMPLO DE APLICACAO - ANEXO 2 BRANCH CEDIL0RIO M0CEDIL8RIO M8CEDIL4RIO M4 LINE CONSTANTS 10.231 0.0840 4 1.25 -39.4 70.97 32.25 15.8 0.0 2 20.231 0.0840 4 1.25 0.0 70.97 32.25 15.8 0.0 2 30.231 0.0840 4 1.25 39.4 70.97 32.25 15.8 0.0 2 00.5 0.58 4 .575 -25.76 100.07 76.12 00.5 0.58 4 .575 25.76 100.07 76.12 BLANK card ending conductor cards of imbedded "LINE CONSTANTS" data case 250.0 5000.0 1 85.75 1 250.0 50.0 1 85.75 1 250.0 .01 1 85.75 9 10 1 BLANK card ending frequency cards of imbedded "LINE CONSTANTS" data case BLANK card ending "LINE CONSTANTS" data cases C 345678901234567890123456789012345678901234567890123456789012345678901234567890 C LECT +1 -4 1 .48D-7 C .30 30 0 1 1 0 0 C .30 30 0 1 1 0 0 .0 .30 30 0 1 3 0 0 .30 30 0 1 3 0 0 .0 $PUNCH C Punched cards of "JMARTI SETUP" that began at 11-Aug-87 14.15.29 C -1CEDIL0RIO M0 2. 0.00 -2 3 C 13 0.41285804706219800000E+03 C 0.117013827128701600E+04 0.109328155224148600E+05 0.565854842148188200E+04 C 0.486728968691983400E+04 0.707738202163730700E+05 0.683747081419827600E+06 C 0.353541780131411400E+07 0.158128637229056500E+08 0.342732557798577400E+08 C 0.393504407219363000E+08 0.341282907181993800E+08 0.204464663801654300E+08 C 0.576218378853537400E+08 C 0.658089815842887500E+01 0.690898889439574300E+02 0.124258964382556600E+03 C 0.223516853729638100E+03 0.270128036573363500E+04 0.269143242842107800E+05 C 0.146979058765937700E+06 0.698197935976008100E+06 0.314038859911038900E+07 C 0.746337177863556100E+07 0.266392944394791100E+08 0.155396621655433500E+08 C 0.219748648076802900E+08 C 14 0.50730990999172660000E-03 C 0.349887778185420700E+00 0.385183799451747400E+00 0.266248332878260700E+00 C 0.101302703627398000E+01 0.108212455213644600E+01 0.239521017601984200E+01 C 0.279893774926274700E+02 0.197810245556054600E+03 0.567095451317796800E+03 C 0.101827034029151500E+05 0.198712743344115800E+06 0.108881451987841900E+04 C -0.162015835275688900E+08 0.159908008795460600E+08 C 0.707841765138503400E+02 0.823999719984897900E+02 0.532720131788135200E+02 C 0.202956475231847800E+03 0.233040781840643500E+03 0.490731499989644900E+03 C 0.235012914821689300E+04 0.605706538137324200E+04 0.493208405106231800E+04 C 0.198963455259852500E+05 0.369358647098982800E+05 0.521960681977429600E+06 C 0.445816728405941500E+05 0.446262545134347200E+05 C -2CEDIL8RIO M8 2. 0.00 -2 3 C 14 0.32434436644879670000E+03 C 0.252272221491383900E+04 -0.783260653525317500E+03 0.843607760308153000E+03 C 0.405537446805999600E+03 0.172880473437216700E+03 0.230195610382546000E+03 C 0.660733376851304400E+02 0.356064973429300400E+04 0.360876600484696400E+05 C 0.458067308721601900E+06 0.155426951919825100E+07 0.349466741634103300E+07 C 0.831361356808634600E+07 0.267191748194155000E+08 C 0.395606796435294500E+01 0.463636540679569400E+01 0.937357763832626400E+01 C 0.153396557146477500E+02 0.203063670799592800E+02 0.331878663859569700E+02 C 0.611409934802217200E+02 0.268325289733025000E+04 0.271015097222632500E+05 C 0.344182439511496500E+06 0.117037498667061200E+07 0.264512776776344300E+07 C 0.633446917806616700E+07 0.205787510722953700E+08 C 17 0.46498754589668860000E-03 C 0.288517470048902000E-01 0.725364888134314100E+01 0.999131789083491600E+01 C 0.153030921256421000E+02 0.244482193358247200E+02 0.400115122032868100E+02 C 0.122877187849543700E+03 0.820855190798686500E+03 0.879480508898586200E+04 C 0.441511425024394600E+05 0.248753815647050200E+06 0.990837040650586800E+06 C -0.264421592112195600E+06 0.262172272279167900E+06 0.502053009085352100E+11 C -0.100926392651793000E+12 0.507198004150070300E+11 C 0.115636081824773900E+02 0.282211316665886800E+04 0.392830379714747600E+04 C 0.586261414284695500E+04 0.960221199316188700E+04 0.143387940303068500E+05 C 0.206931546242506100E+05 0.334707109352508900E+05 0.103328087252337900E+06 C 0.171203461893639500E+06 0.273018287506075900E+06 0.403228025347721000E+06 C 0.703742994384159800E+07 0.704446737378543500E+07 0.635086021579934400E+06 C 0.635721107601514000E+06 0.636356828709114700E+06 C -3CEDIL4RIO M4 2. 0.00 -2 3 C 11 0.28528038941022320000E+03 C 0.105748940794460900E+04 0.844409223042977800E+03 0.989216010327055800E+03 C 0.195811001056642700E+03 0.654649799725072200E+02 0.207639338793363700E+03 C 0.194266409628094900E+03 0.878586181627698400E+02 0.171400400143261400E+03 C 0.585371148745679700E+04 0.290267832276061300E+06 C 0.291079232024637700E+01 0.541078559821623500E+01 0.100447144771793100E+02 C 0.147071822048762200E+02 0.182300039492140500E+02 0.316863800567847800E+02 C 0.503710160254755600E+02 0.800538621073617100E+02 0.148612948128817500E+03 C 0.452161584069224300E+04 0.224665949971083600E+06 C 16 0.46054840037857230000E-03 C 0.169232871886761800E-01 0.286935787436301400E+01 0.103580120757909800E+02 C 0.155815266183155700E+02 0.242534497739954700E+02 0.278423880349024200E+02 C 0.282101800696905700E+02 0.274475140114087100E+02 0.657890370763595900E+03 C 0.118250445284277200E+05 0.251953429527934600E+05 0.114341679104464300E+07 C 0.731159039137183900E+07 0.440975123230174800E+09 0.951854529816820700E+10 C -0.996801324343811200E+10 C 0.702889834006048500E+01 0.118911520927119400E+04 0.425255626568891100E+04 C 0.638485518828208100E+04 0.973498555538559800E+04 0.123767428315422400E+05 C 0.117119756543999000E+05 0.109106508709576700E+05 0.339297808354646600E+05 C 0.229778262786752900E+06 0.319262580516430800E+06 0.179705932460679200E+07 C 0.294640166322813700E+07 0.449642996330517300E+07 0.437920000657641100E+07 C 0.438357920658298600E+07 C 0.58951122 -0.70710678 -0.41265641 C 0.00000000 0.00000000 0.00000000 C 0.55222554 0.00000000 0.81205257 C 0.00000000 0.00000000 0.00000000 C 0.58951122 0.70710678 -0.41265641 C 0.00000000 0.00000000 0.00000000 BLANK card ending "JMARTI SETUP" data cases BEGIN NEW DATA CASE C 2nd of 4 subcases involves an imbedded CABLE CONSTANTS data subcase. C This came from Dr. Ivano Bonfanti of CESI in Milano, Italy, as explained C in the April, 1994, newsletter. This is for part of a GIS substation C (note cable length is 90.1 meters). JMARTI SETUP 1. $ERASE CABLE CONSTANTS 2 1 1 0 1 0 0 0 C CARDS INDICATING NUMBER OF CONDUCTORS PER SC COAXIAL CABLE 2 C GEOMETRICAL AND PHYSICAL DATA CARDS .055 .07 .4 .415 .415 1.78E-8 1. 1. 1. 2.8E-8 1. 1. 1. C CROSS-SECTION LOCATION CARD 3.52 0. C EARTH RESISTIVITY AND FREQUENCY CARDS 300. 50. 90.1 { 1st of 3 is Xform matrix [T] 300. 50. 90.1 { 2nd of 3 is not actually used 300. 10. 6 10 90.1 { 3rd of 3 is frequency scan BLANK card ending frequency cards of imbedded "CABLE CONSTANTS" data case BLANK card ending "CABLE CONSTANTS" data cases 1 .48D-7 .30 30 0 1 3 0 0 .30 30 0 1 3 0 0 .0 C DEFAULT Note preceding cards result in printer plots showing the fitting. C If the user does not want or need such extra output, remove the C preceding 3 cards and replace by a DEFAULT (cols. 1-7) request. C <++++++> Cards punched by support routine on 11-Nov-18 11.00.00 <++++++> C -1 2. 1.00 -2 2 C 8 0.18076681444364581400E+03 C 0.395968340410608926E+04 0.295005854204842545E+05 0.174454160909646075E+06 C 0.869158190521383775E+06 0.364662556399659720E+07 0.145520745582613517E+08 C 0.531131969090856612E+08 0.227672588502988071E+09 C 0.203427156262552586E+03 0.162748596590855437E+04 0.102833044902314323E+05 C 0.546335287528216505E+05 0.245598170841564134E+06 0.104748764596227382E+07 C 0.412423792356411741E+07 0.191322161376811006E+08 C 13 0.30392899026088639700E-06 C 0.539904879139891847E+02 0.402547067289846439E+03 0.123039457481669524E+04 C 0.138908901013311129E+04 0.753649773022243608E+03 0.111542892895427417E+04 C 0.968760873289468691E+04 -0.126043771137276141E+06 0.459787149292387185E+06 C 0.103557188869358611E+07 0.432459583380759974E+07 -0.684191165642414429E+07 C 0.277530558943064958E+08 C 0.325587875143335441E+05 0.235697143057869719E+06 0.352631470503751655E+06 C 0.411189168490573532E+06 0.469188430536555418E+06 0.596470338226857014E+06 C 0.724813346723028459E+06 0.384515778316817433E+07 0.374877567556418246E+07 C 0.870763693318005093E+07 0.192604243620749750E+08 0.499116989461822033E+09 C 0.501590336460134462E+08 C -2 2. 1.00 -2 2 C 4 0.10451137020107520900E+03 C 0.895419253289177704E+02 0.237227186378166011E+03 0.826336359229880146E+03 C 0.107114138211958689E+05 C 0.114628603535565802E+03 0.307890218295652688E+03 0.108228738484974816E+04 C 0.140985532559441509E+05 C 1 0.29902748896168657700E-06 C 0.626888577669461132E+09 C 0.626889581814834476E+09 C 0.00127764 0.99872717 C 0.00307146 -0.00307629 C 0.99872717 -0.99951855 C -0.00307629 0.00119017 $PUNCH BLANK card ending "JMARTI SETUP" data cases BEGIN NEW DATA CASE C 3rd of 4 subcases involves an imbedded CABLE PARAMETERS data subcase. C Data came from Dr. Ivano Bonfanti of CESI in Milano, Italy, as explained C in the July, 1994, newsletter. This is for part of a GIS substation C (note cable length is 600 meters). There was a problem with stability. C Prof. Akihiro Ametani's new cable constants routine will be used for C variety (neither will stabilize the JMARTI result). JMARTI SETUP $ERASE CABLE CONSTANTS CABLE PARAMETERS { Transfer to new (August, 1994) cable constants code C MISCELLANEOUS DATA CARD 2 -1 1 0 1 0 0 0 C CARDS INDICATING NUMBER OF CONDUCTORS PER SC COAXIAL CABLE 2 C GEOMETRICAL AND PHYSICAL DATA CARDS .0206 .02865 .06395 .0689 .0775 1.775E-8 1. 1. 3.5 9.1E-8 1. 1. 2.25 C CROSS-SECTION LOCATION CARD 1.05 0. C EARTH RESISTIVITY AND FREQUENCY CARDS 300. 2.E+4 600. 300. 50. 600. 300. 10. 5 10 600. BLANK card ending frequency cards of imbedded "CABLE PARAMETERS" data case BLANK card ending "CABLE CONSTANTS" data cases 1 3.E-8 .30 30 0 1 1 0 0 .30 30 0 1 1 0 0 C -1 2. 0.00 -2 2 C 10 0.31527489296784590800E+02 C 0.333339478968597576E+03 0.385843242263242076E+03 0.156105194364770319E+04 C 0.265516564851181875E+04 0.744240773394244571E+04 0.146072675345095268E+05 C 0.100782791919086580E+06 0.423953540641733329E+06 0.181153979783380428E+07 C 0.426224308131755703E+07 C 0.127926506879103514E+03 0.191247861270285398E+03 0.140962601658519885E+04 C 0.256898278361350595E+04 0.719411752689923379E+04 0.147776521260217124E+05 C 0.523694703068609670E+05 0.231983450765360467E+06 0.103162566951607960E+07 C 0.268772287968215719E+07 C 14 0.19555725754330640800E-04 C 0.595812176849891984E+01 0.136598153216988170E+02 0.223062390699905855E+02 C 0.218765134842157032E+02 0.310078805673917373E+03 0.136948065507252182E+04 C 0.264828160326721173E+04 0.165371507469560675E+05 0.135218581976466987E+06 C 0.996294408964042087E+05 0.366762134130366892E+07 0.114366933510132088E+07 C 0.294681671932204888E+08 -0.345352346849985198E+08 C 0.225946305527752293E+04 0.654716563550326282E+04 0.858605165251129074E+04 C 0.647118553095189873E+04 0.258164726071371879E+05 0.633236066083839105E+05 C 0.500078550365742049E+05 0.136383733902013104E+06 0.312757035620522686E+06 C 0.431774484935509624E+06 0.845946787270792295E+06 0.146906560954104667E+07 C 0.100731432338616043E+07 0.100832163770954637E+07 C -2 2. 0.00 -2 2 C 25 0.25771619753860701300E+02 C 0.162475266380510348E+05 -0.344572775671657292E+05 0.221567870581859061E+05 C 0.548062829423196334E+02 0.294282733485373968E+03 0.114737235662244826E+04 C 0.111296958728743471E+04 0.616103629405914149E+03 0.593452262716258020E+03 C 0.410058929951589903E+03 0.560072870416973956E+03 0.116370288832337929E+04 C 0.140583167608246482E+04 0.160667392412559639E+04 0.126071408506419517E+04 C 0.135512678818804875E+04 0.154298926737988064E+04 0.200097285266074278E+04 C 0.243118215728368933E+04 0.389766818470841236E+04 0.543470105184340992E+04 C 0.869585056432428246E+04 0.145112660409044493E+05 0.286736334873477863E+05 C 0.911383809777398274E+05 C 0.128563209709157150E+03 0.145112581702906482E+03 0.148169399876459551E+03 C 0.169117065053743460E+03 0.356066420752889315E+03 0.431980837395433992E+04 C 0.409460631393748736E+04 0.424759484612483994E+04 0.650771124366744516E+04 C 0.864646715220207080E+04 0.117531733286592735E+05 0.230153259465291995E+05 C 0.279579746317147329E+05 0.327400282639653088E+05 0.253412598146896344E+05 C 0.290327207150341456E+05 0.343261562459948328E+05 0.423611078024432791E+05 C 0.516735710276618266E+05 0.810350212150682346E+05 0.115590605200089733E+06 C 0.182107115100442927E+06 0.304620180935359910E+06 0.604499969185806928E+06 C 0.192967270147408662E+07 C 7 0.37195367174999995500E-05 C 0.133780513819153044E+03 0.145518656151808569E+04 0.398591319635461878E+04 C 0.538269437655963885E+04 0.728753925137749230E+04 -0.330812244300116029E+06 C 0.341166613783013151E+08 C 0.313417235208839230E+05 0.337996604035161087E+06 0.908977622632142622E+06 C 0.123253835419277894E+07 0.169290700467920979E+07 0.617846954569379092E+08 C 0.346920308936563656E+08 C -0.00000404 1.00000404 C -0.00000626 0.00000000 C 1.00000404 -0.99998632 C 0.00000626 0.00002749 $PUNCH BLANK card ending "JMARTI SETUP" data cases BEGIN NEW DATA CASE C 4th of 4 subcases comes from Dr. Antonio Carlos Siqueira de Lima C of Gerência de Estudos Especiais, Proteção e Controle (GPE/DPP) C of ONS (Operador Nacional do Sistema Elétrico) in Rio de Janeiro, C Brazil. On September 24, 2001, he used the EEUG list server C to request help with some JMARTI SETUP fitting. He wrote: C "ATP fails to produce the JMARTI model. Either it stops with C an error message or enter an infinite loop." In her response on C October 12th, BPA's Dr. Tsu-huei Liu explained that an improvement C to the JMARTI fitter logic was required: "a similar problem would C be encountered when one runs JMARTI SETUP with the cable data of the C 9th subcase of DC-27 (the second pipe-type example). ... While C performing the curve fitting of either the characteristic impedance C or the propagation function (as function of frequency), ATP now will C skip the splitting of a segment to allocate poles and zeros if the C frequency range of the segment is too small ( f1/f2 < 1.0023 )." Of C course, there should be further explanation in some future newsletter. C Unusual roundoff was controlled by a change to MAIN23 on 10 October 2001. C Why such behavior was not seen with other data during years past is not C known. For some reason, the new data printed the JMARTI fitting table C that was one row short of the full ending decade. For example: C Table of Zc vs. ZCEQ for mode 1 follows: C Frequency ZCMAG ZCEQMAG DELMAG ZCPH ZCEQPH DELPH C in Hertz in Ohms in Ohms in percent degrees degrees degrees C 1.000000000E-01 3.557566852E+03 3.204331628E+03 -9.929124001E+00 -4.461019744E+01 -2.083729607E+01 2.377290137E+01 C 1.258925412E-01 3.170699522E+03 3.041719567E+03 -4.067870646E+00 -4.464329838E+01 -2.500977735E+01 1.963352103E+01 C 1.584893192E-01 2.825902051E+03 2.830476983E+03 1.618928111E-01 -4.467639932E+01 -2.937423126E+01 1.530216806E+01 C 1.995262315E-01 2.518599553E+03 2.576915901E+03 2.315427547E+00 -4.470950026E+01 -3.355375069E+01 1.115574958E+01 C 2.511886432E-01 2.244714642E+03 2.299542112E+03 2.442514030E+00 -4.474260120E+01 -3.709485066E+01 7.647750545E+00 C 3.162277660E-01 2.000613323E+03 2.024419985E+03 1.189968182E+00 -4.477570214E+01 -3.962787342E+01 5.147828724E+00 C 3.981071706E-01 1.783056784E+03 1.775021830E+03 -4.506280699E-01 -4.480880309E+01 -4.104536614E+01 3.763436947E+00 C 5.011872336E-01 1.589158413E+03 1.563720658E+03 -1.600706094E+00 -4.484190403E+01 -4.156692498E+01 3.274979049E+00 C 6.309573445E-01 1.416345506E+03 1.390377117E+03 -1.833478431E+00 -4.487500497E+01 -4.162790291E+01 3.247102053E+00 C 7.943282347E-01 1.262325125E+03 1.246766400E+03 -1.232544947E+00 -4.490810591E+01 -4.166770484E+01 3.240401064E+00 C 1.000000000E+00 1.125053678E+03 1.122474287E+03 -2.292681872E-01 -4.494120685E+01 -4.195192005E+01 2.989286801E+00 C ... Etc. through following decades until this termination: C 7.943282347E+01 1.264854975E+02 1.338400942E+02 5.814577075E+00 -4.333035447E+01 -2.514828779E+01 1.818206668E+01 C Plot of Zc in ohms vs. frequency in Hertz, for mode IMODE = 1. C This behavior was corrected by adding near-zero FLZERO to at least C one check for termination of a loop. $ERASE JMARTI SETUP C CABLE CONSTANTS { For use inside of JMARTI, this CC declaration is required C The preceding was required prior to improvement on 26 October 2001. Now C the CABLE CONSTANTS declaration is optional. Illustrate non-use here. C If enabled, the answer should be unchanged as prior to 26 October 2001. CABLE PARAMETERS C IT>IS-->NPC->IE-->KM-->IZF->IYF->NPP->NGR->ID->IYG-> 3 1 3 0 0 0 1 0 C Begin pipe parameter card C RP1---->RP2------>RP3------>RHOPIPE-->MUPIPE--->EPS1PIPE->EPS2PIPE-> C Following data missing input on "RP3": outer radius of pipe outer insulator C .0094641 .010 1.71800E-8 1. 2.85 1.0 .0094641 .010 .0151.71800E-8 1. 2.85 1.0 C Begin cross section location cards C D1----->THETA1--->D2------->THETA2--->D3------->THETA3---> .004 -90. .004 30. .004 -210. C Begin card indicating number of conductors C NP>NP2->NP3->NP4-> 1 1 1 C Begin geometrical & physical data cards C RIN1--->ROUT1---->RIN2----->ROUT2---->SEPP----->SEPG-----> C 0 .002 .002 003 0 .002 .002 { No more data needed for single conductor cable C RHOC--->MUC------>MUi1----->EPSi1----> 3.36500E-8 1. 1. 2.85 C RIN1--->ROUT1---->RIN2----->ROUT2---->SEPP----->SEPG-----> C 0 .002 .002 003 0 .002 .002 { No more data needed for single conductor cable C RHOC--->MUC------>MUi1----->EPSi1----> 3.36500E-8 1. 1. 2.85 C RIN1--->ROUT1---->RIN2----->ROUT2---->SEPP----->SEPG-----> C 0 .002 .002 003 0 .002 .002 { No more data needed for single conductor cable C RHOC--->MUC------>MUi1----->EPSi1----> 3.36500E-8 1. 1. 2.85 C C center->---------> 1.0 0.0 C Begin frequency & earth resistivty cards C RHO----------FREQ-----------IDEC-IPNT-DIST----IPUN C 100. 6.E+3 200. C 100. 60. 200. C 100. 1. 5 10 200. C Preceding 3 cards were the originals. Replace by the more-standard 3: 100. 5000. 200. 100. 60. 200. 100. .1 3 1 200. { Artificially small freq scan C 100. .001 8 10 200. { Nominal card needed to hang C To speed execution, 8 decades have been reduced to 3 and 10 points/decade C have been reduced to 1. But an old program will not hang using this choice. C To demonstrate an infinite loop for versions prior to October of 2001, C restore the original frequency card immediately above. WSM, 15 Oct. 2001. BLANK card ending frequency cards within CABLE PARAMETERS BLANK card ending CABLE PARAMETERS data subcases C DEFAULT C 1 3.E-8 C .30 60 0 1 1 0 0 C .30 60 0 1 1 0 0 1 { Default for GMODE is 2.E-9 mhos/km .30 30 0 1 1 0 0 .30 30 0 1 1 0 0 $PUNCH BLANK CARD ending JMARTI SETUP data subcases BEGIN NEW DATA CASE BLANK