13 View<double****, CLayout, HostType>
pv_grid;
14 View<double****, CLayout, HostType>
pv_ptl;
49 int period_default = (sml.
is_XGCa ? 10 : 1);
52 bool is_on = nlr.
get<
bool>(
"diag_f0_df_on", is_on_default);
53 int period = nlr.
get<
int>(
"diag_f0_df_period", period_default);
56 std::string name =
"diagnosis.fsource";
59 int n_initial_sources = 1;
Definition: diag_f0_df.hpp:39
void port1(const DomainDecomposition< DeviceType > &pol_decomp, const Grid< DeviceType > &grid, const VelocityGrid &vgrid, Plasma &plasma, const VGridDistribution< HostType > &f0_f_h, const VGridDistribution< HostType > &f0_df0g_h, std::string source_name)
Definition: diag_f0_df.cpp:313
void clear_active_sources()
Definition: diag_f0_df.hpp:28
void write(const Simulation< DeviceType > &sml, const Grid< DeviceType > &grid, const MagneticField< DeviceType > &magnetic_field, Plasma &plasma)
Definition: diag_f0_df.cpp:3
int n_active_sources() const
Definition: diag_f0_df.hpp:19
T get(const string ¶m, const T default_val, int val_ind=0)
Definition: NamelistReader.hpp:373
Definition: diag_f0_df.hpp:9
void set_ptl_baseline(const Simulation< DeviceType > &sml, const Grid< DeviceType > &grid, const MagneticField< DeviceType > &magnetic_field, Plasma &plasma, const VelocityGrid &vgrid, const DomainDecomposition< DeviceType > &pol_decomp)
Definition: diag_f0_df.cpp:290
Definition: velocity_grid.hpp:8
static constexpr bool is_XGCa
Equivalent to the preprocessor flag for now.
Definition: sml.hpp:17
subroutine plasma(grid, itr, p, dene_out, deni_out, Te_out, Ti_out, Vparai_out)
Calculate the plasma density, temperature, and parallel velocity for a point in triangle itr using pl...
Definition: neutral_totalf.F90:1224
Definition: NamelistReader.hpp:193
Definition: magnetic_field.hpp:13
Definition: diag_f0_df.hpp:37
void init(NLReader::NamelistReader &nlr, const Simulation< DeviceType > &sml, const Grid< DeviceType > &grid, const Plasma &plasma, int f_source_period)
Definition: diag_f0_df.hpp:45
void init(const std::string &name, int period, int first_step=-1)
Definition: diagnostic.hpp:46
bool is_on
Definition: diagnostic.hpp:20
Definition: diag_f0_df.hpp:40
void use_namelist(const string &namelist)
Definition: NamelistReader.hpp:355
void calculate_ptl_change(const Simulation< DeviceType > &sml, const Grid< DeviceType > &grid, const MagneticField< DeviceType > &magnetic_field, Plasma &plasma, const VelocityGrid &vgrid, const DomainDecomposition< DeviceType > &pol_decomp)
Definition: diag_f0_df.cpp:303
int n_nonadiabatic_species
Number of nonadiabatic species.
Definition: plasma.hpp:95
Definition: diag_f0_df.hpp:36
View< double ****, CLayout, HostType > pv_ptl
Definition: diag_f0_df.hpp:14
bool f0_grid
Whether to use f0 grid.
Definition: sml.hpp:64
Definition: diag_f0_df.hpp:35
Definition: magnetic_field.F90:1
int get_num_cpu_threads()
Definition: globals.hpp:17
Definition: plasma.hpp:14
Definition: diagnostic.hpp:10
View< double ****, CLayout, HostType > pv_grid
Definition: diag_f0_df.hpp:13
int tindex_period
Definition: diag_f0_df.hpp:16
std::vector< std::string > source_names
Definition: diag_f0_df.hpp:11
int nnode
Number of grid nodes.
Definition: grid.hpp:154
void add_active_source(const std::string &new_source_name)
Definition: diag_f0_df.hpp:24
DiagF0DF()
Definition: diag_f0_df.hpp:43
Definition: diag_f0_df.hpp:38