1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
//
// Each problem must have its own version of PhysicalParticleContainer::InitData()
// to initialize the particle data. It must also initialize charge and mass.
//
#include <cmath>
#include <AMReX_BLProfiler.H>
#include <ParticleContainer.H>
#include <WarpXConst.H>
using namespace amrex;
void
PhysicalParticleContainer::InitData()
{
BL_PROFILE("PPC::InitData()");
charge = -PhysConst::q_e;
mass = PhysConst::m_e;
const int lev = 0;
const Geometry& geom = Geom(lev);
const Real* dx = geom.CellSize();
Real weight, ux, uy, uz;
Real particle_xmin, particle_xmax, particle_ymin, particle_ymax, particle_zmin, particle_zmax;
int n_part_per_cell;
{
ParmParse pp("langmuirwave");
n_part_per_cell = 1;
pp.query("num_particles_per_cell", n_part_per_cell);
weight = 1.e25;
pp.query("n_e", weight);
#if BL_SPACEDIM==3
weight *= dx[0]*dx[1]*dx[2]/n_part_per_cell;
#elif BL_SPACEDIM==2
weight *= dx[0]*dx[1]/n_part_per_cell;
#endif
particle_xmin = particle_ymin = particle_zmin = -2.e-5;
particle_xmax = particle_ymax = particle_zmax = 2.e-5;
pp.query("particle_xmin", particle_xmin);
pp.query("particle_xmax", particle_xmax);
pp.query("particle_ymin", particle_ymin);
pp.query("particle_ymax", particle_ymax);
pp.query("particle_zmin", particle_zmin);
pp.query("particle_zmax", particle_zmax);
ux = 0.;
uy = 0.;
uz = 0.;
pp.query("ux", ux);
pp.query("uy", uy);
pp.query("uz", uz);
ux *= PhysConst::c;
uy *= PhysConst::c;
uz *= PhysConst::c;
}
std::array<Real,PIdx::nattribs> attribs;
attribs.fill(0.0);
attribs[PIdx::w ] = weight;
attribs[PIdx::ux] = ux;
attribs[PIdx::uy] = uy;
attribs[PIdx::uz] = uz;
for (MFIter mfi = MakeMFIter(lev); mfi.isValid(); ++mfi)
{
const Box& tile_box = mfi.tilebox();
RealBox tile_real_box { tile_box, dx, geom.ProbLo() };
const int grid_id = mfi.index();
const int tile_id = mfi.LocalTileIndex();
const IntVect& boxlo = tile_box.smallEnd();
for (IntVect iv = boxlo; iv <= tile_box.bigEnd(); tile_box.next(iv))
{
for (int i_part=0; i_part<n_part_per_cell;i_part++)
{
Real particle_shift = (0.5+i_part)/n_part_per_cell;
#if (BL_SPACEDIM == 3)
Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0];
Real y = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1];
Real z = tile_real_box.lo(2) + (iv[2]-boxlo[2] + particle_shift)*dx[2];
#elif (BL_SPACEDIM == 2)
Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0];
Real y = 0.0;
Real z = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1];
#endif
if (x >= particle_xmax || x < particle_xmin ||
y >= particle_ymax || y < particle_ymin ||
z >= particle_zmax || z < particle_zmin )
{
continue;
}
else
{
AddOneParticle(lev, grid_id, tile_id, x, y, z, attribs);
}
}
}
}
}
|
