{
This problem shows some of the possible 'ON' qualifiers for 3D plots.
}
title '3D Test -- Plot Qualifiers'
coordinates
cartesian3
Variables
u
definitions
k = 0.1
heat = 4
equations
div(K*grad(u)) + heat = 0
extrusion
surface z = 0
surface z = 0.8-0.3*(x^2+y^2)
surface z = 1.0-0.3*(x^2+y^2)
boundaries
Region 1 'outer'
layer 2 k = 1
start(-1,-1)
value(u) = 0
line to (1,-1) to (1,1) to (-1,1) to close
Region 2 'plug'
layer 2 k = 1
start 'dot' (0.5,0.5) arc(center=0,0) angle=360
plots
grid(x,y,z) on region 1
grid(y,z) on x=0 on 'plug'
contour(u) on x=0.51 on layer 2
contour(u) on z=0.51 on region 2
contour(u) on surface 2 on region 2 { on paraboloidal layer interface }
vector(grad(u)) on surface 2 on 'outer'
end