python: calculating centroid of multi-model PDB using numpy

Question

I am working with multi-model pdb file arranged in the following format:

# example of 3 model pdb file, where each model separated by "ENDMDL"
# the