In (3), g is an

equivalent conductance. This parameter depends on the value of the element that was discretized and the specific time step [DELTA]t chosen for the simulation; the term h is a history current source, whose magnitude depends on past values of voltages at element terminals, the current through the element, or both, depending on the specific integration method used [9].

Accordingly, for each phase of the load, equivalent conductance ([G.sub.e]) is calculated as;

From the outlines of the proposed power resolution, it can be qualitatively concluded that two power components could be used to detect the harmonic producing loads: These are; scattered conductance power ([D.sub.sc]), which occurs by the difference between nth harmonic balanced conductance and equivalent conductance, and scattered susceptance power ([D.sub.ss]), which occurs by the difference between nth harmonic balanced susceptance and nth harmonic equivalent susceptance.

The equivalent conductance matrix [[G.sub.e]] is obtained by summing the conductance matrices of the two-port components ([[G.sub.e]]=[G]+[[G.sub.L]]) as it is developed in the equation (6), where [[G.sub.L]] is the conductance matrix of the reciprocal two-port longitudinal resistances [21-24]:

This equation verifies the equality of the equivalent conductance ([G.sub.e1k]=[G.sub.e10]).

Finally, 490PT or

equivalent conductance testers, NSN 613001-510-9594, are found in the SATS and the FRS.

where F is a soil-structure dependent formation factor, [delta] is the empirical ratio between

equivalent conductance of clay counter ions to the maximum value of this

equivalent conductance, and [I.sub.n] is the intercept of the linear part of the [[sigma.sub.a]-[[sigma].sub.w] curve.

Cn =

Equivalent conductance of combined paths except i'th.

Then, the

equivalent conductance [G.sub.eq](V, f), the susceptance [B.sub.eq](V, f) and the equivalent capacitance [C.sub.eq](V, f) of the modified simple equivalent circuit, shown in Fig.

ECa itself depends on EC of the soil solution (ECw), water content ([theta], and this dependence is referred to as `tortuosity' or the `formation factor', F), EC of the solid phase charged surfaces (ECs), and the interaction between the electrolyte concentration and ECs (represented by [delta] which is the empirical ratio of the

equivalent conductance of clay counter ions to the maximal value of this

equivalent conductance).