Add constant for invalid surface tokens. #5001

Job	Run time
Python 3.11 (omp=n, mpi=n, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	1h 9m 38s
Python 3.11 (omp=y, mpi=n, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	1h 26m 39s
Python 3.11 (omp=n, mpi=y, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	1h 10m 13s
Python 3.11 (omp=y, mpi=y, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	1h 35m 43s
Python 3.12 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	49m 20s
Python 3.13 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	52m 36s
Python 3.11 (omp=y, mpi=y, dagmc=y, ncrystal=, libmesh=, event= vectfit=)	1h 49m 59s
Python 3.11 (omp=n, mpi=n, dagmc=, ncrystal=y, libmesh=, event= vectfit=)	1h 8m 39s
Python 3.11 (omp=y, mpi=y, dagmc=, ncrystal=, libmesh=y, event= vectfit=)	1h 34m 42s
Python 3.11 (omp=y, mpi=n, dagmc=, ncrystal=, libmesh=y, event= vectfit=)	1h 44m 5s
Python 3.11 (omp=y, mpi=n, dagmc=, ncrystal=, libmesh=, event=y vectfit=)	1h 51m 2s
Python 3.11 (omp=n, mpi=y, dagmc=, ncrystal=, libmesh=, event= vectfit=y)	1h 7m 49s
finish	2s
	16h 20m 27s

Provide feedback