About identifying different particle types in source code

MFiX
4

This looks reasonable to me. Note that you can write double-precision constants in Fortran with a D instead of E (D for “double”), so you can write 1d-12 instead of 0.000000000001d0