import math
import numpy
import ephem
import argparse
from lsl.common import stations, metabundle
from lsl.reader.ldp import LWA1DataFile
from lsl.correlator import fx as fxc
from lsl.astro import unix_to_utcjd, DJD_OFFSET
from lsl.misc import parser as aph
import matplotlib.pyplot as plt
from lsl.misc import telemetry
telemetry.track_script()
def _best_freq_units(freq):
"""Given a numpy array of frequencies in Hz, return a new array with the
frequencies in the best units possible (kHz, MHz, etc.)."""
# Figure out how large the data are
scale = int(math.log10(freq.max()))
if scale >= 9:
divis = 1e9
units = 'GHz'
elif scale >= 6:
divis = 1e6
units = 'MHz'
elif scale >= 3:
def test_script(self):
"""Test telemetry tracking of scripts"""
telemetry.track_script()
