scale_localmesh_1d.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!> module FElib / Mesh / Local 1D
!!
!! @par Description
!!      Module to manage 1D local mesh for element-based methods
!!
!! @author Yuta Kawai, Team SCALE
!<
#include "scaleFElib.h"
module scale_localmesh_1d
 
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_localmesh_base, only: &
    localmeshbase, localmeshbase_init, localmeshbase_final
  use scale_element_base, only: elementbase, elementbase1d
 
  !-----------------------------------------------------------------------------
  implicit none
  private
 
  !-----------------------------------------------------------------------------
  !
  !++ Public type & procedure
  ! 
  !> Derived type representing a local mesh for 1D domain
  type, extends(localmeshbase), public :: localmesh1d
    type(elementbase1d), pointer :: refelem1d !< Pointer to the reference element for 1D
    real(rp) :: xmin !< Minimum x-coordinate of the local mesh
    real(rp) :: xmax !< Maximum x-coordinate of the local mesh
  end type localmesh1d
 
  public :: localmesh1d_init, localmesh1d_final
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
 
contains
  !> Initialize an object to manage a local mesh for 1D domain
!OCL SERIAL
  subroutine localmesh1d_init( this, &
    lcdomID, refElem, myrank )
    
    implicit none
 
    type(localmesh1d), intent(inout) :: this
    integer, intent(in) :: lcdomid
    class(elementbase1d), intent(in), target :: refelem
    integer, intent(in), optional :: myrank
    !-------------------------------------------------
 
    this%refElem1D => refelem
    call localmeshbase_init(this, lcdomid, refelem, 1, myrank)
    !$acc enter data copyin(this)
    return
  end subroutine localmesh1d_init
 
  !> Finalize an object to manage a local mesh for 1D domain
!OCL SERIAL
  subroutine localmesh1d_final( this, is_generated )
    implicit none
 
    type(localmesh1d), intent(inout) :: this
    logical, intent(in) :: is_generated
    !-------------------------------------------------
 
    call localmeshbase_final( this, is_generated )
    !$acc exit data delete(this)
    return
  end subroutine localmesh1d_final
  
end module scale_localmesh_1d