Data_Type_Mesh.f90

Go to the documentation of this file.

module data_type_mesh
!-----------------------------------------------------------------------------------------------------------------------------------
USE ir_precision        ! Integers and reals precision definition.
USE data_type_cell_quad ! Definition of Type_Cell_Quad.
USE data_type_node      ! Definition of Type_Node.
USE data_type_face      ! Definition of Type_Face.
!-----------------------------------------------------------------------------------------------------------------------------------

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

!-----------------------------------------------------------------------------------------------------------------------------------
type, public:: type_mesh
  integer(I4P)::                      n = 0_i4p
  real(R8P)::                         dh = 0._r8p
  real(R8P), allocatable::            x(:)
  real(R8P), allocatable::            y(:)

  integer(I4P)::                      nnode = 0_i4p
  integer(I4P)::                      nface = 0_i4p
  integer(I4P)::                      ncell = 0_i4p
  type(type_node),      allocatable:: node(:)
  type(type_face),      allocatable:: face(:)
  type(type_cell_quad), allocatable:: cell(:)
  contains
    procedure, non_overridable:: initialize ! Procedure for initializing Type_Mesh.
    procedure, non_overridable:: finalize   ! Procedure for finalizing   Type_Mesh.
    procedure, non_overridable:: pprint     ! Procedure for printing Type_Mesh components with a "pretty" format.
endtype type_mesh
!-----------------------------------------------------------------------------------------------------------------------------------
contains

  elemental subroutine initialize(self,N)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_mesh), intent(INOUT):: self   
  integer(I4P),     intent(IN)::    n
  integer(I4P)::                    nnode
  integer(I4P)::                    nface
  integer(I4P)::                    ncell
  integer(I4P)::                    ij
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  self%N = n
  self%dh = 1._r8p/real(n,r8p)
  if (allocated(self%x)) deallocate(self%x) ; allocate(self%x(0:n))
  if (allocated(self%y)) deallocate(self%y) ; allocate(self%y(0:n))
  do ij=0,n
    self%x(ij) = real(ij,r8p)/real(n,r8p)
  enddo
  self%y = self%x

  nnode = (n+1)*(n+1)
  nface = 2*n*(n+1)
  ncell = n*n
  if (allocated(self%node)) deallocate(self%node) ; allocate(self%node(1:nnode))
  if (allocated(self%face)) deallocate(self%face) ; allocate(self%face(1:nface))
  if (allocated(self%cell)) deallocate(self%cell) ; allocate(self%cell(1:ncell))
!  ij = 0
  !do j=0,N
    !do i=0,N
      !ij = ij + 1 ; call xyz%set(x=real(i,R8P)/real(N,R8P),y=real(j,R8P)/real(N,R8P),z=0._R8P) ; call self%node(ig)%initialize(xyz)
    !enddo
  !enddo
  !ij = 0
  !do j=1,N
    !do i=1,N
      !ij = ij + 1 ; call self%face(ij)%initialize(node=[self%node(ij)],N=null(),S=null())
    !enddo
!  enddo
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine initialize

  elemental subroutine finalize(self)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_mesh), intent(INOUT):: self 
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  self%N = 0_i4p
  self%dh = 0._r8p
  if (allocated(self%x)) deallocate(self%x)
  if (allocated(self%y)) deallocate(self%y)
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine finalize

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

  !---------------------------------------------------------------------------------------------------------------------------------
  write(unit,'(A)',iostat=err)' N='//str(n=self%N)
  if (allocated(self%x).AND.allocated(self%y)) then
    do ij=lbound(self%x,dim=1),ubound(self%x,dim=1)
      write(unit,'(A)',iostat=err)' ij='//str(n=ij)//' x='//str(n=self%x(ij))//' y='//str(n=self%y(ij))
    enddo
  endif
  write(unit,'(A)',iostat=err)' dh='//str(n=self%dh)
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endfunction pprint
endmodule data_type_mesh