Hi julia.hartig. The problem has been solved. I modified two places:
1.Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation
For ep_s(IJK,M), M means phase number. But I input particle number before.
2.ERROR des/read_particle_input.f:329
Because I check “user scalar tracking” in GUI, the beginning of particle_input.dat file should be modified like this:
Then everything is fine.