chem_luscheme_siream.f90 File Reference

Routines for atmospheric chemical and aerosols scheme 3. More...

Functions/Subroutines
subroutine	kinetic_siream (nr, rk, temp, xlw, press, azi, att, option_photolysis)
	kinetic_siream More...
subroutine	fexchem_siream (ns, nr, y, rk, zcsourc, convers_factor, chem)
	fexchem_siream More...
subroutine	jacdchemdc_siream (ns, nr, y, convers_factor, convers_factor_jac, rk, jacc)
	jacdchemdc_siream More...
subroutine	rates_siream (ns, nr, rk, y, w)
	rates_siream More...
subroutine	dratedc_siream (ns, nr, rk, y, dw)
	dratedc_siream More...
subroutine	lu_decompose_siream (ns, m)
	lu_decompose_siream More...
subroutine	lu_solve_siream (ns, m, x)
	lu_solve_siream More...

Detailed Description

Routines for atmospheric chemical and aerosols scheme 3.

Remarks

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

Function/Subroutine Documentation

dratedc_siream()

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

dratedc_siream

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

subroutine fexchem_siream	(	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_siream

Computation of the chemical production terms

Parameters

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

subroutine jacdchemdc_siream	(	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_siream

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

subroutine kinetic_siream	(	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_siream

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

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

lu_decompose_siream

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

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

lu_solve_siream

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.

rates_siream()

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

rates_siream

Computation 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]	w	reaction rates