#! /usr/bin/env python
#
# A module for working with HydroTrend.
#
# Mark Piper (mark.piper@colorado.edu)

import numpy as np


def load_series(output_file):
    '''
    Reads a column of text containing HydroTrend output. Returns a numpy array,
    or None on an error.
    '''
    try:
        series = np.loadtxt(output_file, skiprows=2)
    except (IOError, StopIteration):
        pass
    else:
        return(series)

@@ -100,16 +100,44 @@
 o.edu)%0A%0A
+from subprocess import call%0A
 import n
@@ -161,15 +161,8 @@
 load
-_series
 (out
@@ -438,8 +438,174 @@
 series)%0A
+%0A%0Adef call(input_dir, output_dir):%0A    '''%0A    Invokes HydroTrend through the shell.%0A    '''%0A    call(%5B'hydrotrend', '--in-dir', input_dir, '--out-dir', output_dir%5D)%0A

Add call() function, rename load() function