scale_model_mesh_manager.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!> FElib / model framework / mesh manager
!!
!! @par Description
!!          A module for managing mesh used in models
!!
!! @author Yuta Kawai, Team SCALE
!!
!<
!-------------------------------------------------------------------------------
#include "scaleFElib.h"
module scale_model_mesh_manager
  !-----------------------------------------------------------------------------
  !
  !++ Used modules
  !
  use scale_precision
  use scale_io
  use scale_prc, only: &
    prc_abort
 
  use scale_mesh_base, only: meshbase
  use scale_mesh_base1d, only: meshbase1d
  use scale_mesh_base2d, only: meshbase2d
  use scale_mesh_base3d, only: meshbase3d
  use scale_meshfield_base, only: meshfield3d
 
  use scale_sparsemat, only: sparsemat
 
  use scale_model_meshbase_manager, only: &
    modelmeshbase, &
    modelmeshbase1d, modelmeshbase2d, modelmeshbase3d
 
  use scale_model_var_manager, only: modelvarmanager
  
  !-----------------------------------------------------------------------------
  implicit none
  private
 
  !-----------------------------------------------------------------------------
  !
  !++ Public type & procedures
  !
  
  type, extends(modelmeshbase1d), public :: modelmesh1d
  contains
    procedure, public :: modelmesh1d_init
    procedure, public :: modelmesh1d_final
  end type modelmesh1d
 
  type, extends(modelmeshbase2d), public :: modelmesh2d
  contains
    procedure, public :: modelmesh2d_init
    procedure, public :: modelmesh2d_final
  end type modelmesh2d
 
  type, abstract, extends(modelmeshbase3d), public :: modelmesh3d
  contains
    procedure, public :: modelmesh3d_init
    procedure, public :: modelmesh3d_final
    procedure(modelmesh3d_create_communicator), public, deferred :: create_communicator
  end type modelmesh3d
 
  interface
    subroutine modelmesh3d_create_communicator( this, sfield_num, hvfield_num, htensorfield_num, &
      var_manager, field_list, commid, &
      field_list_is, field_list_ie     )
      import modelmesh3d
      import meshbase3d
      import modelvarmanager
      import meshfield3d
      class(modelmesh3d), target, intent(inout) :: this
      integer, intent(in) :: sfield_num
      integer, intent(in) :: hvfield_num
      integer, intent(in) :: htensorfield_num
      class(modelvarmanager), intent(inout) :: var_manager
      class(meshfield3d), intent(in) :: field_list(:)
      integer, intent(out) :: commid
      integer, intent(in), optional :: field_list_is
      integer, intent(in), optional :: field_list_ie
    end subroutine modelmesh3d_create_communicator
  end interface
 
  ! Cascade
  public :: modelmeshbase
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedures & variables
  !
  !------------------
 
contains
 
  !* 1D *************************************************************
 
!OCL SERIAL
  subroutine modelmesh1d_init( this, mesh )
    implicit none
    class(modelmesh1d), target, intent(inout) :: this
    class(meshbase1d), target, intent(in) :: mesh
    !-----------------------------------------------------
 
    call this%ModelMeshBase1D_Init(mesh)
 
    return
  end subroutine modelmesh1d_init
 
!OCL SERIAL
  subroutine modelmesh1d_final( this )
    implicit none
    class(modelmesh1d), target, intent(inout) :: this
 
    !-----------------------------------------------------
 
    call this%ModelMeshBase1D_Final()
 
    return
  end subroutine modelmesh1d_final
 
  !* 2D *************************************************************
 
!OCL SERIAL
  subroutine modelmesh2d_init( this, mesh )
    implicit none
    class(modelmesh2d), target, intent(inout) :: this
    class(meshbase2d), target, intent(in) :: mesh
    !-----------------------------------------------------
 
    call this%ModelMeshBase2D_Init(mesh)
 
    return
  end subroutine modelmesh2d_init
 
!OCL SERIAL
  subroutine modelmesh2d_final( this )
    implicit none
    class(modelmesh2d), target, intent(inout) :: this  
    !-----------------------------------------------------
 
    call this%ModelMeshBase2D_Final()
 
    return
  end subroutine modelmesh2d_final
 
  !* 3D *************************************************************
 
!OCL SERIAL
  subroutine modelmesh3d_init( this, mesh )
    implicit none
    class(modelmesh3d), target, intent(inout) :: this
    class(meshbase3d), target, intent(in) :: mesh
    !-----------------------------------------------------
 
    this%ptr_mesh => mesh
    call this%ModelMeshBase3D_Init(mesh)
 
    return
  end subroutine modelmesh3d_init  
 
!OCL SERIAL
  subroutine modelmesh3d_final( this )
    implicit none
    class(modelmesh3d), target, intent(inout) :: this
    !-----------------------------------------------------
 
    nullify( this%ptr_mesh )
    call this%ModelMeshBase3D_Final()
 
    return
  end subroutine modelmesh3d_final
 
end module scale_model_mesh_manager