add create_facetdomain and add boundary term to spin diffusion