loop_voltage_diag.hpp

Go to the documentation of this file.

 
#ifndef LOOP_VOLTAGE_DIAG_HPP
#define LOOP_VOLTAGE_DIAG_HPP
 
#include "grid.hpp"
#include "diagnostic.hpp"
#include "grid_field.hpp"
 
 
class LoopVolDiagnostics : public Diagnostic {
 
    public:
 
    LoopVolDiagnostics() : Diagnostic() {}
 
    void init(NLReader::NamelistReader& nlr, int f_source_period, const std::string& name){
        nlr.use_namelist("diag_param");
        bool is_on = nlr.get<bool>("diag_current_drive_on",false); 
        int loop_vol_diag_period = f_source_period*nlr.get<int>("diag_current_drive_period",1); 
 
        if (is_on) Diagnostic::init(name, loop_vol_diag_period);
    }
 
    void write(const Grid<DeviceType>& grid,
                     const GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin>& loop_voltage_h,
                     const GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin>& jpara_diff_previous_h,
                     const GridField<HostType,VarType::Scalar,PhiInterpType::None,TorType::OnePlane,KinType::DriftKin>& jpara_diff_integral_h,
                     const View<double*, CLayout, HostType>& jpara_diff_current, int gstep, double time);
};
 
#endif