rate_coeff_neutrals_m Module

Module for calculating fits for rate coefficients k_iz, k_rec and cooling rate coefficients w_iz, p_rec

Uses

- variable_m
- error_handling_grillix_m
- screen_io_m
- precision_grillix_m
- status_codes_grillix_m

Contents

Derived Types

rate_coeff_neutrals_t

Functions

eval_rate_coeff poly8x8_eval poly8_eval

Subroutines

create_rate_coeff

Derived Types

type, public :: rate_coeff_neutrals_t

Class for rate coefficients (ionization & recombination)

Type-Bound Procedures

procedure , public :: create => create_rate_coeff Subroutine
procedure , public :: eval => eval_rate_coeff Function

Functions

public function eval_rate_coeff(self, te_in, ne_in)

Evaluate the normalized rate coefficient from a double polynomial fit

Arguments

Type	Intent	Optional	Attributes		Name
class(rate_coeff_neutrals_t),	intent(in)			::	self	Instance of class
real(kind=GP),	intent(in)			::	te_in	Electron temperature
real(kind=GP),	intent(in)			::	ne_in	Plasma density

Return Value real(kind=gp)

public function poly8x8_eval(pc, x, y)

Returns double polynomial of degree 8x8, with coefficients pc(x,y), evaluated at x and y.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=GP),	intent(in),		dimension(0:8,0:8)	::	pc	polynomial coefficients
real(kind=GP),	intent(in)			::	x	input variables
real(kind=GP),	intent(in)			::	y	input variables

Return Value real(kind=gp)

public function poly8_eval(pc, x)

Returns polynomial of degree 8, with coefficients pc(0:8), evaluated at x.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=GP),	intent(in),		dimension(0:8)	::	pc	polynomial coefficients
real(kind=GP),	intent(in)			::	x	input variable

Return Value real(kind=gp)

Subroutines

public subroutine create_rate_coeff(self, file_name, ne_ref, te_ref, k_ref)

Create rate coefficient object with fit coefficients taken from data file

Arguments

Type	Intent	Optional	Attributes		Name
class(rate_coeff_neutrals_t),	intent(inout)			::	self	Instance of class
character(len=*),	intent(in)			::	file_name	Name of atomic data file to read coefficients from
real(kind=GP),	intent(in)			::	ne_ref	Reference for electron density
real(kind=GP),	intent(in)			::	te_ref	Reference for electron temperature
real(kind=GP),	intent(in)			::	k_ref	Reference for rate coefficient