[freeflow][momentum][fcstaggered] Check dimensions of neumann flux

The dimension is non-standard in the sense that we only have one dof per face but require the problem
for convenience to return a full dimension vector for each coordinate direction.