I've ported a simple finite difference time domain code into Numpy, and it was thoroughly enjoyable to see that it runs at about 1/2 the speed of C. Are others doing similar things? I want to try to port a hybrid FEM-MOM code, but that going to take a lot of effort. Rob. "EM Computation- where God licenses the source code"