Fix errors, code follows equivalent to old fortran for main coils
include aux current, exact agreement with previous fortran code