Computing the fields is pretty easy. It's just a matter of integration. I've already got code for creating the field from a delta function ring - it would just be a matter of putting a bunch of rings in the region of the coil shapes suggested.
I'll also see if I can find some visualization stuff to make it easier to see the results.
darn this is too much fun!!!