bpy.ops.object.mode_set(mode='EDIT')
bpy.ops.mesh.subdivide(number_cuts=xsubds, smoothness=0.0)
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.mode_set(mode='OBJECT')
me.update()
bm.from_mesh(me)

for e in bm.edges:
    if ((e.verts[0].co[1]) != (e.verts[1].co[1])):
        edgeindex = e.index
        obj.data.edges[edgeindex].select = True


bpy.ops.object.mode_set(mode='EDIT')
bpy.ops.mesh.subdivide(number_cuts=ysubds, smoothness=0.0)
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.mode_set(mode='OBJECT')
me.update()
bm.from_mesh(me)

for e in bm.edges:
    if ((e.verts[0].co[2]) != (e.verts[1].co[2])):
        edgeindex = e.index
        obj.data.edges[edgeindex].select = True

bpy.ops.object.mode_set(mode='EDIT')
bpy.ops.mesh.subdivide(number_cuts=zsubds, smoothness=0.0)
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.mode_set(mode='OBJECT')
me.update()
bm.from_mesh(me)

bpy.ops.object.editmode_toggle()
bpy.ops.mesh.subdivide(number_cuts=14, smoothness=0, quadcorner='INNERVERT')
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.editmode_toggle()
bpy.ops.object.editmode_toggle()
bpy.ops.mesh.subdivide(number_cuts=8, smoothness=0, quadcorner='INNERVERT')
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.editmode_toggle()
bpy.ops.object.editmode_toggle()
bpy.ops.mesh.subdivide(number_cuts=10, smoothness=0, quadcorner='INNERVERT')
bpy.ops.mesh.select_all(action='DESELECT')
bpy.ops.object.editmode_toggle()