A small problem in Incompact3d's postprocess

[jushi-123]

Dear incompact3d developers,

I have a question to consult you, the subroutines "init_coarser_mesh_statV" in post.f90 enable to write the whole array of statistics. Now, if I have an array "t11" which is xsize(1)*xsize(2)*xsize(3), but I want to write the part of this array(at every other grid point in x and z directions). So the final array is xsize(1)/2 * xsize(2) * xsize(3)/2, how can this requirement be implemented?

I look forward to your reply, thanks!

[mathrack]

Hello,

The safe way is to define a decomp_info object for the target grid and write your own interpolation subroutine to interpolate the data from the main grid to the target one. Then call decomp_2d_write_one and provide the decomp_info object of the target grid. See https://github.com/2decomp-fft/2decomp-fft/blob/main/src/io.f90#L169

You can also try to tune the call to the subroutines init_coarser_mesh_stat*. At your own risk.