move computation of normal derivatives field to subroutine
add first version of basis-function based construction of us
move compute_normals to subroutine
bugfix: parallelization problems in spectral operators when runnning in 2nd order mode