User:Jon Peli Oleaga/FORTRAN program for calculating representative parameters and operating conditions of AC overhead transmission lines

Software is currently used to study the operation of electric power transmission lines. As an example, this article presents a program developed in the late 1970s to analyze the electrical behaviour of power lines.

The program was used for about a decade in the second half of the 20th century by the European electric utility Iberduero, based in Bilbao (Basque Country). Written by the author of this article, who was also its only user, it ceased to be used when he left the transmission power line department during a company restructuring. Several years after his retirement, he was the only person aware of its existence, capabilities, and accuracy.

Although the program is now obsolete, it can still be used to calculate line parameters, the electric and magnetic fields produced by the line, and the voltage and current imbalances generated under different operating conditions. It can also be used to study overvoltages caused by ferroresonance when one circuit of a double-circuit line is de-energized while connected to a shunt reactor at one end. The resulting overvoltages are not calculated with high accuracy, but the program allows identification of the reactance values for which they may occur.

The program is written in the FORTRAN programming language, which was commonly used for technical calculations at the time. Although minor modifications have been made to allow it to run on more recent computers and compilers, the source code largely retains the structure of the original version

Initially, the software consisted of two separate programs: EREMU, mainly intended to calculate the electric and magnetic fields produced by AC overhead power lines, and PARLI, used to calculate their principal parameters and representative operating conditions.

These programs were later merged into a single program in a way that preserved compatibility with the previously used input data files. This origin is still evident in the source code.

The electric and magnetic field calculations are performed in accordance with the methods described in Transmission Line Reference Book: 345 kV and Above (TLRB)^[1], published by the Electric Power Research Institute in 1975, and L’Effet de Couronne en Tension Alternative (ECTA)^[2], published by Eyrolles in 1976 as part of the Collection de la Direction des Études et Recherches d’Électricité de France.

Calculations of other line parameters, including series impedances, shunt admittances, ABCD parameters, noise profiles, and imbalances under operating conditions, are also carried out following the procedures described in the TLRB.

The analysis of resonant overvoltages is performed by means of successive imbalance calculations with different shunt reactors connected at one end, the other end, or both ends of the line.

The complete program is provided so that interested users can copy it and run it using a free FORTRAN compiler.

One suitable compiler is available from the Silverfrost website, which provides a free version for private use. Using the FTN95 Personal Edition compiler, the program can be compiled and executed without modification.

Input data for running the program

As was common in files of that period, the input data are written in lines with a maximum length of 80 characters.

Although minor changes have been made to improve the documentation of the program, it is still not fully documented and no detailed instructions for its use are available. In addition, the method of entering input data may not be entirely clear.

To assist anyone wishing to run the program, this section provides the data used to generate the results presented in this Wikibook. To run the program initially, copy this data into a file.

Some of the input data may appear unusual, but they were chosen to illustrate how the data are entered rather than to represent realistic cases.

The first two cases are intended to illustrate how the results of the program for the electric and magnetic fields should be interpreted at any location near a power line.

The third case presents some of the line parameters used in the calculation of the electric and magnetic fields and shows that the values of the potential gradients (electric fields) at the conductor surface, as calculated by the program for single-circuit lines, are consistent with those given in Figure 3.2.8 of the TLRB.

The fourth and fifth cases show that the values of the potential gradients (electric fields) at the conductor surface, as calculated by the program for double-circuit lines, are consistent with those given in Figures 3.2.24 and 3.2.25 of the TLRB. In these cases, in addition to the conductor surface gradients used for comparison with the aforementioned figures, the electric and magnetic fields at other locations near the lines are also calculated. This makes it possible to observe the significant influence of the phase arrangement of the second circuit relative to that of the first circuit on the electric and magnetic fields in the surrounding space.

The sixth case compares the acoustic noise profile calculated by the program with the values presented in Chapter 6 of the TLRB.

 ESPLA-N FITXATEGIEN HASIERAN SARTUTAKO HIRU AZALPEN-LERROEN EREDUA.              2026-02-08 The program does not read the comments on the right.
LERRO HAUEK FITXATEGIEN EDUKIA IDENTIFIKATZEKO ERABILTEN DIRA,                   They are only useful for showing how to enter data.
PROGRAMA LAGUNTZAILEAK ERABILIZ.  (DATA FOR TESTING THE LIPARE PROGRAM)
 496890862    0    0    2 371935712                                              Default language in column 25 (1=>Basque or 2=>English)
***EREMU                                                                         ***EREMU => calculate fields
800.EPRI Tansmission line ref.book 345kV & above,pag.47,52  80001320002026.03.08
  5       0.   0.00 0 0 0 0 0  -10.002.00 22  0.002.00  0                        5 conductors, (0. 0. not used im EREMU), 1 0 1 1 0 output (SEE PROGRAM)
 1.2300 1.5000 4  .4572 0 -14.00  18.50  18.50   0.200  1.00-120.0  800.0-120.0  cm  cm  kop  m  zerden1  m  m  m  ohm/km  zerden2  /  kA  gr  kV  gr
 1.2300 1.5000 4  .4572 0    .00  18.50  18.50   0.200  1.00    .0  800.0   0.0  zerden1 0=>previous value dist.between cond.1=>prev.val.bundle diameter
  12300  15000 4   4572 0  14000  18500  18500   0.200  1.00 12000 800000 12000  IN THIS LINE other way to insert similar data to data en previous lines
 0.6000 0.8000 1        0 -11.05  30.70  30.70   0.500                           ohm/m and zerden2 are not necesary in EREMU
   6000   8000 1        0  11050  30700  30700   0.500                           IN THIS LINE other way to insert data that allows more accurate figures.
