## [Paramesh-users] Cylindrical coordinates, guard cell filling

 [Paramesh-users] Cylindrical coordinates, guard cell filling From: Gross, Markus - 2011-11-09 12:30:17 ```Hi, Just looking at paramesh for cylindrical/spherical coordinates I noticed that, in my test setup anyway, the guard cells are not filled correctly at r=0. Here we have to take into account the special geometry, i.e. map in reverse order the adjacent block. However, paramesh seems to happily just copy. Has anybody else noticed that? Coded a workaround or knows which flags to set to avoid this? Trying to be a bit more illustrative, if the row at r=0 is filled with sequential numbers such as G G G G G 1 2 3 4 5 Assuming only one block, with five gridpoints, the guardcells (G) should be filled such that 5 4 3 2 1 1 2 3 4 5 But I only get 1 2 3 4 5 1 2 3 4 5 Of course, a bit more tricky with leafs and more blocks. Many thanks for you time! Regards, Markus ```

 [Paramesh-users] Cylindrical coordinates, guard cell filling From: Gross, Markus - 2011-11-09 12:30:17 ```Hi, Just looking at paramesh for cylindrical/spherical coordinates I noticed that, in my test setup anyway, the guard cells are not filled correctly at r=0. Here we have to take into account the special geometry, i.e. map in reverse order the adjacent block. However, paramesh seems to happily just copy. Has anybody else noticed that? Coded a workaround or knows which flags to set to avoid this? Trying to be a bit more illustrative, if the row at r=0 is filled with sequential numbers such as G G G G G 1 2 3 4 5 Assuming only one block, with five gridpoints, the guardcells (G) should be filled such that 5 4 3 2 1 1 2 3 4 5 But I only get 1 2 3 4 5 1 2 3 4 5 Of course, a bit more tricky with leafs and more blocks. Many thanks for you time! Regards, Markus ```
 Re: [Paramesh-users] Cylindrical coordinates, guard cell filling From: Klaus Weide - 2011-11-09 15:23:20 ```On Wed, 9 Nov 2011, Gross, Markus wrote: > Just looking at paramesh for cylindrical/spherical coordinates I noticed > that, in my test setup anyway, the guard cells are not filled correctly > at r=0. Here we have to take into account the special geometry, i.e. map > in reverse order the adjacent block. However, paramesh seems to happily > just copy. > > Has anybody else noticed that? Coded a workaround or knows which flags > to set to avoid this? Markus, Have you tried turning lsingular_line on (set to T in amr_runtime_parameters) ? That may only work for spherical coordinates along the axis, not sure whether it is supposed to work with cylindrical or polar coordinates. I have never tried it myself. Klaus ```