def cylinder(center, R, H, N=100, ntype='STRUCT'): """Create a cylinder of center C, radius R and hieght H. Usage: cylinder((xc,yc,zc), R, H, N)""" A = Geom.cylinder(center, R, H, N, ntype) if ntype == 'STRUCT': return [ C.convertArrays2ZoneNode('cyl-part1', [A[0]]), C.convertArrays2ZoneNode('cyl-part2', [A[1]]), C.convertArrays2ZoneNode('cyl-part3', [A[2]]), C.convertArrays2ZoneNode('cyl-part4', [A[3]]), C.convertArrays2ZoneNode('cyl-part5', [A[4]]), C.convertArrays2ZoneNode('cyl-part6', [A[5]]), C.convertArrays2ZoneNode('cyl-part7', [A[6]]), C.convertArrays2ZoneNode('cyl-part8', [A[7]]), C.convertArrays2ZoneNode('cyl-part9', [A[8]]), C.convertArrays2ZoneNode('cyl-part10', [A[9]]), C.convertArrays2ZoneNode('cyl-part11', [A[10]]) ] else: return C.convertArrays2ZoneNode('cyl', [A])
# - cylinder (array) - import Geom as D import Converter as C a = D.cylinder((0,0,0), 1., 10.) C.convertArrays2File(a, 'out.plt')