added option in script for choosing between formulae for multipoles