IsothermInterface.f90

Go to the documentation of this file.

module isotherminterfacemodule
  use kindmodule, only: dp, i4b
 
  implicit none
  private
  public :: isothermtype
 
  !> @brief Interface for isotherms.
  !<
  type, abstract :: isothermtype
  contains
    procedure(value_if), deferred :: value !< Evaluate isotherm value at node n
    procedure(derivative_if), deferred :: derivative !< Evaluate derivative d(value)/dc at node n
  end type isothermtype
 
  abstract interface
    !> @brief Evaluate the isotherm at a given node
    !<
    function value_if(this, c, n) result(val)
      ! -- import
      import :: isothermtype, dp, i4b
      ! -- return
      real(dp) :: val !< isotherm value
      ! -- dummy
      class(isothermtype), intent(in) :: this !< isotherm object
      real(dp), dimension(:), intent(in) :: c !< concentration array
      integer(I4B), intent(in) :: n !< node index
    end function value_if
 
    !> @brief Evaluate derivative of the isotherm at a given node
    !<
    function derivative_if(this, c, n) result(derv)
      ! -- import
      import :: isothermtype, dp, i4b
      ! -- return
      real(dp) :: derv !< derivative d(value)/dc evaluated at c
      ! -- dummy
      class(isothermtype), intent(in) :: this !< isotherm object
      real(dp), dimension(:), intent(in) :: c !< concentration array
      integer(I4B), intent(in) :: n !< node index
    end function derivative_if
  end interface
 
end module isotherminterfacemodule