# Print current simulation name
print(sim.Name)

# Remove all sources before adding them
to_del = []
for s in sim.AllSettings:
	if isinstance(s,acoustic.SourceSettings):
		to_del.append(s)
for d in to_del:
	sim.Remove(d)

# Add sources with correct element and parameters
for source_idx in np.arange(len(elements)):
	source_settings = acoustic.SourceSettings()
	source_settings.Amplitude = amp[source_idx]
	source_settings.Phase = phases[source_idx]
	sim.Add(source_settings, rings[source_idx])
	print(rings[source_idx].Name, amp[source_idx], phases[source_idx])

Sorry I'm not providing a full example. Easiest way to do that would be to manually create the whole simulation without the sources (just make all sources as materials), then right click on simulation and select 'To Python...' and then add the lines above to your code

You could do something like this to get the elements as a list:

src_ents = []
entities = model.AllEntities()
for e in entities:
    if 'Element' in e.Name:
        src_ents.append(e)

Careful, I have no idea how the elements will be sorted in that case but you could use a fixed size numerical suffix and then sort them.

If I remember correctly the SEFT tutorial python script shows a similar method to create the SEFT (single element curved transducer) transducer from scratch