Outputting entire 3D field for larger domains

[DHRUBO ROY]

Hi

I have been recently started using basilisk for rectangular domain (nx= 15, ny =30, nz =30). To take the output for the entire 3D field I have used a code to get output using reduction operation.

Using reduction operations limits the number of grid points due to restriction in RAM usage per core. Is there any code to get the entire 3D output for a scalar .

Regards 

Dhrubo