Hi all,
I used the command “fix 100 flow ave/spatial 1 10000 10000 x 1.0 4.0 z 1.0 4.0 vx density/number file spatial_velocity.out” and in the file “spatial_velocity.out” the third line is “# Bin Coord1 Coord2 Ncount z 1.0” but not “# Bin Coord1 Coord2 Ncount vx density/number” as expected, so I found that in fix_ave_spatial.cpp there is an error in the following lines:
if (ndim == 1) fprintf(fp,"# Bin Coord Ncount");
else if (ndim == 2) fprintf(fp,"# Bin Coord1 Coord2 Ncount");
else if (ndim == 3) fprintf(fp,"# Bin Coord1 Coord2 Coord3 Ncount");
for (int i = 0; i < nvalues; i++) fprintf(fp," %s",arg[9+i]);
the last line should be
for (int i = 0; i < nvalues; i++) fprintf(fp," %s",arg[6+ndim*3+i]);
Best regards
Wenjing Zhou