Updated with tower rho calculation
Updated to add region functionality
fixed up the division so numpy doesn't complain