chem_luscheme1.f90 File Reference

Routines for atmospheric chemical scheme 1. More...

Functions/Subroutines
subroutine	kinetic_1 (nr, rk, temp, xlw, press, azi, att, option_photolysis)
	kinetic_1 More...
subroutine	fexchem_1 (ns, nr, y, rk, zcsourc, convers_factor, chem)
	fexchem_1 More...
subroutine	jacdchemdc_1 (ns, nr, y, convers_factor, convers_factor_jac, rk, jacc)
	jacdchemdc_1 More...
subroutine	rates_1 (ns, nr, rk, y, w)
	rates_1 More...
subroutine	dratedc_1 (ns, nr, rk, y, dw)
	dratedc_1 More...
subroutine	lu_decompose_1 (ns, m)
	lu_decompose_1 More...
subroutine	lu_solve_1 (ns, m, x)
	lu_solve_1 More...

Detailed Description

Routines for atmospheric chemical scheme 1.

Remarks

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

Function/Subroutine Documentation

dratedc_1()

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

dratedc_1

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

subroutine fexchem_1	(	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_1

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

subroutine jacdchemdc_1	(	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_1

Computation of the Jacobian matrix for atmospheric chemistry

Parameters

[in]	nr	total number of chemical reactions
[in]	ns	total number of chemical species
[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_1()

subroutine kinetic_1	(	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_1

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

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

lu_decompose_1

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

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

lu_solve_1

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

Parameters

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

rates_1()

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

rates_1

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