Problems with rebind_disk
Hi,
I am trying to do the following: Take an existing disc with a resolution of 1024x3000 (radii x sectors) and rebind it to do an azimuthal average, to end up with a disc with a resolution of 1024x1.
I've used make_grid
to make a 1024x1 grid, and rebind_disk
to make the new disc, then I've used fargoInit
to get a disk0.h5
file to run an integration. I checked the radial surface density profile from disk0.h5
and compared it to the original 1024x3000 file. The two are similar except that there are some unexpected bumps which appear in disk0.h5
, which should not be there. I've tried also with 1024x2 and 1024x4 grids (maybe only 1 sector is not good) but I get similar results (with bumps at different locations).
Why is this happening? Is this a bug?
Thank you,
Gabriele