Program received signal SIGFPE

Dear professors,
In my source.f, I deal with both single primary and multiple primaries. I follow the guide in manual which tells me that I should use a DO cycle until SOEVSV. However, I still encounter the error:
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation. I have uploaded related files(I remove some sources in muon1326.txt in order not to exceed the maximum size). Any help will be appreciated!
Best regards,
Hechong Han
Te-LS_1326.inp (3.3 KB)
Te-LS_1326001.err (714 Bytes)
Te-LS_1326001.log (559 Bytes)
Te-LS_1326001.out (104.1 KB)
usrrnc.f (1.4 KB)
mgdraw.f (7.5 KB)
source.f (7.5 KB)
muon1326.txt (3.9 MB)

Dear @hanhechong

some of the lines in the muon1326.txt contains lines, where some parameters are undefined, like:

-13 nan nan nan 0.105658 0 -34250.000000 9135.286416 -14865.904643 

Removing these lines should fix the crashes.

Cheers,
David

Dear @horvathd
Thank you for your help! I did find some nan in the input and fixed it.
Best regards,
Hechong Han