search

nbia-astro / skeletor

Parallel PIC code written in Python and based on the skeleton codes provided by PICKSC
GNU General Public License v3.0
1 stars 0 forks source link

The locations of the boundaries of the computational domain should be arbitrary #120

Open tobson opened 7 years ago

tobson commented 7 years ago

Right now the computational domain always goes from (0, 0) to (Lx, Ly). In the shearing sheet, for instance, we'd rather have it go from (-Lx/2, -Ly/2) to (Lx/2, Ly/2).

Probably all it takes to implement this are tiny changes in the particle push and deposit routines, specifically here and here.

tberlok commented 7 years ago

I agree, let's make this change sooner than later.

We would also need to modify particles.py here and the grid x and y arrays but I think that's about it. We will not have to modify the particle boundaries as they are written in the grid distance units which we can keep as [0, nx) x [0, ny).

I will be back in business on Wednesday.