So I implemented the projection method as in Chorin (1969) for the solution of the Navier-Stokes Equations with the Boussinesq Approximation, using a 2nd Order accurate FD scheme. The primitive variables are arranged using a MAC scheme. Basically, what I want is to be able to generate a file (or files) to see the solution on a rendering software (a.k.a, Tecplot, paraView, etc). How can I arrange the data? I mean, from a pseudo-code kind of view.