swap density pressure

Author: LasNikas

src/schemes/boundary/open_boundary/mirroring.jl

@@ -136,13 +136,6 @@ function extrapolate_values!(system, v_open_boundary, v_fluid, u_open_boundary,
             end
         end
 
-        if prescribed_density
-            density[particle] = reference_value(reference_density, density[particle],
-                                                particle_coords, t)
-        else
-            inverse_state_equation!(density, state_equation, pressure, particle)
-        end
-
         if prescribed_pressure
             pressure[particle] = reference_value(reference_pressure, pressure[particle],
                                                  particle_coords, t)
@@ -152,6 +145,13 @@ function extrapolate_values!(system, v_open_boundary, v_fluid, u_open_boundary,
 
             pressure[particle] = f_d[1] + dot(pos_diff, df_d)
         end
+
+        if prescribed_density
+            density[particle] = reference_value(reference_density, density[particle],
+                                                particle_coords, t)
+        else
+            inverse_state_equation!(density, state_equation, pressure, particle)
+        end
     end
 
     if !(prescribed_velocity) && boundary_zone.average_inflow_velocity