Thank you very much. :) I'll try this tomorrow and post the results.. It's definitely neater than mymod = __import__('make_sys_file') and getattrib(mymod, relevant_make_sys_func) with make_sys_file.py containing 100 different make_sys functions.
The application is about an interaction with a big external spin of value J, and the number of lattices is essentially 2J+1. I want to be able to see the effect of J vs. transmission etc.
Also, I've got to ask: why do you want to do this? I'm struggling to
think of an application where you would have ~100 separate channels like
that that are all separately measureable. This is probably due to a lack
of imagination on my part ;)