Energia elektrikoa garraiatzeko aireko linea

Wikipedia(e)tik
Hona jo: nabigazioa, Bilatu
Aireko linea elektrikoa

Aireko linea elektrikoak energia elektrikoa garraiatzeko edo/eta banatzeko erabilitako egiturak dira. Euskarri batzuetatik eskegita dauden eroale elektrikoz hornituta egoten dira (eroale bat edo gehiagoz, askotan hiru edo lauz). Haietan isolamendu elektrikoa gehien bat airearen bitartez ziurtatzen denez, aireko linea elektrikoak energia elektrikoa kopuru handitan garraiatzeko egitura merkeenak izaten dira.

Lineen euskarriak zurez (zur naturalez edo laminatuz), altzairuz (sareta egiturez edo hodi moduko zurtoinez), hormigoiz, aluminioz eta, batzuetan, indartutako plastikoz egiten dira, Txikiak eta sinpleak direnean zutoin deitzen zaie eta handiak eta konplexuak direnean dorre.

Lineen eroale biluziak (isolamendu bereziaz estali gabeak) eskuarki aluminioz egiten dira, aluminio soilez edo altzairuzko edo beira- eta karbono-zuntzez egindako material konposatuzko erdigunez indartutako aluminioz. Hala ere, tentsio ertaineko banaketa-linea eta bezeroekiko behe-tentsioko lotura-linea batzuetako eroaleak oraindino kobrezkoak dira.

Energia garraiatzeko aireko linea elektrikoak diseinatzerakoan helburuen arteko garrantzitsuenetarikoak hurrengoak dira[1]:

  • linearekin ukitze arriskutsuak ekiditeko tentsiopean dauden eroaleen eta lurraren arteko distantziak egokiak izango direla ziurtatzea

eta

  • ekaitzek, izotz kargek, lurrikarek eta litezkeen bestelako kalte-eragileek sor ditzaketen indarrak gainditzeko eta horrelako egoeretan eroaleei eusteko gauza izango diren euskarriak erabiltzen direla ziurtatzea

Tentsioaren araberako sailkapena[aldatu | aldatu iturburu kodea]

Tentsio handi eta ertaineko linea elektrikoak Łomża-n, Polonian

Industria elektrikoak euren tentsioaren arabera sailkatzen ditu aireko linea elektrikoak, gehien erabilitako sailak hurrengoak dira:

  • Behe-tentsiokoak edo tentsio txikikoak (LV hizkiez identifikatuta ingelesezko dokumentu teknikoetan) – eroaleen arteko1000 V-eko baino tentsio txikiagokoak, bizitoki- edo merkataritza-bezero txikiak enpresa elektrikoarekin lotzeko erabiliak.
  • Tentsio ertainekoak (MV) – 1000 V (1 kV) eta 33 kV arteko tentsioekin, hirigune eta landa-gunetan energia banatzeko erabiliak, banaketa sarea osatuz.
  • Goi-tentsiokoak edo tentsio handikoak (HV) – 100 kV baino tentsio txikiagokoak, gehien bat, azpi-garraioan (energia kopuru dezenteak urrunegi ez dauden tokien arteko garraioan) erabiltzen dira, eta 115 kV eta 138 kV-eko tentsioen antzekoak azpi-garraioan, energia kopuru handixeagoak garraiatzeko elkarrengandik urrunegi ez dauden tokien artean eta, baita ere, bezero oso handiak sare elektrikora lotzeko.
  • Estra goi-tentsiokoak edo tentsio estra handikoak (EHV) – 230 kV eta 800 kV artekoak, energia kopuru oso handiak urrun dauden tokien artean garraiatzeko erabiliak, garraio sarea osatuz.
  • Ultra goi-tentsiokoak edo tentsio ultra handikoak (UHV) – 800 kV baino gehiagokoak.

Gaur egun toki askotan erabiltzen dira eroaleen arteko 765.000 volt-eko tentsiorainoko aireko linea elektrikoak, baina ez sartaldeko Europan, hemen, orokorrean, zentralen eta kontsumo-guneen arteko distantziak beste toki batzuetan baino laburragoak direlako.

Egiturak[aldatu | aldatu iturburu kodea]

Sakontzeko, irakurri: «Energia garraiatzeko linea elektrikoen euskarri»

Aireko linea elektrikoen egiturak mota askotakoak izan daitezke, linea motaren arabera. Egiturak lurrean zuzenean jarritako zurezko zutoinez egindakoak bezain sinpleak izan daitezke; batzuetan eroaleei eusteko zeharraga bat edo gehiagorekin eta beste batzuetan “besobakoak”, zutoinari zuzenean atxikitako isolagailutan lotutako eroalerekin. Hirigunetan, batzuetan, altzairuz egindako hodi itxurako euskarriak erabiltzen dira. Tentsio handiko lineen eroaleak sarritan dorreak deitutako altzairuzko sareta-euskarrietan jarri ohi dira. Batzuetan, urrun dauden tokietan, helikopteroz garraiatutako aluminiozko euskarriak erabiltzen dira. Betoiz egindako zurtoinak ere erabiltzen dira.[1] Indartutako plastikoz egindako zurtoinak ere badaude erabilgarri, baina haien kostu handiagatik gutxitan erabiltzen dira.

Egiturak diseinatzerakoan eroaleek sortutako indar guztiak kontuan izan behar dira.[1] Eroaleen pisuaz gain, haizeaz, izotz metaketaren ondorioz eta bibrazioez sortutako gainkargak ere kontuan hartu behar dira. Euren muturretatik eskegita dauden eroale malguek katenaria baten itxura hartzen dute, eta aireko linea elektrikoen diseinua hein handi batean forma horren propietateetan oinarritzen da.[1] Linea elektrikoek lerro zuzenak jarraitzen dutenean, eroaleen tentsio mekanikoak orekatuta egoteagatik dorreek normalean ez dute jasaten eroaleen tentsio mekanikoengatiko luzerako indar handiak eta arinagoak izan daitezke; normalean eroaleen (eta haiei atxikita egon daitekeen izotzaren) pisuagatik eta haien aurka haizea jotzeagatik sortutako indarrak jasan behar dute gehien bat eta.

Energia garraiatzeko aireko linea handien proiektuetan zenbait motako euskarri edo dorreak erabiltzen dira:

  • lerrokadura- edo esekidura-euskarriak, lineek bide zuzenak jarraitzen dituztenean esekitako isolagailuekin erabili ahal diren egitura arinekoak;
  • angelu-euskarriak, lineek norabidea aldatzen dutenean euskarriaren albo bietan isolagailuak ainguratuta erabiltzekoak;
  • aingura-euskarriak edo muturreko euskarriak, sendoenak izan behar direnak, eroaleak euskarriaren albo batean bakarrik ainguratuta egoteagatik luzerako indar handienak jasan behar dituztelako;
  • tarte luzeak gainditzeko euskarri bereziak eta abar... .

Linea batzuk, euskarrien albo bietako eroaleen tentsio mekanikoen orekaz baliatuz, nahiko malguak diren euskarriz diseinatzen dira. Egitura zurrunagoak erabili behar dira eroale bat edo gehiago apurtzen bada ere, zutunik manten daitezen nahi bada. Aireko linea askotan tarteka horrelako egitura zurrunagoko euskarriak jartzen dira, eroale bat edo gehiago apurtu izanagatik bata bestearen ondoan jausitako euskarri arrunten kopurua mugatzeko.[1]

Euskarrien oinarriak handiak eta garestiak izan daitezke, gehien bat lurraren ezaugarriak egokiak ez direnean, hezegune askotan legez. Egiturak egonkortu daitezke nahikotxo eroaleek sortutako indarrei aurre egiten dieten lokarriak erabiliz.

Energia garraiatzeko linea elektrikoek eta haien euskarriek ikus-kutsadura sor dezakete. Hori ekiditeko batzuetan lineak lurperatzen dira baina lurperatze hori oso garestia da eta, beraz, ez da sarritan egiten.

Aga bakarreko zurezko euskarri baterako, zurezko aga bat jartzen da lurrean eta hari itsatsita zeharraga batzuk, isolagailuak zeharragetan jarriz. "H" motako zurezko euskarri baterako, zurezko bi aga jartzen dira lurrean eta haiek baino harantzago luzatzen den haien arteko zeharraga luze bat, isolagailuak zeharraga horren muturretan eta erdian jartzen direla normalean. Altzairuzko sareta-dorreen egiturak orokorrean bi formakoak izaten dira: Isolagailu kateak piramide-enbor itxurako oinarri baten gaineko zati zuzen batetik zabaltzen diren hiru edo sei beso antzeko hegal-habeetara lotuta dituztenek eta kate horiek piramide-enbor itxurako oinarri gainean jarrita dauden albo azpiegitura biren gainean kokatzen den habe batera lotuta dituztenek. Hala ere , beste itxurak dituzten beste egitura batzuk ere erabiltzen dira.

Batzuetan lurrera jarritako “babes-eroale” edo “lurreko eroale” bat edo gehiago jartzen da euskarrien goi partean eroale aktiboei tximisten aurkako babesa emateko. Gaur egun, noizbait, barnean zuntz optikoak dituzten babes eroale optikoak erabiltzen dira eta eroale aktiboak babesteko eta, aldi berean, komunikazioak bideratzeko balio dute.

Zirkuituak[aldatu | aldatu iturburu kodea]

Zirkuitu bakarreko garraio linea batek zirkuitu bakar bat osatzeko behar diren eroaleak ditu. Beraz, hiru faseko sistema bateko lineetan euskarriek, gutxienez, isolatutako hiru eroaleri eutsi behar diete.

Zirkuitu bikoitzeko garraio linea batek zirkuitu bi ditu, eta hiru faseko sistemetan euskarriek, gutxienez, isolatutako sei eroaleri eutsi behar diete.

Lineak tentsio handikoak izateagatik faseetan eroale-sortak erabiltzen ba dira, euskarriek isolatutako eroale gehiagori eusten diete. Bestalde, Lineek zirkuituak babesteko babes-eroaleak badituzte, isolaturik ez dauden beste eroale batzuei ere eutsi behar diete, normalean bati edo biri.

Trakziotarako erabiltzen diren fase bakarreko korronte alternoko garraio lineek lau eroale izaten dituzte, bi zirkuitutarako. Normalean, zirkuitu biak tentsio berekoak izaten dira.

Tentsio handiko korronte zuzeneko sistemetan (HVDC literatura tekniko ingelesean, high voltage direct current) normalean linea bakoitzak eroale bi izaten ditu, baina inoiz edo behin polo bakar bateko eroalea baino ez da jartzen euskarri batzuetan.

Toki batzuetan, Alemanian adibidez, 100 kV baino tentsio handiagoko garraio-linea gehienak bi, lau edo, kasu bakan batzuetan, sei zirkuitukoak izaten dira, igarotze-eskubideak bakanak direlako. Batzuetan zirkuitu guztiak jartzen dira egiturak jartzen direnean baina beste batzuetan zirkuitu batzuk jarri barik uzten dira, geroago, beharrezkoak izango direnean, jartzeko. Zirkuitu bat baino gehiagoko lineen desabantaila nagusia mantentze-lana zailtzen dela da, tentsiopean dagoen zirkuitu baten ondoan lan egitea edo zirkuitu biak batera deskonektatzea beharrezkoa delako. Bestalde, akats baten ondorioz bi zirkuituak gal daitezke.

Isolagailuak eta isolagailu-kateak[aldatu | aldatu iturburu kodea]

Tentsio ertaineko lineak zeramikazko isolagailuekin Kalifornian
Tentsio handiko lineetan eskegitako isolagailu-kateak erabiltzen dira.

Isolagailuek eroaleei eutsi eta ohiko tentsioa eta tximisten eta maniobren ondorioz sortutako gaintentsioak jasan behar dituzte. Isolagailuak, orokorki, mota bikoak izaten dira: eroaleei euskarriaren egituraren gainetik eusten dieten isolagailu zurrunak eta, eskegita, eroaleei euskarriaren egituraren azpitik eusten dieten esekidura-isolagailuak; goi-tentsioko lineetan erabiltzen diren isolagailu-kateak bigarren motako isolagailu batzuk bata bestearen atzean jarriz osatzen dira. Isolagailu-kateen asmakuntza berebiziko garrantzia izan zuen tentsio handiak erabili ahal izateko. Kate horiek osatzeko erabiltzen diren isolagailuak normalean beira-tenplatuzko edo portzelanazko disko antzeko bati zementuz goitik buloi baten burua sartzeko barrunbea duen txano metaliko bat eta behetik buloi burudun metaliko bat lotuz egiten dira.

Lau isolagailuz egindako isolagailu-katea

Hemeretzigarren mendeko amaieran telegrafo lineetan erabiltzen ziren isolagailu zurrunen antzeko isolagailuek ezin dute 69.000 volt baino gehiagoko tentsioak jasan eta, horren ondorioz,. 33 kV arte (69 kV Estatu Batuetan) bi motakoak erabiltzen badira ere[1], tentsio horretatik gora eskuarki isolagailu-kateak baino ez dira erabiltzen eroaleei eusteko.

Isolagailuak normalean portzelanaz edo beira tenplatuz egiten dira, beira-zuntzez indartutako polimeroez egindakoen kopurua handituz badoa ere; garestiagoak izan arren erabiltzen diren tentsio gero eta handiagoekin eta klima baldintza gogorretan erabiltzeko ezaugarri egokiagoak izateagatik. [2] Txinak dagoeneko garatu ditu 1.100 kV-rainoko lan-tentsioak jasan ahal dituzten isolagailu polimerikoak eta Indiak 1.200 kV-rainoko lan-tentsioak jasateko gauza izango direnak garatzeko bidean dago.

Isolagailu-kateak isolagailu batzuk beste batzuei lotuz egiten dira; tentsio mailarekin erabil behar den isolagailu kopurua handitzen dela. Erabiliko den isolagailu kopurua linearen lan-tentsioaren eta jasan beharreko lainoa, kutsadura edo gatz zipriztinduen antzeko ingurumen faktoreetara lotuta dago. Egoerak ezegokiak direnetan ihes-korronterako bide luzeagoko kate luzeagoak erabili behar dira. Isolagailuak eta isolagailu-kateak eroaleen baoko pisuak eta haizearen eta batutako izotzaren ondorioz sortutako indarrei aurre egiteko bezain sendoak izan behar dira.[3]

Portzelanazko isolagailu batzuek erdieroale den akabera esmalteztatua izaten dute eta korronte txiki bat (miliampere gutxi batzuk) zeharkatzen du haien azala. Horrek haien gainazala apur bat berotu eta lainoaren eta hauts metaketaren ondorio txar batzuk txikitzen laguntzen du. Erdieroalea den geruza horrek isolagailu-katearen isolagailuen artean tentsioa modu egokiagoan banatzeko ere balio du.

Isolagailu polimerikoak, naturaz, hidrofoboak dira, eta giroa hezea denean portaera egokiagoa dute. Gainera, ikerketa batzuek erakutsi dute isolagailu polimerikoetan ihes-korronteak ekiditeko beharrezkoa den ihes-bidearen luzera beirazkoetan edo portzelanazkoetan baino txikiagoa dela. Bestalde, isolagailu polimerikoen masa beirazko edo portzelanazko kateena baino %30-tik %50-era txikiagoa da (bereziki tentsio handietan). Kutsadura eta giro hezearen aurreko euren portaera hobea garrantzi handikoa da haien erabilera gero eta handiagorako.

Oso tentsio handitarako isolagailu-kateak, 200 kV baino handiagotarakoak, eraztun ekipotentzialez hornituta egon daitezke. Horrela eremu elektrikoaren banaketa isolagailuen inguruan hobetzen da eta gaintentsioak sortzen direnean deskarga disruptiboak agertzea zailtzen du.

Eroaleak[aldatu | aldatu iturburu kodea]

ACSR motako eroale bat

Altzairuz indartutako aluminiozko eroaleak (ACSR ingelesezko literatura teknikoan, aluminum-conductor steel-reinforced ) gehien bat tentsio ertain eta handiko lineetan erabiltzen dira baina banakako bezeroentzako lineetan ere erabiltzen dira. Erresistentzia berbera duten kobrezko eroaleen pisu erdia inguru dutelako eta merkeagoak direlako erabiltzen dira aluminiozko eroaleak. Diametroa, ordea, kobrezkoen baino handiagoa izan behar dute, aluminioak konduktibitate txikiagoa izateagatik.[1] Hala ere, oraindino erabiltzen dira kobrezko eroaleak kasu batzuetan, batez ere behe-tentsioekin eta lurreratzetan. Beste motatako eroaleak ere, erabiltzen dira; hala nola AAC (all-aluminum conductor), AAAC (all-aluminum-alloy conductor) eta ACAR (aluminum conductor alloy reinforced) deituraz izendatuak.

Eroale handiagorekin, haien erresistentzia elektrikoa txikiago izanda, energia gutxiago galtzen da baina garestiagoak dira. Kelvin-en legea deitutako optimizaziorako arauaren arabera, linea bateko eroale tamaina optimoa da urtean eroalean galdutako energiaren kostua eta eroalea jartzeko egin beharreko inbertsioaren urteko korrituak berdinak direnean. Hala ere, optimizazioaren kontua asko zailtzen da urtean zeharreko lineen erabilera aldakorra, jartze kostu aldakorra, eroaleen tamaina diskretua eta horren antzeko beste faktore batzuk kontuan hartuz gero.[1]

Eroalea luzera unitateko pisu konstantea duen objektua malgua izanda, haren muturretik eskegita dagoenean hartzen duen forma geometrikoa katenaria batekoa da. Eroalearen gezia (eroalea euskarrietan lotuta dagoen lotunetatik igarotzen den lerro zuzena eta zuzen horri paraleloa den katenariaren zuzen ukitzailearen arteko distantzia bertikala edo, beste definizio ez bateragarriaren arabera, katenariaren puntu garaien eta baxuenen arteko distantzia bertikala) tenperaturarekin eta izotz geruzarenaren moduko beste karga gehigarrirekin aldatzen da. Segurtasun ziorengatik eroaleen eta lurraren arteko distantzia beti izan behar da linearen tentsioaren funtzioa den minimo bat baino handiagoa. Eroaleen tenperatura eroaletik igarotzen den korronteak sortutako beroarekin igotzen denez, batzuetan energia garraiatzeko linearen ahalmena handitu ahal da dilatazio-koefiziente termiko txikiagoa edo funtzionamendurako tenperatura onargarri handiagoa duten eroaleak aukeratuz.

Ohiko ACSR eroale bat eta ACCC berri bat

Gezi txikiagoak izaten dituzten mota horretako eroale batzuk ingelesezko literatura teknikoan ACCC akronimoaz ezagutuak dira (aluminum conductor composite core). Eroaleak indartzeko ACSR-etan erabiltzen diren altzairuzko harien ordez ACCC-etan erabilitako beirazko edo karbonozko zuntzez egindako erdigune indargarriekin eroaleen beroagatiko hedapen-koefizientea 10 bider txikiagoa da. Erdigune konposatua eroalea ez bada ere, altzairua baino arinago eta sendoagoa da eta horrek eroaleetan %28 aluminio gehiago jartzea ahalbidetzen du, haien diametroa edo pisua handitu barik (trapezio itxurako hari konpaktuak erabiliz). Gehitutako aluminio kopuruak, diametro eta pisu bereko beste eroale motekin konparatuta, linearen-energia galerak %25etik %40ra gutxitzera laguntzen du, linearen intentsitatearen arabera. ACSR edo AAC motakoekin konparatuta, ACCC motako eroaleen beroagatiko gezi txikiagoek linearen intentsitatea, gehitzea ahalbidetzen dute eta, horren ondorioz, energia garraiatzeko ahalmena ere.

Aireko linea elektrikoek askotan izaten dituzte haien eroaleei lotutako gailu batzuk; esate baterako, sorten eroaleak bereizita mantentzeko bereizgailuak, zio ezberdinengatik eroaleetan ager daitezkeen bibrazioak ekiditeko edo gutxitzeko bibrazioen aurkako motelgailuak, toki eta egoera berezietan isolagailu-kateak euskarriei larregi hurbil ez daitezen pisuak, hegaztiak eroaleez konturatu daitezen ohargailuak edo Nazioarteko Abiazio Zibilaren Erakundearen aholkuei[4] jarraituz jartzen diren hegazkin eta helikopteroentzako ohargailuak.

Eroale-sortak[aldatu | aldatu iturburu kodea]

Bereizgailu bat eroale-sorta batean

Tentsio oso handiekin eroale-sortak erabiltzen dira, Koroa-efektuaren ondorioz izaten diren galerak txikiagotzeko. Eroale-sortak tarteka bereizgailuz konektatuta dauden eroale paraleloren multzoak dira, sarritan konfigurazio zilindrikotan kokatuta. Erabilitako eroaleen tamainaren arabera, 220 kV lineetan ez dira erabiltzen eroale-sortak edo bi eroaleko sortak erabiltzen dira normalean, 380 kV lineetan normalean bi, hiru edo lau eroaleko sortak erabiltzen dira. 765 kV-eko lineetan ohikoak dira sei eroaleko sortak. Bereizgailuak, sorta bateko eroaleak geometrikoki bananduta baina elektrikoki lotuta mantentzen dituzten osagaiak, haizeak sortzen dituen indarrei eta zirkuitulaburrak gertatzen direnean sortzen diren indar magnetikoei aurre egiteko gauza izan behar dira.

Tentsio estra handiekin eremu elektrikoaren gradientea eroale bakar baten azalean airea ionizatzeko bezain handia da eta, horren ondorioz, energia galtzen da zarata entzungarria eta interferentzia elektromagnetikoak sortuz. Eroale-sorta baten inguruan sortzen den eremu elektrikoa askoz handiagoa izango zen eroale bakar baten inguruan sortutakoaren antzekoa da; beraz, gradienteak txikiagoak dira eta indar handiko eremuez agertzen diren arazoak txikiagotu ahal dira.

Gainera, eroale-sortak erabiliz azal-efektuaren ondorioz korronte alternoarekin gertatzen den eroale handien garraiatzeko ahalmenaren txikiagotzea ekidin ahal da. Bestalde, azalera handiagoa izateagatik, sortetako eroaleak eraginkortasun handiagoarekin hozten dira eta, eroaleen erresistentzia tenperatura baxuagorekin txikiagoa denez, galerak gehiago txikiagotzen dira.

Koroa galerak eta Azal-efektuaren ondorio kaltegarriak txikituz gain, eroale-sortak erabiliz lineen induktantzia ere txikiagoa da, eta linea induktantzia txikiak guztiz desiragarriak dira, induktantziak txikiagotuz erreaktantzia induktiboak ere txikiagotzen direlako, eta eurekin batera tentsio jauskerak linearen zehar ere.

