Radio Noise Analysis Program¶

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:

AL/ACSR data: Aluminum Electrical Conductor Handbook, 1982.
Alumoweld data: Copperweld Steel Company, Southern Division, Fayetteville, TN.
Steel data: Mechanical Data from ASTM Standards, ASTM Standard A475.
CU data: Westinghouse T&D Reference, 1964, ABB Power Systems Inc.

EARTH RETURN EFFECTS ARE INCLUDED.

{rho: EARTH RESISTIVITY IN OHM-METERS}

{epsrel: RELATIVE EARTH PERMITTIVITY }

THE PROGRAM CALCULATES NOISE AT GROUND LEVEL, Y=0, FOR A RANGE OF X COORDINATES.

For example, if xstart is 0, xstop is 10 and the number of data points is 3, then the program will calculate noise at ground level at the following distances away from the distribution line: 0 feet, 5 feet, 10 feet.

{xstart: xstart (feet)}

{xstop: xstop (feet)}

{n_datapoints: NUMBER OF DATA POINTS}

In [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'

rho = 100
epsrel = 2.8
xstart = 0
xstop = 125
n_datapoints = 51

target_folder = '/tmp'

In [3]:

from os.path import join
output_text_path = join(target_folder, 'rnoise.log')
target_text_path = join(target_folder, 'rnoise.txt')

In [4]:

from macros import transmogrify_conductors
conductor_text_path = transmogrify_conductors(join(
    target_folder, 'conductors.dat'), conductor_table_path)

In [5]:

from macros import transmogrify_pylons
pylon_text_path = transmogrify_pylons(join(
    target_folder, 'pylons.dat',
), phase_table_path, shield_table_path, span_table_path)

In [6]:

from subprocess import run

executable = 'rnoise.out'
script = 'rnoise.f'

output_text = run([
    'gfortran',
    '-o',
    executable,
    script,])

In [7]:

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(rho),
    str(epsrel),
    str(xstart),
    str(xstop),
    str(n_datapoints),
], stderr=STDOUT)
with open(output_text_path, 'wb') as output_file:
    output_file.write(output_text)

In [8]:

print('output_text_path = ' + output_text_path)
print('target_text_path = ' + target_text_path)

output_text_path = /tmp/rnoise.log
target_text_path = /tmp/rnoise.txt

Pay Notebook Creator: Salah Ahmed	0
Set Container: Numerical CPU with TINY Memory for 10 Minutes	0
Total	0

American Public Power Association tools

Radio Noise Analysis Program¶