***PARLI
765.TWO CIRCUIT E62SB LINE, ABC ABC Super Bundle, 1-456456  62100765202014-06-27
 0 2 2 1    50.                                                                  0. def.lang., 2 tasks, 2 circuits, MERGE CIRCUITS ,Frecuency
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 kalkulatzen 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
                                                                      These two final lines are required if the final calculation is of type "***PARLI".

Electric and magnetic fields

Calculations

Electric field on the surface of conductors

The electric field on the surface of conductors is calculated as señaled on pages 7 to 15 of ECTA^[2] and 109 to 111 of TLRF^[1].

Taking into acount that the ground surface can be considered an equipotential surface at zero potential, if the ground surface is represented as a plane (flat surface) at zero potential, the potential at any point above that surface created by the electric charges in $n$ conductors is calculated considering these $n$ charges and other $n$ opposite charges situadeted simetricaly to the original ones with respect to the ground plane; as shown in the figure F1.

With that disposition the potential of an uncharged condutor $j$ dew to the charge of the conductor $i$ would be:

$V_{ji}={-{\frac {q_{i}}{2\pi \epsilon _{0}}}}\ln {D_{ji}}-{\frac {-q_{i}}{2\pi \epsilon _{0}}}\ln {D'_{ji}}={\frac {q_{i}}{2\pi \epsilon _{0}}}\ln {\frac {D'_{ji}}{D_{ji}}}$

At the surface of the $i$ conductor of radius $r_{i}$ and placed at an height $h_{i}$ above the gound, taking into account that $2h_{i}>>r_{i}$ , would be:

 $V_{ii}={\frac {q_{i}}{2\pi \epsilon _{0}}}\ln {\frac {2h_{i}}{r_{i}}}$

where:

$q_{i}$ = uniform linear charge density of $i$ conductor.

$h_{i}$ = $i$ conductor height.

$r_{i}$ = $i$ conductor radius.

$D'_{ji}$ = direct distance from conductor $j$ to image conductor $i$ .

$D_{ji}$ = direct distance from conductor $j$ to conductor $i$

$\epsilon$ = dielectric constant = ${\frac {1}{36\pi }}10^{-9}$ farads/meter for air.

So

Considering the $n$ conductors:

$V_{1}=P_{11}q_{1}+..+P_{1i}q_{i}+..+P_{1n}q_{n}$

..............................................................

$V_{i}=P_{i1}q_{1}+..+P_{ii}q_{i}+..+P_{in}q_{n}$

..............................................................

$V_{n}=P_{n1}q_{1}+..+P_{ni}q_{i}+..+P_{nn}q_{n}$

in matrix notation

${\begin{Bmatrix}V_{1}\\..\\V_{i}\\..\\V_{n}\\\end{Bmatrix}}={\begin{bmatrix}P_{11}&..&P_{1i}&..&P_{1n}\\..&..&..&..&..\\P_{i1}&..&P_{ii}&..&P_{in}\\..&..&..&..&..\\P_{n1}&..&P_{ni}&..&P_{nn}\\\end{bmatrix}}{\begin{Bmatrix}q_{1}\\..\\q_{i}\\..\\q_{n}\\\end{Bmatrix}}$

and more succinctly for a system of $m+n$ conductors

${\begin{Bmatrix}V_{m+n}\\\end{Bmatrix}}={\begin{bmatrix}P_{m+n,m+n}\\\end{bmatrix}}{\begin{Bmatrix}q_{m+n}\\\end{Bmatrix}}$

Where ${\begin{Bmatrix}V_{m+n}\\\end{Bmatrix}}$ is the column vector of electric potentials of the $m+n$ conductors, and ${\begin{Bmatrix}q_{m+n}\\\end{Bmatrix}}$ is the column vector of charge per unit length of these conductors. The $(m+n).(m+n)$ square matrix ${\begin{bmatrix}P_{m+n,m+n}\\\end{bmatrix}}$ is known as the matrix of Maxwell potential coefficients of the system.

Thus, if the Maxwell potential coefficients of an electric line are known, the conductor potentials can be expressed as a function of their linear charge densities. However, in practical applications, the conductor potentials are typically known, whereas the quantities required to compute the resulting electric fields are the linear charge densities.

To obtain these, both sides of the previous equation are premultiplied by the inverse of the Maxwell potential coefficient matrix.

This yields:

${\begin{bmatrix}P_{m+n,m+n}\\\end{bmatrix}}^{-1}.{\begin{Bmatrix}V_{m+n}\\\end{Bmatrix}}={\begin{bmatrix}C_{m+n,m+n}\\\end{bmatrix}}.{\begin{Bmatrix}V_{m+n}\\\end{Bmatrix}}={\begin{bmatrix}P_{m+n,m+n}\\\end{bmatrix}}^{-1}.{\begin{bmatrix}P_{m+n,m+n}\\\end{bmatrix}}.{\begin{Bmatrix}q_{m+n}\\\end{Bmatrix}}={\begin{bmatrix}I_{m+n}\\\end{bmatrix}}.{\begin{Bmatrix}q_{m+n}\\\end{Bmatrix}}={\begin{Bmatrix}q_{m+n}\\\end{Bmatrix}}$

where ${\begin{bmatrix}I_{m+n}\\\end{bmatrix}}$ is the $(m+n)\times (m+n)$ identity matrix, and ${\begin{bmatrix}C_{m+n,m+n}\\\end{bmatrix}}$ , which is the inverse of the Maxwell potential coefficient matrix, is called the capacitance matrix of the conductor system. In this matrix, $C_{ii}$ represents the self-capacitance of conductor $i$ , while $C_{ij}$ represents the mutual capacitance between conductors $i$ and $j$ .

