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 30 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 60 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
 
    if (c(n) > 0.0_dp) then
      derv = this%a(n) * this%Kf(n) * c(n)**(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 43 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
 
    if (c(n) > 0.0_dp) then
      val = this%Kf(n) * c(n)**this%a(n)
    else
      val = 0.0_dp
    end if