current_drive.hpp

Go to the documentation of this file.

#ifndef CURRENT_DRIVE_HPP
#define CURRENT_DRIVE_HPP
 
#include "toroidal_average.hpp"
#include "moments.hpp"
#include "loop_voltage_diag.hpp"
#include "source.hpp"
 
class CurrentDrive : public Source{
    public: // temporary, for checkpointing
    GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin> jpara_diff_previous_h;
    GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin> jpara_diff_integral_h;
    View<double*,CLayout,HostType> jpar0_h; // Equilibrium current \f$j_\parallel=\hat{\boldsymbol{b}}\cdot(\nabla\times\boldsymbol{B})/\mu_0\f$
 
    double loop_voltage_psimax; 
 
    bool loop_voltage_fsa;      
 
    double current_drive_p;     
    double current_drive_i;     
    double current_drive_d;     
 
    public:
 
    CurrentDrive(){}
 
    CurrentDrive(NLReader::NamelistReader& nlr, const Grid<DeviceType>& grid, const MagneticField<DeviceType>& magnetic_field);
 
    void apply(const Grid<DeviceType>& grid, Plasma& plasma, Smoothing& smoothing,
                     const Moments& f0_moments, const DomainDecomposition<DeviceType>& pol_decomp,
                     const MagneticField<DeviceType>& magnetic_field,
                     const GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin>& loop_voltage_h,
                     LoopVolDiagnostics& loop_vol_diag, int gstep, double time);
};
 
#endif