Extrudes a 2D mesh to 3D. If the 2D mesh has triangular elements the 3D elements will be prisms. If the 2D mesh has quadrilaterals then the 3D elements will be hexahedral.
Version : 1.0
Author : George Kourakos
Department of Land Air and Water
University of California Davis
[p MSH]=extrude_mesh(p_2d, MSH_2d, top, bot, t, el_order)
p_2d: Coordinates of 2D mesh nodes
MSH_2d: Mesh structure generated by Gmsh.
top: top elevation of the 2D mesh nodes
bot: bottom elevation of the 2D mesh nodes
t: parametric variable that defines how many layers will be generated and their vertical distribution between the top and bottom elevation. The t values must be t(1) = 0 (corresponds to bottom elevation) and t(end) = 1 (coresponds to top elevation) A quick way to plot the layer distribution is plot(ones(length(t),1),t,'.-')
el_order:element order 'linear' or 'quadratic'. At the moment incomplete elements are not supported
p: coordinates of 3D mesh nodes
MSH 3D mesh structure.
see also: Centerfor2points