a, b, c... n eta a', b', c'... m' eroalez osatutako joaneko eta itzulerarako A eta B bi eroale-sortaz egindako zirkuitu monofasiko baten induktantzia kalkulatzeko garrantzi handia dute hurrengo bi kontzeptuek:

  • Batezbesteko geometrikoa erabiliz kalkulatutako eroale-sorta bien arteko batez besteko distantzia; hemen batezbesteko geometrikozko distantzia edo, espainieraren ohiturei jarraituz, batez besteko distantzia geometriko deituko eta BGD edo GMD akronimoez ordezkatuko dena (GMD ingelesezko literaturan erabilitako akronimoa da, “geometric mean distance”-tik sortua) eta
  • Batezbesteko geometrikoa erabiliz kalkulatutako eroale-sorta baten eroaleen arteko batez besteko distantzia; hemen eroale-sortaren batezbesteko geometrikozko erradio edo eroale-sortaren batez besteko erradio geometriko deituko eta BGE edo GMR akronimoez ordezkatuko dena (GMR ingelesezko literaturan erabilitako akronimoa da, “geometric mean radius”-etik sortua).

Bi eroale-sortaren arteko GMD-a kalkulatzeko lehenengo eroale-sortaren eroale guztietatik bigarren eroale-sortaren eroale guztietara dauden distantzia guztien batezbesteko geometrikoa lortu behar da, hau da:

GMD_{AB} = \sqrt[nm']{(D_{aa'}D_{ab'}...D_{am'})(D_{ba'}D_{bb'}...D_{bm'})...(D_{na'}D_{nb'}...D_{nm'})}

Eta eroale-sorten GMR-a kalkulatzeko eroale-sorta bakoitzaren eroale guztietatik eroale-sorta berberaren eroale guztietara dauden distantzia guztien batezbesteko geometrikoa:

GMR_{A} = \sqrt[nn]{(D_{aa}D_{ab}...D_{an})(D_{ba}D_{bb}...D_{bn})...(D_{na}D_{nb}...D_{nn})}

GMR_{B} = \sqrt[m'm']{(D_{a'a'}D_{a'b'}...D_{a'm'})(D_{b'a'}D_{b'b'}...D_{b'm'})...(D_{m'a'}D_{m'b'}...D_{m'm'})}

non eroale batetik berbererako distantzia eroalearen beraren GMR-a den. RE erradioko E eroale zilindriko eta homogeneo batentzat GMR hori hurrengoa da:

D_{EE} = GMR_{E} = R_{E}.e^{-{1 \over 4}}

e logaritmo nepertarren oinarria dela; baina beste eroale mota batzuentzat normalean eroaleen fabrikatzaileek argitaratutako tauletatik ateratzen da; eroaleen osaketa zehatza eta dagozkion ezaugarri induktiboen araberakoa delako, eta horiek analitikoki zehaztea erraza ez delako, batez ere eroaleak konposatuak direnean. GMRE-ren balioak normalean 6 mm eta 18 mm artean egoten dira.[5]

A eta B eroaleen induktantziak hurrengoak dira:

L_A= {\mu _o \over {2 \pi }}\ln {GMD_{AB}\over GMR_A} = 2 .10^{-7}\ln {GMD_{AB}\over GMR_A} [6][7]          H/m

eta

L_B= {\mu _o \over {2 \pi }}\ln {GMD_{AB}\over GMR_B} = 2 .10^{-7}\ln {GMD_{AB}\over GMR_B}           H/m

eta zirkuitu osoarena:

L_{ZIR}=L_{A}+L_{B}

Ikusten denez, eroale-sortez egindako zirkuituen induktantziak kalkulatzeko, erabilitako eroaleen ezaugarriez gain sorta osatzen duten eroaleen kopurua eta haien kokapenak sortan ere, kontuan hartu behar dira; eroale-sortaren GMR-a kalkulatzeko beharrezkoak direlako.[8]

bi eroaleko sorta
berdin banandutako hiru eroaleko sorta
berdin banandutako lau eroaleko sorta

GMR-en kalkulua errazten da eroale-sorta bakoitzaren eroale guztiak berdinak direnean eta era uniformean banatzen direnean zirkunferentzia batean. Kasu horretan, eroaleen kopuruaren arabera eta haien arteko DB distantzia eta GMRE batezbesteko geometrikozko erradioaren funtzioan eroale-sorten batezbesteko geometrikozko erradioak hurrengoak dira:

Variable labels for bundled conductor calculations

eroale bakar baten GMR-a:      GMR_1 = GMR_E

2 eroaleko sortaren GMR-a:      GMR_2 = \sqrt{GMR_E.D_B}

3 eroaleko sortaren GMR-a:      GMR_3 = \sqrt[3]{GMR_E . (D_B)^2}

4 eroaleko sortaren GMR-a:      GMR_4 = \sqrt[4]{GMR_E . (D_B)^2.D_B\sqrt{2}}[9]

n eroaleko sortaren GMR-a:      GMR_n = { D_B \over{2.sin ( { \pi \over{n}  } ) } }.\sqrt[n]{n.GMR_E.2.sin( { \pi \over{n}  } )\over {D_B}}          n bat baino gehiago denean. [10]

Normalean zirkuitu monofasiko baten joaneko eta itzulerako eroale-sorten arteko distantziak eroale-sorta bakoitzeko eroaleen artekoak baino askoz handiagoak dira, eta, horren ondorioz, sorten arteko GMD-a bi sorten erdiguneen arteko D distantziaz ordeztu ahal da. Bestalde, joaneko eta itzulerako eroale-sortak berdinak izaten dira. Kasu horretan zirkuitu monofasikoaren induktantzia hurrengoa da:

L_{ZIR} = 2 .10^{-7}\ln { D \over GMR} + 2 .10^{-7}\ln { D \over GMR} =  4 .10^{-7}\ln { D \over GMR}           H/m


eta erreaktantzia:


X_L = 2 .\pi .f.L_{ZIR}


50 Hz-ko sistemetan


X_L = 1,256637.10^{-4}\ln { D \over GMR}           Ω/m


Sistema polifasikoetan (trifasikoak barne) fase bakoitzeko beste faseekiko elkar-induktantziak esplizituki kontuan hartzen ez badira, aurrerago aurkezten den konputagailu-programan egiten den antzera, ezin da aztertu era egokian zirkuituen funtzionamendua; hala ere kasu bakan batzuetan lortu ahal da, zehaztasun gutxiagorekin bada ere, funtzionamendua aztertzeko gutxi gorabeherako parametroak erabiliz. Esate baterako linea trifasiko baten A, B eta C bere faseekin transposizio perfektua egiten dela suposatzen bada eta haren faseen arteko distantziak DAB, DBC eta DCA eta haien batezbesteko erradio geometrikoa GMR badira, fase bakoitzeko induktantzia hurrengoa dela suposatu ahal da[11]:


L_F = 2.10^{-7}\ln {\sqrt[3] {D_{AB}.D_{BC}.D_{CA}}\over GMR }           H/m


Hurrengo bi irudietan eroale-sorten tamainen eta eroale-sorten arteko distantzien eraginak induktantzietan agertzen dira.


Desabantailak

Eroale-sortak erabiliz haizeak eragindako jasan beharreko indarrak handiagoak dira, hala ere, elektrikoki haien baliokidea izango zen eroale bakarrarekin baino txikiagoak.

Bestalde, eroale sortak jartzea eroale bakunak jartzea baino zailago eta garestiago da.

Babes-eroaleak[aldatu | aldatu iturburu kodea]

Aireko linea elektrikoek sarritan babes-eroaleak dituzte. Babes-eroaleak tximista batek fase-eroale bat zuzenean jotzeko probabilitatea txikitzeko normalean euskarrien gainean lurrera lotuta egoten diren eroaleak dira. Babes-eroaleak, gainera, neutroa lurrera lotuta duten sistemetan akats-korronteak bideratzeko balio duten lurrarekin paraleloan jarritako bideak dira. Tentsio oso handiko lineek babes-eroale bi izan ditzakete., goi partean eta bidearen erdigunetik aldenduta, euskarrien egituretan haiek lotzeko espresuki jarritako zatitan. Lineak tximisten aurka babestuz azpieztazioetako tresneriaren diseinua errazten da, haren isolamenduak tentsio txikiagoak jasan behar izango ditu eta.

Garraiotarako lineen babes-eroaleak komunikazioa eta energia-sistemaren kontrola ahalbidetzeko ere erabili ahal dira, ingelesez OPGW akronimoz (optical ground wire-tik eratorria) izendatutako zuntz optikoz hornituta dauden babes-eroaleak erabiliz. Lehengo Sobiet Batasuneko tentsio oso handiko aireko linea batzuetan, babes-eroaleak PLCko sistemekin (Power Line Communication / Komunikazioa energia-lineen bidez) erabiltzen ziren eta isolagailu bidez lotzen ziren euskarrietara.

HVDC-ko bihurgailu-estazio batzuetan, babes-eroalea urruneko lur-elektrodo batera konektatzeko ere erabiltzen da. Honek HVDC sisteman lurra eroale legez erabiltzea ahalbidetzen du. Babes-eroalea tximistentzako zeharbidea errazten duten gailuz hornitutako isolagailu txikitan jartzen da Isolamendu horrek euskarrien korrosio elektrokimikoa eragozten du.

Tentsio ertaineko banaketa lineek ere ukan ditzakete babes-eroale bat edo gehiago, edo eroale aktiboen azpian kokatutako eta lurrera lotutako eroale bat, ibilgailu edo makina handiek eroale aktiboak ukitzea oztopatzeko edo/eta eroale neutro legez erabiltzeko izar moduan konektatutako sistemetan.

Isolatutako eroaleak[aldatu | aldatu iturburu kodea]

Aireko lineetan gutxitan erabiltzen dira kable berez isolatuak goi-tentsiotan, inoiz distantzia laburretarako (km bat baino laburragoak). Berez isolatutako eroaleak (atal honetan isolatutako eroaleak deituta) zuzenean lotu daitezke euskarrien egituretara, inolako isolagailuren beharrik gabe. Eroale biluzizko airez isolatutako lineak normalean isolatutako eroalez egindakoak baino merkeagoak izaten dira.

Arruntagoa da “estalitako” eroaleak erabiltzea. Eroale biluziak legez erabiltzen dira baina basa-animalientzat seguruagoak dira, eroaleak ukitzeagatik hego luzedun hegaztiak elektrokutatuta suertatzearen probabilitatea txikiagotzen da eta lineen arrisku orokorra apur bat txikitzen da. Horrelako lineak sarritan ikusten dira baso sarriak dauden tokietan, bertan zuhaitzek linea ukitzea probabilitate handiko jazoera izateagatik. Desabantaila bakarra kostua da, eroale isolatuak biluziak baino garestiagoak dira eta. Konpainia elektriko askok eroaleak elkarrengandik edo beste egituretatik hurbilago dauden tokietan erabiltzen dituzte estaliko eroaleak, esate baterako, lurpeko linea batetik aireko linea batera igarotzerakoan, ebakigailuen ondoan eta antzeko egoeratan.

Dena den, gaur egun behe-tentsioko aireko linea asko isolatutako eroalez egindako txirikordatutako kablez egiten dira.

Garraio-linea konpaktuak[aldatu | aldatu iturburu kodea]

Garraiotarako aireko linea konpaktu batek aireko linea arrunt batek baino igarobide txikiagoa behar du, eroaleak elkarrengandik hurbilago daudelako. Eroaleak elkarrengana larregi hurbildu ez daitezen euskarriak elkarrengandik hurbilago jarri, gurutzeta isolatzaileak erabili edo euskarrien artean eroaleak isolagailuz aldendu ahal dira. Lehenengo motakoak eraikitzeko eta mantentzeko errazagoak dira, eroaleen arteko isolagailuak jartzea eta mantentzea erraza ez delako. Linea konpakturen adibideak:

Batzuetan garraio linea konpaktuak sor daitezke martxan dauden tentsio bateko batzuk beste tentsio handiagoz erabiltzeko birdiseinatuz, horrela baimendutako igarobideaz garraiatu ahal den energia kopurua handituz.[12]

Behe-tentsioa[aldatu | aldatu iturburu kodea]

Behe-tentsioko linea bat Old Coulsdon-en, Surrey-n

Behe-tentsioko aireko lineek eroale biluziak beirazko edo portzelanazko isolagailu gainean edo isolatutako eroalez txirikordatutako kableak izan ditzakete. Zirkuitu bakarreko linea trifasikoen eroale-kopurua lautik (hiru faseenak eta neutroa eta lurraren betebeharra betetzen duen laugarren bat) seirainokoa (hiru faseenak, banandutako neutroaren eta lurrarenak gehi argiztapena elikatzeko beste bat) izan daiteke.

Trenentzako energia[aldatu | aldatu iturburu kodea]

Artikulu nagusia: Aireko lineak

Aireko lineak edo aireko hariak erabiltzen dira trolebusei, tranbiei eta trenei energia helarazteko ere. Tranbien eta trenentzako aireko lineak eroaleak trenbideko errailen gainean jarrita egon daitezen diseinatzen dira, ezelango oztopo barik eroale horiek eta errailen artean. Trakzioko aireko lineak elikatzeko goi-tentsioko saretik hartzen da energia, tarte erregularretan tartekatzen diren azpiestazioen bidez. Kasu batzuetan maiztasun txikiko korronte alternoa (KA) erabiltzen da eta trakzio elektrikorako sare berezi baten bidez banatzen.

Energia garraiatzeko aireko lineen ingurua[aldatu | aldatu iturburu kodea]

Aireko linea baten azpiko eremuaren erabilera mugatua dago, ezer ez delako gehiegi urreratu behar tentsiopean dauden eroaleetara. Eroaleen eta euskarrien inguruan metatutako izotzak elkartzen direnean ere, arriskua sortzen da. Irrati-uhinak hartzeko arazoak sor daitezke linea elektrikoen ondoan, bai lineak jatorrizko seinaleak oztopatzeagatik bai eroaleetan, isolagailuetan eta beste gailu laguntzaileen gune zorrotzetan gertatzen diren deskarga partzialek sortutako irrati-uhinengatik.

Aireko lineen inguruan haien eroaleekin kontaktuan jartzea ahalbidetzen duten jarduerak arriskutsuak bihurtzen dira; ahal dela, haien inguruan ez da komeni, esate baterako, kometak edo globoak hegaldaraztea, eskailerak erabiltzea edo makineriaz jardutea.

Aireportuen ondoko aireko garraio- eta banaketa-lineak sarritan irudikatzen dira mapetan, eta lineetan euretan plastikozko islagailu ikusgarriak edo/eta argiak jartzen dira, jakinarazteko pilotuei inguruan eroaleak daudela.

Energia garraiatzeko aireko lineen eraikuntzak nabarmenki eragin dezake ingurumenean, batez ere eremu basatietan. Horrelako proiektuen ingurumen-azterketetan batzuetan kontuan hartzen dira zuhaitzak eta zuhaixkak moztearen efektuak, animalia migratzaileen migrazio-ibilbideen aldaketak, gizakien eta animalien balizko sarrerak korridoreetan, ibaiak igarotzean arrainen habitatetan sortutako aldaketak eta antzerako beste efektu batzuk.

Historia[aldatu | aldatu iturburu kodea]

Pultsu elektrikoen lehenengo transmisioa distantzia gogoangarrian Stephen Gray fisikariak lortu zuen 1729ko uztailaren 14an, Zetazko haritatik eskegitako kalamuzko kordel hezeak erabiliz (garai hartan eroale metalikoen erresistentzia elektriko txikiaz baliatzen ez zirelako).

Aireko linea elektrikoen lehen erabilera praktikoa telegrafo elektrikoarekin izan zen. 1837rako baziren 20 km-rainoko sistema telegrafiko komertzial esperimentalak. Energia elektrikoaren garraioa kopuru handitan 1882an lortu zen, Miesbach eta Munich-en arteko goi-tentsioko linea elektrikoa (60 km-koa) martxan jarri zenean. Hiru faseko korronte alterno ko lehenengo aireko linea 1891ean eraiki zen, Lauffen eta Frankfurten artean, Frankfurteko Nazioarteko Erakusketa Elektrikoa zela eta.

Lehenengo 110 kV-eko aireko linea 1912an jarri zen martxan eta 220 kV-eko lehenengoa 1923an. RWE AG enpresa alemanak izan zen hogeigarren mendeko hirugarren hamarkadan 220 kV-ekoa eraiki zuena eta 1926an Rhin ibaia zeharkatu zuen tentsio horretako linea batekin, 138 m garai ziren euskarriak erabiliz, Voerde-ko dorreak.

American Electric Power enpresak jarri zuen martxan 1953an, Estatu Batuetan, 345 kV-eko lehenengo linea. Suedian 1952an jarri zen martxan lehenengo 400 kV-eko linea, 160 km-koa zen eta hegoaldeko hiri jendetsuagoen eta iparraldeko zentral hidrauliko handienen artean jarri zen. Alemanian 1957koa da 380 kV-eko lehenengo linea. Urte berean Messinako itsasartea igarotzen zuen aireko linea jarri zen martxan Italian; lehenengo 150 kV-eko zirkuitu bat baino ez zuen baina 1971etik 1994arte, eroaleak kendu ziren arte, 220 kV-eko bi zirkuitu izan zituen; erabilitako dorreak, Messinako dorreak, Elbaren gaineko 1. aireko linearenentzat eredu izan ziren eta Elbaren gaineko 2. aireko linea eraiki arte munduko altuenak izan ziren; energia garraiatzeko erabiltzen ez badira ere, oraindino daude zutunik, monumentu historiko izendatuta. 1967tik aurrera Errusia, Kanada eta AEBetan 765 kV-eko aireko lineak eraiki ziren; gaur egun Brasil, Hegoafrika, India eta Txinan ere badaude lan-tentsio horretako aireko lineak. Kazakhstanen 1982an Ekibastuz-eko zentralen eta Kokshetau-ren arteko 1.150 kV-eko korronte alternoko linea trifasiko bat eraiki zuten. Geroago beste nazio batzuetan ere egin dira ikerketak 1.000 kV baino tentsio handiagoko linearekin.


Azterketa matematikoa[aldatu | aldatu iturburu kodea]

Energia elektrikoa garraiatzeko aireko linea bat transmisio-linea bat da eta transmisio-lineak aztertzeko erabiltzen diren ereduak erabili ahal dira bera aztertzeko, baina potentziako sistemen maiztasunekin luzera tipiko lineetan sinplifikazio erabilgarri asko egin daitezke. Potentzia sistemak aztertzerakoan, batzuetan lineen luzera osoan zehar barreiatuta dauden serieko erresistentzia eta induktantzia eta paraleloko ihes-konduktantzia eta kapazitantzia egokiak diren inpedantzia eta admitantzia bildu edo sare sinplez ordezka daitezke.[6][9].

Linea labur eta luzera ertainekoen modeloak[aldatu | aldatu iturburu kodea]

Energia elektrikoa garraiatzeko aireko linea laburren (80 km baino laburragoen) zirkuitu baliokideak erresistentzia bat eta erreaktantzia induktibo bat seriean jarriz osatu ahal dira, admitantzia kontuan hartu barik.

Linea labur baten zirkuitu baliokidea

Zirkuitu baliokide horren egoera iraunkorreko funtzionamenduaren ekuazioak hurrengoak dira:

V_S =V_R + Z.I_R = V_R + I_R.(z.l) = V_R + I_R.(R + j \omega L) = V_R + I_R.(r + j x).l

eta

I_S = I_R

non

  • V_S tentsioa linearen mutur elikatzailean
  • I_S intentsitatea linearen mutur elikatzailean
  • V_R tentsioa linearen mutur elikatuan
  • I_R intentsitatea linearen mutur elikatuan
  • Z linearen serie inpedantzia guztia,
  • z linearen luzera unitateko serie inpedantzia eta
  • l linearen luzera eta
  • R linearen serie erresistentzia guztia
  • L linearen serie induktantzia guztia
  • \omega \ = 2. \pi .f linearen aldagaiak irudikatzeko erabilitako fasoreen abiadura angeluarra
  • f sistemaren maiztasuna
  • r linearen luzera unitateko serie erresistentzia
  • x linearen luzera unitateko erreaktantzia

diren.

Linea luzetxoagoen (80 eta240 km artekoen) zirkuitu baliokideetan shunt kapazitantzia eta hari dagokion suszeptantzia kapazitiboa ere kontuan hartu behar dira (tentsio oso garaiko lineetan baita konduktantzia, ihes korronteak direla eta isolagailuetan izaten diren galerak kontuan hartzeko). Linearen mutur bietan kapazitantzia horri dagokion admitantzia kapazitiboaren erdia jartzen bada  \pi erako zirkuitu baliokidea eratzen da; eta linearen erdian admitantzia hori osoa serie inpedantziaren bi erdien artean kokatzen bada, T erako zirkuitu baliokidea.

Luzera ertaineko linea baten \pi erako zirkuitu baliokidea

Zirkuitu honetan goian aipatutako aldagai eta parametroez gain hurrengoak ere kontuan hartu behar dira:

Y = y.l = G + jB = G + j \omega C = (g + j \omega c).l

non

  • Y linearen shunt admitantzia guztia eta
  • y linearen luzera unitateko shunt admitantzia
  • G linearen shunt konduktantzia osoa
  • B linearen shunt suszeptantzia osoa
  • C linearen shunt kapazitantzia osoa
  • g linearen luzera unitateko shunt konduktantzia
  • c linearen luzera unitateko shunt kapazitantzia

diren. Orokorrean linea elektrikoak ate biko sare baten baliokidetzat har daitezke; ate bi horietako tentsioak eta intentsitateak ABCD parametroen bidez hurrengo eran erlazionatuta daudela:.[13][6][9]:

V_S = A.V_R + B.I_R           eta

I_S =  C.V_R + D.I_R           (izan kontuan ekuazio hauetan agertzen diren parametro eta aldagai guztiak konplexuak direla)


Linea laburren kasuan:

 A = D = 1,       B = Z       eta       C = 0[6][9]


eta luzera ertaineko linetan

 A = D = {Z.Y \over 2} + 1,       B = Z       eta       C = Y.\left ({ 1 + {{Z.Y} \over 4}} \right )[6][9]

Lineak luzeak direnean ezin dira inpedantzia eta admitantzia bilduz osatutako sare sinplez ordeztu, kalkuluak zehaztasun nahikoz egingo badira; baina ate biko modeloa erabili ahal da. Kasu horietan modeloan erabili beharreko parametroak hurrengoak dira:

 A = D = Ch( \gamma .l),       B = Z_C.Sh( \gamma .l)       eta       C = {{Sh( \gamma .l)} \over {Z_C}}[6][9]

Non

 \gamma = \sqrt {z.y} = \sqrt {(r+jx)(g+jb)}       eta


 Z_C = \sqrt {z \over y} = \sqrt {{(r+jx)} \over {(g+jb)}}

