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"
#include "diag_f0_df.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, const VGridDistribution<DeviceType>& f0_f, VGridDistribution<DeviceType>& f0_df0g_cumul, double dt, int gstep, DiagF0DF& diag_f0_df);
};
 
#endif