chem_luscheme2.f90 File Reference

Routines for atmospheric chemical scheme 2. More...

Functions/Subroutines
subroutine	kinetic_2 (nr, rk, temp, xlw, press, azi, att, option_photolysis)
	kinetic_2 More...
subroutine	fexchem_2 (ns, nr, y, rk, zcsourc, convers_factor, chem)
	fexchem_2 More...
subroutine	jacdchemdc_2 (ns, nr, y, convers_factor, convers_factor_jac, rk, jacc)
	jacdchemdc_2 More...
subroutine	rates_2 (ns, nr, rk, y, w)
	rates_2 More...
subroutine	dratedc_2 (ns, nr, rk, y, dw)
	dratedc_2 More...
subroutine	lu_decompose_2 (ns, m)
	lu_decompose_2 More...
subroutine	lu_solve_2 (ns, m, x)
	lu_solve_2 More...

Detailed Description

Routines for atmospheric chemical scheme 2.

Remarks

These routines are automatically generated by SPACK See CEREA: http://cerea.enpc.fr/polyphemus

Function/Subroutine Documentation

dratedc_2()

subroutine dratedc_2	(	integer	ns,
		integer	nr,
		double precision, dimension(nr)	rk,
		double precision, dimension(ns)	y,
		double precision, dimension(nr,ns)	dw
	)

dratedc_2

Computation of derivatives of reaction rates

Parameters

[in]	nr	total number of chemical reactions
[in]	ns	total number of chemical species
[in]	rk	kinetic rates
[in]	y	concentrations vector
[out]	dw	derivatives of reaction rates

fexchem_2()

subroutine fexchem_2	(	integer	ns,
		integer	nr,
		double precision, dimension(ns)	y,
		double precision, dimension(nr)	rk,
		double precision, dimension(ns)	zcsourc,
		double precision, dimension(ns)	convers_factor,
		double precision, dimension(ns)	chem
	)

fexchem_2

Computation of the chemical production terms

Parameters

[in]	ns	total number of chemical species
[in]	nr	total number of chemical reactions
[in]	y	concentrations vector
[in]	rk	kinetic rates
[in]	zcsourc	source term
[in]	convers_factor	conversion factors
[out]	chem	chemical production terms for every species

jacdchemdc_2()

subroutine jacdchemdc_2	(	integer	ns,
		integer	nr,
		double precision, dimension(ns)	y,
		double precision, dimension(ns)	convers_factor,
		double precision, dimension(ns,ns)	convers_factor_jac,
		double precision, dimension(nr)	rk,
		double precision, dimension(ns,ns)	jacc
	)

jacdchemdc_2

Computation of the Jacobian matrix for atmospheric chemistry

Parameters

[in]	ns	total number of chemical species
[in]	nr	total number of chemical reactions
[in]	y	concentrations vector
[in]	convers_factor	conversion factors of mug/m3 to molecules/cm3
[in]	convers_factor_jac	conversion factors for the Jacobian matrix (Wmol(i)/Wmol(j))
[in]	rk	kinetic rates
[out]	jacc	Jacobian matrix

kinetic_2()

subroutine kinetic_2	(	integer	nr,
		double precision, dimension(nr)	rk,
		double precision	temp,
		double precision	xlw,
		double precision	press,
		double precision	azi,
		double precision	att,
		integer	option_photolysis
	)

kinetic_2

Computation of kinetic rates for atmospheric chemistry

Parameters

[in]	nr	total number of chemical reactions
[in]	option_photolysis	flag to activate or not photolysis reactions
[in]	azi	solar zenith angle
[in]	att	atmospheric attenuation variable
[in]	temp	temperature
[in]	press	pressure
[in]	xlw	water massic fraction
[out]	rk(nr)	kinetic rates

lu_decompose_2()

subroutine lu_decompose_2	(	integer	ns,
		double precision, dimension(ns,ns)	m
	)

lu_decompose_2

Computation of LU factorization of matrix m

Parameters

[in]	ns	matrix row number from the chemical species number
[in,out]	m	on entry, an invertible matrix. On exit, an LU factorization of m

lu_solve_2()

subroutine lu_solve_2	(	integer	ns,
		double precision, dimension(ns,ns)	m,
		double precision, dimension(ns)	x
	)

lu_solve_2

Resolution of MY=X where M is an LU factorization computed by lu_decompose_2

Parameters

[in]	ns	matrix row number from the chemical species number
[in]	m	an LU factorization computed by lu_decompose_2
[in,out]	x	on entry, the right-hand side of the equation.

rates_2()

subroutine rates_2	(	integer	ns,
		integer	nr,
		double precision, dimension(nr)	rk,
		double precision, dimension(ns)	y,
		double precision, dimension(nr)	w
	)

rates_2

Computation of reaction rates

Parameters

[in]	ns	total number of chemical species
[in]	nr	total number of chemical reactions
[in]	rk	kinetic rates
[in]	y	concentrations vector
[out]	w	reaction rates