Add a mean-field option to gaussian_bfield