electronicFields = models.ElectronicFields('elec')
elec_model = models.ElectronicModelA(geomFields,electronicFields,meshes)
fgs = mesh0.getBoundaryGroups()
testPoint = VecD3()
testPoint[0] = 0.0
testPoint[1] = 0.0
testPoint[2] = 0.0
#import ddd
octreeCells = fvmbaseExt.Octree()
octreeCells.Impl(mesh0, geomFields)
######### search top surface ####################
bcID = 6
octreeTop = fvmbaseExt.Octree()
octreeTop.Create(mesh0, geomFields, bcID)
nearestPoint = elec_model.findClosestPoint(testPoint, octreeTop)
for fg in fgs:
if fg.id == bcID:
if atype == 'tangent':
metricsCalculator.setTangentCoords(0, 7, 1)
flowFields = models.FlowFields('flow')
fmodel = models.FlowModelA(geomFields, flowFields, meshes)
reader.importFlowBCs(fmodel)
fmodel.init()
import time
t0 = time.time()
solid = fvmbaseExt.MPM()
octree = fvmbaseExt.Octree()
option = 2
fvmbaseExt.CellMark_Impl(mesh0, geomFields, fileBase, octree, solid, option)
cells = mesh0.getCells()
nCells = cells.getSelfCount()
print nCells
t1 = time.time()
if outfile != '/dev/stdout':
print '\nsolution time = %f' % (t1 - t0)