parametro kpnplexuak linearen hedapen konstantea eta inpedantzia adierazgarria deitutako parametroak diren.

Kalkuluak konputagailuz egiteko programa[aldatu | aldatu iturburu kodea]

Linea elektrikoen funtzionamendua zehaztasun handiagoz aztertzeko gaur egun programa informatikoak erabiltzen dira. Hemen, eredu gisa, hogeigarren mendeko bederatzigarren hamarkadako hasieran linea elektrikoen funtzionamendu elektrikoa aztertzeko garatutako programa bat aurkezten dugu.

Programa zaharkituta badago ere, sistema elektriko osoa aztertzerakoan linea ordezkatzeko parametroak, lineak sortutako eremu elektriko eta magnetikoak eta garraio egoeren arabera tentsioetan eta korronteetan sortutako desorekak kalkulatu ahal dira. Gainera, erreaktantziez konpentsatuta dagoen bi zirkuituko linea bateko zirkuitu bat konektatuta ez dagoenean ferroerresonantziaz sor daitezkeen gaintentsioak aztertzeko ere erabili ahal da; sortzen diren gaintentsioak ez du zehaztasunez kalkulatzen baina erreaktantzien zein balioekin sortzen diren bai.

Programa, garai haietan kalkulu tekniko gehienak egiteko erabiltzen zen FORTRAN mintzairan idatzita dago eta, ordenagailu berriagotan eta konpilatzaile berriagorekin erabili ahal izateko aldaketa txiki batzuk egin bazaizkio ere, garai haietako FORTRANez dago idatzita. Bestalde, berez, hasieran programa baten ordez bi programa ezberdin izan ziren. Geroago programa bakar batean batu ziren; sarrerako datuak sartzeko zituzten moduak gehiegi aldatu barik, ordura arte erabilitako fitxategien datuak larregi aldatu behar ez izateko. Jatorri hori nabaria da iturri-kodean.

Luzetxoa bada ere, programa osoa aurkezten da, nahi izanez gero, gura duenak programa kopiatu eta, doako FORTRANerako konpilatzaile baten laguntzaz, erabili ahal dezan.

Silverfrosten web orritik jaitsi ahal den eta partikularrentzat doakoa den FTN95 Personal Edition konpilatzailea erabiliz gero, programa inolako aldaketa barik erabili ahal da.

Creative Commons-en Attribution-ShareAlike 3.0 lizentziapean partekatzen den programaren iturri-kodea hurrengoa da:

C  Lineen parametro elektromagnetikoak eta egoerak kalkulatzeko eta
C  ferroerresonantzia azterketak egiteko PARLI programa eta
C  aireko linea elektrikoek sortzen dituzten eremu elektromagnetikoak
C  kalkulatzeko EREMU programa batuta.
C  Programak Jon Peli Oleaga Olabarriak garatu zituen
C  hogeigarren mendeko bederatzigarren hamarkadaren hasieran.
C  Hamarkada horren amaieran batu ziren eta, geroago, konpilatzaile
C  berriagorekin erabili ahal izateko beharrezko aldaketak egin zaizkie.
C Sarrerako datuak oraindino 80 zutabetan egiten dira, behinolako fitxetan legez
C  Bigarren hizkuntza modu naif eta ez sistematikoan sartu da.
C  Programa Creative Commons-en
C  Attribution-ShareAlike 3.0 lizentziapean partekatzen da
C  http://creativecommons.org/licenses/by-sa/3.0/
C
C  The original PARLI program to calculate electromagnetic parameters
C  and operating regimes of overhead transmission lines
C  (it can also be used to conduct ferroresonance studies)
C  and EREMU program to calculate electromagnetic fields created
C  by overhead transmission lines
C  were created by Jon Peli Oleaga Olabarria
C  in the first eighties of the twentieth century.
C  They were merged in the last eighties and later 
C  they have been changed just to compile them with newer compilers
C  and to add a second language; changing only the strictly necessary.
C  Input data still are inserted in 80 columns as with the old cards.
C  The incorporation of a second language has been done 
C  in a rather naif and unsystematic way.
C  The work is licensed under the Creative Commons 
C  Attribution-ShareAlike 3.0 License.
C  http://creativecommons.org/licenses/by-sa/3.0/
C
      PROGRAM LIPARE
      CHARACTER FITX*14,FILE*40,CHK*4
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      DOUBLE PRECISION EXPE
      COMMON /IZENBN/ IZK,EXPE,NOZ,N,KT,TL,MOZ,NFL
      IZK0=1
      IRT=5
      IRF=4
      IDT=6
      IDF=7
      IDU=9
      FITX='LIPARE.DAT'
      CHK=' '
      NOZ=0
      MOZ=0
      NFL=0
      CALL HIZKUN
      OPEN(UNIT=IDF,FILE='LIPARE.RES',STATUS='UNKNOWN')
      OPEN(UNIT=IDU,FILE='LIPARE.RPU',STATUS='UNKNOWN')
    1 IF(CHK.NE.' ')GOTO 12
      PRINT *,'SAR ERABILIKO DEN DATU-FITXATEGIAREN IZENA:'
      PRINT *,'<RETURN> SAKATUZ GERO, FITXATEGIA = '//FITX
      PRINT *,'(  AMAITZEKO SAKA *<RETURN> BAKARRIK  )'
      READ(IRT,'(A)')FILE
      IF(FILE.EQ.'*')GOTO 9999
      IF(FILE.EQ.' ')FILE=FITX
C-VAX      OPEN(UNIT=IRF,STATUS='OLD',FILE=FILE,READONLY)
C-PC      OPEN(UNIT=IRF,STATUS='OLD',FILE=FILE,MODE='READ')
C-PC-F95-C20041023
      OPEN(UNIT=IRF,STATUS='OLD',FILE=FILE,ACTION='READ')
      READ(IRF,'(//)')
C --- ZAINDU SUBRUTINEA, PROGRAMEA ERABILTEKO BAIMENA AZTERTUTEKO
C     CALL ZAINDU(2,IRF) !"1993-02-03 eta 2003-10-06"ko "parli4.for"ean
      CALL ZAINDU(0,IRF,IZK0)
      IF(IZK0.LT.1 .OR. IZK0.GT.2) IZK0=1  ! 20140626 oraingoz hizkuntza bi bakarrik
      READ(IRF,'(A)')CHK
                              BACKSPACE IRF
   12 IZK=IZK0
      IF    (CHK.EQ.'***E'.OR.CHK.EQ.'***e')THEN
        CHK=' '
        CALL EREPRO(CHK)
      ELSEIF(CHK.EQ.'***P'.OR.CHK.EQ.'***p')THEN
        CHK=' '
        CALL LINPAR(CHK)
      ELSE
        PRINT'(A)',' ',
     +' Lehengo PARLI eta EREMU programak LIPARE barri honetan'//
     +' batu ziranetik,',' datu fitxategietan euretatik'//
     +' zein erabili gura dan esan behar da,',
     +' datu fitxategiaren bostgarren lerroaren lehenengo tokietan',
     +' ***EREMU ala ***PARLI iminiz.',' ',
     +' Aldatu egizuz lehengo fitxategiak.',' '
        GOTO 9999
      ENDIF
      GOTO 1
 9999 END
CJFO      PROGRAM LINPAR
C     Aireko linea elektrikoen parametro esanguratsuenak
C     eta haien funtzionatzeko moduak kalkulatzeko programa
      SUBROUTINE LINPAR(CHK)
      PARAMETER (ND1=10,ND2=12)
C-20131026AN ND1 8TIK 10ERA ALDATUTA ZIRKUITU 2 4 LUR KABLEREKIN ONARTZEKO
      COMPLEX AD(ND2,ND2),BC(ND2,ND2),RI(ND2,ND2),RF(ND2,ND2)
     &,Z(ND1,ND1),Y(ND1,ND1),W(ND1,ND1)
      COMPLEX ER(ND2),ES(ND2),PR,PS,PRT,PST,PTR,PTS,ERO
      REAL DX(ND1),DY(ND1),R(ND1),G(ND1),RMG(ND1),RAD(ND1),D(ND1)
      INTEGER IN(ND1),ID(ND1),KR(ND1),LOC(ND2)
      CHARACTER EGUNA*10,IZENA*56,TT1*27,TT2*27,SII*44,SAK*44,SI*29
     &,SA*29,KG*24,OS*18,ZAT*5,EUSK*5,EROA*15,AZALP*5,CHK*4
      DOUBLE PRECISION EXPE
      COMMON /IZENBN/ IZK,EXPE,NOZ,N,KT,TL,MOZ,NFL
      COMMON /IZENBC/ EGUNA,IZENA,EUSK,EROA
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      PBAS=100.
    2 CALL HIRUAT(IRF,CHK,*9999,*6,*9998)
    6 READ(IRF,5011,END=9998)TL,IZENA,EXPE,EGUNA
 5011 FORMAT(F4.0,A56,F10.2,A10)
      READ(IRF,5012)IZK2,NT,NC,KP,FNY
      IF (IZK2.NE.0) IZK=IZK2
      IF (IZK.LT.1 .OR. IZK.GT.2) IZK=1
C-20131027 5012 FORMAT(4I5,F10.5)
 5012 FORMAT(4I2,F7.3)
                print *,'LINPAR  ',izena
C     IZK = ERABILIKO DEN HIZKUNTZA
C     NT  = LINEAN AZTERTUKO DIREN ZATI GEHI EGOEREN KOPURUA
C     NC  = LINEAREN ZIRKUITU KOPURUA
C     KP  = 1=> EROALE BATZUK BESTE BATZUEKIN PARALELOAN JARTZEN DIRELA
C           EROALE BAKOTZAREN 54. ZUTABEAN JARRI ZEIN BESTE EROALEREKIN.
C           ( PARALELOAN JARRI ETA GERO EZIN DIRA AZTERTU EGOERAK,
C             BAKARRIk ATARA AHAL DIRA ALDAEZIN OROKORRAK )
C     FNY = LINEAREN FREKUENTZIA
      IF(IZK-1)11,10,11
   10 ZAT='ZATI   '
      OS='OSAGAI SIMETRIKOAK'
      KG='   ALDAEZIN  OROKORRAK  '
      TT1='  (EGIAZKO KONFIGURAZIOA)  '
      TT2=' (KONFIGURAZIO BALIOKIDEA) '
      SI='  SERIE INPEDANTZEAK (OHM)   '
      SA='   SHUNT ADMITANTZEAK (MHO)  '
      SII=' SERIE INPEDANTZE INDUKTIBOAK (OHM:KM)  '
      SAK=' SHUNT ADMITANTZE KAPAZITIBOAK (MHO:KM) '
      GO TO 12
   11 ZAT='SECTION'						!7
      OS=' SIMETRICAL COMP. '			!18
      KG='A,B,C,D  TRANSFER MATRIX'		!24
      TT1=' ( ACTUAL CONFIGURATION )'	!27
      TT2='(EQUIVALENT CONFIGURATION)'	!27
      SI='   SERIES IMPEDANCES (OHM)   '	!29
      SA='   SHUNT ADMITANCES (MHO)    '	!29
      SII='  SERIES INDUCTIVE IMPEDANCES (OHM:KM)  '	!40
      SAK='  SHUNT CAPACITIVE ADMITANCES (MHO:KM)  '	!40
 
   12 CONTINUE
      N3=NC*3
      N2=N3*2
      DO 20 I=1,N2
      DO 20 J=1,N2
      RF(I,J)=(0.,0.)
      IF(I.EQ.J)RF(I,I)=(1.,0.)
   20 CONTINUE
      XLUZ=0.
      DO 1000 N=1,NT
C-20131027      READ(IRF,5021)KT,K1,K2,K3,NCO,NST,ELT,RHO,EUSK,EROA
      READ(IRF,5021)KT,NCO,K1,K2,K3,NST,ELT,RHO,EUSK,EROA
C-20131027 5021 FORMAT(6I5,2F10.5,A5,A15)
 5021 FORMAT(I1,I3,4I1,F7.3,F5.2,33X,A5,A15)
C     KT =0 => TENTSIO ETA INTENTSITATEEN AZTERKETA,
C        =1 => LINEAREN SERIE INPEDANTZE ETA SHUNT ADMITANTZEEN KALKULUA
C        =2 => KONEKTATUTAKO SERIE INPEDANTZEEN DATUAK
C        =3 => KONEKTATUTAKO SHUNT ADMITANTZEEN DATUAK
C        =4 => LINEAREN ERRESONANTZIA FERROMAGNETIKOA AZTERTUKO DA
C     K1 <1 => EZ DIRA ATARATZEN SERIE INPEDANTZEAK ETA SHUNT ADMITANTZEAK
C        =1 => KONFIGURAZIO BALIOKIDEAREN SI ETA SA-AK IDAZTEN DIRA
C        >1 => KONFIGURAZIO ORIGINAL ETA BALIOKIDEAREN SI ETA SA-AK IDAZTEN DIRA
C        (KT=4)2=> ERREAKTAN.ALBO BATEAN, 2=>ER.BESTE ALBOAN, 3=>ER.ALBO BIETAN
C     K2 <1 => EZ DIRA ATARATZEN ALDAEZIN OROKORTUAK
C        >0 => ALDAEZIN OROKORTUAK ATARATZEN DIRA (ZATIARENAK)
C        (KT=4)0=>, 1=>
C     K3 <1 => EZ DIRA ATARATZEN ALDAEZIN OROKORTUAK
C        =1 => FASEEN ALDAEZIN OROKORTUAK ATARATZEN DIRA (ZATI GUZTIENAK)
C        >1 => FASEEN ETA KONPOSATU SIMETRIKOEN ALDAEZIN OROKORTUAK ATARATZEN DIRA
C        (KT=4)0=>IRTEERA GRAFIKO LOGARITMIKOAN, 1=>IRTEERA GRAFIKO LINEALEAN
C     NCO = EROALEEN KOPURUA, EDO KT=4 BADA PUNTU KOPURUA GRAFIKOAN
C     NST = AZPITARTEEN KOPURUA (KALKULUAK ZEHAZTASUN HANDIAGOZ EGITEKO)
C     ELT = TARTEAREN LUZERA KM-TAN
C     RHO = LINEA TARTE HONETAKO LURRAREN ERRESISTIBITATEA OHM-M-TAN
C     EUSK = ERABILITAKO EUSKARRIEN IZENA (IRTEERAN IDAZTEKO BAKARRIK)
C     EROA = ERABILITAKO EROALEEN IZENA (IRTEERAN IDAZTEKO BAKARRIK)
      IF(KT.LE.0)GO TO 300
      IF(KT.GE.2)GO TO 200
C
C --- LINE BATEN SERIE INPEDANTZE ETA SHUNT ADMITANTZEEN KALKULUA
      IF(K1.LT.0)K1=0     !C-20131026
      IF(K2.LT.0)K2=0     !C-20131026
      IF(K3.LT.0)K3=0     !C-20131026
      XLUZ=XLUZ+ELT
      ELT=ELT/NST
      DO 100 I=1,NCO
      READ(IRF,5100)RMG(I),RAD(I),IN(I),D(I),ID(I),DX(I),DY(I)
     +,R(I),KR(I)
 5100 FORMAT(2F7.4,I2,F7.4,I2,2F7.3,7X,F7.5,I2)
C     RMG = EROALEAREN "BATEZBESTEKO ERRADIO GEOMETRIKOA" CM-TAN
C     RAD = EROALEAREN ERRADIOA CM-TAN
C     IN  = EROALE SORTAREN AZPIEROALE KOPURUA
C     D   = AZPIEROALEEN ARTEKO DISTANTZIA EDO EROALE SORTAREN DIAMETROA M-TAN
C     ID  =0=> AURREKO DATUA AZPIEROALEN ARTEKO DISTANTZIA DA
C         =1=> AURREKO DATUA EROALE SORTAREN DIAMETROA DA
C     DX  = EROALEAREN X KOORDENATUA M-TAN
C     DY  = EROALEAREN y KOORDENATUA M-TAN
C     R   = EROALEEN ERRESISTENTZIA OHM/KM-TAN
C     KR  =
      IF(R(I).GT.0.)GO TO 100
      J=I-1
      RMG(I)=RMG(J)
      RAD(I)=RAD(J)
      IN(I)=IN(J)
      ID(I)=ID(J)
      D(I)=D(J)
      R(I)=R(J)
  100 CONTINUE
C
C --- SARRERAKO DATUAK ATARATZEN DIRA
      CALL IZENBU
      IF(N.EQ.1)LOZ=NOZ
      IF(IZK-1)102,101,102
  101 WRITE(IDF,6101)NST,ELT,NST*ELT,RHO,FNY
 6101 FORMAT(1H0,23X,'SARRERAKO DATUAK'/1H0,23X,'LINE-LUZERA',I3,'*',
     1F6.2,'=',F7.2,' KM'/1H ,23X,'LUR-ERRESISTIBITATEA',F9.2,' OHM.M'/
     21H ,23X,'MAIZTASUNA', 8X,F11.2,' HZ'/ 1H0,23X,28('*'),
     3' EROALEEN EZAUGARRIAK '          ,29('*') )
      GO TO 103
  102 WRITE(IDF,6102)NST,ELT,NST*ELT,RHO,FNY
 6102 FORMAT(1H0,23X,'INPUT DATA'/1H0,23X,'LINE LONG. ',I3,'*',
     1F6.2,'=',F7.2,' KM'/1H ,23X,'GROUND RESISTIVITY',F11.2,' OHM.M'/
     21H ,23X,'FREQUENCY ', 8X,F11.2,' HZ'/ 1H0,23X,27('*'),
     3'  CONDUCTORS´ FEATURES  ',27('*') )
  103 WRITE(IDF,6103)(I,RMG(I),RAD(I),IN(I),D(I),ID(I),DX(I),DY(I),R(I),
     1I=1,NCO)
 6103 FORMAT(1H ,26X,'..   GMR(CM)   RAD(CM)     IN   D(M)  ID      DX(M
     1)   DY(M)      R(OHM:KM)'//(1H ,26X,I2,2F10.5,I7,F7.2,I4,3X,2F8.2,
     24X,F11.7))
      CALL ZY(ND1,NCO,FNY,RHO,Z,Y,DX,DY,RAD,RMG,IN,D,ID,R,G,KR,W,KP)
      IF(K1.LE.1)GO TO 112
C
C --- BENETAKO EROALE GUZTIEI JAGOKEZAN BALIOAK ATARATZEN DIRA URTEERAN
      CALL IZENBU
      WRITE(IDF,6111)SII,TT1
 6111 FORMAT(1H0,39X,A40,/,46X,A27 )
      CALL ATARAC(Z,NCO,ND1)
      IF(NCO.GT.6)CALL IZENBU
      WRITE(IDF,6111)SAK,TT1
      CALL ATARAC(Y,NCO,ND1)
C
C --- BENETAKO EROALEEN BALIOKIDEEI JAGOKEZAN BALIOAK ATARATZEN DIRA
  112 CALL LABURT(Z,NCO,N3,ND1,KR,KP,ID)
      N2=N3*2
      IF(K1.LT.1)GO TO 113
      CALL IZENBU
      WRITE(IDF,6111)SII,TT2
      CALL ATARAC(Z,N3,ND1)
      IF(N3.GT.6)CALL IZENBU
      WRITE(IDF,6111)SAK,TT2
      CALL ATARAC(W,N3,ND1)
C
C --- LINEKO ZATI AU BENETAKO LINE BAT DA
  113 CONTINUE
      DO 115 I=1,N3
      DO 115 J=1,N3
      Z(I,J)=ELT*Z(I,J)
  115 Y(I,J)=ELT*W(I,J)
      GO TO 211
C
C --- LINEKO ZATI AU SERIE INPEDANTZE ALA SHUNT ADMITANTZE BAT DA
  200 DO 201 I=1,N3
      DO 201 J=1,N3
  201 Z(I,J)=CMPLX(0.,0.)
      IF(KT.LT.4)GO TO 205
C
C --- ERRESONANTZE KASU BAT AZTERTZEN DA
      ELT=ELT*FNY*6.2831853
      RHO=RHO*FNY*6.2831853
      CALL RES(ND2,NC,K1,K2,K3,ELT,RHO,NCO,RF,RI,BC,AD,ES,ER,LOC)
      GO TO 1000
C
C     ******************************************************************
C     ASTERISKOZ OSATUTAKO LERROEN ARTEAN PROGRAMA EZ DAGO PRESTATUTA
C     SARTZEKO FASE ARTEKO SERIE ELKAR INPEDANTZEAK EDO ADMITANTZEAK.
C     HORRETARAKO PROGRAMAREN ZATI HAU ALDATU BEHAR DA
  205 READ(IRF,5205)(Z(I,I),I=1,N3)
 5205 FORMAT(12F6.3)
C     ******************************************************************
      IF(K1.LE.0)GO TO 211
      CALL IZENBU
      IF(N.EQ.1)LOZ=NOZ
      IF(KT.EQ.2)WRITE(IDF,6210)ZAT,N,SI
      IF(KT.EQ.3)WRITE(IDF,6210)ZAT,N,SA
 6210 FORMAT(1H0, 5X,A7,I3,33X,A29   )
      CALL ATARAC(Z,N3,ND1)
C
C --- ZATIAREN ALDAEZIN GENERALIZATUAK KALKULATZEKO DO 220
  211 DO 220 I1=1,N3
      DO 220 J1=1,N3
      I2=I1+N3
      J2=J1+N3
      AD(I1,J1)=(0.,0.)
      AD(I2,J2)=(0.,0.)
      IF(I1.NE.J1)GO TO 219
      AD(I1,J1)=(1.,0.)
      AD(I2,J2)=(1.,0.)
  219 AD(I1,J2)=(0.,0.)
      IF(KT.EQ.2)AD(I1,J2)=Z(I1,J1)
      AD(I2,J1)=(0.,0.)
      IF(KT.EQ.3)AD(I2,J1)=Z(I1,J1)
      IF(KT.NE.1)GO TO 220
      AD(I2,J1)=Y(I1,J1)/2.
      BC(I1,J1)=AD(I1,J1)
      BC(I2,J2)=AD(I2,J2)
      BC(I1,J2)=Z(I1,J1)
      BC(I2,J1)=AD(I1,J2)
  220 CONTINUE
