Hi Sumaiya,
1. I don't know about the 58GB thing, sounds a bit larger than I expect but you could be right. When you visualize it in ParaView, consider visualizing the points only without generating gpyph; or if you would like glyph, keep the mode being "All Points", but reduce the Theta and Phi resolution, to maybe 3 or something like that, and that will greatly reduce the memory used. However, if the output size keeps growing, at some point you'll have to use some sort of MPI or partitioned view.
2. Sorry for the lack of comments. Part3 is still about preparing a GRC-3 particle bed. It requires finishing Part2. It makes several copies of the results from Part2, put them side by side to make a larger material patch, then let the gravity do the work. After settling, it compresses the resultant material bed so the top is a bit more even/flat. Then it saves the settled material to a file. Note that this file is very large, representing a 4m × 2m soil bin which would allow a rover running on it. If you don't need a test environment this big, you can safely ignore Part3.
Thank you,
Ruochun