print("Input data = "+species+" from " + data_dir)
print("Total no. of particles =" +str(N_part_total))
print("No. of particles used =" +str(N_part))
#Partition the total input file with N particles into mpi_size parts:
mpi_borders = np.linspace(0, N_part, mpi_size + 1).astype('int')
my_border_low = mpi_borders[mpi_rank]
my_border_high = mpi_borders[mpi_rank+1]
#Redefining N_part as number of particles analysed per processor
N_part = my_border_high - my_border_low
#Read Trajectory Data
time, U, E, B, Energy = Read_Data.Read_Trajectory_Data(my_border_low, my_border_high, data_dir, species)
#Comment: U is used to denote the proper relativistic velocity, which is U = gamma*V
N = len(time)
#Process Trajectory Data using chosen wavelet and edge mode
if mpi_rank==0:
print('Processing Data for particles...')
#Transfrom the data and compute new quantities:
time, V, E, b, V_par, V_perp, vpar, vperp, Acc_mag, Acc_par, Acc_perp1, Acc_perp2 = Process_Data.Process_raw(time, U, E, B, Omega_c, q, m)
N_pad = len(vpar[0,:])
ref=int(np.ceil(N/4))
#vpar_Global=Diffusion(vpar,ref,N)
