Initial implementation.
Further steps towards implementation, but it doesn't work. Schmidt functions seem OK for orders 1 and 2 then slowly become worse and worse. But using the order 2 model still doesn't get the magnetic field values I would expect.