chem_luscheme3.f90 File Reference

Routines for atmospheric chemical scheme 3. More...

Functions/Subroutines
subroutine	kinetic_3 (nr, rk, temp, xlw, press, azi, att, option_photolysis)
	kinetic_3 More...
subroutine	fexchem_3 (ns, nr, y, rk, zcsourc, convers_factor, chem)
	fexchem_3 More...
subroutine	jacdchemdc_3 (ns, nr, y, convers_factor, convers_factor_jac, rk, jacc)
	jacdchemdc_3 More...
subroutine	rates_3 (ns, nr, rk, y, w)
	rates_3 More...
subroutine	dratedc_3 (ns, nr, rk, y, dw)
	dratedc_3 More...
subroutine	lu_decompose_3 (ns, m)
	lu_decompose_3 More...
subroutine	lu_solve_3 (ns, m, x)
	lu_solve_3 More...

Detailed Description

Routines for atmospheric chemical scheme 3.

Remarks

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

Function/Subroutine Documentation

dratedc_3()

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

dratedc_3

Computation of derivatives 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]	dw	derivatives of reaction rates

fexchem_3()

subroutine fexchem_3	(	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_3

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_3()

subroutine jacdchemdc_3	(	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_3

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_3()

subroutine kinetic_3	(	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_3

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_3()

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

lu_decompose_3

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_3()

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

lu_solve_3

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

Parameters

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

rates_3()

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

rates_3

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