compiletime load "module_ppp.so";
case = BATCH.file;
 
X0 = 115;   X1 = 635;    # Cylinder elongation in X direction
Y0 = 167.5; Z0 = 247;    # Center Coordinates
R0 = 70;                 # Cylinder radius
Delta = 5;               # Search radius
 
NL   = 60;               # Point resolution in longitudinal direction
NPHI = 75;               # Point resolution in angular direction
 
PPP obj = { };
obj.read("$case/field/B_0.00us.pos", "B");
 
obj.cylinder(X0, Y0, Z0,   X1, Y0, Z0,    R0, Delta,  NL, NPHI, "interpolate");
 
obj.write_pos_scalar_binary("extract/$case-B-radial.pos", "B radial [mT]", B_r*1000);
obj.write_pos_scalar_binary("extract/$case-B-parallel.pos", "B parallel [mT]", sqrt(B_x^2 + B_phi^2)*1000);