What needs to be done for Simulation

1 separate the plot module after the simulation end at t_max or turn_max

2 use OpenMP for parallel computation with omp.h head file

3 calculate the emittance for x-x' and center of x-direction