# - reorderAll (pyTree) - import Converter.PyTree as C import Generator.PyTree as G import Transform.PyTree as T ni = 30 nj = 40 nk = 1 m1 = G.cart((0, 0, 0), (10. / (ni - 1), 10. / (nj - 1), 1), (ni, nj, nk)) m1[0] = 'cart1' m2 = T.rotate(m1, (0.2, 0.2, 0.), (0., 0., 1.), 15.) m2 = T.reorder(m2, (-1, 2, 3)) m2[0] = 'cart2' t = C.newPyTree(['Base', 2]) t[2][1][2] += [m1, m2] t = T.reorderAll(t, 1) C.convertPyTree2File(t, "out.cgns")
# - selectCells (pyTree) - import Converter.PyTree as C import Generator.PyTree as G import Post.PyTree as P import KCore.test as test def F(x, y, z): if (x + 2 * y + z > 20.): return True else: return False # CAS 1D a = G.cart((0, 0, 0), (1, 1, 1), (30, 1, 1)) a = C.addVars(a, 'Density') a = C.addVars(a, 'centers:cellN') a = P.selectCells(a, F, ['CoordinateX', 'CoordinateY', 'CoordinateZ']) test.testT(a, 1) # CAS 2D a = G.cart((0, 0, 0), (1, 1, 1), (11, 11, 1)) a = C.addVars(a, 'Density') a = C.addVars(a, 'centers:cellN') a = C.addBC2Zone(a, 'wall', 'BCWall', 'imin') a = P.selectCells(a, F, ['CoordinateX', 'CoordinateY', 'CoordinateZ']) test.testT(a, 2) # CAS 3D a = G.cart((0, 0, 0), (1, 1, 1), (11, 11, 11)) a = C.addVars(a, 'Density') a = C.addVars(a, 'centers:cellN')
# - rotate (PyTree) - import Generator.PyTree as G import Transform.PyTree as T import Converter.PyTree as C a = G.cart((0, 0, 0), (1, 1, 1), (10, 10, 2)) # Rotate with an axis and an angle b = T.rotate(a, (0., 0., 0.), (0., 0., 1.), 30.) b[0] = 'cartRot1' # Rotate with two axis c = T.rotate(a, (0., 0., 0.), ((1., 0., 0.), (0, 1, 0), (0, 0, 1)), ((1, 1, 0), (1, -1, 0), (0, 0, 1))) c[0] = 'cartRot2' # Rotate with three angles c = T.rotate(a, (0., 0., 0.), (0, 0, 90)) c[0] = 'cartRot3' C.convertPyTree2File([a, b, c], 'out.cgns')
# - surfaceWalk (pyTree) import Converter.PyTree as C import Geom.PyTree as D import Transform.PyTree as T import Generator.PyTree as G # User definition of parametric curve def f(t, u): x = t+u y = t*t+1+u*u z = u return x, y, z # Array definition of geometry a = D.surface(f) c = D.circle((1.2, 1.7, 0.6), 0.1) c = T.rotate(c, (1.2, 1.7, 0.6), (0, 1, 0), 90.) c = T.reorder(c, (-1, 2, 3)) c = T.projectOrtho(c, [a]) h = G.cart((0., 0., 0.), (0.01, 1, 1), (15, 1, 1)) r = G.surfaceWalk([a], c, h, niter=100) C.convertPyTree2File(r, "out.cgns")