Functions/Subroutines
subroutine	main_plasma
	Main plasma. More...

subroutine	current (GEOMETRY, PROFILES, TRANSPORT, SOURCES, EVOLUTION, CONTROL, j_boun, ifail, failstring)
	CURRENT TRANSPORT EQUATION. More...

subroutine	ion_density
	PARTICLE TRANSPORT EQUATIONS. More...

subroutine	electron_density (GEOMETRY, PROFILES, TRANSPORT, SOURCES, EVOLUTION, CONTROL, HYPER_DIFF, ifail, failstring)

subroutine	quasi_neutrality (GEOMETRY, PROFILES, IMPURITY, CONTROL, ifail, failstring)
	QUASI NEUTRALITY. More...

subroutine	temperatures
	HEAT TRANSPORT EQUATIONS. More...

subroutine	rotation (GEOMETRY, PROFILES, TRANSPORT, SOURCES, EVOLUTION, CONTROL, ifail, failstring)
	ROTATION TRANSPORT EQUATIONS. More...

subroutine	derivn (N, X, Y, DY1)
	These subroutines calculate first and second derivatives, DY1 and DY2, of function Y respect to argument X. More...

subroutine	integr (N, X, Y, INTY)
	This subroutine calculates integral of function Y(X)*X from X=0 until X=X(N) More...

subroutine	integr2 (N, X, Y, INTY)

subroutine	f_axis (n, r, f)
	This subroutine finds f(r_1=0) from f(r_2), f(r_3) and f(r_4) More...

subroutine	f_par_axis (n, r, f)
	This subroutine finds f(r_1=0) from f(r_2), f(r_3) d/dr f(r_1)=0. More...

Function/Subroutine Documentation

subroutine current	(	type (magnetic_geometry)	GEOMETRY,
		type (plasma_profiles)	PROFILES,
		type (transport_coefficients)	TRANSPORT,
		type (sources_and_sinks)	SOURCES,
		type (time_evolution)	EVOLUTION,
		type (run_control)	CONTROL,
		integer	j_boun,
		integer, intent(inout)	ifail,
		character(len=500)	failstring
	)

CURRENT TRANSPORT EQUATION.

This subroutine solves current diffusion equation and provides the flux function, density of parallel current, density of toroidal current, safety factor power density due to Ohmic heating and parallel electric field

