compiletime load "module_ppp.so";

case = "sample2";

PPP obj={};

for (time=100; time<=200; time=time+20) {
  obj.read("$case/density_Ar_$(time).00ms.pos", "d");
  obj.read("$case/velocity_Ar_$(time).00ms.pos", "v");
}

obj.xy_plane(0, 10, "nearest_neighbour");
obj.write_pos_vector_binary("flux-cutplane.pos", "Ar flux [1e18/(m²s)]", 1e-18*d*v_x, 1e-18*d*v_y, 1e-18*d*v_z);