Hi Daniel,
Thank you for spending so much time on replying my question. I try to
simplify my problem. Are the following right statements?
-------------------------------------------------------
mesh3D = Grid3D(dx, dy, dz, nx, ny, nz)
mesh2D = Grid2D(dx, dy, nx, ny)
var3D = CellVariable(mesh = mesh3D, value = random.random(nx * ny *nz))
var2D = CellVariable(mesh = mesh2D)
x3D, y3D, z3D =mesh3D.getCellCenters()
x2D, y2D = mesh2D.getCellCenters()
var2D.setValue(max([z3D, where=(var3D>0.5)]))
----------------------------------------------------
The expansion into 3D is exactly what I need, if this var2D is good, .
To put data from 2D mesh into 3D mesh:
var3D.setValue(var2D(var3D.getMesh().getCellCenters()[:2], order=1))
Best Regards,
Kundar