Once the matrix ${\begin{bmatrix}C_{m+n,m+n}\\\end{bmatrix}}$ is known, the conductor charges can be calculated as a function of the conductor potentials. From these charges, the electric potential and the electric field at any point can also be determined.

Electric field above ground

The expressions in Figure F1 are valid for instantaneous values. In the case of sinusoidal scalar or vectorial periodic quantities of fixed direction (such as electric charges, electric potentials, horizontal and vertical components of the electric field), other expressions can be established based on maximum and rms values, simply taking into account that:

Auxiliary formulas to be inserted later into a figure.

( provisional while the figures are being prepared )

$v_{pi}={-{\frac {q_{i}}{2\pi \epsilon _{0}}}}\ln {D_{pi}}-{\frac {-q_{i}}{2\pi \epsilon _{0}}}\ln {D'_{pi}}={\frac {q_{i}}{2\pi \epsilon _{0}}}\ln {\frac {D'_{pi}}{D_{pi}}}$

$v_{p}=\sum _{i=1}^{n}{\frac {q_{i}}{2\pi \epsilon _{0}}}\ln {\frac {D'_{pi}}{D_{pi}}}$

$e_{xpi}={\frac {q_{i}}{2\pi \epsilon _{0}}}{\frac {1}{\left|D_{pi}\right|}}\cos {\alpha _{pi}}+{\frac {-q_{i}}{2\pi \epsilon _{0}}}{\frac {1}{\left|D'_{pi}\right|}}\cos {\alpha '_{pi}}$

$e_{ypi}={\frac {q_{i}}{2\pi \epsilon _{0}}}{\frac {1}{\left|D_{pi}\right|}}\sin {\alpha _{pi}}+{\frac {-q_{i}}{2\pi \epsilon _{0}}}{\frac {1}{\left|D'_{pi}\right|}}\sin {\alpha '_{pi}}$

$e_{xp}=\sum _{i=1}^{n}{\frac {q_{i}}{2\pi \epsilon _{0}}}\left({\frac {1}{\left|D_{pi}\right|}}{\frac {(x_{p}-x_{i})}{\left|D_{pi}\right|}}-{\frac {1}{\left|D'_{pi}\right|}}{\frac {(x_{p}-x_{i})}{\left|D'_{pi}\right|}}\right)$

$e_{yp}=\sum _{i=1}^{n}{\frac {q_{i}}{2\pi \epsilon _{0}}}\left({\frac {1}{\left|D_{pi}\right|}}{\frac {(h_{p}-h_{i})}{\left|D_{pi}\right|}}-{\frac {1}{\left|D'_{pi}\right|}}{\frac {(h_{p}+h_{i})}{\left|D'_{pi}\right|}}\right)$

Magnetic field above ground

Program results for that data

This section shows how the program presents the calculated electric and magnetic fields, how the data provided for these fields at locations other than the conductor surface should be interpreted, and the influence of phase arrangement in double-circuit lines on the resulting fields.

The adjacent figure illustrates the meaning of the data presented by the program on pages three and four for the electric and magnetic fields at a point located 12 m from the line centre and 4 m above ground level.

The most relevant values are those shown in the last four columns of these pages, and they are presented in static form in the figure. The program is not designed to calculate directly the time-varying values also shown in the figure. To obtain these values without modifying the program, it is necessary to repeat the calculations 36 times, each time changing the phase angles of the input voltages and currents by 10 degrees.

Theoretically, the three phases A, B, and C of an AC three-phase system can be assigned to the conductors of a single-circuit transmission line in six different ways: ABC, BCA, CAB, ACB, CBA, and BAC. In practice, however, the first three configurations are equivalent, since under steady-state conditions the voltages and currents in the conductors behave identically in each case. The same equivalence applies to the last three configurations.

Output for the lines L01 and L02 of the shown input data file:

1     ENPRESELEK,E.A.                                                                                       PAGE                 1
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L01 H32 CX-DIPPER for Wikibooks figure                 EXPEDIENT   7500132.01
 
 
        ***************************  INPUT DATA, LINE CHARACTERISTICS  ***************************
 
        CON.NUM.  FNCY     RHO          PC AN CE CM PL             XA    XG   NX     YA    YG   NY
          5         0.     0.00          0  0  1  1  0          -14.00  1.00   4    4.00  2.00   0
 
        EROA  BGE  ERR  ZEN SOR.D ZER      X       Y1      Y2      INT KA  FAS G    TEN KV  PHAS G
 
          1  1.37 1.76   4  45.72  0    -15.00   15.00   15.00      1.00 -140.00    750.00 -120.00
          2  1.37 1.76   4  45.72  0      0.00   15.00   15.00      1.00  -20.00    750.00    0.00
          3  1.37 1.76   4  45.72  0     15.00   15.00   15.00      1.00  100.00    750.00  120.00
          4  0.60 0.80   1   0.00  0    -10.00   24.00   24.00      0.00    0.00      0.00    0.00
          5  0.60 0.80   1   0.00  0     10.00   24.00   24.00      0.00    0.00      0.00    0.00
1     ENPRESELEK,E.A.                                                                                       PAGE                 2
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L01 H32 CX-DIPPER for Wikibooks figure                 EXPEDIENT   7500132.01
 
 
 
                                                 ORIGINAL CONFIGURATION
 
                                             ELECTRIC FIELD ON THE CONDUCTORS
 
 
                           CONDUC.    CHARGE Q/M         GR.MV KV/CM        GRMAX KV/CM
 
                             1      0.56294989E-05       14.40179E+00       16.75257E+00
 
                             2      0.59636313E-05       15.25659E+00       17.74690E+00
 
                             3      0.56294989E-05       14.40179E+00       16.75257E+00
 
                             4      0.47318039E-06       10.64656E+00       10.64656E+00
 
                             5      0.47318034E-06       10.64656E+00       10.64656E+00
1     ENPRESELEK,E.A.                                                                                       PAGE                 3
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L01 H32 CX-DIPPER for Wikibooks figure                 EXPEDIENT   7500132.01
 
                                                   ORIGINAL CONFIGURATION
 
                                            ELECTRIC FIELD  4.00 M ABOVE GROUND
 
 
  DIST.   HORIZONTAL FIELD    VERTICAL FIELD      MAXIMUM FIELD      MINIMUM FIELD      PRESENT FIELD    FIELD 90 G LATER
           MOD.KV/M  ARG.G    MOD.KV/M  ARG.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G
 
  -14.00  2.207E+00 -154.86  9.612E+00   75.18  1.375E+01  -81.35  2.366E+00    8.65  4.481E+00  129.10  1.321E+01  -84.24
  -13.00  2.774E+00 -149.91  9.147E+00   78.23  1.321E+01  -77.99  2.861E+00   12.01  4.299E+00  142.13  1.282E+01  -81.17
  -12.00  3.289E+00 -146.97  8.564E+00   82.26  1.252E+01  -74.77  3.408E+00   15.23  4.227E+00  157.31  1.227E+01  -78.07
  -11.00  3.722E+00 -145.11  7.918E+00   87.56  1.171E+01  -71.90  4.003E+00   18.10  4.344E+00  173.70  1.159E+01  -74.94
  -10.00  4.055E+00 -143.83  7.278E+00   94.41  1.083E+01  -69.81  4.632E+00   20.19  4.696E+00 -170.30  1.081E+01  -71.75
1     ENPRESELEK,E.A.                                                                                       PAGE                 4
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L01 H32 CX-DIPPER for Wikibooks figure                 EXPEDIENT   7500132.01
 
                                                   ORIGINAL CONFIGURATION
 
                                            MAGNETIC INDUCTION  4.00 M ABOVE GROUND
 
 
   DIST.  INDUC.HORIZONTAL    INDUC.VERTICAL      MAXIM.INDUC.       MINIM.INDUC.      PRESENT INDUC.    INDUC. 90 G LAT.
         MOD.TESLA   ARG.G  MOD.TESLA   ARG.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G
 
  -14.00  1.401E-05   56.72  9.384E-06   24.98  2.308E-05   32.08  5.992E-06  -57.92  1.622E-05   47.89  1.749E-05 -161.31
  -13.00  1.336E-05   59.52  1.104E-05   27.39  2.359E-05   38.60  6.650E-06  -51.40  1.685E-05   55.34  1.780E-05 -156.21
  -12.00  1.250E-05   63.32  1.253E-05   29.22  2.393E-05   45.10  7.341E-06  -44.90  1.739E-05   62.84  1.801E-05 -151.29
  -11.00  1.152E-05   68.42  1.382E-05   30.73  2.413E-05   51.51  8.067E-06  -38.49  1.784E-05   70.37  1.815E-05 -146.61
  -10.00  1.054E-05   75.15  1.486E-05   32.12  2.420E-05   57.80  8.830E-06  -32.20  1.820E-05   77.89  1.823E-05 -142.20
1     ENPRESELEK,E.A.                                                                                       PAGE                 5
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L02 H32 CX-DIPPER for Wfigure, two phases interchanged EXPEDIENT   7500132.02
 
 
        ***************************  INPUT DATA, LINE CHARACTERISTICS  ***************************
 
        CON.NUM.  FNCY     RHO          PC AN CE CM PL             XA    XG   NX     YA    YG   NY
          5         0.     0.00          0  0  1  1  0          -14.00  1.00   4    4.00  2.00   0
 
        EROA  BGE  ERR  ZEN SOR.D ZER      X       Y1      Y2      INT KA  FAS G    TEN KV  PHAS G
 
          1  1.37 1.76   4  45.72  0    -15.00   15.00   15.00      1.00  100.00    750.00  120.00
          2  1.37 1.76   4  45.72  0      0.00   15.00   15.00      1.00  -20.00    750.00    0.00
          3  1.37 1.76   4  45.72  0     15.00   15.00   15.00      1.00 -140.00    750.00 -120.00
          4  0.60 0.80   1   0.00  0    -10.00   24.00   24.00      0.00    0.00      0.00    0.00
          5  0.60 0.80   1   0.00  0     10.00   24.00   24.00      0.00    0.00      0.00    0.00
1     ENPRESELEK,E.A.                                                                                       PAGE                 6
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L02 H32 CX-DIPPER for Wfigure, two phases interchanged EXPEDIENT   7500132.02
 
 
 
                                                 ORIGINAL CONFIGURATION
 
                                             ELECTRIC FIELD ON THE CONDUCTORS
 
 
                           CONDUC.    CHARGE Q/M         GR.MV KV/CM        GRMAX KV/CM
 
                             1      0.56294989E-05       14.40179E+00       16.75257E+00
 
                             2      0.59636313E-05       15.25659E+00       17.74690E+00
 
                             3      0.56294989E-05       14.40179E+00       16.75257E+00
 
                             4      0.47318039E-06       10.64656E+00       10.64656E+00
 
                             5      0.47318034E-06       10.64656E+00       10.64656E+00
1     ENPRESELEK,E.A.                                                                                       PAGE                 7
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L02 H32 CX-DIPPER for Wfigure, two phases interchanged EXPEDIENT   7500132.02
 
                                                   ORIGINAL CONFIGURATION
 
                                            ELECTRIC FIELD  4.00 M ABOVE GROUND
 
 
  DIST.   HORIZONTAL FIELD    VERTICAL FIELD      MAXIMUM FIELD      MINIMUM FIELD      PRESENT FIELD    FIELD 90 G LATER
           MOD.KV/M  ARG.G    MOD.KV/M  ARG.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G   MAG.KV/M   DIR.G
 
  -14.00  2.207E+00  154.86  9.612E+00  -75.18  1.375E+01  -81.35  2.366E+00    8.65  4.481E+00  129.10  1.321E+01   95.76
  -13.00  2.774E+00  149.91  9.147E+00  -78.23  1.321E+01  -77.99  2.861E+00   12.01  4.299E+00  142.13  1.282E+01   98.83
  -12.00  3.289E+00  146.97  8.564E+00  -82.26  1.252E+01  -74.77  3.408E+00   15.23  4.227E+00  157.31  1.227E+01  101.93
  -11.00  3.722E+00  145.11  7.918E+00  -87.56  1.171E+01  -71.90  4.003E+00   18.10  4.344E+00  173.70  1.159E+01  105.06
  -10.00  4.055E+00  143.83  7.278E+00  -94.41  1.083E+01  -69.81  4.632E+00   20.19  4.696E+00 -170.30  1.081E+01  108.25
1     ENPRESELEK,E.A.                                                                                       PAGE                 8
      LINES DIVISION                                                                                        DATE        2026-02-19
                       TRANSMISSION LINE OF  750. KV L02 H32 CX-DIPPER for Wfigure, two phases interchanged EXPEDIENT   7500132.02
 
                                                   ORIGINAL CONFIGURATION
 
                                            MAGNETIC INDUCTION  4.00 M ABOVE GROUND
 
 
   DIST.  INDUC.HORIZONTAL    INDUC.VERTICAL      MAXIM.INDUC.       MINIM.INDUC.      PRESENT INDUC.    INDUC. 90 G LAT.
         MOD.TESLA   ARG.G  MOD.TESLA   ARG.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G  MAG.TESLA   DIR.G
 
  -14.00  1.401E-05  -96.72  9.384E-06  -64.98  2.308E-05   32.08  5.992E-06  -57.92  6.072E-06  112.43  2.306E-05   31.43
  -13.00  1.336E-05  -99.52  1.104E-05  -67.39  2.359E-05   38.60  6.650E-06  -51.40  6.766E-06  117.52  2.355E-05   37.71
  -12.00  1.250E-05 -103.32  1.253E-05  -69.22  2.393E-05   45.10  7.341E-06  -44.90  7.494E-06  122.92  2.389E-05   43.93
  -11.00  1.152E-05 -108.42  1.382E-05  -70.73  2.413E-05   51.51  8.067E-06  -38.49  8.253E-06  128.59  2.407E-05   50.04
  -10.00  1.054E-05 -115.15  1.486E-05  -72.12  2.420E-05   57.80  8.830E-06  -32.20  9.041E-06  134.47  2.412E-05   56.00

An analysis of the results shown on pages 1 to 8 of the program output for lines L01 and L02 indicates that the electric and magnetic fields generated by these two phase configurations are identical. The only difference is that the direction of rotation of the fields outside the conductors is reversed.

_ The voltage gradients (electric fields) at the conductor surfaces do not rotate and are identical in both cases.

_ The maximum and minimum values of the electric and magnetic fields at all points outside the conductors are also identical in both cases.

_ At all locations where the minimum field value is nonzero, the direction of field rotation is opposite in the two configurations.

To interpret these results, it is necessary to consider that the absolute difference in angular position of the fields over a 5 ms interval must be less than 180°. Accordingly, defining FAP1 and FAP2 as the angular positions of the fields at instants 1 and 2, respectively, the rotation over 5 ms is given by:

   • FAP2 − FAP1, if |FAP2 − FAP1| < 180° 
   • FAP2 − FAP1 + 360°, if |FAP2 − FAP1| ≥ 180° and FAP2 − FAP1 < 0 
   • FAP2 − FAP1 − 360°, if |FAP2 − FAP1| ≥ 180° and FAP2 − FAP1 > 0

Note: In the analyzed cases, the angular velocities of the rotating fields may appear to differ. However, this is due to the fact that the angular velocity is not constant: it is higher when the field magnitude is lower and lower when the field magnitude is higher. Additionally, the initial field values (corresponding to the instantaneous voltage and current conditions specified in the input data) are different in each case.

Conversely, in transmission lines with more than one circuit, the phase arrangement of the additional circuits relative to that of the first circuit has a significant impact on the electric field values, as illustrated in Figure 3. This figure shows the maximum electric field at conductor surfaces and ground level for two of the six possible phase arrangements of the second circuit in a double-circuit 800 kV AC transmission line.

Audible noise profiles

Power line parameters

Resonance evaluations

Operating conditions

On pages 51, 52, 58, and 61 of transmission lines L10, L11, and L12 are illustrated the effects of using shunt reactors connected to an end of a long transmission line or a weak system to prevent overvoltages under light-load conditions, and series capacitors to compensate for the inductive reactance of loaded transmission lines, as described in the Voltage control and reactive power management Wikipedia article.

The output of the program for those lines and two operating conditions (moderate inductive load and no load) is the following:

1     ENPRESELEK,E.A.                                                                                       PAGE                44
      LINES DIVISION                                                                                        DATE        2013-10-24
                       TRANSMISSION LINE OF  765. KV L09 SINGLE CIRCUIT H32 1x1 sections                    EXPEDIENT  75013211.00
0     SECTION  1       TRANSMISSION LINE BUILT WITH 71HA STRUCTURES AND CX-XXX CONDUCTORS
0                       INPUT DATA
0                       LINE LONG.   1*300.00= 300.00 KM
                        GROUND RESISTIVITY      75.00 OHM.M
                        FREQUENCY               50.00 HZ
0                       ***************************  CONDUCTORS� FEATURES  ***************************
                           ..   GMR(CM)   RAD(CM)     IN   D(M)  ID      DX(M)   DY(M)      R(OHM:KM)
 
                            1   1.37000   1.75900      4   0.46   0     -15.24   15.24      0.2000000
                            2   1.37000   1.75900      4   0.46   0       0.00   15.24      0.2000000
                            3   1.37000   1.75900      4   0.46   0      15.24   15.24      0.2000000
                            4   0.60000   0.80000      1   0.00   0     -10.00   24.00      0.5000000
                            5   0.60000   0.80000      1   0.00   0      10.00   24.00      0.5000000
1     ENPRESELEK,E.A.                                                                                       PAGE                45
      LINES DIVISION                                                                                        DATE        2013-10-24
                       TRANSMISSION LINE OF  765. KV L09 SINGLE CIRCUIT H32 1x1 sections                    EXPEDIENT  75013211.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      746.450        0.585      188.145     -167.872        750.000        0.471      200.000       40.000     42.393   0 0
  1   2      741.212        0.637      204.246     -180.220        750.000        0.471      200.000       40.000     39.758   0 0
  1   3      733.405        0.640      219.038     -159.847        750.000        0.471      200.000       40.000     45.890   0 0
0 1                                    611.429     -507.939                                  600.000      120.000
 
 
0 0                                    611.429     -507.939                                  600.000      120.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                46
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L10 SINGLE CIRCUIT H32 1x6 sections                    EXPEDIENT  75013216.00
0     SECTION  1       TRANSMISSION LINE BUILT WITH 71HA STRUCTURES AND CX-XXX CONDUCTORS
0                       INPUT DATA
0                       LINE LONG.   6* 50.00= 300.00 KM
                        GROUND RESISTIVITY      75.00 OHM.M
                        FREQUENCY               50.00 HZ
0                       ***************************  CONDUCTORS� FEATURES  ***************************
                           ..   GMR(CM)   RAD(CM)     IN   D(M)  ID      DX(M)   DY(M)      R(OHM:KM)
 
                            1   1.37000   1.75900      4   0.46   0     -15.24   15.24      0.2000000
                            2   1.37000   1.75900      4   0.46   0       0.00   15.24      0.2000000
                            3   1.37000   1.75900      4   0.46   0      15.24   15.24      0.2000000
                            4   0.60000   0.80000      1   0.00   0     -10.00   24.00      0.5000000
                            5   0.60000   0.80000      1   0.00   0      10.00   24.00      0.5000000
1     ENPRESELEK,E.A.                                                                                       PAGE                47
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L10 SINGLE CIRCUIT H32 1x6 sections                    EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      745.885        0.589      188.280     -169.966        750.000        0.471      200.000       40.000     41.764   0 0
  1   2      740.790        0.640      204.507     -182.246        750.000        0.471      200.000       40.000     39.212   0 0
  1   3      733.096        0.644      219.417     -161.946        750.000        0.471      200.000       40.000     45.204   0 0
0 1                                    612.204     -514.158                                  600.000      120.000
 
 
0 0                                    612.204     -514.158                                  600.000      120.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                48
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L10 SINGLE CIRCUIT H32 1x6 sections                    EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      750.000        0.549      -16.272     -236.953        792.853        0.000        0.000        0.000     25.210   1 0
  1   2      750.000        0.579        1.646     -250.835        788.678        0.000        0.000        0.000     22.729   2 0
  1   3      750.000        0.548       19.506     -236.479        791.493        0.000        0.000        0.000     24.332   3 0
0 1                                      4.880     -724.267                                    0.000        0.000
 
 
0 0                                      4.880     -724.267                                    0.000        0.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                49
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L11 1C H32 1x6 sections + receiving end shunt reactor  EXPEDIENT  75013216.00
0       SECTI  1                                    SHUNT ADMITANCES (MHO)
 
                   1                2                3
 
        1    0.0000000E+00    0.0000000E+00    0.0000000E+00
            -0.4000000E-03    0.0000000E+00    0.0000000E+00
 
        2    0.0000000E+00    0.0000000E+00    0.0000000E+00
             0.0000000E+00   -0.4000000E-03    0.0000000E+00
 
        3    0.0000000E+00    0.0000000E+00    0.0000000E+00
             0.0000000E+00    0.0000000E+00   -0.4000000E-03
 
1     ENPRESELEK,E.A.                                                                                       PAGE                50
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L11 1C H32 1x6 sections + receiving end shunt reactor  EXPEDIENT  75013216.00
0     SECTION  2       TRANSMISSION LINE BUILT WITH 71HA STRUCTURES AND CX-XXX CONDUCTORS
0                       INPUT DATA
0                       LINE LONG.   6* 50.00= 300.00 KM
                        GROUND RESISTIVITY      75.00 OHM.M
                        FREQUENCY               50.00 HZ
0                       ***************************  CONDUCTORS� FEATURES  ***************************
                           ..   GMR(CM)   RAD(CM)     IN   D(M)  ID      DX(M)   DY(M)      R(OHM:KM)
 
                            1   1.37000   1.75900      4   0.46   0     -15.24   15.24      0.2000000
                            2   1.37000   1.75900      4   0.46   0       0.00   15.24      0.2000000
                            3   1.37000   1.75900      4   0.46   0      15.24   15.24      0.2000000
                            4   0.60000   0.80000      1   0.00   0     -10.00   24.00      0.5000000
                            5   0.60000   0.80000      1   0.00   0      10.00   24.00      0.5000000
1     ENPRESELEK,E.A.                                                                                       PAGE                51
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L11 1C H32 1x6 sections + receiving end shunt reactor  EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      770.663        0.482      187.027     -105.234        750.000        0.471      200.000       40.000     38.759   0 0
  1   2      763.706        0.534      204.028     -117.360        750.000        0.471      200.000       40.000     36.518   0 0
  1   3      758.906        0.549      219.753      -97.347        750.000        0.471      200.000       40.000     42.684   0 0
0 1                                    610.808     -319.941                                  600.000      120.000
 
 
0 0                                    610.808     -319.941                                  600.000      120.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                52
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L11 1C H32 1x6 sections + receiving end shunt reactor  EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      750.000        0.363      -17.257     -156.154        765.313        0.000        0.000        0.000      8.846   1 0
  1   2      750.000        0.394        0.569     -170.438        763.694        0.000        0.000        0.000      8.023   2 0
  1   3      750.000        0.363       18.297     -155.914        763.484        0.000        0.000        0.000      8.337   3 0
0 1                                      1.608     -482.506                                    0.000        0.000
 
 
0 0                                      1.608     -482.506                                    0.000        0.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                53
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L12 1C H32 1x6 sections + receiving end series capacit EXPEDIENT  75013216.00
0       SECTI  1                                    SERIES IMPEDANCES (OHM)
 
                   1                2                3
 
        1    0.0000000E+00    0.0000000E+00    0.0000000E+00
            -0.4000000E+03    0.0000000E+00    0.0000000E+00
 
        2    0.0000000E+00    0.0000000E+00    0.0000000E+00
             0.0000000E+00   -0.4000000E+03    0.0000000E+00
 
        3    0.0000000E+00    0.0000000E+00    0.0000000E+00
             0.0000000E+00    0.0000000E+00   -0.4000000E+03
 
1     ENPRESELEK,E.A.                                                                                       PAGE                54
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L12 1C H32 1x6 sections + receiving end series capacit EXPEDIENT  75013216.00
0     SECTION  2       TRANSMISSION LINE BUILT WITH 71HA STRUCTURES AND CX-XXX CONDUCTORS
0                       INPUT DATA
0                       LINE LONG.   6* 50.00= 300.00 KM
                        GROUND RESISTIVITY      75.00 OHM.M
                        FREQUENCY               50.00 HZ
0                       ***************************  CONDUCTORS� FEATURES  ***************************
                           ..   GMR(CM)   RAD(CM)     IN   D(M)  ID      DX(M)   DY(M)      R(OHM:KM)
 
                            1   1.37000   1.75900      4   0.46   0     -15.24   15.24      0.2000000
                            2   1.37000   1.75900      4   0.46   0       0.00   15.24      0.2000000
                            3   1.37000   1.75900      4   0.46   0      15.24   15.24      0.2000000
                            4   0.60000   0.80000      1   0.00   0     -10.00   24.00      0.5000000
                            5   0.60000   0.80000      1   0.00   0      10.00   24.00      0.5000000
1     ENPRESELEK,E.A.                                                                                       PAGE                55
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L12 1C H32 1x6 sections + receiving end series capacit EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      723.186        0.744      191.235     -244.626        750.000        0.471      200.000       40.000    138.759   0 0
  1   2      720.314        0.794      206.213     -257.834        750.000        0.471      200.000       40.000    143.290   0 0
  1   3      709.199        0.788      219.800     -236.299        750.000        0.471      200.000       40.000    139.391   0 0
0 1                                    617.249     -738.760                                  600.000      120.000
 
 
0 0                                    617.249     -738.760                                  600.000      120.000
 
 
1     ENPRESELEK,E.A.                                                                                       PAGE                56
      LINES DIVISION                                                                                        DATE        2014-06-27
                       TRANSMISSION LINE OF  765. KV L12 1C H32 1x6 sections + receiving end series capacit EXPEDIENT  75013216.00
0         ***************      SENDING END     **************   *****************  RECEIVING END   ****************
 *** ***   EQUIV. LINE       LINE       ACTIVE      REACTIVE     EQUIV. LINE       LINE       ACTIVE      REACTIVE   VOLTAGE   ***
 CIR PHA     VOLTAGE      INTENSITY      POWER        POWER        VOLTAGE      INTENSITY      POWER        POWER      DROP    E I
 *** ***       (KV)         (KA)         (MW)        (MVAR)          (KV)         (KA)         (MW)        (MVAR)      (KV)    ***
 
 
 
  1   1      750.000        0.549      -16.272     -236.953        792.853        0.000        0.000        0.000     25.210   1 0
  1   2      750.000        0.579        1.646     -250.835        788.678        0.000        0.000        0.000     22.729   2 0
  1   3      750.000        0.548       19.506     -236.479        791.493        0.000        0.000        0.000     24.332   3 0
0 1                                      4.880     -724.267                                    0.000        0.000
 
 
0 0                                      4.880     -724.267                                    0.000        0.000

Source file

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  1970eko hamarkadaren amaieran.
C  1980ko hamarkadaren amaieran batu ziren eta, geroago, konpilatzaile
C  berriagorekin erabili ahal izateko beharrezko aldaketak egin zaizkie.
C  Sarrerako datuak oraindino 80 zutabetan sartzen dira, 
C  behinolako fitxetan legez.
C  Bigarren hizkuntza modu naif eta ez sistematikoan sartu da.
C  Letra xehez idatzitako lerroak programaren dokumentazioa
C  apur bat hobetzeko sartu dira geroago.
C  Programa Creative Commons-en
C  Attribution-ShareAlike 4.0 lizentziapean partekatzen da
C  http://creativecommons.org/licenses/by-sa/4.0/
C
C  The original PARLI program to calculate electromagnetic parameters
C  and operating conditions of AC overhead transmission lines
C  (it can also be used to conduct ferroresonance studies)
C  and the EREMU program to calculate electromagnetic fields
C  created by AC overhead transmission lines
C  were created by Jon Peli Oleaga Olabarria
C  at the end of the 1970s.
C  They were merged in the last 1980s, 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 lines in lowercase were added later 
C  to slightly improve the program documentation.
C  The work is licensed under the Creative Commons 
C  Attribution-ShareAlike 4.0 License.
C  http://creativecommons.org/licenses/by-sa/4.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  )'
      PRINT *,' '
      PRINT *,'Enter the name of the input data file:'
      PRINT *,'<RETURN> implies default file name, = '//FITX
      PRINT *,'( To exit the program enter *<RETURN> )'
      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,IZK0) !2026-03-31 gehitu IZK0
      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,IZK0) !2026-03-31 gehitu IZK0
      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
      IZK=IZK0  !2026-03-31
      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 EPRI's TLRB345(1975)pag.116
      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   ! Linearen zati honen azpizatien inp. eta adm.
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)! 2026-03-28 lehen 12F6.3 (datuak S eta ohm-etan) 
      IF(KT.EQ.3)THEN
      DO 206 I=1,N3    ! inpedanziak eta admitantziak sartu ordez
  206 Z(I,I)= 1/Z(I,I) ! beti erreaktoreen inpedantziak sartzeko 
      ENDIF            ! Ohiko erreaktoreakaz datu zehatzagoak sartzeko 
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   !zatiaren azpizati honen ABCD matrizea OOJP-20260312
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 !DO honetan barneko azpizatiak batzen dira OOJP
      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)
      !=>hasi sartzen zatiak hartzaile muturretik OOJP-20260315
      GO TO 245
  243 IF(K2.GE.1)CALL ATARAC(AD,N2,ND2)
  244 CALL BITRIC(AD,N2,N2,RF,N2,N2,BC,ND2)
      !=>hasi sartzen zatiak hartzaile muturretik OOJP-20260315
  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.1,3I1))        !OOJP-20260219
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
     1,MAXP1(2)*35,MAXP2(2)*19,KAP1(2)*27,KAP2(2)*19
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2,MAXP1,MAXP2,KAP1,KAP2
      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
