def axisym(t, center, axis, angle=360., Ntheta=360, rmod=None):
"""Create an axisymmetric mesh given an azimuthal surface mesh.
Usage: axisym(t, (xo,yo,zo), (nx,ny,nz), teta, Nteta, rmod)"""
# Attention en centres, les coord. des centres ne sont pas justes! mais
# elles ne sont pas utilisees dans la fonction
if rmod is not None:
rmod = C.getFields(Internal.__GridCoordinates__, rmod)[0]
return C.TZAGC(t, 'both', 'both', Geom.axisym, Geom.axisym, center, axis,
angle, Ntheta, rmod, center, axis, angle, Ntheta - 1, rmod)
def orthoDrive(t, line, mode=0):
"""Generate a surface mesh by using 1D array (defining a mesh)
and following orthogonally the curve defined in line.
Usage: orthoDrive(t, line)"""
al = C.getFields(Internal.__GridCoordinates__, line)
if len(al) == 1: al = al[0]
al2 = Converter.node2Center(al)
# Attention les coord. des centres ne sont pas justes! mais
# elles ne sont pas utilisees dans la fonction
return C.TZAGC(t, 'both', 'both', True, Geom.orthoDrive, Geom.orthoDrive,
al, mode, al2, mode)
