[Nek5000-users] 3D cylinder mesh converting error

[nek5000-users at lists.mcs.anl.gov]

Hi! nek,
I met a trouble in using the .msh to .rea convert program mshconvert. When I tried to convert a 3D cylindrical mesh, I input the command  as follow:

convert('cylRB.msh',curves={4:{'type':'C','radius':0.25 ,'depth': 4,'circle_center':(0,0,0)}},bcs={3:'W',4:'W'},temperature={3:'t',4:'I'})

and then an error appeared on the screen :

/home/xbl/case/mshconvert.pyc in read_curved_sides(curves)
    769             cell = max(face[2])
    770             nds = face[1]
--> 771             cell_face_m.append((cell, face_map[cell][(nds[0], nds[1])], face_number))
    772             if curve['type'] == 'C':
    773                 for jj in range(curve['depth']):

KeyError: (6741, 97)

'4' in the command means that the curved surface in the mesh is zone 4,says the lateral surface of the cylinder.

I succeeded in converting a 2D-circle mesh before, so I want to know how to input correct command for 3D case.

Does anyone meet that kind of error ever? Any help is appreciated.


regards,
Xu




--
Xu Bolun,
University of Science and Technology of China,
Hefei,Anhui,China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161019/f1958870/attachment.html>