3D_Surf_Export

{

   This problem shows data export on extrusion surfaces in 3D.

   It is a modification of 3D_SPHERE.PDE.

}

title '3D Export Test - Sphere'

coordinates

    cartesian3

Variables

    u

definitions

    k = 0.1                    { conductivity }

    heat =6*k                  { internal heat source }

   u0 = exp(-x^2-y^2)

equations

    div(K*grad(u)) + heat   = 0

extrusion

   surface z = -sqrt(1-(x^2+y^2))      { the bottom hemisphere }

   surface z = sqrt(1-(x^2+y^2))       { the top hemisphere }

boundaries

    surface 1 value(u) = u0     { fixed value on sphere surfaces }

    surface 2 value(u) = u0

    Region 1

       start(1,0) arc(center=0,0) angle=360

plots

    grid(x,y,z)

    contour(u) on x=0          { YZ plane through diameter }

       export

    contour(u) on z=0.5        { XY plane above center }

               export format "#x#b#y#b#z#b#1"

    contour(u) on surface 2            { top surface }

                export format "#x#b#y#b#z#b#1"

end