Implementation of mass properties

ya, I think you would want to cover all possibilities: subs & equip, subs & no equip, no subs & equip

if len(subsystems) > 0:
   ...

for name in el:
    ...

moments = [self.Mx, self.My, self.Mz]
self.Cg = [moment/self.dry_mass for moment in moments]