Routines for user defined atmospheric chemical scheme. More...
Functions/Subroutines | |
subroutine | kinetic_4 (nr, rk, temp, xlw, press, azi, att, option_photolysis) |
kinetic_4 | |
subroutine | fexchem_4 (ns, nr, y, rk, zcsourc, convers_factor, chem) |
fexchem_4 | |
subroutine | ssh_jacdchemdc (ns, nr, y, convers_factor, convers_factor_jac, rk, jacc) |
ssh_jacdchemdc | |
subroutine | ssh_lu_decompose (ns, m) |
ssh_lu_decompose | |
subroutine | ssh_lu_solve (ns, m, x) |
ssh_lu_solve | |
subroutine | ssh_dimensions (ns, nr, nr_photolysis) |
ssh_dimensions | |
subroutine | ssh_kinetic (ns, nbin_aer, nr, iheter, icld, rk, temp, xlw, press, azi, att, lwctmp, granulo, wetdiam, dsf_aero, ispeclost, wmol, lwcmin, option_photolysis) |
ssh_kinetic | |
subroutine | ssh_fexchem (ns, nr, nemis, y, rk, zcsourc, convers_factor, chem) |
ssh_fexchem | |
subroutine | hetrxn (ns, nbin_aer, temp, press, icld, lwctmp, wetdiam, granulo, rk156, rk157, rk158, rk159, dsf_aero, ispeclost, wmol, lwcmin) |
hetrxn |
Routines for user defined atmospheric chemical scheme.
See CEREA: http://cerea.enpc.fr/polyphemus https://sshaerosol.wordpress.com/ https://github.com/sshaerosol/ssh-aerosol/
subroutine fexchem_4 | ( | 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_4
Computes the chemical production terms
[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 |
subroutine hetrxn | ( | integer | ns, |
integer | nbin_aer, | ||
double precision | temp, | ||
double precision | press, | ||
integer | icld, | ||
double precision | lwctmp, | ||
double precision, dimension(nbin_aer) | wetdiam, | ||
double precision, dimension(nbin_aer) | granulo, | ||
double precision | rk156, | ||
double precision | rk157, | ||
double precision | rk158, | ||
double precision | rk159, | ||
double precision, dimension(nbin_aer) | dsf_aero, | ||
integer, dimension(4) | ispeclost, | ||
double precision, dimension(ns) | wmol, | ||
double precision | lwcmin ) |
hetrxn
Dummy function for compatibility with SPACK
subroutine kinetic_4 | ( | 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_4
Computation of kinetic rates for atmospheric chemistry
[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 |
subroutine ssh_dimensions | ( | integer | ns, |
integer | nr, | ||
integer | nr_photolysis ) |
ssh_dimensions
Rountine provided by SPACK. Return number of species / reactions
subroutine ssh_fexchem | ( | integer | ns, |
integer | nr, | ||
integer | nemis, | ||
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 ) |
ssh_fexchem
Routine provided by SPACK. Computes the chemical production terms
[in] | nr | total number of chemical reactions |
[in] | ns | total number of chemical species |
[in] | nemis | flag to activate source terms |
[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 |
subroutine ssh_jacdchemdc | ( | 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 ) |
ssh_jacdchemdc
Routine provided by SPACK. Computes the Jacobian matrix for chemistry
[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 |
subroutine ssh_kinetic | ( | integer | ns, |
integer | nbin_aer, | ||
integer | nr, | ||
integer | iheter, | ||
integer | icld, | ||
double precision, dimension(nr) | rk, | ||
double precision | temp, | ||
double precision | xlw, | ||
double precision | press, | ||
double precision | azi, | ||
double precision | att, | ||
double precision | lwctmp, | ||
double precision, dimension(nbin_aer) | granulo, | ||
double precision, dimension(nbin_aer) | wetdiam, | ||
double precision, dimension(nbin_aer) | dsf_aero, | ||
integer, dimension(4) | ispeclost, | ||
double precision, dimension(ns) | wmol, | ||
double precision | lwcmin, | ||
integer | option_photolysis ) |
ssh_kinetic
Routine provided by SPACK. Computes kinetic rates for the gas-phase.
subroutine ssh_lu_decompose | ( | integer | ns, |
double precision, dimension(ns,ns) | m ) |
ssh_lu_decompose
Routine provided by SPACK. Computes LU factorization of matrix m
[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 |
subroutine ssh_lu_solve | ( | integer | ns, |
double precision, dimension(ns,ns) | m, | ||
double precision, dimension(ns) | x ) |
ssh_lu_solve
Resolution of MY=X where M is an LU factorization computed by lu_decompose. Routine provided by SPACK.
[in] | ns | matrix row number from the chemical species number |
[in] | m | an LU factorization computed by lu_decompose |
[in,out] | x | on entry, the right-hand side of the equation |