pellet.hpp

Go to the documentation of this file.

#ifndef PELLET_HPP
#define PELLET_HPP

#include "grid.hpp"
#include "magnetic_field.hpp"
#include "domain_decomposition.hpp"
#include "velocity_grid.hpp"
#include "plasma.hpp"
#include "vgrid_distribution.hpp"
#include "source.hpp"

class Pellet : public Source{
    double pellet_etemp;
    double pellet_itemp;
    double pellet_n_atoms;
    int ion_species;
    double rate_abl_dv; // Ablation rate per volume unit
    View<bool*,CLayout,HostType> use_node_h;

    double n_atoms_ionized; //<! This value accumulates

    public:

    Pellet(){}

    Pellet(NLReader::NamelistReader& nlr, const Grid<DeviceType>& grid, Plasma& plasma, const MagneticField<DeviceType>& magnetic_field);

    void apply(const Grid<DeviceType>& grid, const DomainDecomposition<DeviceType>& pol_decomp, Plasma& plasma, const VelocityGrid& vgrid, VGridDistribution<HostType>& df0g_h, double dt);
};

#endif