How to take the atmospheric neutron spectrum as the input, that is, the incident particles. The abscissa is energy and the ordinate is fluence rate. There are no reference examples. To write their own related functions, how to write. Can both GPS and pariclegun be implemented?