diff --git a/Examples/FlowAroundSphere.py b/Examples/FlowAroundSphere.py
index 0a9ffa1ad9b601dda7e02ae0e8e5b6898aac9b83..6a5afdc171015ded3fec41f0c60f625cbb81df1c 100644
--- a/Examples/FlowAroundSphere.py
+++ b/Examples/FlowAroundSphere.py
@@ -66,7 +66,7 @@ d3d = Discretization([Nx, Ny, Nz])
 # ====== Sphere inside the domain ======
 RADIUS = 0.5
 pos = [0., 0., 0.]
-from hysop.domain.subsets.sphere import Sphere, HemiSphere
+from hysop.domain.subsets import Sphere, HemiSphere
 #sphere = Sphere(origin=pos, radius=RADIUS, parent=box)
 sphere = HemiSphere(origin=pos, radius=RADIUS, parent=box)
 
@@ -168,7 +168,7 @@ op['vort_absorption'] = AbsorptionBC(velo, vorti, discretization=topofft,
 op['vort_absorption'].discretize()
 
 # =====  Penalization of the vorticity on a sphere inside the domain =====
-from hysop.operator.penalize_vorticity import PenalizeVorticity
+from hysop.operator.penalization import PenalizeVorticity
 op['penalVort'] = PenalizeVorticity(velocity=velo, vorticity=vorti,
                                     discretization=topo_with_ghosts,
                                     obstacles=[sphere], coeff=1e8,
@@ -201,7 +201,7 @@ io_ener = IOParams('energy_enstrophy')
 monitors['energy'] = EnergyEnstrophy(velo, vorti, discretization=topofft,
                                      io_params=io_ener, is_normalized=False)
 
-from hysop.domain.subsets.control_box import ControlBox
+from hysop.domain.control_box import ControlBox
 from hysop.operator.drag_and_lift import MomentumForces, NocaForces
 ref_step = topo_with_ghosts.mesh.space_step
 cbpos = npw.zeros(dim)