Data_Type_Cavity.f90

Go to the documentation of this file.

module data_type_cavity
!-----------------------------------------------------------------------------------------------------------------------------------
USE ir_precision           ! Integers and reals precision definition.
USE data_type_conservative ! Definition of Type_Conservative.
USE data_type_mesh         ! Definition of Type_Mesh.
!-----------------------------------------------------------------------------------------------------------------------------------

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

!-----------------------------------------------------------------------------------------------------------------------------------
type, public:: type_cavity
  real(R8P)::                            re   = 500._r8p
  real(R8P)::                            beta = 0.6_r8p
  type(type_mesh)::                      mesh
  type(type_conservative), allocatable:: cons(:,:)
  contains
    procedure, non_overridable:: initialize ! Procedure for initializing Type_Cavity.
    procedure, non_overridable:: finalize   ! Procedure for finalizing   Type_Cavity.
    procedure, non_overridable:: pprint     ! Procedure for printing Type_Cavity components with a "pretty" format.
endtype type_cavity
!-----------------------------------------------------------------------------------------------------------------------------------
contains

  elemental subroutine initialize(self,Re,beta,N)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_cavity), intent(INOUT):: self 
  real(R8P),          intent(IN)::    re
  real(R8P),          intent(IN)::    beta
  integer(I4P),       intent(IN)::    n
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  self%Re = re
  self%beta = beta
  call self%mesh%initialize(n=n)
  if (allocated(self%cons)) deallocate(self%cons) ; allocate(self%cons(0:n,0:n)) ! call self%cons%set(s=,v=)
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine initialize

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

  !---------------------------------------------------------------------------------------------------------------------------------
  self%Re = 1._r8p
  self%beta = 0.6_r8p
  call self%mesh%finalize
  if (allocated(self%cons)) deallocate(self%cons)
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine finalize

  function pprint(self,unit) result(err)
  !---------------------------------------------------------------------------------------------------------------------------------
  implicit none
  class(type_cavity), intent(IN):: self 
  integer(I4P),       intent(IN):: unit
  integer(I4P)::                   err
  integer(I4P)::                   i,j
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  write(unit,'(A)',iostat=err)' Re='//str(n=self%Re)
  write(unit,'(A)',iostat=err)' beta='//str(n=self%beta)
  write(unit,'(A)',iostat=err)' Mesh data'
  err = self%mesh%pprint(unit)
  write(unit,'(A)',iostat=err)' Conservative variables data'
  do j=0,self%mesh%N
    do i=0,self%mesh%N
      err = self%cons(i,j)%pprint(unit)
    enddo
  enddo
  return
  !---------------------------------------------------------------------------------------------------------------------------------
  endfunction pprint
endmodule data_type_cavity