C
      IF(KT.GT.1)GO TO 244
      IF(K2.LE.0)GO TO 225
      CALL IZENBU
      WRITE(IDF,6224)KG
 6224 FORMAT(1H0,49X,A24   )
  225 CALL BITRIC(BC,N2,N2,AD,N2,N2,RI,ND2)
      CALL BITRIC(AD,N2,N2,RI,N2,N2,BC,ND2)
      DO 230 I=1,N2
      DO 230 J=1,N2
      RI(I,J)=(0.,0.)
  230 IF(I.EQ.J)RI(I,J)=(1.,0.)
      KRL=0
      DO 240 I=1,NST
      IF(KRL.GT.0)GO TO 239
      CALL BITRIC(BC,N2,N2,RI,N2,N2,AD,ND2)
      KRL=1
      GO TO 240
  239 CALL BITRIC(BC,N2,N2,AD,N2,N2,RI,ND2)
      KRL=0
  240 CONTINUE
      IF(KRL.GT.0)GO TO 243
      IF(K2.GE.1)CALL ATARAC(RI,N2,ND2)
      CALL BITRIC(RI,N2,N2,RF,N2,N2,BC,ND2)
      GO TO 245
  243 IF(K2.GE.1)CALL ATARAC(AD,N2,ND2)
  244 CALL BITRIC(AD,N2,N2,RF,N2,N2,BC,ND2)
  245 DO 250 I=1,N2
      DO 250 J=1,N2
  250 RF(I,J)=BC(I,J)
      IF(K3.EQ.0)GO TO 1000
      KT=3
      CALL IZENBU
      WRITE(IDF,6224)KG
      CALL ATARAC(RF,N2,ND2)
      IF(K3.LT.2)GO TO 1000
C
C --- LINEKO ALDAEZIN GENERALIZATUEN OSAGAI SIMETRIKOAK ATARATZEN DIRA
      CALL IZENBU
      WRITE(IDF,6224)KG
      WRITE(IDF,6255)OS
 6255 FORMAT(1H ,51X,1H(,A18,1H) )
      CALL OSASYM(2,NC,ND2,RF,RI,AD,BC)
      CALL ATARAC(RI,N2,ND2)
      IF(K3.LT.2) GOTO 1000
C
C --- BESTE FITXATEGI BAT EMOITZA LABURTUAKAZ IDAZTEN DA
      IF(MOD(NFL,60).EQ.0)THEN
       MOZ=MOZ+1
       IF(IZK.EQ.2)WRITE(IDU,6260)PBAS,MOZ
 6260 FORMAT(1H1,'POWER BASE =',F5.0,' MVA',103X,'PAGE',I3/
     1 1X,49('-'),'LINE CHARACTERISTICS IN P.U.     ',50('-')/
     3 1X,17('-'),'-LINE',18('-'),'   LONG.',7X,'  SERIES IMPEDANCE ',
     4 11X,'   SHUNT ADMITANCE    ',4X,'COMENT.',1X,' PAGES   BASE'/
     5 43X,'(KM)',2(6X,'DIRECT C.',5X,'HOMOPOLAR C.',1X),16X,'(KV)')
       IF(IZK.NE.2)WRITE(IDU,6261)PBAS,MOZ
 6261 FORMAT(1H1,'POTENTZIA BASE =',F5.0,' MVA',95X,'ORRIALDE',I3/
     1 1X,52('-'),'LINEEN EZAUGARRIAK  P.U.TAN',53('-')/
     3 1X,17('-'),'LINEA',18('-'),'  LUZER.',7X,' SERIE  INPEDANTZEA',
     4 11X,'   SHUNT  ADMITANTZEA ',4X,'ARGIBI.',1X,'ORRIAL.  BASE'/
     5 43X,'(KM)',2(6X,'O.ZUZENAK',6X,'O.HOMOPOLAR',1X),16X,'(KV)')
       NFL=5
      ENDIF
      AZALP='     '
      IF(NC.GT.1)THEN
       AZALP='CIR.1'
       IF(KP.EQ.1)AZALP='CIR.P'
      ENDIF
      TBAS=1200.
      IF(TL.LT.899.)TBAS=700.
      IF(TL.LT.499.)TBAS=400.
      IF(TL.LT.299.)TBAS=230.
      IF(TL.LT.199.)TBAS=132.
      IF(TL.LT.99.)TBAS=66.
      IF(TL.LT.55.)TBAS=45.
      IF(TL.LT.37.)TBAS=30.
      IF(TL.LT.25.)TBAS=20.
      IF(TL.LT.17.)TBAS=15.
      IF(TL.LT.14.)TBAS=13.2
      IF(TL.LT.8.)TBAS=5.
      BASE=TBAS**2/PBAS
      IF(MOD(NFL,5).EQ.0)THEN
       WRITE(IDU,'(1X)')
       NFL=NFL+1
      ENDIF
      WRITE(IDU,6265)IZENA,XLUZ
     1,RI(1,N3+1)/BASE,RI(3,N3+3)/BASE,RI(N3+1,1)*BASE,RI(N3+3,3)*BASE
     2,AZALP,LOZ,NOZ,TBAS
 6265 FORMAT(1X,A40,F8.2,8F8.4,1X,A5,1X,I3,'-',I3,F6.1)
      NFL=NFL+1
      GO TO 1000
