maxwell – FRTС project

Все что относится к распределению Максвелла

Uses

- constants
- kind_module

Variables

Type	Visibility	Attributes		Name		Initial
integer,	public,	parameter	::	i0	=	1002
real(kind=wp),	public		::	v_grid(i0,100)			сетка обычных скоростей
real(kind=wp),	public		::	vij(i0,100)
real(kind=wp),	public		::	fij0(i0,100,2)
real(kind=wp),	public		::	fij(i0,100,2)
real(kind=wp),	public		::	dfij(i0,100,2)
real(kind=wp),	public		::	dij(i0,100,2)
logical,	public		::	flag_d0			бывший d0
integer,	public		::	jindex			common/dddql/ d0,jindex,kindex
integer,	public		::	kindex			common/dddql/ d0,jindex,kindex

Functions

public function currlhcd(v, f) result(curs)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	v(:)
real(kind=wp),	intent(in)			::	f(:)

Return Value real(kind=wp)

public function create_vt_grid(vclt) result(vt_grid)

создание сетки тепловых скоростей

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	vclt

Return Value real(kind=wp), (i0)

public function funmaxwell(v, dfunmaxwell)

распределение Максвелла $f(v) = \frac{1}{\sqrt{2\pi}} \exp(-\frac{1}{2} v^2 ))$ и его производная $dfmaxw = - v \cdot f(v)$

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp)				::	v
real(kind=wp)				::	dfunmaxwell

Return Value real(kind=wp)

public function fmaxw_classic(v, alfa2, dfmaxw)

распределение Максвелла с альфа-частицами $f(v) = \frac{1}{\sqrt{2\pi}} \exp(-\frac{1}{2} v^2 (1.0 + \frac{1}{2} \cdot alfa_2 \cdot v^2))$ и его производная $dfmaxw = - v \cdot (1.0 + alfa_2 \cdot v^2) \cdot f(v)$

Arguments

Type		Name
real(kind=wp)	::	v
real(kind=wp)	::	alfa2
real(kind=wp)	::	dfmaxw

Return Value real(kind=wp)

public function fmaxw_ext(v, alfa2, dfmaxw)

$alfa = \sqrt{alfa_2}$ $api = 2 \cdot alfa \cdot \exp({-\frac{1}{4 alfa_2}})$ $b = 2 - erf(0.5/alfa) + api$ $f = psiq(v, alfa_2)$ $fmaxw = \frac{f+api}{b \sqrt{2\pi}}$

Arguments

Type		Name
real(kind=wp)	::	v
real(kind=wp)	::	alfa2
real(kind=wp)	::	dfmaxw

Return Value real(kind=wp)

public function fmaxw(v, alfa2, dfmaxw)

Arguments

Type		Name
real(kind=wp)	::	v
real(kind=wp)	::	alfa2
real(kind=wp)	::	dfmaxw

Return Value real(kind=wp)

public function psiq(v, alfa2)

$psiq=exp(ksiV**2)*erfcc(ksiV)*exp(-0.25/alfa2)$

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp)				::	v
real(kind=wp)				::	alfa2

Return Value real(kind=wp)

public function erfcc(x)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp)				::	x

Return Value real(kind=wp)

Subroutines

public subroutine init_vi(vclt, vi)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	vclt
real(kind=wp),	intent(out)			::	vi(i0)

public subroutine init_fmaxw_classic(vclt, enorm, fi, dfi)

Arguments

Type	Intent		Name
real(kind=wp),	intent(in)	::	vclt
real(kind=wp),	intent(in)	::	enorm
real(kind=wp),	intent(out)	::	fi(i0)
real(kind=wp),	intent(out)	::	dfi(i0)

public subroutine init_fmaxw_ext(vclt, enorm, fi, dfi)

Arguments

Type	Intent		Name
real(kind=wp),	intent(in)	::	vclt
real(kind=wp),	intent(in)	::	enorm
real(kind=wp),	intent(out)	::	fi(i0)
real(kind=wp),	intent(out)	::	dfi(i0)

maxwell Module

Contents

Variables

Functions

Subroutines

Uses

Variables

Functions

public function currlhcd(v, f) result(curs)

Arguments

Return Value real(kind=wp)

public function create_vt_grid(vclt) result(vt_grid)

Arguments

Return Value real(kind=wp), (i0)

public function funmaxwell(v, dfunmaxwell)

Arguments

Return Value real(kind=wp)

public function fmaxw_classic(v, alfa2, dfmaxw)

Arguments

Return Value real(kind=wp)

public function fmaxw_ext(v, alfa2, dfmaxw)

Arguments

Return Value real(kind=wp)

public function fmaxw(v, alfa2, dfmaxw)

Arguments

Return Value real(kind=wp)

public function psiq(v, alfa2)

Arguments

Return Value real(kind=wp)

public function erfcc(x)

Arguments

Return Value real(kind=wp)

Subroutines

public subroutine init_vi(vclt, vi)

Arguments

public subroutine init_fmaxw_classic(vclt, enorm, fi, dfi)

Arguments

public subroutine init_fmaxw_ext(vclt, enorm, fi, dfi)

Arguments