Given the geometric architecture of a transmission pylon, estimate the electricity loss along the transmission line.
This program was originally written in Fortran in 1989 by David R. Brown and the American Public Power Association.
The following sources provided data on the electrical properties of different conductors:
EFFECTS OF EARTH RESISTANCE ARE INCLUDED.
{earth_resisitivity_in_ohm_meters: EARTH RESISTIVITY IN OHM-METERS}
IF SHIELDS ARE SEGMENTED (INSULATED) TO PREVENT 60 HZ SHIELD WIRE CURRENTS, SHIELD WIRE CURRENTS SHOULD BE SET TO ZERO.
{with_shield_currents: shield_currents ? TO INCLUDE SHIELD CURRENTS, TYPE A 1. TO SET SHIELD CURRENTS TO ZERO, TYPE A 2.}
# Click the Blue Plane to preview this as a CrossCompute tool
conductor_table_path = 'conductors.csv'
phase_table_path = 'phases.csv'
shield_table_path = 'shields.csv'
span_table_path = 'spans.csv'
earth_resisitivity_in_ohm_meters = 100
with_shield_currents = 1
target_folder = '/tmp'
from os.path import join
output_text_path = join(target_folder, 'tloss.log')
target_text_path = join(target_folder, 'tloss.txt')
from macros import transmogrify_conductors
conductor_text_path = transmogrify_conductors(join(
target_folder, 'conductors.dat'), conductor_table_path)
from macros import transmogrify_pylons
pylon_text_path = transmogrify_pylons(join(
target_folder, 'pylons.dat',
), phase_table_path, shield_table_path, span_table_path)
from subprocess import run
executable = 'tloss.out'
script = 'tloss.f'
output_text = run([
'gfortran',
'-o',
executable,
script,])
from os.path import expanduser
from subprocess import check_output, STDOUT
output_text = check_output([
join('.', executable),
conductor_text_path,
pylon_text_path,
expanduser(target_text_path),
str(earth_resisitivity_in_ohm_meters),
str(with_shield_currents),
], stderr=STDOUT)
with open(output_text_path, 'wb') as output_file:
output_file.write(output_text)
print('output_text_path = ' + output_text_path)
print('target_text_path = ' + target_text_path)