sphdem variational compiles but omega (and kappa?) not being calculated correctly