$\sigma _\| \left(\frac{\partial }{\partial t}- \frac{\rho{\dot{B}}_0}{2B_0}\cdot\frac{\partial}{\partial\rho} \right) \Psi = \frac{F^2}{\mu_0B_0\rho } \frac{\partial}{\partial\rho} \left[\frac{V'}{4\pi^2} \left\langle\left|\frac{\nabla\rho}{R}\right|^2 \right\rangle \frac{1}{F}\frac{\partial\Psi}{\partial\rho} \right]- \frac{V'}{2\pi\rho} \left(j_\mathrm{ni,exp}+j_\mathrm{ni,imp}\cdot\Psi \right)$

These subroutines define generic numerical coefficients and boundary conditions, required by standardized interface to numerical solver

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 205 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine derivn	(	integer	N,
		real (r8), dimension(n)	X,
		real (r8), dimension(n)	Y,
		real (r8), dimension(n)	DY1
	)

These subroutines calculate first and second derivatives, DY1 and DY2, of function Y respect to argument X.

These subroutines have been extracted from RITM code and consist of derivation and integration routines

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 3856 of file main_plasma.f90.

Here is the caller graph for this function:

subroutine electron_density	(	type (magnetic_geometry)	GEOMETRY,
		type (plasma_profiles)	PROFILES,
		type (transport_coefficients)	TRANSPORT,
		type (sources_and_sinks)	SOURCES,
		type (time_evolution)	EVOLUTION,
		type (run_control)	CONTROL,
		real (r8), dimension(2)	HYPER_DIFF,
		integer, intent(inout)	ifail,
		character(len=500)	failstring
	)

Definition at line 1342 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine f_axis	(	integer	n,
		real *8, dimension(n)	r,
		real *8, dimension(n)	f
	)

This subroutine finds f(r_1=0) from f(r_2), f(r_3) and f(r_4)

Author: ??

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 3978 of file main_plasma.f90.

Here is the caller graph for this function:

subroutine f_par_axis	(	integer	n,
		real(r8), dimension(n)	r,
		real(r8), dimension(n)	f
	)

This subroutine finds f(r_1=0) from f(r_2), f(r_3) d/dr f(r_1)=0.

Author: David Coster

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 4019 of file main_plasma.f90.

Here is the caller graph for this function:

subroutine integr	(	integer	N,
		real (r8), dimension(n)	X,
		real (r8), dimension(n)	Y,
		real (r8), dimension(n)	INTY
	)

This subroutine calculates integral of function Y(X)*X from X=0 until X=X(N)

These subroutines have been extracted from RITM code and consist of derivation and integration routines

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 3916 of file main_plasma.f90.

Here is the caller graph for this function:

subroutine integr2	(	integer	N,
		real (r8), dimension(n)	X,
		real (r8), dimension(n)	Y,
		real (r8), dimension(n)	INTY
	)

Definition at line 3944 of file main_plasma.f90.

Here is the caller graph for this function:

subroutine ion_density ( )

PARTICLE TRANSPORT EQUATIONS.

This subroutine solves ion particle transport equations for ion components from 1 to NION, and provides: density and flux of ion components from 1 to NION

$\left(\frac{\partial }{\partial t}- \frac{{\dot{B}}_0}{2B_0}\cdot \frac{\partial}{\partial\rho}\rho \right) \left(V'n_i \right)+ \frac{\partial }{\partial\rho }\Gamma_i= V' \left(S_{i\mathrm{,exp}}-S_{i,\mathrm{imp}}\cdot n_i \right)$

These subroutines define generic numerical coefficients and boundary conditions, required by standardized interface to numerical solver

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 811 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine main_plasma ( )

Main plasma.

This routine finds the solution for the set of transport equations describing the main plasma, with given sources and transport coefficients for all components.

Equations are derived following Hinton&Hazeltine, Rev. Mod. Phys. vol. 48 (1976), pp.239-308

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 16 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine quasi_neutrality	(	type (magnetic_geometry)	GEOMETRY,
		type (plasma_profiles)	PROFILES,
		type (impurity_profiles)	IMPURITY,
		type (run_control)	CONTROL,
		integer, intent(inout)	ifail,
		character(len=500)	failstring
	)

QUASI NEUTRALITY.

This subroutine calculates electron density, electron flux plasma effective charge and convective contribution to electron energy transport from density and flux of background ions (all ion components computed by the ETS) and impurity ions (all ion components computed by separate impurity routine) using quasi-neutrality condition

These subroutines define generic numerical coefficients and boundary conditions, required by standardized interface to numerical solver

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 1838 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine rotation	(	type (magnetic_geometry)	GEOMETRY,
		type (plasma_profiles)	PROFILES,
		type (transport_coefficients)	TRANSPORT,
		type (sources_and_sinks)	SOURCES,
		type (time_evolution)	EVOLUTION,
		type (run_control)	CONTROL,
		integer, intent(inout)	ifail,
		character(len=500)	failstring
	)

ROTATION TRANSPORT EQUATIONS.

This subroutine solves the momentum transport equations for ion components fron 1 to NION, and provides: ion toroidal rotation velocity, ion angular velocity, ion momentum (total and individual per ion component), ion momentum flux (total and individual per ion component)

$\left( \frac{\partial}{\partial t}- \frac{{\dot{B}}_0}{2B_0}\cdot \frac{\partial }{\partial\rho }\rho \right) \left(V' \left\langle R\right\rangle m_in_iu_{i,\varphi} \right)+ \frac{\partial }{\partial\rho}\Phi_i= V' \left( U_{i,\varphi\mathrm{,exp}}- U_{i,\varphi,\mathrm{imp}}\cdot u_{i,\varphi}+ U_{\mathrm{zi},\varphi} \right)$

These subroutines define generic numerical coefficients and boundary conditions, required by standardized interface to numerical solver

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 3281 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine temperatures ( )

HEAT TRANSPORT EQUATIONS.

This subroutine solves transport equations for ion components from 1 to NION and electrons, and provides: temperatures, heat fluxes and its convective and conductive components

$\frac{3}{2} \left(\frac{\partial}{\partial t}- \frac{{\dot{B}}_0}{2B_0}\cdot \frac{\partial}{\partial\rho}\rho \right) \left(n_iT_iV'^{\frac{5}{3}} \right)+ V'^{\frac{2}{3}} \frac{\partial}{\partial\rho} \left(q_i+T_i\gamma_i \right)= V'^{\frac{5}{3}} \left[Q_{i\mathrm{,exp}}-Q_{i,\mathrm{imp}}\cdot T_i+Q_\mathrm{ei}+Q_\mathrm{zi}+Q_\mathit{\gamma i} \right]$

$\frac{3}{2} \left(\frac{\partial }{\partial t}- \frac{{\dot{B}}_0}{2B_0}\cdot \frac{\partial }{\partial\rho}\rho \right) \left(n_eT_eV'^{\frac{5}{3}}\right)+ V'^{\frac{2}{3}}\frac{\partial}{\partial\rho} \left(q_e+T_e\gamma_e \right)= V'^{\frac{5}{3}} \left[Q_{e\mathrm{,exp}}-Q_{e,{\mathrm{imp}}}\cdot T_e+Q_{\mathrm{ie}}-Q_{\mathit{\gamma i}} \right]$

These subroutines define generic numerical coefficients and boundary conditions, required by standardized interface to numerical solver

Author: D.Kalupin, R.Stankiewicz

Version: "$Id: main_plasma.f90 2093 2019-10-02 06:48:06Z g2diy $"

Definition at line 2290 of file main_plasma.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

Functions/Subroutines

Function/Subroutine Documentation