Compute Non-standard Statistics

Allow computation of statistics not directly implemented by the program
or the np.ma package.

Added two options:
add_stats = [‘stat1’, ‘stat2’]   # List of statistics to be computed
add_func=‘myfunc’    # Module where the functions that compute these
statistics are defined
Reading that web page actually gave me the idea of using the list as a mutable object...see the problem I had below