Cleaned up some tests
powermean() can now take a vector of qs, so qD() and qDZ() don't need to handle that separately.
Fixed some comments, and corrected powermeans used to do ecosystem A and ecosystem A bar calculations.
Added in sub-community and ecosystem gamma diversity (γ̄ and γ, Ḡ and G), and some simple matching tests.