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/absorption/absorption_Arplus_10.00us.pos", "A");
obj.read("$case/cell_surfaces.pos", "S");

obj.create_scalar("AxS", A*S);
obj.cylinder(X0, Y0, Z0,   X1, Y0, Z0,    R0, Delta,  NL, NPHI, "average");
obj.write_pos_scalar_binary("extract/$case-abs-cyl-Arplus.pos", 
                            "Ar+ absorption [1/m^2 s]", 
                            (S>0) ? AxS / S  ! 0);