From: Kyunghoon L. <aer...@gm...> - 2013-03-28 06:07:42
|
Dear all, I'd like to make a slight change to transient_ex1 such that half of the domain has a low diffusivity and the other half of the domain has a high diffusivity. I wonder how I can get x coordinates of a quadrature. I tried putting if (qface_points[qp](0) < 0.0) { // subdivice the domain w.r.t x = 0.0 diff_coeff = 1.0; // high diffusivity } else { diff_coeff = 0.5; // low diffusivity } in the quadrature for loop, but it seemed not working. Any adivce will be appreicated. Best regards, K. Lee. |