Phase 0 works for PIC in Asciifile but not DEM

MFIX-Exa Bug report
1

This isn’t really a bug report but I’m going to report it anyway. I was running my tenCate single particle terminal velocity case that’s basically tests/DEM01-*. It was for PIC; ran fine. I switched it over to DEM and it was crashing b/c it couldn’t calculate an ave dp to calculate kn, b/c it didn’t have any particles, b/c it deleted my one particle. It took me too long to figure out that it was being deleted because the first character on the first data line was 0 but it’s 1 in the tests. I guess I just figured we indexed phase (state?) from 0. ( RIP FORTRAN :headstone:). Anyway… that’s it. Not really a bug, just strange 0 works for PIC but not for DEM.

1
0  0.25  0.13  0.25  7.5e-03  1120.0  0.0  0.0  0.0
2

You are correct in that the type indexing starting at 1 (and not 0), traces back to the humble FORTRAN origins of the code. This could probably be improved so that the processing routine throws an error when an invalid index is read rather than failing silently.

3

pour_one_out