quick fix

src/schemes/boundary/open_boundary/system.jl

@@ -263,24 +263,23 @@ end
 @inline function convert_particle!(system::OpenBoundarySPHSystem, fluid_system,
                                    boundary_zone::BoundaryZone{BidirectionalFlow}, particle,
                                    v, u, v_fluid, u_fluid)
-    particle_position = particle_coords - boundary_zone.zone_origin
+    particle_position = current_coords(u, system, particle) - boundary_zone.zone_origin
 
     # Check if particle is in- or outside the fluid domain.
     # `plane_normal` is always pointing in the fluid domain.
-    if dot(particle_position, boundary_zone.plane_normal)
-
-        # Activate a new particle in simulation domain
-        transfer_particle!(fluid_system, system, particle, v_fluid, u_fluid, v, u)
-
-        # Reset position of boundary particle
-        for dim in 1:ndims(system)
-            u[dim, particle] = boundary_zone.plane_normal[dim]
-        end
+    if signbit(dot(particle_position, boundary_zone.plane_normal))
+        deactivate_particle!(system, particle, u)
 
         return system
     end
 
-    deactivate_particle!(system, particle, u)
+    # Activate a new particle in simulation domain
+    transfer_particle!(fluid_system, system, particle, v_fluid, u_fluid, v, u)
+
+    # Reset position of boundary particle
+    for dim in 1:ndims(system)
+        u[dim, particle] = boundary_zone.plane_normal[dim]
+    end
 
     return system
 end

src/visualization/write2vtk.jl

@@ -329,7 +329,7 @@ function write2vtk!(vtk, v, u, t, system::OpenBoundarySPHSystem; write_meta_data
                        for particle in active_particles(system)]
 
     if write_meta_data
-        vtk["boundary_zone"] = first(typeof(boundary_zone).parameters)
+        vtk["boundary_zone"] = type2string(first(typeof(system.boundary_zone).parameters))
         vtk["width"] = round(system.boundary_zone.zone_width, digits=3)
         vtk["velocity_function"] = type2string(system.reference_velocity)
         vtk["pressure_function"] = type2string(system.reference_pressure)