langmuirisothermmodule Module Reference

Data Types
interface	langmuirisothermtype
	Langmuir isotherm implementation of IsothermType. More...

Functions/Subroutines
type(langmuirisothermtype) function	constructor (Kl, Sbar)
	Constructor for Langmuir isotherm. More...
real(dp) function	value (this, c, n)
	Evaluate the isotherm at a given node. More...
real(dp) function	derivative (this, c, n)
	Evaluate derivative of the isotherm at a given node. More...

Function/Subroutine Documentation

constructor()

type(langmuirisothermtype) function langmuirisothermmodule::constructor ( real(dp), dimension(:), intent(in), pointer  Kl, real(dp), dimension(:), intent(in), pointer  Sbar  )

private

Parameters

[in]	kl	Langmuir constant
[in]	sbar	Total concentration of sorption sites

Definition at line 29 of file LangmuirIsotherm.f90.

    ! -- return
    type(LangmuirIsothermType) :: isotherm
    ! -- dummy
    real(DP), pointer, dimension(:), intent(in) :: Kl !< Langmuir constant
    real(DP), pointer, dimension(:), intent(in) :: Sbar !< Total concentration of sorption sites
    ! -- local
 
    isotherm%Kl => kl
    isotherm%Sbar => sbar
 

Here is the caller graph for this function:

derivative()

real(dp) function langmuirisothermmodule::derivative ( class(langmuirisothermtype), intent(in)  this, real(dp), dimension(:), intent(in)  c, integer(i4b), intent(in)  n  )

private

Returns

derivative d(value)/dc evaluated at c

Parameters

[in]	c	concentration array
[in]	n	node index

Definition at line 64 of file LangmuirIsotherm.f90.

    ! -- return
    real(DP) :: derv !< derivative d(value)/dc evaluated at c
    ! -- dummy
    class(LangmuirIsothermType), intent(in) :: this
    real(DP), dimension(:), intent(in) :: c !< concentration array
    integer(I4B), intent(in) :: n !< node index
    ! -- local
    real(DP) :: denom
 
    if (c(n) > 0.0_dp) then
      denom = (1.0_dp + this%Kl(n) * c(n))**2.0_dp
      derv = (this%Sbar(n) * this%Kl(n)) / denom
    else
      derv = 0.0_dp
    end if

value()

real(dp) function langmuirisothermmodule::value ( class(langmuirisothermtype), intent(in)  this, real(dp), dimension(:), intent(in)  c, integer(i4b), intent(in)  n  )

private

Returns

isotherm value

Parameters

[in]	c	concentration array
[in]	n	node index

Definition at line 44 of file LangmuirIsotherm.f90.

    ! -- return
    real(DP) :: val !< isotherm value
    ! -- dummy
    class(LangmuirIsothermType), intent(in) :: this
    real(DP), dimension(:), intent(in) :: c !< concentration array
    integer(I4B), intent(in) :: n !< node index
    ! -- local
    real(DP) :: denom
 
    if (c(n) > 0.0_dp) then
      denom = 1.0_dp + this%Kl(n) * c(n)
      val = (this%Sbar(n) * this%Kl(n) * c(n)) / denom
    else
      val = 0.0_dp
    end if