C
C --- LINEKO TENTSIO ETA INTENTSITATEEN AZTERKETA
  300 CALL IRAKUR(NC,ER,LOC,ND2)
      DO 310 I=1,N2
      DO 310 J=1,N2
  310 RI(I,J)=RF(I,J)
      CALL EXCHAN(RI,N2,LOC,ND2,*1000,2)
      CALL BITRIC(RI,N2,N2,ER,N2,1,ES,ND2)
      CALL IZENBU
      DO 320 I=1,N2
      IF(LOC(I).EQ.0)GO TO 320
      ERO=ER(I)
      ER(I)=ES(LOC(I))
      ES(LOC(I))=ERO
  320 CONTINUE
      PTR=(0.,0.)
      PTS=(0.,0.)
      DO 330 I=1,NC
      I1=3*(I-1)
      PRT=(0.,0.)
      PST=(0.,0.)
      DO 325 J=1,3
      I1=I1+1
      I2=I1+N3
      PR=ER(I1)*CONJG(ER(I2))
      PS=ES(I1)*CONJG(ES(I2))
      TJ=CABS(ES(I1)-ER(I1))
      PRT=PRT+PR
      PST=PST+PS
  325 WRITE(IDF,6325)I,J,CABS(ES(I1))*SQRT(3.),CABS(ES(I2)),PS,
     1CABS(ER(I1))*SQRT(3.),CABS(ER(I2)),PR,TJ,LOC(I1),LOC(I2)
 6325 FORMAT(1H ,I2,I4,2(4F13.3,2X),F9.3,I4,I2)
      PTR=PTR+PRT
      PTS=PTS+PST
  330 WRITE(IDF,6330)I,PST,PRT
 6330 FORMAT(1H0,I2,4X,2(26X,2F13.3,2X)//)
      WRITE(IDF,6330)0,PTS,PTR
 1000 CONTINUE
      GOTO 2
 9998 CLOSE(UNIT=IRF)
      CHK=' '
 9999 RETURN
      END
C---------------------------------
      SUBROUTINE ATARAC(ZC,NEU,ND)
      COMPLEX ZC(ND,ND)
      CHARACTER AX*20,AY*45
      INTEGER ASI,BUK,ZEN
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      LZ=1
      ASI=1
      BUK=6
   11 IF(NEU.LE.BUK)GO TO 21
      WRITE(IDF,911)(K,K=ASI,BUK)
  911 FORMAT(1H / 3X,6(15X,I2))
      WRITE(IDF,912)(K,(REAL(ZC(K,I)),I=ASI,BUK),(AIMAG(ZC(K,I))
     1,I=ASI,BUK),K=LZ,NEU)
  912 FORMAT(1H /(7X,I2,6(3X,E14.7),/,9X,6(3X,E14.7)/))
      ASI=ASI+6
      BUK=BUK+6
      GO TO 11
   21 ZEN=NEU-ASI+1
      WRITE(AX,921)ZEN
  921 FORMAT(9H(1H / 3X,,I1,9H(15X,I2)))
      WRITE(IDF,AX)(K,K=ASI,NEU)
      WRITE(AY,922)ZEN,ZEN
  922 FORMAT('(1H /(7X,I2,',I1,'(3X,E14.7)/9X,',I1,'(3X,E14.7)/))')
      WRITE(IDF,AY)(K,(REAL(ZC(K,I)),I=ASI,NEU),(AIMAG(ZC(K,I))
     1,I=ASI,NEU),K=LZ,NEU)
      RETURN
      END
      SUBROUTINE BITRIC(A,NFA,NCA,B,NFB,NCB,C,ND)
C     A ETA B MATRIZEAK BIDERKATZEKO AZPIERRUTINA
C-20041122      COMPLEX A(ND, 1),B(ND, 1),C(ND, 1)
      COMPLEX A(ND, *),B(ND, *),C(ND, *)  !C-20041122
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      IF(NCA.NE.NFB)GO TO 88
       DO 200 J=1,NCB
        DO 200 I=1,NFA
        C(I,J)=CMPLX(0.,0.)
         DO 200 K=1,NCA
  200    C(I,J)=C(I,J)+A(I,K)*B(K,J)
   99 RETURN
   88 WRITE(IDF,888)
  888 FORMAT(6X/ 6X,'SR. BITRIC - MULTZOEN DIMENTSIOAK EZ DIRA EGOKIAK')
      GO TO 99
      END
      SUBROUTINE EXCHAN(A,KA,LOC,ND,*,KON)
C     VS = A x VR ekuazioetan VS eta VR bektoreetako bi balio
C     bata bestearen tokian jartzeko azpierrutina.
C     Egingo diren bikote aldaketak LOC(K)-z definitzen dira.
C     LOC(K) guztiak =K badira EXCHAN eta INTRIX antzekoak dira.
C     (ez guztiz berdinak)
C20041023      COMPLEX A(ND, 1),DIV,FAC
      COMPLEX A(ND, ND),DIV,FAC
      INTEGER LOC(ND)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      XL=1.E-20
      IF(KON-1)3,1,2
    1 DO 100 I=1,KA
  100 LOC(I)=I
    2 DO 200 K=1,KA
      IF(LOC(K).EQ.0)GO TO 200
      DIV=-A(LOC(K),K)
      IF(CABS(DIV).LT.XL)GO TO 1001
      A(LOC(K),K)=(1.,0.)
      DO 210 J=1,KA
  210 A(LOC(K),J)=A(LOC(K),J)/DIV
      DO 220 I=1,KA
      IF(I.EQ.LOC(K))GO TO 220
      FAC=A(I,K)
      A(I,K)=(0.,0.)
      DO 221 J=1,KA
  221 A(I,J)=A(I,J)+FAC*A(LOC(K),J)
  220 CONTINUE
      DO 230 I=1,KA
  230 A(I,K)=-A(I,K)
  200 CONTINUE
    3 RETURN
 1001 WRITE(IDF,999)LOC(K),K
  999 FORMAT(1H //11X,1H(,I2,1H,,I2,22H) ARDATZA TXIKIEGIA DA)
      RETURN 1
      END
      SUBROUTINE LABURT(A,NI,NF,ND,KS,KQ,KR)
C20131023      COMPLEX A(ND, 1)
      COMPLEX A(ND,*)
      INTEGER KR(ND),KS(ND)
      NF=NI
      DO 1100 I=1,NI
 1100 KR(I)=KS(I)
      IF(KQ.EQ.0)GO TO 1
      DO 1300 I=1,NI
      IF(KR(I).LE.I.OR.KR(I).GT.NI)GO TO 1300
      K=KR(I)
      DO 310 J=1,NI
  310 A(I,J)=A(I,J)-A(K,J)
      DO 320 J=1,NI
  320 A(J,I)=A(J,I)-A(J,K)
      KR(I)=0
 1300 CONTINUE
    1 DO 1600 L=1,NI
      M=NI-L+1
      IF(KR(M).GT.0)GO TO 1600
      DO 630 I=1,NF
      DO 630 J=1,NF
      IF(I.EQ.M.OR.J.EQ.M)GO TO 630
      A(I,J)=A(I,J)-A(I,M)/A(M,M)*A(M,J)
  630 CONTINUE
      IF(M.EQ.NF)GO TO 601
      DO 610 I=M,NF-1
      DO 610 J=1,NF
  610 A(J,I)=A(J,I+1)
      DO 620 I=M,NF-1
      DO 620 J=1,NF-1
  620 A(I,J)=A(I+1,J)
  601 NF=NF-1
 1600 CONTINUE
      RETURN
      END
      SUBROUTINE IRAKUR(NC,ER,LOC,ND)
      COMPLEX ER(ND)
      REAL VM(3),VF(3),AM(3),AF(3)
      INTEGER IV(3),NCV(3),NFV(3),IA(3),NCA(3),NFA(3),LOC(ND)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      DO 2100 I=1,NC
      READ(IRF,106)NCT,(VM(J),VF(J),IV(J),NCV(J),NFV(J),J=1,3)
     +                ,(AM(J),AF(J),IA(J),NCA(J),NFA(J),J=1,3)
  106 FORMAT(I2,6(2F5.2,3I1))
C     NCT = ZIRKUITUA
C     VM  = TENTSIOAREN (INTENTSITATEAREN) MODULUA
C     VF  = TENTSIOAREN (INTENTSITATEAREN) FASEA
C     IV  =0=>MUTUR ELIKATUAREN TENTSIOA
C         =1=>MUTUR ELIKATZILEAREN TENTSIOA
C         =2=>MUTUR ELIKATZAILEAREN INTENTSITATEA
C     NCV = ZIRKUITUA
C     NFV = FASEA
C     AM  = INTENTSITATEAREN (TENTSIOAREN) MODULUA
C     AF  = INTENTSITATEAREN (TENTSIOAREN) FASEA
C     IA  =0=>MUTUR ELIKATUAREN INTENTSITATEA
C         =1=>MUTUR ELIKATZAILEAREN TENTSIOA
C         =2=>MUTUR ELIKATZAILEAREN INTENTSITATEA
C         =3=>POTENTZIA MUTUR ELIKATUAN (EZ BETI)
C     NCA = ZIRKUITUA
C     NFA = FASEA
      DO 2110 J=1,3
      I1=3*(NCT-1)+J
      I2=I1+NC*3
      LOC(I1)=0
      IF(IV(J).EQ.1)LOC(I1)=3*(NCV(J)-1)+NFV(J)
      IF(IV(J).EQ.2)LOC(I1)=3*(NCV(J)-1)+NFV(J)+NC*3
      LOC(I2)=0
      IF(IA(J).EQ.1)LOC(I2)=3*(NCA(J)-1)+NFA(J)
      IF(IA(J).EQ.2)LOC(I2)=3*(NCA(J)-1)+NFA(J)+NC*3
      VF(J)=VF(J)/57.295779
      IF(IV(J).NE.2)VM(J)=VM(J)/SQRT(3.)
      ER(I1)=CMPLX(VM(J)*COS(VF(J)),VM(J)*SIN(VF(J)))
      IF(IA(J).EQ.3)GO TO 2111
      AF(J)=AF(J)/57.295779
      IF(IA(J).EQ.1)AM(J)=AM(J)/SQRT(3.)
      GO TO 2110
 2111 AM(J)=AM(J)/VM(J)/COS(ATAN(AF(J)))
      AF(J)=VF(J)-ATAN(AF(J))
 2110 ER(I2)=CMPLX(AM(J)*COS(AF(J)),AM(J)*SIN(AF(J)))
 2100 CONTINUE
      RETURN
      END
      SUBROUTINE OSASYM(NB,NC,ND,Z123,ZDIH,SD,SI)
      COMPLEX A0,A1,A2,Z123(ND,ND),ZDIH(ND,ND),SD(ND,ND),SI(ND,ND)
      NTT=3*NB*NC
      A0=(1.,0.)
      A1=(-0.5, 0.8660254)
      A2=(-0.5,-0.8660254)
      DO 100 I=1,NTT
      DO 100 J=1,NTT
  100 SD(I,J)=(0.,0.)
      DO 200 I=1,NTT-2,3
      SD(I  ,I  )=A0
      SD(I  ,I+1)=A0
      SD(I  ,I+2)=A0
      SD(I+1,I  )=A2
      SD(I+1,I+1)=A1
      SD(I+1,I+2)=A0
      SD(I+2,I  )=A1
      SD(I+2,I+1)=A2
  200 SD(I+2,I+2)=A0
      N=ND
      CALL BITRIC(Z123,NTT,NTT,SD,NTT,NTT,SI,N)
      DO 300 I=1,NTT-2,3
      SD(I  ,I  )=A0/3.
      SD(I  ,I+1)=A1/3.
      SD(I  ,I+2)=A2/3.
      SD(I+1,I  )=A0/3.
      SD(I+1,I+1)=A2/3.
      SD(I+1,I+2)=A1/3.
      SD(I+2,I  )=A0/3.
      SD(I+2,I+1)=A0/3.
  300 SD(I+2,I+2)=A0/3.
      CALL BITRIC(SD,NTT,NTT,SI,NTT,NTT,ZDIH,N)
      RETURN
      END
      SUBROUTINE REGRA(Y,NG,XG,X ,K3)
C20131024      REAL Y( 6,1)
      REAL Y( 6,236)
C20041023      CHARACTER*1 CL(121),PUN/'.'/,ZUR/' '/
C20041023      CHARACTER*1 E1/'1'/,E2/'2'/,E3/'3'/,E4/'4'/,E5/'5'/,E6/'6'/
      CHARACTER*1 CL(121),PUN,ZUR
      CHARACTER*1 E1,E2,E3,E4,E5,E6
      CHARACTER MEZUA1(4)*38,MEZUA2(4)*35,MEZUA3(4)*5
      DOUBLE PRECISION EXPE
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      COMMON /IZENBN/ IZK,EXPE,NOZ,N,KT,TL,MOZ,NFL
      MEZUA1(1)='FASEEN ETA NEUTROAREN ARTEKO TENTSIOAK'
      MEZUA1(2)='PHASE NEUTRAL VOLTAGES                '
      MEZUA2(1)='FASEKO TENTSIOEN BALIORIK HANDIENA:'
      MEZUA2(2)='      PHASE VOLTAGE, MAXIMUM VALUE:'
      MEZUA3(1)=' FASE'
      MEZUA3(2)='PHASE'
      PUN='.'
      ZUR=' '
      E1='1'
      E2='2'
      E3='3'
      E4='4'
      E5='5'
      E6='6'
      IF(K3.GT.0)GO TO 100
      WRITE(IDF,105)MEZUA1(IZK)
  105 FORMAT(1H1,10X,A38,77X,'V:VMAX',/1H0
     1,2X,'XL(OHM)',1X,3('1',11X,'2',6X,'3    4   5  6  7 8 9 '),'1')
      WRITE(IDF,103)
  103 FORMAT(11X,121('.'))
      GO TO 200
  100 WRITE(IDF,205)MEZUA1(IZK)
  205 FORMAT(1H1,10X,A38,77X,'V:VMAX',/ 1H0,2X,'XL(OHM)',1X,
     1'0.00',25X,'0.25',26X,'0.50',26X,'0.75',24X,'1.00')
      WRITE(IDF,103)
  200 YMAX=ABS(Y(1,1))
      DO 1 I=1,6
      DO 1 IJ=1,NG
    1 IF(YMAX.LT.ABS(Y(I,IJ))) YMAX=ABS(Y(I,IJ))
      DO 10 I=1,121
   10 CL(I)=ZUR
      DO 20 I=1,NG
      X=X+XG
      CL(1)=PUN
      IF(K3.GT.0)GO TO 21
      CL(41)=PUN
      CL(81)=PUN
      CL(121)=PUN
      J1=40.*ALOG10(Y(1,I)/YMAX)+121.5
      J2=40.*ALOG10(Y(2,I)/YMAX)+121.5
      J3=40.*ALOG10(Y(3,I)/YMAX)+121.5
      J4=40.*ALOG10(Y(4,I)/YMAX)+121.5
      J5=40.*ALOG10(Y(5,I)/YMAX)+121.5
      J6=40.*ALOG10(Y(6,I)/YMAX)+121.5
      IF(J1.GE.1)CL(J1)=E1
      IF(J2.GE.1)CL(J2)=E2
      IF(J3.GE.1)CL(J3)=E3
      IF(J4.GE.1)CL(J4)=E4
      IF(J5.GE.1)CL(J5)=E5
      IF(J6.GE.1)CL(J6)=E6
      WRITE(IDF,104)X,(CL(J),J=1,121)
  104 FORMAT(1H , F9.3,1X,121A1)
      GO TO 22
   21 CL(61)=PUN
      CL(31)=PUN
      CL(91)=PUN
      CL(121)=PUN
      J1=120.*Y(1,I)/YMAX+1.5
      J2=120.*Y(2,I)/YMAX+1.5
      J3=120.*Y(3,I)/YMAX+1.5
      J4=120.*Y(4,I)/YMAX+1.5
      J5=120.*Y(5,I)/YMAX+1.5
      J6=120.*Y(6,I)/YMAX+1.5
      CL(J1)=E1
      CL(J2)=E2
      CL(J3)=E3
      CL(J4)=E4
      CL(J5)=E5
      CL(J6)=E6
      WRITE(IDF,104)X,(CL(J),J=1,121)
C  204 FORMAT(1H ,F10.3,4X,121A1)
   22 IF (J1.GT.0) CL(J1)=ZUR
      IF (J2.GT.0) CL(J2)=ZUR
      IF (J3.GT.0) CL(J3)=ZUR
      IF (J4.GT.0) CL(J4)=ZUR
      IF (J5.GT.0) CL(J5)=ZUR
      IF (J6.GT.0) CL(J6)=ZUR
   20 CONTINUE
      IF(K3.EQ.0)WRITE(IDF,103)
      IF(K3.GT.0)WRITE(IDF,103)
      WRITE(IDF,108)(MEZUA3(IZK),I=1,3),MEZUA2(IZK),YMAX
C 108 FORMAT(1H // 15X,'1144 FASE 1',5X,'2255 FASE 2',5X,'3366 FASE 3',
C    1 7X,'FASEKO TENTSIOEN BALIORIK HANDIENA:',F13.3,' KV')
  108 FORMAT(1H // 15X,'1144 ',A5,' 1',5X,'2255 ',A5,' 2',5X,'3366 ',A5
     1,' 3',7X,A35,F13.3,' KV')
      RETURN
      END
      SUBROUTINE RES(ND,NC,N1,N2,N3,XA,XG,NG,RF,RI,BC,AD,ES,ER,LOC)
      COMPLEX RF(ND,ND),RI(ND,ND),BC(ND,ND),AD(ND,ND),ES(ND),ER(ND)
      REAL VOLT(6,236)
      INTEGER LOC(ND)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      NT=ND
      N6=NC*6
      N5=NC*3-3
      N4=N6-3
      IF(NG.GT.236)NG=236
      CALL IRAKUR(NC,ER,LOC,NT)
      WRITE(IDF,199)(J,J=1,N6)
  199 FORMAT(1H1// 3X,7HXL(OHM),3X,12(3X,2HES,I2,3X)// )
      XA=XA-XG
      XH=XA
      DO 1000 I=1,NG
      XA=XA+XG
      DO 1100 I1=1,NC*3
      DO 1100 J1=1,NC*3
      I2=I1+NC*3
      J2=J1+NC*3
      BC(I1,J1)=(0.,0.)
      BC(I2,J2)=(0.,0.)
      BC(I1,J2)=(0.,0.)
      BC(I2,J1)=(0.,0.)
      IF(I1.NE.J1) GO TO 1100
      BC(I1,J1)=(1.,0.)
      BC(I2,J2)=(1.,0.)
      IF(N2.EQ.0.OR.J1.GT.N5)BC(I2,J1)=CMPLX(0.,-1./XA)
 1100 CONTINUE
C     N1 ESATEN DU ERREAKTANTZIAK NON DAUDEN KOKATUTA
      IF(N1.GE.3)CALL BITRIC(RF,N6,N6,BC,N6,N6,RI,NT)
      IF(N1.EQ.2)CALL BITRIC(BC,N6,N6,RF,N6,N6,RI,NT)
      IF(N1.NE.4)GO TO 1001
      CALL BITRIC(BC,N6,N6,RI,N6,N6,AD,NT)
      DO 1200 I1=1,N6
      DO 1200 J1=1,N6
 1200 RI(I1,J1)=AD(I1,J1)
 1001 CALL EXCHAN(RI,N6,LOC,NT,*999,2)
      CALL BITRIC(RI,N6,N6,ER,N6,1,ES,NT)
      WRITE(IDF,198)XA,(CABS(ES(J)),J=1,N6)
  198 FORMAT(1H ,F10.3,1X,12(1PE10.2))
      DO 1300 I1=1,3
      VOLT(I1,I)=CABS(ES(I1+N5))
 1300 VOLT(I1+3,I)=CABS(ES(I1+N4))
 1000 CONTINUE
C     N3 ESATEN DU IRTEERA GRAFIKOA NOLAKOA DEN
      CALL REGRA(VOLT,NG,XG,XH,N3)
  999 RETURN
      END
      SUBROUTINE ZY(L,N,FNCY,RHO,Z,Y,DX,DY,RAD,RMG,IN,D,ID,R,G,KR,W,KW)
C***** SUBROUTINE ONEN BIDEZ LINE BATEN (R+JX) INPEDANTZEAK ETA   ******
C***** (G+JB) ADMITANTZEAK KALKULATZEN DIRA                       ******
C***** "Transmission Line Reference Book 345 kV and Above"        ******
C***** liburuaren 95. orrialdean esandakoaren arabera             ******
C***** Electric Power Research Institute, Palo Alto CA94304, 1975 ******
      REAL FNCY,DX(L),DY(L),RAD(L),RMG(L),D(L),R(L),G(L),K,MUCY
      INTEGER IN(L),ID(L),KR(L)
      COMPLEX Z(L,L),Y(L,L),W(L,L)
      FK(A,F,RLOC)=2.80E-3*A*SQRT(F/RLOC)
      FP(K,T,P)=P/8.-K*COS(T)/3./SQRT(2.)+K**2 *(COS(2.*T)*(0.6728+ALOG(
     12./K))+T*SIN(2.*T))/16.+K**3 *COS(3.*T)/45./SQRT(2.)-P*K**4 *COS(4
     2.*T)/1536.
      FQ(K,T,P)=ALOG(2./K)/2.-0.0386+K*COS(T)/3./SQRT(2.)-P*K**2 *COS(2.
     1*T)/64.+K**3 *COS(3.*T)/45./SQRT(2.)-K**4 *(T*SIN(4.*T)+COS(4.*T)*
     2(ALOG(2./K)+1.0895))/384.
      PI=3.1415926
      C=3.E8
      OMEG=2.*PI*FNCY
      MUCY=2.E-4*OMEG
C *** MUCY=MU ZATI (BI BIDER PI) BIDER OMEGA (MU=4*PI*E-7, MUCY=MU*FNCY)
C *** MUCY=MU ENTRE DOS PI POR OMEGA         (MU=4*PI*E-7, MUCY=MU*FNCY)
      EPCY=2.E-10*C**2/OMEG
C *** EPCY=BAT ZATI (BI BIDER PI EPSILON OMEGA)(EPS=1/MU/C**2=E-9/36/PI)
C *** EPCY=UNO ENTRE DOS PI EPSILON OMEGA      (EPS=1/MU/C**2=E-9/36/PI)
        DO 1 I=1,N
          G(I)=0.  !C-20041122
        R(I)=R(I)/IN(I)
        RMG(I)=RMG(I)/100.
        RAD(I)=RAD(I)/100.
        RG=GMR(IN(I),RMG(I),D(I),ID(I))
        RR=GMR(IN(I),RAD(I),D(I),ID(I))
        DPMM=2.*DY(I)
        TETA=0.
        K=FK(DPMM,FNCY,RHO)
        P=FP(K,TETA,PI)*2.*MUCY
        Q=FQ(K,TETA,PI)*2.*MUCY
        X=ALOG(DPMM/RG)*MUCY
        Z(I,I)=CMPLX(R(I)+P,X+Q)
        X=ALOG(DPMM/RR)*EPCY
        Y(I,I)=CMPLX(0.,-X)
        W(I,I)=Y(I,I)
        J1=I+1
        IF(N-J1)1,15,15
   15     DO 16 J=J1,N
          DPMN=DY(I)+DY(J)
          DMN =DY(I)-DY(J)
          HMN =DX(I)-DX(J)
          DPMN=SQRT(DPMN**2+HMN**2)
          DMN =SQRT(DMN **2+HMN**2)
          TETA=ASIN(HMN/DPMN)
          K=FK(DPMN,FNCY,RHO)
          P=FP(K,TETA,PI)*2.*MUCY
          Q=FQ(K,TETA,PI)*2.*MUCY
          X=ALOG(DPMN/DMN)*MUCY
          Z(I,J)=CMPLX(P,X+Q)
          Z(J,I)=Z(I,J)
          X=ALOG(DPMN/DMN)*EPCY
          Y(I,J)=CMPLX(0.,-X)
          W(I,J)=Y(I,J)
          W(J,I)=W(I,J)
   16     Y(J,I)=Y(I,J)
    1   CONTINUE
      M=L
      CALL LABURT(W,N,N3,M,KR,KW,ID)
      CALL EXCHAN(Y,N,ID,M,*99,1)
      CALL EXCHAN(W,N3,ID,M,*99,1)
        DO 2 I=1,N
        W(I,I)=CMPLX(G(I),0.)+W(I,I)
    2   Y(I,I)=CMPLX(G(I),0.)+Y(I,I)
   99 RETURN
      END
C --- *** --- *** --- *** --- *** --- *** --- *** --- *** --- *** ---
      SUBROUTINE EREPRO(CHK)
C     AIREKO LINEA ELEKTRIKOEK SORTUTAKO
C     EREMU ELEKTROMAGNETIKOAK KALKULATZEKO PROGRAMA.
C     PARALELOAN JARRITAKO 10 LINEA ARTE KONTUAN HARTU AHAL DITU.
C$      FORTY   MAP,XREF,FDS                                             EREMU
CEREMU
      PARAMETER (N1=32)
      DIMENSION V(N1),PHI(N1),RA(N1),RMG(N1),D(N1),IDL(N1),IN(N1),X(N1),
     &Y1(N1),Y2(N1),PC(N1,N1),VR(N1),VI(N1),XR(N1),XI(N1),FI(N1),JCK(N1)
      REAL FNCY
      INTEGER POPT
      CHARACTER  IZENA*56,EGUNA*10,  EUSK*5,EROA*15,  CHK*4
      DOUBLE PRECISION EXPE
      EQUIVALENCE(RMG(1),PC(1,N1))
      COMMON/BI/PC,VR,VI,RA,D,IN,IDL
      COMMON/IRU/XA,XG,NX,YA,YG,NY,X,Y1
      COMMON /IZENBN/ IZK,EXPE,NOZ,NNN,KTKTKT,TL,MOZ,NFL
      COMMON /IZENBC/ EGUNA,IZENA,EUSK,EROA
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      CHARACTER     IZBR21(2)*7, IZBR22(2)*28,IZBR23(2)*15,IZBR24(2)*33
     4,IZBR31(2)*22,IZBR32(2)*18,IZBR33(2)*54,IZBR34(2)*16,IZBR35(2)*8
     5,IZBR36(2)*54,IZBR37(2)*16
      COMMON /HZENBU/IZBR21,IZBR22,IZBR23,IZBR24
     4       ,IZBR31,IZBR32,IZBR33,IZBR34,IZBR35,IZBR36,IZBR37
      CHARACTER IZBE21(2)*90,IZBE22(2)*46
      COMMON /HZEPRO/IZBE21,IZBE22
   99 CALL HIRUAT(IRF,CHK,*98,*9999,*9998)
   98 READ(IRF,'(F4.0,A56,F10.2,A10)',END=9998)TL,IZENA,EXPE,EGUNA
      READ(IRF,599)L2,FNCY,RHO,POPT,NA,IE,IM,L,XA,XG,NX,YA,YG,NY
  599 FORMAT(I3,F9.0,F7.2,5I2,2X,2(F6.2,F4.2,I3))
                         print *,'EREMU   ',izena
C     L2   = EROALEEN KOPURUA
C     FNCY = LINEAREN FREKUENTZIA
C     RHO  = LINE TARTE HONEN LURRAREN ERRESISTIBITATEA
C     POPT <1 => KONF. ORIGINALAREN MAXWELL-EN POTENTZIAL KOEF. ATARATZEN DIRA
C          =1 => KONF. ORIGINALAREN MAXWELL-EN POTENTZIAL KOEF. ATARATZEN DIRA
C          =2 => KONF. BALIOKIDEAREN M-EN POT. KOEF. ATARATZEN DIRA
C          >2 => BI KONFIGURAZIOEN M-EN POT. KOEFIZIENTEAK ATARATZEN DIRA
C     NA   =1 => ZARATA AKUSTIKOARI BURUZKO DATUAK ATARATZEN DIRA
C     IE   >0 => EREMU ELEKTRIKOAREN EMAITZAK ATARATZEN DIRA
C     IM   >0 => EREMU MAGNETIKOAREN EMAITZAK ATARATZEN DIRA
C     L    =
C     XA   = EREMUAK KALKULATZEKO HASIERAKO X KOORDENADA (M-TAN)
C     XG   = EREMUAK KALKULATZEKO X KOORDENADARENTZAKO GEHIGARRIA (M-TAN)
C     NX   = ZENBAT BIDER GEHITUKO DEN XG ( <101 )
C     YA   = EREMUAK KALKULATZEKO HASIERAKO Y KOORDENADA (M-TAN)
C     YG   = EREMUAK KALKULATZEKO Y KOORDENADARENTZAKO GEHIGARRIA (M-TAN)
C     NY   = ZENBAT BIDER GEHITUKO DEN YG ( <11 )
      READ(IRF,598)(RMG(I),RA(I),IN(I),D(I),IDL(I),X(I),Y1(I),Y2(I)
     &,XI(I),FI(I),V(I),PHI(I),I=1,L2)
  598 FORMAT(2F7.4,I2,F7.4,I2,3F7.3,9X,F5.3,F6.2,F7.3,F6.2)
C     RMG = EROALEAREN "BATEZBESTEKO ERRADIO GEOMETRIKOA" CM-TAN
C     RA  = EROALEAREN ERRADIOA CM-TAN
C     IN  = EROALE SORTAREN AZPIEROALE KOPURUA
C     D   = AZPIEROALEEN ARTEKO DISTANTZIA EDO EROALE SORTAREN DIAMETROA M-TAN
C     IDL =0 => AURREKO DATUA AZPIEROALEN ARTEKO DISTANTZIA DA
C         =1 => AURREKO DATUA EROALE SORTAREN DIAMETROA DA
C     X   = EROALEAREN X KOORDENATUA M-TAN
C     Y1  = EROALEAREN Y KOORDENATUA M-TAN EUSKARRIEN ONDOAN
C     Y2  = EROALEAREN Y KOORDENATUA M-TAN EUSKARRIEN ARTEKO TARTEAREN ERDIAN
C     XI  = INTENTSITATEAREN BALIOA KA-TAN
C     FI  = INTENTSITATEAREN FASEA GRADUTAN
C     V   = TENTSIOAREN BALIOA KV-TAN
C     PHI = TENTSIOAREN FASEA GRADUTAN
                  KTKTKT=20
      CALL IZENBU
C     IZBE21='SARRERAKO DATUAK, AZTERTUTAKO LINEAREN EZAUGARRIAK'
      WRITE(IDF,696)IZBE21(IZK),IZBE22(IZK)
     1,L2,FNCY,RHO,POPT,NA,IE,IM,L,XA,XG,NX,YA,YG,NY
  696 FORMAT(//8X,A//8X,A,13X,'XA    XG   NX     YA    YG   NY'/
     & 9X,I2,F11.0,F9.2,8X,5I3, 8X,2(F8.2,F6.2,I4))
 
      WRITE(IDF,697)(I,RMG(I),RA(I),IN(I),100.*D(I),IDL(I),X(I),Y1(I),
     1Y2(I),XI(I),FI(I),V(I),PHI(I),I=1,L2)
C 697 FORMAT(1H /8X,'COND  GMR  RAD  ZEN BUN D ZER      X       Y1     '
C    1,' Y2      INT KA  FAS G    TEN KV  PHAS G'//
C    2,(9X,I2,F6.2,F5.2,I4,F7.2,I3,2X,3F8.2,2(F10.2,F8.2)))
  697 FORMAT(1H /8X,'EROA  BGE  ERR  ZEN SOR.D ZER      X       Y1     '
     1,' Y2      INT KA  FAS G    TEN KV  PHAS G'//
     2,(9X,I2,F6.2,F5.2,I4,F7.2,I3,2X,3F8.2,2(F10.2,F8.2)))
      IF(RHO)95,95,97
   97 IF(FNCY)95,95,96
   96 DK=660.*SQRT(RHO/FNCY)
      GO TO 94
   95 DK=800.
   94 CONTINUE
          DO 1 I=1,L2
    1     Y1(I)=Y1(I)+(Y2(I)-Y1(I))*2./3.
        DO 2 I=1,L2
        PHI(I)=PHI(I)*0.017453293
        V(I)=1000.*V(I)/1.7320508
        VR(I)=V(I)*COS(PHI(I))
        VI(I)=V(I)*SIN(PHI(I))
        FI(I)=FI(I)*0.017453293
        XI(I)=1000.*XI(I)
        XR(I)=XI(I)*COS(FI(I))
        XI(I)=XI(I)*SIN(FI(I))
        RA (I)=RA (I)/100.
    2   CONTINUE
      KOPT=0
      IF(POPT.LT.0)POPT=0     !C-20131026
      IF(POPT.GT.3)POPT=3     !C-20131026
      POPT=POPT+1
  100 CALL PMAXW(L2,X,Y1,IN,RA,D,IDL,PC)
      GO TO(203,202,201,200),POPT
  200 IF(KOPT)202,202,201
  201 CALL GITRIX(PC,L2,V,L4)
      KOPT=1                 !C-20131026
      L2=L4
  202 CALL IDATZM(PC,L2,KOPT,1,IZK)
  203 CALL INTRIX(PC,L2,JCK,N1,*99)
      IF(POPT-1)204,205,204
  204 CALL IDATZM(PC,L2,KOPT,2,IZK)
  205 CALL IDATZG(L2,KOPT,IZK)
      IF(IE.GT.0.OR.IM.GT.0.OR.L.GT.0)
     &CALL EREMU(L2,KOPT,XR,XI,DK,IE,IM,L)
      IF(NA.GT.0)CALL ANOISE(L2,KOPT,V,PHI,RA,IN,NA,IZK)
      IF(KOPT) 99,208, 99
  208 KOPT=1
      IF(POPT-3)99,99,100
 9998 CLOSE(UNIT=IRF)
      CHK=' '
 9999 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             ANOISE
      SUBROUTINE ANOISE(N,KOPT,V,PHI,RA,IN,AN,IZK)
C  *********************************************************************
C  *  SUBRUTINA HONETAN LINEEK SORTUTAKO PRESIO AKUSTIKOAREN ZEHARKAKO *
C  *  PROFILA KALKULATZEN DA,   ''E.P.R.I.''K 1975EAN ARGITARATUTAKO   *
C  *  '' TRANSMISSION  LINE  REFERENCE  BOOK  345  KV  AND  ABOVE ''   *
C  *  LIBURUAREN 6. KAPITULUAN (AUDIBLE NOISE) ESANDAKOAREN ARABERA    *
C  *********************************************************************
      PARAMETER (N1=32)
      REAL V(N1),RA(N1),PHI(N1),X(N1),Y(N1),GRS(N1),Z(N1),GRM(N1)
      INTEGER AN,IN(N1)
      COMPLEX P
      REAL LONG,JHR,JWC,NM2HR,NM2WC,NM2PT,LANDA,LANDA1,LANDA2
      CHARACTER ANOI11(2)*4, ANOI12(2)*72,ANOI21(2)*47,ANOI22(2)*54
     1,ANOI23(2)*72,ANOI24(2)*26,ANOI25(2)*25,ANOI26(2)*23
      COMMON /HZOISE/ ANOI11,ANOI12,ANOI21,ANOI22,ANOI23,ANOI24,ANOI25
     1,ANOI26
      EQUIVALENCE(GRM(1),Z(1))
      COMMON/BAT/A1(N1),A2(N1)
      COMMON/IRU/XA,XG,NX,YA,YG,NY,X,Y
      COMMON/LAU/GRM
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      CHARACTER*25 KONFOR(2),KONFBA(2),POTGR1(2)*32,POTGR2(2)*60
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2
      DB(XNM2)=20.*ALOG10(XNM2/2.E-5)
      AB(AX)=10.**(AX/10.)
      IF(AN.EQ.0) GO TO 99
      ATE=-0.0075
      ATE=-0.0075
      LANDA=3.43
      DC=20.5
      XK1=0.
      XK2=1.
      LONG=5000.
      L1=100
      L2=10
      XM=0.01
C     YM=0.01
      ZM=10.
      KONTR=0
      PI=3.14159265
      IF(NX.GT.L1.OR.NY.GT.L2.OR.XG.LT.XM.OR.YG.LT.XM)GO TO 98
        DO 100 I=1,N
        IF(IN(I)-2)101,102,103
  101   XKN=5.6
        GO TO 104
  102   XKN=1.8
        GO TO 104
  103   XKN=1.0
  104   A1(I)=46.4-665.0/GRM(I)
        A1(I)=AB(A1(I))
        A1(I)=IN(I)**2*(200.*RA(I)/3.8)**4.4*A1(I)*XKN
        A2(I)=53.5-505.5/GRM(I)
        C1=10.6-0.41/(2.*RA(I))
        A2(I)=A2(I)+C1
        C1=47.4-1000./(IN(I)+15)
        C2=24.1-390.0/(IN(I)+10)
        C=(C2-C1)/2.28
        C=C1+C*(200.*RA(I)-2.35)
        A2(I)=A2(I)+C
        A2(I)=AB(A2(I))
        GRS(I)=(12.5*200.*RA(I)-4.57)/(200.*RA(I)-1.07)
        IF(IN(I)-4)106,106,105
  105   GRS(I)=GRS(I)/(1.+0.027*(IN(I)-4))
  106   GRS(I)=GRM(I)/GRS(I)
        U=10.*(GRS(I)-0.8)
        Z(I)=(63.4*U**2+1.87*U**3-1.15*U**4)/1000.
  100   CONTINUE
      Y2=YA-YG
        DO 200 K=1,NY+1
        CALL IZENBU
        Y2=Y2+YG
        IF(Y2)200,201,201
  201   IF(KOPT)202,202,203
C 997 FORMAT(1H /51X,22HKONFIGURAZIO ORIGINALA)
  202   WRITE(IDF,'(1H /51X,A)')KONFOR(IZK)
        GO TO 204
C 996 FORMAT(1H /51X,23HKONFIGURAZIO BALIOKIDEA)
  203   WRITE(IDF,'(1H /51X,A)')KONFBA(IZK)
  204     DO 300 I=1,N
          IF(ABS(V(I))-ZM)301,301,300
  301     R=RA(I)*100.
          G=GRS(I)
          IF(R.GE.1.0.AND.R.LE.3.0.AND.G.LE.1.4.AND.G.GE.0.8)GO TO 300
C 998 FORMAT(/I3,' EROALEA (LUR KABLEA) EZ DA KONTUAN HARTZEN'
C    1,' EROALE HEZEKO KASUAN ')
          WRITE(IDF,'(/A,I3,A)')TRIM(ANOI11(IZK)),I,TRIM(ANOI12(IZK))
  300     CONTINUE
 
C 995 FORMAT(/16X,'PRESIO AKUSTIKOAREN ZEHARKAKO PROFILA, LURRETIK',F6.2
C    1,' M-RA N/M2 ETA 20 MICRON/M2-REN GAINEKO DB-TAN' // 2X,5HDIST.,8X
C    2,'*******  BANDA ZABALEKO ZARATA ALEATORIOAREN MAILA  *******',12X
C    3,'100 HZ-EKO ZARATAREN MAILA' / 15X,'** EURITE GOGORRAK ***',14X,
C    4'**** EROALE HEZEAK ****', 13X,'**** EURITE GOGORRAK ****' /
C    5 18X,'N/M2',9X,5HDB(A),18X,4HN/M2,9X,5HDB(A),19X,4HN/M2,9X,2HDB)
        WRITE(IDF,995)ANOI21(IZK),Y2,ANOI22(IZK),ANOI23(IZK),ANOI24(IZK)
     1, ANOI25(IZK),ANOI26(IZK),ANOI25(IZK)
  995 FORMAT(/16X,A,F6.2,A//2X,A,12X,A/15X,A,11X,A,13X,A/
     5 19X,'N/M2',9X,5HDB(A),18X,4HN/M2,9X,5HDB(A),19X,4HN/M2,10X,2HDB)
 
        ILK=0     !Datuekin idatzitako lerro kopurua
        X2=XA-XG
         DO 400 J=1,NX+1
         JHR=0.
         JWC=0.
         P=(0.,0.)
         X2=X2+XG
         IF(MOD(ILK,6).EQ.0)THEN
          WRITE(IDF,'(1X)')
          ILK=ILK+1
         ENDIF
         ILK=ILK+1   !Idatzitako lerroen kopurua
          DO 500 I=1,N
          XD=X2-X(I)
          YH1=Y2-Y(I)
          YH2=Y2+Y(I)
          IF(ABS(XD)-XM)501,502,502
  501     IF(ABS(YH1)-XM)403,502,502
  502     R=RA(I)*100.
          G=GRS(I)
          R1=SQRT(XD*XD+YH1*YH1)
          R2=SQRT(XD*XD+YH2*YH2)
          IF(LONG-1000.)505,505,506
  505     HM1=(ATAN2(LONG,2.*R1)*EXP(ATE*R1)/R1 + XK1*ATAN2(LONG,2.*R2)*
     1    EXP(ATE*R2)/R2)/(2.*PI)
          GO TO 507
  506     HM1=(EXP(ATE*R1)/R1+XK1*EXP(ATE*R2)/R2)/4.
  507     JHR=JHR+A1(I)*HM1
          COEF=DC*A2(I)/SQRT(2.)/PI
          LANDA1=0.-PHI(I)-2.*PI*R1/LANDA
          LANDA2=0.-PHI(I)-2.*PI*R2/LANDA
          P1=COEF*(COS(LANDA1)/R1+XK2*COS(LANDA2)/R2)
          P2=COEF*(SIN(LANDA1)/R1+XK2*SIN(LANDA2)/R2)
          P=P+CMPLX(P1,P2)
          IF(R.GE.1..AND.R.LE.3..AND.G.LE.1.4.AND.G.GE..8) GO TO 504
          IF(ABS(V(I))-ZM)500,500,503
  503     KONTR=1
  504     JWC=JWC+A1(I)*HM1*Z(I)
  500     CONTINUE
         NM2HR=DC*SQRT(JHR)*1.E-3
         NM2WC=DC*SQRT(JWC)*1.E-3
         NM2PT=CABS(P)*1.E-3
         DBAHR=DB(NM2HR)
         DBAWC=DB(NM2WC)
         DBAPT=DB(NM2PT)
         IF(KONTR)402,401,402
  401 WRITE(IDF,994)X2,NM2HR,DBAHR,NM2WC,DBAWC,NM2PT,DBAPT
  994 FORMAT(1H , F7.2,8X,1PE9.2,0PF12.2,15X,1PE9.2,0PF12.2,16X,1PE9.2
     1,0PF12.2,1X)
         GO TO 400
  402 WRITE(IDF,993)X2,NM2HR,DBAHR,NM2PT,DBAPT
C 993 FORMAT(1H / F7.2,8X,1PE9.2,0PF12.2,10X,
C    1'NOT ENOUGH STATISTICAL DATA',13X,1PE9.2,0PF12.2 )
  993 FORMAT(1H , F7.2,8X,1PE9.2,0PF12.2,10X,
     1'DATU ESTATISTIKO GUTXIEGI',15X,1PE9.2,0PF12.2 )
         GO TO 400
  403    WRITE(IDF,992)X2,X2,Y2
C 992 FORMAT(1H / F7.2,25X,' RESULTS FOR X=',F6.2,' E Y=',F6.2,'
C    1 ARE NOT SIGNIFICANT')
  992 FORMAT(1H , F7.2,25X,' X=',F6.2,' ETA Y=',F6.2,
     1' TOKIRAKO EMAITZAK EZ DIRA ESANGURATSUAK')
  400    CONTINUE
  200   CONTINUE
      GO TO 99
   98 WRITE(IDF,999)
  999 FORMAT(1H1///10HANOISE'N, ,'MUGAEI JAGOKEZAN (DAGOKIEZAN, DAGOZKIE
     1N) BALDINTZAK EZ DIRA BETETZEN')
   99 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             ATARA
      SUBROUTINE ATARA (MULTZO,NEURRI)
      PARAMETER (N1=32)
      REAL MULTZO(N1,N1)
      CHARACTER *20 AX
      CHARACTER *25 AY
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      INTEGER ASIERA,AMAIA
      ASIERA=1
      AMAIA=6
   33 IF(NEURRI-AMAIA)1,1,2
    2 WRITE(IDF,999)(K,K=ASIERA,AMAIA)
  999 FORMAT(1H //3X,6(15X,I2))
      WRITE(IDF,997)(K,(MULTZO(K,I),I=ASIERA,AMAIA),K=1,NEURRI)
  997 FORMAT(1H /(7X,I2,6(3X,E14.7)))
      ASIERA=ASIERA+6
      AMAIA=AMAIA+6
      GO TO 33
    1 NBARRI=NEURRI-ASIERA+1
      WRITE(AX,998)NBARRI
  998 FORMAT(9H(1H //3X,,I1,9H(15X,I2)))
      WRITE(IDF,AX)(K,K=ASIERA,NEURRI)
      WRITE(AY,996)NBARRI
  996 FORMAT(12H(1H /(7X,I2,,I1,12H(3X,E14.7))))
      WRITE(IDF,AY)(K,(MULTZO(K,I),I=ASIERA,NEURRI),K=1,NEURRI)
      RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             BITRIX
      SUBROUTINE BITRIX(A,NFA,NCA,B,NFB,NCB,P)
      PARAMETER (N1=32)
      DIMENSION A(N1,N1),B(N1,1),P(N1,1)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
  888 FORMAT(1H ,5X,'MULTZO ONEIK EZIN DIRA BIDERTU')
      IF(NCA-NFB)3,1,3
    1 CONTINUE
        DO 2 J=1,NCB
          DO 21 I=1,NFA
          P(I,J)=0.
            DO 211 K=1,NCA
  211       P(I,J)=P(I,J)+A(I,K)*B(K,J)
   21     CONTINUE
    2   CONTINUE
      RETURN
    3 WRITE(IDT,888)
      STOP 1002
      END
C$      FORTY   MAP,XREF,FDS                                             EREMU
      SUBROUTINE EREMU(N,KOPT,XIR,XII,DEK,IELEC,IMAG,L)
C *** ******************************************************************
C     SUBRUTINA HONETAN LINEA BATEK SORTUTAKO EREMU ELEKTRIKO ETA
C     MAGNETIKOA KALKULATZEN DIRA,   D.W. DENO-REN 'TRANSMISSION
C     LINE FIELDS' (IEEE TRANSACTIONS VOL.PAS-95 NO.5 SEP/OCT 1976) 
C     ARTIKULUAREN ETA WG 36-01(EFFECTS OF CORONA AND FIELD) LAN TALDEAK
C     SC 36(INTERFERENCE) KOMITEAREN IZENEAN AURKEZTUTAKO
C     36-07(CIGRE 1978) TXOSTENAREN ARABERA
C *** ******************************************************************
      PARAMETER (N1=32)
      DOUBLE PRECISION EXPE
      COMMON/BAT/QR,QI
      COMMON/IRU/XASI,XGEI,NXG,YASI,YGEI,NYG,X,Y
      COMMON /IZENBN/ IZK,EXPE,NOZ,NNNN,KT,TL,MOZ,NFL
      REAL X(N1),Y(N1),QR(N1),QI(N1),XIR(N1),XII(N1)
      REAL MUPI
      COMPLEX POTEN
      CHARACTER*25 KONFOR(2),KONFBA(2)
      COMMON /HIZKON/ KONFOR,KONFBA
      CHARACTER*25  POTEL1(2)*48,POTEL2(2)*50
     1,ERMAG1(2)*28,ERMAG2(2)*15,ERMAG3(2)*122,ERMAG4(2)*19,ERMAG5(2)*19
     2,ERELE1(2)*25,ERELE2(2)*15,ERELE3(2)*121,ERELE4(2)*19,ERELE5(2)*19
     3,MEZU21(2)*15,MEZU22(2)*7, MEZU23(2)*36, MEZU24(2)*74
      COMMON /HZREMU/ POTEL1,POTEL2
     1       ,ERMAG1,ERMAG2,ERMAG3,ERMAG4,ERMAG5
     2       ,ERELE1,ERELE2,ERELE3,ERELE4,ERELE5
     3       ,MEZU21,MEZU22,MEZU23,MEZU24
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      FUNK(R,S,T,U,A)=SQRT((R*SIN(A)+S*COS(A))**2 +(T*SIN(A)+U*COS(A))**
     12 )
      BAT(A,B)=A/(A*A+B*B)
 
      LIM1=100
      LIM2=10
      XMUGA=0.01
      YMUGA=0.01
C20041023      ZMG=1.E-38
      ZMG=1.E-37
      RD=57.295779
      EPPI=18.E6
C     EPPI=BAT ZATI (BI BIDER PI BIDER EPSILON)
C     =1./(2.*PI*EPS)=18.E9(S.I.) C=3.E8 BADA
      MUPI=2.E-7
C     MUPI=MU ZATI (BI BIDER PI)=MU/(2.*PI)=2.E-7(S.I.)
      IF(IELEC)11,11,12
   11 IF(IMAG)99,99,12
   12 IF(NXG-LIM1)14,14,13
   13 WRITE(IDF,613)
  613 FORMAT(1H1///10X,'BIDEA ZABALEGIA DA, NXG 100 BAINO AUNDIAGOA DA')
      GO TO 99
   14 IF(NYG-LIM2)16,16,15
   15 WRITE(IDF,615)
  615 FORMAT(1H1///10X,'BIDEA GARAIEGIA DA, NYG 10 BAINO AUNDIAGOA DA')
      GO TO 99
   16 IF(XGEI-0.01)17,18,18
   17 WRITE(IDF,617)
  617 FORMAT(1H1///10X,'XGEI 0.01 BAINO TXIKIAGOA DA')
      GO TO 99
   18 IF(YGEI-0.01)19,20,20
   19 WRITE(IDF,619)
  619 FORMAT(1H1///10X,'YGEI 0.01 BAINO TXIKIAGOA DA')
      GO TO 99
   20 CONTINUE
       DO 100 KB=1,2
       IF(KB-1)101,101,102
  101  IF(IELEC)103,103,102
  102  YA=YASI-YGEI
C*200   DO 200 EREMU ELEKTRIKO EDO INDUKZIO MAGNETIKOAREN KALKULUA
C       AZTERTU GURA DIREN SESTRA GUZTIETAN
        DO 200 K=1,NYG+1
        CALL IZENBU
        ILK=0     !Datuekin idatzitako lerroen kopurua
        IF(KOPT)201,201,203
  201   WRITE(IDF,995)KONFOR(IZK)     !Konfigurazio originala
  995 FORMAT(1H /51X,A)
        GO TO 204
  203   WRITE(IDF,995)KONFBA(IZK)     !Konfigurazio baliokidea
  204   YA=YA+YGEI
        IF(KB-1)208,208,206
  206   WRITE(IDF,988)TRIM(ERMAG1(IZK)),YA,ERMAG2(IZK)
     1  ,ERMAG3(IZK),(ERMAG4(IZK),I=1,2),(ERMAG5(IZK),I=1,4)
  988   FORMAT(1H /44X,A,F6.2,A///3X,A/8X,2(A),4(A))
        GO TO 210
  208   WRITE(IDF,989)TRIM(ERELE1(IZK)),YA,ERELE2(IZK)
     1  ,ERELE3(IZK),(ERELE4(IZK),I=1,2),(ERELE5(IZK),I=1,4)
  989   FORMAT(1H /44X,A,F6.2,A///2X,A/1H ,8X,2(A),4(A))
  210   XA=XASI-XGEI
C*500    DO 500 E.ELEK. EDO I.MAG.-AREN KALKULUA SESTRA BAKOITZEAN
         DO 500 I=1,NXG+1
         GXR=0.
         GXI=0.
         GYR=0.
         GYI=0.
         XA=XA+XGEI
C*600     DO 600 E.ELEK. EDO I.MAG.-AREN KALKULUA
C         SESTRA BATEKO PUNTUETAN
          DO 600 J=1,N
          XD=XA-X(J)
          YH1=YA-Y(J)
          IF(KB-1)602,602,601
  601     YH2=YA+DEK
          GO TO 603
  602     YH2=YA+Y(J)
  603     CONTINUE
          IF(ABS(XD).GE.XMUGA)GO TO 606
          IF(ABS(YH1).LT.XMUGA.OR.ABS(YH2).LT.XMUGA)GO TO 504
  606     S11=BAT(XD,YH1)
          S12=BAT(XD,YH2)
          S21=BAT(YH1,XD)
          S22=BAT(YH2,XD)
          IF(KB-1)608,608,607
  607     XM=(S21-S22)*MUPI
          YM=(S12-S11)*MUPI
          GXR=GXR+XIR(J)*XM
          GXI=GXI+XII(J)*XM
          GYR=GYR+XIR(J)*YM
          GYI=GYI+XII(J)*YM
          GO TO 600
  608     XM=(S11-S12)*EPPI
          YM=(S21-S22)*EPPI
          GXR=GXR+QR(J)*XM
          GXI=GXI+QI(J)*XM
          GYR=GYR+QR(J)*YM
          GYI=GYI+QI(J)*YM
  600     CONTINUE
C        ELIPSE BAKOITZAREN ARDATZ NAGUSIEN KALKULUA
         IF(ABS(YA).LE.YMUGA.AND.KB.EQ.1)GO TO 591
         ECUA=GYR*GXR+GYI*GXI
         IF(ABS(ECUA).LE.ZMG)GO TO 599
         ECUB=GXI*GXI-GYI*GYI+GXR*GXR-GYR*GYR
         ECUC=-ECUA
         TANA1=(-ECUB+SQRT(ECUB**2 -4.*ECUA*ECUC))/(2.*ECUA)
         TANA2=(-ECUB-SQRT(ECUB**2 -4.*ECUA*ECUC))/(2.*ECUA)
         A1=ATAN(TANA1)
         A2=ATAN(TANA2)
  592    FA1=FUNK(GYI,GXI,GYR,GXR,A1)
         FA2=FUNK(GYI,GXI,GYR,GXR,A2)
         A11=A1*57.295779
         A22=A2*57.295779
         IF(FA1-FA2)502,501,501
  501    F1=FA1*SQRT(2.)
         F2=FA2*SQRT(2.)
         A1=A11
         A2=A22
         GO TO 503
  502    F1=FA2*SQRT(2.)
         F2=FA1*SQRT(2.)
         A1=A22
         A2=A11
C        ALDIUNEKO BALIOEN KALKULUA
  503    GXR2=-GXI
         GYR2=-GYI
         G1M=SQRT(GXR *GXR +GYR *GYR )*SQRT(2.)
         G2M=SQRT(GXR2*GXR2+GYR2*GYR2)*SQRT(2.)
         G1D=0.
         IF(ABS(GYR ).GE.ZMG.OR.ABS(GXR ).GE.ZMG)G1D=ATAN2(GYR ,GXR )*RD
         G2D=0.
         IF(ABS(GYR2).GE.ZMG.OR.ABS(GXR2).GE.ZMG)G2D=ATAN2(GYR2,GXR2)*RD
         GXABS=SQRT(GXR*GXR+GXI*GXI)
         GYABS=SQRT(GYR*GYR+GYI*GYI)
         GXFAS=0.
         IF(ABS(GXR).GE.ZMG.OR.ABS(GXI).GE.ZMG)GXFAS=ATAN2(GXI,GXR)*RD
         GYFAS=ATAN2(GYI,GYR)*57.295779
         IF(MOD(ILK,6).EQ.0)THEN
           WRITE(IDF,'(1X)')
          ILK=ILK+1
         ENDIF
         ILK=ILK+1   !Idatzitako lerroen kopurua
         WRITE(IDF,987)XA,GXABS,GXFAS,GYABS,GYFAS,F1,A1,F2,A2,G1M,G1D,
     1   G2M,G2D
  987 FORMAT(1H ,F7.2,6(2X,1PE9.2,0PF8.2))
         GO TO 500
C MEZU23='-RAKO EMAITZAK EZ DIRA ESANGURATSUAK'
  504    WRITE(IDF,985)XA,TRIM(MEZU21(IZK)),XA,TRIM(MEZU22(IZK)),YA,
     1   MEZU23(IZK)
  985 FORMAT(1H ,F7.2,2X,A,F6.2,A,F6.2,A)
         ILK=ILK+1
         GO TO 500
C MEZU24='BEGIRATU IA EROALEEN TENTSIOAK EDO INTENTSITATEAK SARTU DIREN'
  599    WRITE(IDF,984)MEZU24(IZK)
  984 FORMAT(1H ///10X,A///)
 
         IF(KB-1)103,103,99
  591    GXR=0.
         GXI=0.
         A1=90./57.295779
         A2= 0.
         GO TO 592
  500    CONTINUE
  200   CONTINUE
  103  IF(IMAG)99,99,100
  100  CONTINUE
   99 IF(L.LT.1.OR.NYG.GT.LIM1.OR.NXG.GT.LIM2)GO TO 98
      XA=XASI-XGEI
      DO 700 K=1,NXG+1
      IF(KOPT.GT.0)GO TO 701
      WRITE(IDF,995)KONFOR(IZK)     !Konf. originala
      GO TO 702
  701 WRITE(IDF,995)KONFBA(IZK)     !Konf. baliokidea
  702 XA=XA+XGEI
      WRITE(IDF,983)XA
C 983 FORMAT(1H /34X,'POTENCIAL ELECTRICO A DIFERENTES ALTURAS PARA X=',
C    &F6.2,' M'///36X,'ALTURA    VA.EF.POT.     A.FAS.POT.     VA.AC.POT
C    &.'/37X,'(M)',9X,'(KV)',11X,'(GS)',11X,'(KV)')
  983 FORMAT(1H /34X,'POTENTZIAL ELEKTRIKOAK GARIAERA EZBERDINETAN (X=',
     &F6.2,' M)'
     &///36X,'GOIERA    POT.BA.EF.     POT.FAS.A.     POT.UN.BA.'
     &/37X,'(M)',9X,'(KV)',11X,'(GS)',11X,'(KV)')
      YA=YASI-YGEI
      DO 700 I=1,NYG+1
      YA=YA+YGEI
      IF(YA.LT.0.)GO TO 700
      POTEN=(0.,0.)
      DO 800 J=1,N
      XD=XA-X(J)
      YH1=YA-Y(J)
      YH2=YA+Y(J)
      IF(ABS(XD).LT.XMUGA.AND.ABS(YH1).LT.XMUGA)GO TO 703
      YH1=SQRT(XD*XD+YH1*YH1)
      YH2=SQRT(XD*XD+YH2*YH2)
      XD=EPPI*QR(J)*ALOG(YH2/YH1)
      YH1=EPPI*QI(J)*ALOG(YH2/YH1)
  800 POTEN=POTEN+CMPLX(XD,YH1)
      XD=CABS(POTEN)
C     HURRENGO AGINDUA 1999.07.14.AN GEHITU ZEN, EMOITZAK EGOKITZEKO
      XD=XD*1000.
      IF(XD.LT.ZMG)GO TO 700
      YH1=ATAN2(AIMAG(POTEN),REAL(POTEN))
      YH2=XD*COS(YH1)*SQRT(2.)
      YH1=YH1*57.295779
      WRITE(IDF,982)YA,XD,YH1,YH2
  982 FORMAT(1H /35X,F6.2,3(6X,1PE9.2))
      GO TO 700
  703 WRITE(IDF,985)YA,XA,YA
  700 CONTINUE
   98 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             GITRIX
      SUBROUTINE GITRIX(Z,L2,V,L4)
      PARAMETER (N1=32)
      DIMENSION Z(N1,N1),V(N1)
      L3=L2+1
      L4=L2
        DO 1 L=1,L2
        L1=L3-L
        IF(V(L1)-10.E-6)2,2,4
    2   L4=L1-1
          DO 11 I=1,L4
            DO 111 K=1,L4
  111       Z(I,K)=Z(I,K)-Z(I,L1)/Z(L1,L1)*Z(L1,K)
   11     CONTINUE
    1   CONTINUE
    4 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             GMR
      FUNCTION GMR(IN,RA,D,IDL)
C     GMR = EROALE SORTA BATEN BATEZBESTEKO GEOMETRIKOKO ERRADIOA
C     IN = EROALE KOPURUA
C     RA = EROALE BATEN ERRADIOA
C     D = SORTAREN DIAMETROA EDO EROALEN ARTEKO DISTANTZIA
C     IDL = "D" ZER DEN ADIERAZTEKO ALDAGAIA
      IF(IN-2)1,22,2
    2 IF(IDL)21,21,22
   21 A=D/(2.*SIN(3.1415926/IN))
      GO TO 3
   22 A=D/2.
    3 GMR=(IN*RA*A**(IN-1))**(1./IN)
      GO TO 4
    1 GMR=RA
    4 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             IDATZG
      SUBROUTINE IDATZG(N,KOPT,IZK)
C     EROALEETAKO BATEZ BESTEKO POTENTZIAL GRADIENTEAK ETA
C     POTENTZIAL GRADIENTE HANDIENAK IDAZTEKO SUBRUTINA
      PARAMETER (N1=32)
      COMMON/BAT/QR,QI/BI/MULTZO,VR,VI,RA,D,IN,IDL
      COMMON/LAU/XGM(N1)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      CHARACTER*25 KONFOR(2),KONFBA(2),POTGR1(2)*32,POTGR2(2)*60
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2
      REAL MULTZO(N1,N1),VR(N1),VI(N1),QR(N1),QI(N1),RA(N1),D(N1)
      INTEGER IN(N1),IDL(N1)
      M=N
      CALL BITRIX(MULTZO,M,M,VR,M,1,QR)
      CALL BITRIX(MULTZO,M,M,VI,M,1,QI)
      CALL IZENBU
      IF(KOPT)1,1,2
    1 WRITE(IDF,601)KONFOR(IZK)
  601 FORMAT(1H ///49X,A)
      GO TO 3
    2 WRITE(IDF,602)KONFBA(IZK)
  602 FORMAT(1H ///49X,A)
C     POTGR1='POTENTZIAL GRADIENTEA EROALEETAN'
    3 WRITE(IDF,603)POTGR1(IZK),POTGR2(IZK)
  603 FORMAT(1H /45X,A32///27X,A60)
        DO 4 I=1,N
        Q=QR(I)*QR(I)+QI(I)*QI(I)
        Q=SQRT(Q)
        XGR=18.E4*Q/(RA(I)*IN(I))
        IF(IN(I)-1)2001,45,44
   45   XGRM=XGR
        GO TO 43
   44   IF(IDL(I))2001,41,42
   41   XGRM=XGR*(1+(IN(I)-1)*RA(I)*2./D(I)*SIN(3.1415927/IN(I)))
        GO TO 43
   42   XGRM=XGR*(1+(IN(I)-1)*RA(I)*2./D(I))
   43   XGM(I)=XGRM
        WRITE(IDF,643)I,Q,XGR,XGRM
  643 FORMAT(1H /(28X,I2,5X,E15.8,4X,2(2PE15.6,4X)))
    4   CONTINUE
 2001 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             IDATZM
      SUBROUTINE IDATZM(PC,L2,KOPT,K,IZK)
      PARAMETER (N1=32)
      DIMENSION PC(N1,N1)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      CHARACTER*25 KONFOR(2),KONFBA(2),POTGR1(2)*32,POTGR2(2)*60
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2
      IF(KOPT)1,1,2
    1 WRITE(IDF,601)KONFOR(IZK)
  601 FORMAT(1H ///30X,A)
      GO TO 3
    2 WRITE(IDF,602)KONFBA(IZK)
  602 FORMAT(1H ///30X,A)
    3 IF(K-1)7,4,5
    4 WRITE(IDF,604)
  604 FORMAT(1H /25X,36HMAXWELL-EN POTENTZIAL KOEFIZIENTEAK, /
     &32X,19H( METRO FARADIOKO ))
      GO TO 6
    5 WRITE(IDF,605)
  605 FORMAT(1H /36X,12HKAPAZITATEAK,/32X,19H( FARADIO METROKO ))
    6 L3=L2
      CALL ATARA(PC,L3)
    7 RETURN
      END
C$      FORTY   MAP,XREF,FDS                                             INTRIX
      SUBROUTINE INTRIX(A,KA,LCK,ND,*)
C     A MATRIZEA ALDERANTZIKATZEKO AZPIERRUTINA
C20041023      DIMENSION A(ND, 1),LCK(ND)
      DIMENSION A(ND, ND),LCK(ND)
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
        DO 1 I=1,KA
    1   LCK(I)=0
        DO 2 I=1,KA
        AMAX=0.
          DO 21 K=1,KA
          IF(AMAX-ABS(A(K,K)))211,21,21
  211     IF(LCK(K))212,212,21
  212     L=K
          AMAX=ABS(A(K,K))
   21     CONTINUE
        IF(ABS(AMAX)-10.E-6)1001,1001,22
   22   LCK(L)=1
        DIV=A(L,L)
        A(L,L)=1.
          DO 23 J=1,KA
   23     A(L,J)=A(L,J)/DIV
          DO 24 J=1,KA
          IF(L-J)241,24,241
  241     EAJ=A(J,L)
          A(J,L)=0.
            DO 242 K=1,KA
  242       A(J,K)=A(J,K)-A(L,K)*EAJ
   24     CONTINUE
    2   CONTINUE
      RETURN
 1001 WRITE(IDT,999)I,L,L
  999 FORMAT(1H /9X,'INTRIX,',I3,1H(,I2,1H,,I2,') ARDATZA TXIKIEGIA DA')
      RETURN 1
      END
C$      FORTY   MAP,XREF,FDS                                             PMAXW
      SUBROUTINE PMAXW(L2,X,Y,IN,RA,D,IDL,P)
C     MAXWELL-EN "P" POTENTZIAL KOEFIZIENTEAK KALKULATZEKO SUBRUTINA
      PARAMETER (N1=32)
      DIMENSION X(N1),Y(N1),IN(N1),RA(N1),D(N1),IDL(N1),P(N1,N1)
        DO 1 I=1,L2
        X1=2.*Y(I)
        X3=GMR(IN(I),RA(I),D(I),IDL(I))
        P(I,I)=18.E9*ALOG(X1/X3)
        J1=I+1
        IF(L2-J1)1,3,3
    3     DO 2 J=J1,L2
          X1=X(I)-X(J)
          X1=X1*X1
          X2=Y(I)+Y(J)
          X3=Y(I)-Y(J)
          X2=SQRT(X1+X2*X2)
          X3=SQRT(X1+X3*X3)
          P(I,J)=18.E9*ALOG(X2/X3)
          P(J,I)=P(I,J)
    2     CONTINUE
    1   CONTINUE
      RETURN
      END
      SUBROUTINE IZENBU
      CHARACTER EGUNA*10,IZENA*56,EUSK*5,EROA*15
      DOUBLE PRECISION EXPE
      CHARACTER     IZBR21(2)*7, IZBR22(2)*28,IZBR23(2)*15,IZBR24(2)*33
     4,IZBR31(2)*22,IZBR32(2)*18,IZBR33(2)*54,IZBR34(2)*16,IZBR35(2)*8
     5,IZBR36(2)*54,IZBR37(2)*16
      COMMON /HZENBU/IZBR21,IZBR22,IZBR23,IZBR24
     4       ,IZBR31,IZBR32,IZBR33,IZBR34,IZBR35,IZBR36,IZBR37
      COMMON /IZENBN/ IZK,EXPE,NOZ,N,KT,TL,MOZ,NFL
      COMMON /IZENBC/ EGUNA,IZENA,EUSK,EROA
      COMMON /IRKIDT/ IRT,IDT,IRF,IDF,IDU
      NOZ=NOZ+1
      IF(IZK-1)30,20,30
   20 WRITE(IDF,620) NOZ,EGUNA,IZENA,TL,EXPE
  620 FORMAT(1H1,5X,'ENPRESELEK,E.A.',87X,'ORRIALDE',I14/
     1           6X,'LINETAKO SAILA',88X,'DATA ',7X,A10/
     2 23X,A56,F5.0,'KV-EKO LINE ELEKTRIKOA',2X,'ESPEDIENTE',F12.2 )
      GO TO 40
   30 WRITE(IDF,630) NOZ,EGUNA,TL,IZENA,EXPE
  630 FORMAT(1H1,5X,'ENPRESELEK,E.A.',87X,'PAGE',I18/
     1           6X,'LINES DIVISION',88X,'DATE',8X,A10/
     223X,'TRANSMISSION LINE OF ',F5.0,' KV ',A54 ,1X,'EXPEDIENT',F13.2)
 
C     IZBR23=' EUSKARRIZ ETA ',IZBR24=' EROALEZ EGINDAKO LINE ELEKTRIKOA'
   40 IF(KT.EQ.1)WRITE(IDF,621)TRIM(IZBR21(IZK)),N
     1,TRIM(IZBR22(IZK))//' '//TRIM(EUSK)//TRIM(IZBR23(IZK))//
     2' '//TRIM(EROA)//TRIM(IZBR24(IZK))
  621 FORMAT(1H0, 5X,A,I3,T24,A)
 
C     IZBR31='  MUTUR ELIKATZAILEA  ', IZBR32='  MUTUR ELIKATUA  '
      IF(KT.EQ.0)WRITE(IDF,622)IZBR31(IZK),IZBR32(IZK)
     1,            (IZBR33(IZK),I=1,2),IZBR34(IZK)
     2,IZBR35(IZK),(IZBR36(IZK),I=1,2),IZBR37(IZK)
  622 FORMAT(1H0, 9X,15('*'),A22,14('*'),3X,17('*'),A18,16('*')/
     1  2(' ***'),2(A54),1X,A/ A,2(A54),1X,A/
     2 2(' ***'),2(7X,'(KV)', 9X,'(KA)', 9X,'(MW)', 8X,'(MVAR)',3X)
     3,3X,'(KV)',4X,'***'  ///)
      RETURN
      END
 
      SUBROUTINE HIRUAT(IR,CAUX,*,*,*)
      CHARACTER CAUX*4
CJF
      READ (IR,'(A)',END=3) CAUX
      IF (CAUX.EQ.'***E'.OR.CAUX.EQ.'***e') RETURN1
      IF (CAUX.EQ.'***P'.OR.CAUX.EQ.'***P') RETURN2
      BACKSPACE IR
      RETURN
    3 RETURN3
      END
C--+------------------------------------------------------
      SUBROUTINE ZAINDU(ERA,IRF,IZK0)
C  Behinolako programaren erabilera-baimenak kontrolatzeko azpierrutinea
C  gaur egun kasu batzutan erabiliko den hizkuntza zein den jakiteko
C  baino ez da erabiltzen bertsio laburrean
      INTEGER ERA
      READ(IRF,'(20X,I5)') IZK0
      RETURN
      END
C--+------------------------------------------------
      SUBROUTINE HIZKUN
C     Hizkuntzen araberako kateak definitzeko;
C     bi hizkuntzatarako dimentsionatuta
      CHARACTER CHL*132
      CHARACTER*25 KONFOR(2),KONFBA(2),POTGR1(2)*32,POTGR2(2)*60
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2
      CHARACTER     POTEL1(2)*48,POTEL2(2)*50
     1,ERMAG1(2)*28,ERMAG2(2)*15,ERMAG3(2)*122,ERMAG4(2)*19,ERMAG5(2)*19
     2,ERELE1(2)*25,ERELE2(2)*15,ERELE3(2)*121,ERELE4(2)*19,ERELE5(2)*19
     3,MEZU21(2)*15,MEZU22(2)*7, MEZU23(2)*36, MEZU24(2)*74
      COMMON /HZREMU/ POTEL1,POTEL2
     1       ,ERMAG1,ERMAG2,ERMAG3,ERMAG4,ERMAG5
     2       ,ERELE1,ERELE2,ERELE3,ERELE4,ERELE5
     3       ,MEZU21,MEZU22,MEZU23,MEZU24
      CHARACTER     IZBR21(2)*7, IZBR22(2)*28,IZBR23(2)*15,IZBR24(2)*33
     1,IZBR31(2)*22,IZBR32(2)*18,IZBR33(2)*54,IZBR34(2)*16,IZBR35(2)*8
     2,IZBR36(2)*54,IZBR37(2)*16
      COMMON /HZENBU/IZBR21,IZBR22,IZBR23,IZBR24
     1       ,IZBR31,IZBR32,IZBR33,IZBR34,IZBR35,IZBR36,IZBR37
      CHARACTER IZBE21(2)*90,IZBE22(2)*46
      COMMON /HZEPRO/IZBE21,IZBE22
      CHARACTER ANOI11(2)*4, ANOI12(2)*72,ANOI21(2)*47,ANOI22(2)*54
     1,ANOI23(2)*72,ANOI24(2)*26,ANOI25(2)*25,ANOI26(2)*23
      COMMON /HZOISE/ ANOI11,ANOI12,ANOI21,ANOI22,ANOI23,ANOI24,ANOI25
     1,ANOI26
      KONFOR(1)='KONFIGURAZIO ORIGINALA'
      KONFOR(2)='ORIGINAL CONFIGURATION'
      KONFBA(1)='KONFIGURAZIO BALIOKIDEA'
      KONFBA(2)='EQUIVALENT CONFIGURATION'
      POTGR1(1)='POTENTZIAL GRADIENTEA EROALEETAN'
      POTGR1(2)='ELECTRIC FIELD ON THE CONDUCTORS'
      CHL='EROALE     KARGA Q/M          BB-GR KV/CM'
      POTGR2(1)=TRIM(CHL)//'        GRMAX KV/CM'
      CHL='CONDUC.    CHARGE Q/M         GR.MV KV/CM'
      POTGR2(2)=TRIM(CHL)//'        GRMAX KV/CM'
      ERMAG1(1)='INDUKZIO MAGNETIKOA LURRETIK'
      ERMAG1(2)='MAGNETIC INDUCTION'
      ERMAG2(1)=' M-RA'
      ERMAG2(2)=' M ABOVE GROUND'
      CHL='DIST.  INDUK.HORIZONTALA   INDUK.BERTIKALA'
      CHL=TRIM(CHL)//'     INDUK.MAXIMOA      INDUK.MINIMOA'
      CHL=TRIM(CHL)//'      UNEKO INDUK.   INDUK. 90 G GEROAGO'
      ERMAG3(1)=TRIM(CHL)
      CHL='DIST.  INDUC.HORIZONTAL    INDUC.VERTICAL'
      CHL=TRIM(CHL)//'      MAXIM.INDUC.       MINIM.INDUC.'
      CHL=TRIM(CHL)//'      PRESENT INDUC.    INDUC. 90 G LAT.'
      ERMAG3(2)=TRIM(CHL)
      ERMAG4(1)=' MOD.TESLA   ARG.G '
      ERMAG4(2)=' MOD.TESLA   ARG.G '
      ERMAG5(1)=' MAG.TESLA   DIR.G '
      ERMAG5(2)=' MAG.TESLA   DIR.G '
      ERELE1(1)='EREMU ELEKTRIKOA LURRETIK'
      ERELE1(2)='ELECTRIC FIELD'
      ERELE2(1)=' M-RA'
      ERELE2(2)=' M ABOVE GROUND'
      CHL='DIST.   EREMU HORIZONTALA   EREMU BERTIKALA'
      CHL=TRIM(CHL)//'     EREMU MAXIMOA      EREMU MINIMOA'
      CHL=TRIM(CHL)//'      UNEKO EREMUA    EREMUA 90 G GEROAGO'
      ERELE3(1)=TRIM(CHL)
      CHL='DIST.   HORIZONTAL FIELD    VERTICAL FIELD'
      CHL=TRIM(CHL)//'      MAXIMUM FIELD      MINIMUM FIELD'
      CHL=TRIM(CHL)//'      PRESENT FIELD    FIELD 90 G LATER'
      ERELE3(2)=TRIM(CHL)
      ERELE4(1)='  MOD.KV/M  ARG.G  '
      ERELE4(2)='  MOD.KV/M  ARG.G  '
      ERELE5(1)=' MAG.KV/M   DIR.G  '
      ERELE5(2)=' MAG.KV/M   DIR.G  '
 
      MEZU21(1)=' X='
      MEZU21(2)=' RESULTS FOR X='
      MEZU22(1)=' ETA Y='
      MEZU22(2)=' AND Y='
      MEZU23(1)='-RAKO EMAITZAK EZ DIRA ESANGURATSUAK'
      MEZU23(2)=' ARE NOT SIGNIFICATIV'
      CHL='BEGIRATU IA EROALEEN TENTSIOAK EDO INTENTSITATEAK SARTU'
      MEZU24(1)=TRIM(CHL)//' DIREN'
      CHL='CHECK IF THE VOLTAGES AND CURRENTS OF ALL THE CONDUCTORS'
      MEZU24(2)=TRIM(CHL)//' HAVE BEEN DEFINED'
 
      POTEL1(1)='POTENTZIAL ELEKTRIKOAK GARIAERA EZBERDINETAN (X='
      POTEL1(2)='ELECTRICAL POTENTIAL AT DIFFERENT HEIGHTS (X='
      POTEL2(1)='GOIERA    POT.BA.EF.     POT.FAS.A.     POT.UN.BA.'
      POTEL2(2)='HEIGHT    POT.VA.EF.     POT.PHA.A.     POT.UN.BA.'
      IZBR21(1)='ZATI'
      IZBR21(2)='SECTION'
      IZBR22(1)=' '
      IZBR22(2)='TRANSMISSION LINE BUILT WITH'
      IZBR23(1)=' EUSKARRIZ ETA'
      IZBR23(2)=' STRUCTURES AND'
      IZBR24(1)=' EROALEZ EGINDAKO LINE ELEKTRIKOA'
      IZBR24(2)=' CONDUCTORS'
      IZBR31(1)='  MUTUR ELIKATZAILEA  '
      IZBR31(2)='      SENDING END     '
      IZBR32(1)='  MUTUR ELIKATUA  '
      IZBR32(2)='  RECEIVING END   '
      IZBR33(1)='     TENTSIO      LINETAKO    POTENTZIA    POTENTZIA '
      IZBR33(2)='   EQUIV. LINE       LINE       ACTIVE      REACTIVE '
      IZBR34(1)='TENTSIO-  ***'
      IZBR34(2)='VOLTAGE   ***'
      IZBR35(1)=' ZIR FAS'
      IZBR35(2)=' CIR PHA'
      IZBR36(1)='  KONPOSATU B. INTENTSITATEA  ERAGINKOR. BERRERAGINK.'
      IZBR36(2)='     VOLTAGE      INTENSITY      POWER        POWER  '
      IZBR37(1)='JAUSKERA  E I'
      IZBR37(2)='  DROP    E I'
      CHL='******************  SARRERAKO DATUAK, AZTERTUTAKO LI'
      IZBE21(1)=TRIM(CHL)//'NEAREN EZAUGARRIAK  ******************'
      CHL='***************************  INPUT DATA, LINE CHARAC'
      IZBE21(2)=TRIM(CHL)//'TERISTICS  ***************************'
      IZBE22(1)='ER.KOP.   FNCY     RHO          PC AN CE CM PL'
      IZBE22(2)='CON.NUM.  FNCY     RHO          PC AN CE CM PL'
 
      ANOI11(1)=' '
      ANOI11(2)=' THE'
      CHL=' EROALEA (LUR KABLEA) EZ DA KONTUAN HARTZEN'
      ANOI12(1)=TRIM(CHL)//' EROALE HEZEKO KASUAN '
      CHL=' CONDUCTOR  (GROUND WIRE) NOT CONSIDERED'
      ANOI12(2)=TRIM(CHL)//' FOR THE WET CONDUCTOR CONDITION'
      ANOI21(1)='PRESIO AKUSTIKOAREN ZEHARKAKO PROFILA, LURRETIK'
      ANOI21(2)='SOUND PRESSURE, LATERAL PROFILE'
      ANOI22(1)=' M-RA N/M2 ETA 20 MICRON/M2-REN GAINEKO DB-TAN'
      ANOI22(2)=' M ABOVE THE GROUND, IN N/M2 AND DB ABOVE 20 MICRON/M2'
      CHL=' DIST.       *******  BANDA ZABALEKO ZARATA ALEATORIOAREN'
      ANOI23(1)=TRIM(CHL)//' MAILA  *******'
      CHL=' DIST.       **************  BROADBAND RANDOM NOISE LEVEL'
      ANOI23(2)=TRIM(CHL)//'  *************'
      ANOI24(1)='100 HZ-EKO ZARATAREN MAILA'
      ANOI24(2)=' 100 HZ HUM  NOISE LEVEL  '
      ANOI25(1)='**** EURITE GOGORRAK ****'
      ANOI25(2)='****** HEAVY RAIN *****'
      ANOI26(1)='**** EROALE HEZEAK ****'
      ANOI26(2)='**** WET CONDUCTOR ****'
 
      RETURN
      END

Programa erabiltzeko datuak[aldatu | aldatu iturburu kodea]

Programa hobeto dokumentatzeko aldaketa txikitxo batzuk egin bazaizkio ere, ez dago guztiz ondo dokumentatua eta ez daude bera erabiltzeko argibiderik. Horren ordez, aurrean esandako zioengatik datuak sartzeko modua guztiz argi ez egoteagatik, programarekin saioren bat egin nahi duenari laguntzearren atal honetan programarekin saioak egiteko erabili izan diren datuak ematen dira. Saioak egin ahal izateko kopiatu datu horiek fitxategi batean.

Datuak, garai haietako fitxetan bezala, gehienez 80 karaktereko zabalera duten lerrotan sartzen dira,

Emaitzak ere, garai haietan bezala aurkezten ditu programak; garaiko konputagailuen inprimagailuetan erabilgarri zeuden ikurrez bakarrik testu alfanumerikoak osatuz. Ferroerresonantzia aztertzerakoan lortutako emaitzak baino ez dira ematen grafikoen bidez, inprimagailuen ikurrez sortutako grafikoen bidez. Emaitzak ”LIPARE.RES” eta “LIPARE.RPU” fitxategietan ematen dira gaur egun; ez zuzenean inprimagailura bidaliz lehen bezala. Lehenengo fitxategian lineak aztertzeko baliagarriak diren emaitzak ematen dira, eta bigarrenean sistema elektriko oso bateko zama-ikerketak egiteko lineen ordez erabili behar diren parametroak.

Datu batzuek bitxiak dirudite, baina datuak nola sartzen diren hobeto ikusteko aukeratu dira ez kasu “errealak” aztertzeko. Saiotarako erabil ahal diren datuak hurrengoak dira:

ESPLA-N FITXATEGIEN HASIERAN SARTUTAKO HIRU AZALPEN-LERROEN EREDUA.              2013-10-24 Programak ez ditu irakurtzen eskuineko iruzkinak.
LERRO HAUEK FITXATEGIEN EDUKIA IDENTIFIKATZEKO ERABILTEN DIRA,                   Datuak nola sartu behar diren ikusteko baino ez dute balio.
PROGRAMA LAGUNTZAILEAK ERABILIZ.  (PARLI PROGRAMAREKIN PROBAK EGITEKO DATUAK)
 496890862    0    0    1 371935712                                              25 zutabean erabiliko den hizkuntza (1=>euskara ala 2=>ingelesa)
***EREMU
765.SAIOTARAKO E11 LINEA                                              1990.11.24
  3       .0    .00 1 0 1 1 0   -7.621.27 22  0.007.62  2                        2013-10-24, eroale kop.,aterako diren dauentzat BEGIRA PROGRAMAN
  1.370  1.759 4  .4572 0 -15.24  15.24  15.24         4.000-120.0  765.0-120.1  Hurrengo linearen datuak (E12),
  1.370  1.759 4  .4572 0   0.00  15.24  15.24         4.000   0.0  765.0   0.0  errazago irakurteko moduan jarrita.
  1.370  1.759 4  .4572 0  15.24  15.24  15.24         4.000 120.0  765.0 120.0  cm  cm  kop  m  zerden1  m  m  m  ohm/km  zerden2  /  kA  gr  kV  gr
765.SAIOTARAKO E12 LINEA                                    98765432101990.11.24
  3                 1 0 1 0 0    -762 127 22   000 381  0
  13700  17590 4   4572 0 -15240  15240  15240   0.201  4000-12000 765000-12000  cm  cm  kop  m  zerden1  m  m  m  ohm/km  zerden2  /  kA  gr  kV  gr
  13700  17590 4   4572 0    000  15240  15240   0.200  4000     0 765000     0  zerden1 =0=>aurr. dat. eroalen arteko dist. =1=>aurr.dat. sortaren diametroa
  13700  17590 4   4572 0  15240  15240  15240   0.200  4000 12000 765000 12000  ohm/m eta zerden2 ez dira beharrezkoak EREMUrako
765.SAIOTARAKO E13 LINEA  (datu ez egokiekin)               10765012001990.11.24
  3                 0 0 0 0 0    -762 127 22   000 381  0                        2010-04-25
  13700  17590 4   4572 0 -1524   1524   1524           4   -120   765   -120    Lehen datu hauek eta goiko linearenak (E12)
  13700  17590 4   4572 0   000   1524   1524           4      0   765      0    berdin irakurten ziren
  13700  17590 4   4572 0  1524   1524   1524           4    120   765    120    baina orain ez. (orain hauek ez dira egokiak)
765.SAIOTARAKO E14 LINEA                                    98765432101990.11.24
  5                 0 0 1 1 0   -7.621.27 22  0.003.81  0
 1.3700 1.7590 4  .4572 0 -15.24  15.24  15.24   0.200  4.00-120.0  765.0-120.0
 1.3700 1.7590 4  .4572 0    .00  15.24  15.24   0.200  4.00    .0  765.0    .0
 1.3700 1.7590 4  .4572 0  15.24  15.24  15.24   0.200  4.00 120.0  765.0 120.0
 0.6000 0.8000 1        0 -10.00  24.00  24.00   0.500
 0.6000 0.8000 1        0  10.00  24.00  24.00   0.500
765.SAIOTARAKO E11 LINEA                                              1990.11.24
  3       .0    .00 1 0 1 1 0   -7.621.27 22  0.007.62  2                        2013-10-24
  1.370  1.759 4  .4572 0 -15.24  15.24  15.24         4.000-120.0  000.0-120.1  Aurreko E11 linearen datuak errepikatuta,
  1.370  1.759 4  .4572 0   0.00  15.24  15.24         4.000   0.0  000.0   0.0  ikusteko emaitzak aldatzen diren ala ez.
  1.370  1.759 4  .4572 0  15.24  15.24  15.24         4.000 120.0  000.0 120.0 
***PARLI
765.SAIOTARAKO P11 LINEA                                              2013-10-24
 1 1 1 0    50.                                                                   1. hizkuntza(euskara), 1 azterketa, 1 zirkuitu, ez batu zirk., frekuentzia 
1  32121    10.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX                  1=>linea datuak, 3 eroale,
  1.370  1.759 4  .4572 0 -15.24  15.24  15.24  0.201 1  .400-120.0 765.00-120.0  cm  cm  kop  m  zerden1  m  m  m  ohm/km  zerden2  /  kA  gr  kV  gr
  1.370  1.759 4  .4572 0   0.00  15.24  15.24  0.200 2  .400   0.0 765.00   0.0  PARLIn ez dira erabiltzen azken lau zutabeetako datuak
  1.370  1.759 4  .4572 0  15.24  15.24  15.24  0.200 3  4.00 120.0 765.00 120.0  EREMUrako diren erregistrotan bakarrik dira beharrezkoak
765.SAIOTARAKO P12 LINEA                                              2013-10-24
 2 1 1 0    50.                                                                   2. hizkuntza(ingelesa), 1 azterketa, 1 zirkuitu, ez batu zirk., frekuentzia 
1  30021    10.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX                  2013-10-24    
  13700  17590 4   4572 0 -15240  15240  15240  0.201 1  4000-12000 765000-12000  Aurreko linearen datuak (P11),
  13700  17590 4   4572 0   0000  15240  15240  0.200 2  4000   000 765000   000  dezimal gehiago satu ahal izateko moduan jarrita.
  13700  17590 4   4572 0  15240  15240  15240  0.200 3  4000 12000 765000 12000  Datu gutxiago ataratzea eskatzen da.
765.SAIOTARAKO P13 LINEA  (datu ez egokiekin)                         2010-04-25
 0 1 1 0    50.                                                                   0=>hizkuntza lehenetsia, 1 azterketa, 1 zirkuitu, ez batu z., frekuentzia 
1  30001    10.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX                  2010-04-25
  13700  17590 4   4572 0 -1524   1524   1524  0.200 1  4   -120   765   -120     Lehen datu hauek eta goiko linearenak (P12)
  13700  17590 4   4572 0   000   1524   1524  0.200 2  4      0   765      0     berdin irakurtzen ziren
  13700  17590 4   4572 0  1524   1524   1524  0.200 3  4    120   765    120     baina orain ez. (orain hauek ez dira egokiak)
***EREMU
1145EPRI Tansmission line ref. book 345 kV & above, 214 orr.          2013-10-25
  5       0.   0.00 0 1 0 0 1   -4.002.00100  0.001.50  1                         2013-10-25
        2.0000 8 1.0000 1 -18.50  22.00  22.00         4.000 -120. 1145.0 -120.   EPRIren Transmission line refernce book 345 kV & above
        2.0000 8 1.0000 1   0.00  22.00  22.00         4.000    0. 1145.0    0.   liburuaren 214. orrialdean kalkulatutako
        2.0000 8 1.0000 1  18.50  22.00  22.00         4.000  120. 1145.0  120.   zarata akustikoaren mailarekin konparatzeko
        1.2500 1          -15.25  35.00  35.00
 0.0000 1.2500 1           15.25  35.00  35.00
***PARLI
750.P21a LINEA, TRANS.LINE REF. BOOK 345 kV & ABOVE 91 orr.           2013-10-23
 2 5 2 0    60.                                                                  2. hizkuntza, 5 azterketa, 2 zirkuitu, ez bau zir., Frekuentzia
1  60025  193.1 100. 0000   -7.621.27 22  0.003.81   72HA CX-DIPPER              -1- 1 => Linearen datuak
  1.3990 1.759 4   4572 0 -38.10  21.00  00000.05281 1  4000-120.0 750.00-120.0  babes eroaleak kontuan hartu barik
  1.3990 1.759 4   4572 0 -22.86  21.00      0.05281 2  4000   0.0 750.00   0.0
  1.3990 1.759 4   4572 0  -7.62  21.00      0.05281 3  4000 120.0 750.00 120.0
  1.3990 1.759 4   4572 0   7.62  21.00      0.05281 4  4000 120.0 750.00 120.0
  1.3990 1.759 4   4572 0  22.86  21.00      0.05281 5  4000-120.0 750.00-120.0
  1.3990 1.759 4   4572 0  38.10  21.00      0.05281 6  4000   0.0 750.00   0.0
0                                                                               -2- 0 => aztertu egoera bat
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  1. zirkuituaren egoera, mutur elikatuaren tentsioak eta potentziak emanez
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  2. zirkuituaren egoera, mutur elikatuaren tentsioak eta potentziak emanez
42352010    2.3 0.01                                                            -3- 4 => ferroerresonantziaren azterketa, aztertu 235 balio, erreak. albo batean
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  lanean dagoen zirkuitoaren egoera
 2          221          222          223
42353010    .01 0.02                                                            -4- 4 => ferroerresonantziaren azterketa, aztertu 235 balio, erreak. beste alboan
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  lanean dagoen zirkuitoaren egoera
 2          221          222          223
42354010    4.6 0.02                                                            -5- 4 => ferroerresonantziaren azterketa, aztertu 235 balio, erreak, albo bietan
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  lanean dagoen zirkuitoaren egoera
 2          221          222          223
750.P21b LINEA, TRANS.LINE REF. BOOK 345 kV & ABOVE 91 orr.           2013-10-23
 2 5 2 0    60.                                                                  2. hizkuntza, 5 azterketa, 2 zirkuitu, 0=>ez batu zirk., Frekuentzia
1 100025  193.1 100. 0000   -7.621.27 22  0.003.81   72HA CX-DIPPER              -1- 1 => Linearen datuak
  1.3990 1.759 4   4572 0 -38.10  21.00  00000.05281 1  4000-120.0 750.00-120.0  babes eroeleak kontuan hartuta 
  1.3990 1.759 4   4572 0 -22.86  21.00      0.05281 2  4000   0.0 750.00   0.0
  1.3990 1.759 4   4572 0  -7.62  21.00      0.05281 3  4000 120.0 750.00 120.0
  1.3990 1.759 4   4572 0   7.62  21.00      0.05281 4  4000 120.0 750.00 120.0  cm  cm  kop  m  zerden1  m  m  m  ohm/km  zerden2  /  kA  gr  kV  gr
  1.3990 1.759 4   4572 0  22.86  21.00      0.05281 5  4000-120.0 750.00-120.0  zerden1 =0=>aurr. dat. eroalen arteko dist. =1=>aurr.dat. sortaren diametroa
  1.3990 1.759 4   4572 0  38.10  21.00      0.05281 6  4000   0.0 750.00   0.0  zerden2 =0=>babes-eroalea >0=>fase-eroalea
  0.5    0.555 1          -30.48  41.40      1.42916 0
  0.5    0.555 1          -15.24  41.40      1.42916 0
  0.5    0.555 1           15.24  41.40      1.42916 0
  0.5    0.555 1           30.48  41.40      1.42916 0
0                                                                               -2- 0 => aztertu egoera bat
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  mutur elikatuaren datuak
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  potentziak emanez bi zirkuituetan
42352010    2.3 0.01                                                            -3- 4 => ferroerresonantziaren azterketa, 235 bal., erreaktantziak albo batean
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  lanean dagoen zirkuitoaren egoera
 2          221          222          223
42353010    .01 0.02                                                            -4- 4 => ferroerresonantziaren azterketa, 235 bal., erreaktantziak beste alboan
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3  
 2          221          222          223
42354010    4.6 0.02                                                            -5- 4 => ferroerresonantziaren azterketa, 235 bal., erreaktantziak albo bietan
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2          221          222          223
380.P22 LINEA CN ALCANTARA - ET MORATA (COMP. EST. CESI-HE)           2013.10.23
 014 2 0    50.                                                                  0. hizkuntza, 14 azterketa, 2 zirkuitu, 0=>ez batu zirk., Frekuentzia
1  82223  309.5 100. 0000   -7.621.27 22  0.003.81   32TA DX-CARDINAL            -1- 1 => Linearen datuak
 1.1742 1.5210 2  .4000    -5.45  23.81  23.81 .0597 1  .400-120.0 380.00-120.0  kasu honetan lineak duen tarte ezberdin bakarrarenak
 1.1742 1.5210 2  .4000    -5.85  14.51  14.51 .0597 2  .400   0.0 380.00   0.0
 1.1742 1.5210 2  .4000   -12.65  14.51  14.51 .0597 3  .400 120.0 380.00 120.0
 1.1742 1.5210 2  .4000    12.65  14.51  14.51 .0597 4  .400 120.0 380.00 120.0  0. hizkunta => 4. lerroaren 25. zutabean sartutakoa
 1.1742 1.5210 2  .4000     5.45  23.81  23.81 .0597 5  .400-120.0 380.00-120.0
 1.1742 1.5210 2   4000     5.85  14.51  14.51 .0597 6  .400   0.0 380.00   0.0
 0.6052 0.7880 1  .4000    -7.10  31.97  31.97 .4040 0  .000   0.0   0.00   0.0
 0.6052 0.7880 1  .4000     7.10  31970  31970 .4040 0  .000   0.0   0.00   0.0
0                                                                                -2- 0 => aztertu egoera bat
 1 381.-120.    382. 120.    383.   0.    101.  .2 3   102.  .2 3   103.  .2 3   Mutur elikatuaren datuak
 2 384.-120.    385. 120.    386.   0.    104.  .2 3   105.  .2 3   106.  .2 3   tentsioak(kV,gr) eta potentziak (MW,tan(fi)) emanez
0                                                                                -3-
 1 381. 120.    382.   0.    383.-120.    .401 120.    .402   0.    .403-120.    Mutur elikatuaren datuak
 2 384. 120.    385.   0.    386.-120.    .404 120.    .405   0.    .406-120.    tentsioak (kV,gr) eta intentsitateak (kA,gr) emanez
0                                                                                -4-
 1 381. 120.011 382.   0.012 383.-120.013 .401 120.011 .402   0.012 .403-120.013 Mutur elikatuaren datuak
 2 384. 120.021 385.   0.022 386.-120.023 .404 120.021 .405   0.022 .406-120.023 modu ezohikoan
0                                                                                -5-
 1 381. 120.111 382.   0.112 383.-120.113 .401 120.211 .402   0.212 .403-120.213 Mutur elikatzailearen datuak
 2 384. 120.121 385.   0.122 386.-120.123 .404 120.221 .405   0.222 .406-120.223
0                                                                                -6-
 1 381. 120.111 382.   0.112 383.-120.113 .401 115.211 .402  -5.212 .403-125.213 Mutur elikatzailearen datuak
 2 384. 120.121 385.   0.122 386.-120.123 .404 115.221 .405  -5.222 .406-125.223
0                                                                                -7-
 1 381. 120.    382.   0.    383.-120.    387. 120.111 388.   0.112 389.-120.113 Tentsioak bi muturretan
 2 384. 120.    385.   0.    386.-120.    390. 120.121 391.   0.122 392.-120.123 (ez ohikoa)
0                                                                                -8-
 1 381. 120.011 382.   0.012 383.-120.013 387. 120.123 388.   0.112 389.-120.113 Tentsioak bi muturretan
 2 384. 120.021 385.   0.022 386.-120.023 390. 120.121 391.   0.122 392.-120.111 (ez ohikoa) (beste modu batera)
0                                                                                -9-
 1 .407 120.211 .408   0.212 .409-120.213 .401 120.    .402   0.    .403-120.    Intentsitateak bi muturretan
 2 .410 120.221 .411   0.222 .412-120.223 .404 120.    .405   0.    .406-120.    (ez ohikoa)
0                                                                                -10-
 1 .401 120.211 .402   0.212 .403-120.213 .407 120.211 .408   0.212 .409-120.213 Intentsitateak bi muturretan
 2 .404 120.221 .405   0.222 .406-120.223 .410 120.221 .411   0.222 .412-120.223 (ez ohikoa) (beste modu batera)
0                                                                                -11-
 1 381. 120.111 382.   0.112 383.-120.113 .401 120.011 .402   0.012 .403-120.013 Tentsioak eta intentsitateak
 2 384. 120.021 385.   0.022 386.-120.023 .404 120.221 .405   0.222 .406-120.223 bi muturretan (ez ohikoa)
0                                                                                -12-
 1 381. 120.123 382.   0.112 383.-120.113 .401 120.011 .402   0.012 .403-120.013 SAIOAK SAIOAK SAIOAK
 2 384. 120.021 385.   0.022 386.-120.011 .404 120.221 .405   0.222 .406-120.211 
42362010     .1 .025                                                             -13- 4 => Ferrorresonantziarako ikerketa. 236 balio ikertzen dira
 1 380.-120.    380. 120.    380.   0.    150.  .2 300 150.  .2 300 150.  .2 300 4xxxyz1x => Emaitzak grafiko lineal batean aurkezten dira 
 2          221          222          223                                        4xxx2xyz => Erreaktantziak albo batean bakarrik                                           
42362100     .1 .025                                                             -14- (I1,I3,4I1,F7.3,F5.2,33X,A5,A15)/4,balio kopurua<237,
 1 380.-120.    380. 120.    380.   0.    150.  .2 300 150.  .2 300 150.  .2 300 4xxxyz0x => Emaitzak grafiko logaritmiko batean aurkezten dira
 2          221          222          223                                        
***PARLI
765.SAIOTARAKO P23 LINEA, zati bakar batean kalkulatuta     20765011002013-10-24  2014-06-27
 0 2 1 0    50.
1  50001   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P23 LINEA, sei zatitan kalkulatuta           20765011002014-06-27
 0 2 1 0    50.
1  52026   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P23 LINEA, 3 x 2 zatitan transposizio barik  20765011002014-06-27
 0 4 1 0    50.
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P23 LINEA, 3 x 2 zatitan transposizioarekin  20765011002014-06-27
 0 4 1 0    50.
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572     0000  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    15240  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -15240  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572     0000  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
1  50002   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572     0000  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    15240  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572   -15240  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  24000  24000 0.500
  06000  08000 1           10000  24000  24000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, hiru zatitan kalkulatuta          20765011002014-06-27
 1 2 2 0    50.
1  80023   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, hiru zatitan bananduta            20765011002014-06-27  P24 linea 3 zatitan bananduta
 0 4 2 0    50.                                                                   kalkuluak zehaztasun handiagoz egiteko
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80021   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, 1EZ  transposizioekin             20765011002014-06-27  P24 linea 1EZ transposizioekin
 0 4 2 0    50.                                                                   Balioak nola ekilbratzen diren ikusteko
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -12750  20240  20240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  25240  25240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  25240  25240 0.200 2  4000   000 765000   000
  13700  17590 4   4572   -12750  20240  20240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, 1EZ-2EZ  transposizioekin         20765011002014-06-27  P24 linea 1EZ-2EZ transposizioekin
 0 4 2 0    50.                                                                   Balioak nola ekilbratzen diren ikusteko
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -12750  20240  20240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  25240  25240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  15240  15240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572     7750  25240  25240 0.200 5  4000   000 765000   000
  13700  17590 4   4572    12750  20240  20240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  25240  25240 0.200 2  4000   000 765000   000
  13700  17590 4   4572   -12750  20240  20240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572    12750  20240  20240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572     7750  15240  15240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  25240  25240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, 1EZ-2ES  transposizioekin         20765011002014-06-27  P24 linea 1EZ-2ES transposizioekin
 0 4 2 0    50.                                                                   Balioak nola ekilbratzen diren ikusteko
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572   -12750  20240  20240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  15240  15240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  25240  25240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572    12750  20240  20240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572     7750  15240  15240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  25240  25240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
1  80001   100.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  15240  15240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    -7750  25240  25240 0.200 2  4000   000 765000   000
  13700  17590 4   4572   -12750  20240  20240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  15240  15240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572     7750  25240  25240 0.200 5  4000   000 765000   000
  13700  17590 4   4572    12750  20240  20240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P25 LINEA, ABC CBA Erreaktantze Txikia       20765011002014-06-27
 0 2 2 0    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, ,Frekuentzia
1  80023   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  15240  15240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  25240  25240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 0-123456     20765011002014-06-27
 0 2 2 0    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 0-123123     20765011002014-06-27
 0 2 2 0    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 0-456123     20765011002014-06-27
 0 2 2 0    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572   -12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572    -7750  15240  15240 0.200 6  4000 12000 765000 12000
  13700  17590 4   4572     7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 0-000000     20765011002014-06-27
 0 2 2 0    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 0  4000-12000 765000-12000  DATU EZEGOKIAK 
  13700  17590 4   4572   -12750  20240  20240 0.200 0  4000   000 765000   000  faseetan 54. zutabeko balioa 0 baino handiagoa
  13700  17590 4   4572    -7750  15240  15240 0.200 0  4000 12000 765000 12000  izan behar da
  13700  17590 4   4572     7750  25240  25240 0.200 0  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 0  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 0  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 1-123456     20765011002014-06-27
 0 2 2 1    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, JARRI PARALELOAN ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000  Nahiz eta zirkuituak paraleloan jartzeko eskatu
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000  ez dira paraleloan jartzen 54, zutabeko balioak
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000  horretarako egokiak ez direlako
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 1-123123     20765011002014-06-27
 0 2 2 1    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, jarri paraleloan ,Frekuentzia
1  80023   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 1  4000-12000 765000-12000  Nahiz eta zirkuituak paraleloan jartzeko eskatu
  13700  17590 4   4572   -12750  20240  20240 0.200 2  4000   000 765000   000  ez dira paraleloan jartzen 54, zutabeko balioak
  13700  17590 4   4572    -7750  15240  15240 0.200 3  4000 12000 765000 12000  horretarako egokiak ez direlako
  13700  17590 4   4572     7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 2  4000   000 765000   000
  13700  17590 4   4572     7750  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 1-456123     20765011002014-06-27
 0 2 2 1    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, jarri paraleloan ,Frekuentzia
1  80023   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 4  4000-12000 765000-12000  zirkuituak paraleloan jartzen dira
  13700  17590 4   4572   -12750  20240  20240 0.200 5  4000   000 765000   000  eta zirkuitu konposatuari dagozkion
  13700  17590 4   4572    -7750  15240  15240 0.200 6  4000 12000 765000 12000  aldaezin generalizatuak atera ahal dira
  13700  17590 4   4572     7750  25240  25240 0.200 1  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 2  4000   000 765000   000  BAINA
  13700  17590 4   4572     7750  15240  15240 0.200 3  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3   Zirkuitu konposatuaren aldaezinak lortu eta gero
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3   ezin dira aztertu linearen egoerak
765.SAIOTARAKO P24 LINEA, ABC ABC Super Sorta, 1-456456     20765011002014-06-27
 0 2 2 1    50.                                                                  0. hizkuntza, 2 azterketa, 2 zirkuitu, jarri paraleloan ,Frekuentzia
1  80003   300.  75. 0000   -7.621.27 22  0.003.81   71HA CX-XXX
  13700  17590 4   4572    -7750  25240  25240 0.200 4  4000-12000 765000-12000  zirkuituak paraleloan jartzen dira
  13700  17590 4   4572   -12750  20240  20240 0.200 5  4000   000 765000   000  eta zirkuitu konposatuari dagozkion
  13700  17590 4   4572    -7750  15240  15240 0.200 6  4000 12000 765000 12000  aldaezin generalizatuak atera ahal dira
  13700  17590 4   4572     7750  25240  25240 0.200 4  4000-12000 765000-12000
  13700  17590 4   4572    12750  20240  20240 0.200 5  4000   000 765000   000  BAINA
  13700  17590 4   4572     7750  15240  15240 0.200 6  4000 12000 765000 12000
  06000  08000 1          -10000  35000  35000 0.500
  06000  08000 1           10000  35000  35000 0.500
0
 1 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3   Zirkuitu konposatuaren aldaezinak lortu eta gero
 2 750.-120.    750. 120.    750.   0.    200.  .2 3   200.  .2 3   200.  .2 3   ezin dira aztertu linearen egoerak
          amaierako lerro hau beharrezkoa da.

Begira hurrengo artikuluak ere[aldatu | aldatu iturburu kodea]


Erreferentziak[aldatu | aldatu iturburu kodea]

  1. a b c d e f g h Donald G. Fink eta H. Wayne Beaty, Standard Handbook for Electrical Engineers, 11. argitaraldia, McGraw-Hill, New York, 1978, ISBN 0-07-020974-X, 14. atalburua Overhead Power Transmission
  2. NGK-Locke isolagailu polimerikoen fabrikatzailea
  3. Advanced Rubber Products - Suspension Insulators
  4.   , http://www.avaids.com/icao.pdf .
  5.   Transmission Line Reference Book 345 kV and Above, Electric Power Research Institute, 1975, 76-82. orrialdeak .
  6. a b c d e f   Zorrozua, Miguel Angel (2002), Energia elektrikoa garraiatzeko lineak, Euskal Herriko Unibertsitatea, Zarautz, ISBN 84-8373-534-2 .
  7. Zorrozua, 56. orr.
  8.   Glover, Duncan; Sarma; Overbye (2012), Power System Analysis and Design (5.a. argitaraldia), Cengage Learning, 828. orrialdeak, ISBN 978-1-111-42577-7 .
  9. a b c d e f   Etxegarai, Agurtzane; Larrabe, Zigor, Goi-tentsioko lineak eta kableak, Euskal Herriko Unibertsitatea, ISBN 978-84-9860-668-3, https://addi.ehu.es/bitstream/10810/9339/1/goi-tentsioko-lineak.pdf .
  10.   Gary, Claude; Moreau, Marcel (1976), L'effet de couronne en tension alternative, Eyrolles, 11-12. orrialdeak .
  11.   Zorrozua, Miguel Angel (2002), Energia elektrikoa garraiatzeko lineak, Euskal Herriko Unibertsitatea, 57-59. orrialdeak, ISBN 84-8373-534-2 .
  12. Beaty, H. Wayne; Fink, Donald G. , Standard Handbook for Electrical Engineers (15th Edition) McGraw-Hill, 2007 978-0-07-144146-9 pages 14-105 through 14-106
  13. J. Glover, M. Sarma, and T. Overbye, Power System Analysis and Design, Fifth Edition, Cengage Learning, Connecticut, 2012, ISBN 978-1-111-42577-7, Chapter 5 Transmission Lines: Steady-State Operation

Beste irakurgai batzuk[aldatu | aldatu iturburu kodea]

  • William D. Stevenson, Jr. Elements of Power System Analysis Third Edition, McGraw-Hill, New York (1975) ISBN 0-07-061285-4 (bigarren edizioaren itzulpena behinola ikasliburu bezala erabiltzen zen Bilboko Ingeniaritza Goi Eskola Teknikoan)