Data_Type_Face.f90

Go to the documentation of this file.

module data_type_face
!-----------------------------------------------------------------------------------------------------------------------------------
USE ir_precision     ! Integers and reals precision definition.
USE data_type_node   ! Definition of Type_Node.
USE data_type_vector ! Definition of Type_Vector.
!-----------------------------------------------------------------------------------------------------------------------------------

!-----------------------------------------------------------------------------------------------------------------------------------
private
!-----------------------------------------------------------------------------------------------------------------------------------

!-----------------------------------------------------------------------------------------------------------------------------------
type, public:: type_face
  type(type_node_ptr)::        node(1:2)
  type(type_vector), pointer:: n
  real(R8P),         pointer:: s => null()
  contains
    procedure, non_overridable:: initialize ! Procedure for initializing Type_Face.
    procedure, non_overridable:: finalize   ! Procedure for finalizing   Type_Face.
    procedure, non_overridable:: pprint     ! Procedure for printing Type_Face components with a "pretty" format.
endtype type_face
type, public:: type_face_ptr
  type(type_face), pointer:: p => null()
endtype type_face_ptr
!-----------------------------------------------------------------------------------------------------------------------------------
contains

  subroutine initialize(self,node1,node2,N,S)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_face),           intent(INOUT):: self  
  type(type_node),   pointer, intent(IN)::    node1
  type(type_node),   pointer, intent(IN)::    node2
  type(type_vector), pointer, intent(IN)::    n
  real(R8P),         pointer, intent(IN)::    s
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  self%node(1)%p => node1
  self%node(2)%p => node2
  self%N         => n
  self%S         => s
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine initialize

  subroutine finalize(self)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_face), intent(INOUT):: self 
  integer(I4P)::                    n
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  do n=lbound(self%node,dim=1),ubound(self%node,dim=1)
    if (associated(self%node(n)%p)) deallocate(self%node(n)%p) ; self%node(n)%p => null()
  enddo
  if (associated(self%N)) deallocate(self%N) ; self%N => null()
  if (associated(self%S)) deallocate(self%S) ; self%S => null()
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine finalize

  function pprint(self,unit) result(err)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_face), intent(IN):: self 
  integer(I4P),     intent(IN):: unit
  integer(I4P)::                 err
  integer(I4P)::                 n
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  write(unit,'(A)',iostat=err)' Nodes data:'
  do n=lbound(self%node,dim=1),ubound(self%node,dim=1)
    if (associated(self%node(n)%p)) then
      write(unit,'(A)',iostat=err)' n='//trim(str(n=n))
      err = self%node(n)%p%pprint(unit=unit)
    endif
  enddo
  if (associated(self%N)) then
    write(unit,'(A)',iostat=err)' Face outward unit normal:'
    write(unit,'(A)',iostat=err)' x='//trim(str(n=self%N%x))//' y='//trim(str(n=self%N%y))//' z='//trim(str(n=self%N%z))
  endif
  write(unit,'(A)',iostat=err)' Face area='//trim(str(n=self%S))
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endfunction pprint
endmodule data_type_face