freundlichisothermmodule Module Reference

Data Types
interface	freundlichisothermtype
	Freundlich isotherm implementation of IsothermType. More...

Functions/Subroutines
type(freundlichisothermtype) function	constructor (Kf, a)
	Constructor for Freundlich 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(freundlichisothermtype) function freundlichisothermmodule::constructor ( real(dp), dimension(:), intent(in), pointer  Kf, real(dp), dimension(:), intent(in), pointer  a  )

private

Definition at line 39 of file FreundlichIsotherm.f90.

    type(FreundlichIsothermType) :: isotherm
    ! -- dummy
    real(DP), pointer, dimension(:), intent(in) :: Kf
    real(DP), pointer, dimension(:), intent(in) :: a
    ! -- local
    isotherm%Kf => kf
    isotherm%a => a
 

Here is the caller graph for this function:

derivative()

real(dp) function freundlichisothermmodule::derivative ( class(freundlichisothermtype), 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 73 of file FreundlichIsotherm.f90.

    ! -- return
    real(DP) :: derv !< derivative d(value)/dc evaluated at c
    ! -- dummy
    class(FreundlichIsothermType), intent(in) :: this
    real(DP), dimension(:), intent(in) :: c !< concentration array
    integer(I4B), intent(in) :: n !< node index
    real(DP), parameter :: K = 10.0_dp !< constant to limit derivative at c=0
    real(DP) :: eps !< small concentration offset
 
    eps = (k / (this%a(n) * this%Kf(n)))**(1.0_dp / (this%a(n) - 1.0_dp))
 
    if (c(n) > dzero) then
      derv = this%a(n) * this%Kf(n) * ((c(n) + eps)**(this%a(n) - 1.0_dp))
    else
      derv = 0.0_dp
    end if

value()

real(dp) function freundlichisothermmodule::value ( class(freundlichisothermtype), 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 52 of file FreundlichIsotherm.f90.

    ! -- return
    real(DP) :: val !< isotherm value
    ! -- dummy
    class(FreundlichIsothermType), intent(in) :: this
    real(DP), dimension(:), intent(in) :: c !< concentration array
    integer(I4B), intent(in) :: n !< node index
    real(DP), parameter :: K = 10.0_dp !< constant to limit derivative at c=0
    real(DP) :: eps !< small concentration offset
 
    eps = (k / (this%a(n) * this%Kf(n)))**(1.0_dp / (this%a(n) - 1.0_dp))
 
    if (c(n) > dzero) then
      val = this%Kf(n) * (c(n) + eps)**this%a(n) - this%Kf(n) * eps**this%a(n)
    else
      val = 0.0_dp
    end if