This is a simple 5 deep nested loop over three spatial and two component indices with nearest neighbor spatial structure and "full" matrices for components |
This problem is defined with two sets of Parameters |
1)Large Grid but not many components: |
nc = 5; nx = 1000; ny = 1000; nz = 1000 |
2)Smaller Grid but more components as if many species: Repeat the above calculation with the parameters |
nc = 150; nx = 100; ny = 100; nz = 100 |