Hello, would it be useful to define the CarpetRegrid2 variable as positions[num_centres] instead of positions[10] ? In this way a little output (and a little memory) would be spared. Luca