Adding massprojection operator to dgx subdomain
moving flux projection to a function and fixing fault finding bug