Modified this so it only includes functions that use ArcPy.  Some script similar to this should be written to get basin latitude (in decimal degrees) and basin area (km2).