C     COMMON /HIZKON/ KONFOR,KONFBA   OOJP-20260208
      CHARACTER*25 KONFOR(2),KONFBA(2),POTGR1(2)*32,POTGR2(2)*60
     1,MAXP1(2)*35,MAXP2(2)*19,KAP1(2)*27,KAP2(2)*19
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2,MAXP1,MAXP2,KAP1,KAP2
      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.3,0PF8.2)) !OOJP-20260216 lehen 1PE9.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
     1,MAXP1(2)*35,MAXP2(2)*19,KAP1(2)*27,KAP2(2)*19
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2,MAXP1,MAXP2,KAP1,KAP2
      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
     1,MAXP1(2)*35,MAXP2(2)*19,KAP1(2)*27,KAP2(2)*19
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2,MAXP1,MAXP2,KAP1,KAP2
      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 
C   4 WRITE(IDF,604)
C 604 FORMAT(1H /25X,36HMAXWELL-EN POTENTZIAL KOEFIZIENTEAK, /
C    &32X,19H( METRO FARADIOKO ))
    4 WRITE(IDF,604)MAXP1(IZK),MAXP2(IZK)
  604 FORMAT(1H /25X,A,/32X,A)
      GO TO 6
C    5 WRITE(IDF,605)
C 605 FORMAT(1H /36X,12HKAPAZITATEAK,/32X,19H( FARADIO METROKO ))
    5 WRITE(IDF,605)KAP1(IZK),KAP2(IZK)
  605 FORMAT(1H /29X,A,/32X,A)
    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
     1,MAXP1(2)*35,MAXP2(2)*19,KAP1(2)*27,KAP2(2)*19
      COMMON /HZDATZ/ KONFOR,KONFBA,POTGR1,POTGR2,MAXP1,MAXP2,KAP1,KAP2
      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'
      MAXP1(1)='MAXWELL-EN POTENTZIAL KOEFIZIENTEAK'
      MAXP1(2)=" MAXWELL'S POTENTIAL COEFFICIENTS " 
      MAXP2(1)='( METRO FARADEKO )'
      MAXP2(2)='( METRE PER FARAD )'
      KAP1(1)='KAPAZITANTZIA-KOEFIZIENTEAK'
      KAP1(2)='CAPACITANCE COEFFICIENTS' !EPRI'S TLRB345kV(1975) pag.94
      KAP2(1)=' ( FARAD METROKO )'
      KAP2(2)='( FARAD PER METRE )'
      
      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

↑ ^a ^b Chaston, A.N.; LaForest, J.J.; Ramirez, A.R.; and, others (1975). Transmission Line Reference Book 345 kV and Above (1st ed.). p. 393.
↑ ^a ^b Gary, Claude; Moreau, Marcel (1977). L’Effet de Couronne en Tension Alternative (1st ed.). Eyrolles. p. 464.

[TLRB345kV-1] Chaston, A.N.; LaForest, J.J.; Ramirez, A.R.; and, others (1975). Transmission Line Reference Book 345 kV and Above (1st ed.). p. 393.

[ECTA-EDF-2] Gary, Claude; Moreau, Marcel (1977). L’Effet de Couronne en Tension Alternative (1st ed.). Eyrolles. p. 464.

[1]

[2]