Revision history [back]

Now perform the spacing calculations.

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

Similarly compute for D14 ... D81, D23 ... D82, D34 ... D83, D45 ... D84, D56 ... D85, D67 ... D87, D78, D87.

Then compute to the vertical distances.

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

Similarly compute for d14 ... d81, d23 ... d82, d34 ... d83, d45 ... d84, d56 ... d85, d67 ... d87, d78, d87.

Calculate the Potential Coefficients.

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Similarly compute for P13 ... P81, P23 ... P82, P34 ... P83, P45 ... P84, P56 ... P85, P67 ... P87, P78, P87.

Calculate the phase impedance.

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);

Similarly compute for Z13 ... Z81, Z23 ... Z82, Z34 ... Z83, Z45 ... Z84, Z56 ... Z85, Z67 ... Z87, Z78, Z87.

Now perform the spacing calculations.calculations. Start off with the horizontal distances.

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

Similarly compute for D14 ... D81, D23 ... D82, D34 ... D83, D45 ... D84, D56 ... D85, D67 ... D87, D78, D87.

Then compute to the vertical distances.

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

Similarly compute for d14 ... d81, d23 ... d82, d34 ... d83, d45 ... d84, d56 ... d85, d67 ... d87, d78, d87.

Calculate the Potential Coefficients.

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Similarly compute for P13 ... P81, P23 ... P82, P34 ... P83, P45 ... P84, P56 ... P85, P67 ... P87, P78, P87.

Calculate the phase impedance.

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);

Similarly compute for Z13 ... Z81, Z23 ... Z82, Z34 ... Z83, Z45 ... Z84, Z56 ... Z85, Z67 ... Z87, Z78, Z87.

Now perform the spacing calculations. Start off with the horizontal distances.relocated to single answer

D12 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D21 = np.sqrt((x1-xm2)**2+(y1av-ym2)**2); # in meters
D13 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters
D31 = np.sqrt((x1-xm3)**2+(y1av-ym3)**2); # in meters

d12 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d21 = np.sqrt((x1-x2)**2+(y1av-y2av)**2); # in meters
d13 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters
d31 = np.sqrt((x1-x3)**2+(y1av-y3av)**2); # in meters

P11 = 17.975109*math.log(2*y1av/GMrceq); # in km/microF
P22 = 17.975109*math.log(2*y2av/GMrceq); # in km/microF
P33 = 17.975109*math.log(2*y3av/GMrceq); # in km/microF
P44 = 17.975109*math.log(2*y4av/GMrceq); # in km/microF
P55 = 17.975109*math.log(2*y5av/GMrceq); # in km/microF
P66 = 17.975109*math.log(2*y6av/GMrceq); # in km/microF
P77 = 17.975109*math.log(2*y7av/GMre); # in km/microF
P88 = 17.975109*math.log(2*y8av/GMre); # in km/microF
P12 = 17.975109*math.log(D12/d12); # in km/microF
P21 = 17.975109*math.log(D12/d12); # in km/microF

Z11 = (Rc/nb)+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/(4*nb))+math.log(Derc/GMRceq));
Z22 = Z11; Z33 = Z11; Z44 = Z11; Z55 = Z11; Z66 = Z11;
Z77 = Re+(np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*((1/4)+math.log(Derc/GMRe));
Z88 = Z77;
Z12 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);
Z21 = (np.pi**2)*(10**-4)*f+4j*np.pi*(10**-4)*f*math.log(Derc/d12);