Fixed theta interpolation for stokes operator

[abarret]

There was a bug in the interpolation of theta to cell nodes. I've also included some more input files that test different options of configurations.