scale_localmesh_2d.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!> module FElib / Mesh / Local 2D
!!
!! @par Description
!!      Module to manage 2D local mesh for element-based methods
!!
!! @author Yuta Kawai, Team SCALE
!<
#include "scaleFElib.h"
module scale_localmesh_2d
 
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_localmesh_base, only: &
    localmeshbase, localmeshbase_init, localmeshbase_final
  use scale_element_base, only: elementbase, elementbase2d
 
  !-----------------------------------------------------------------------------
  implicit none
  private
 
  !-----------------------------------------------------------------------------
  !
  !++ Public type & procedure
  ! 
  !> Derived type representing a local mesh for 2D domain
  type, extends(localmeshbase), public :: localmesh2d
    type(elementbase2d), pointer :: refelem2d !< Pointer to the reference element for 2D
    
    real(rp) :: xmin !< Minimum x-coordinate of the local mesh
    real(rp) :: xmax !< Maximum x-coordinate of the local mesh
    real(rp) :: ymin !< Minimum y-coordinate of the local mesh
    real(rp) :: ymax !< Maximum y-coordinate of the local mesh
 
    integer :: nex !< Number of elements in the x-direction
    integer :: ney !< Number of elements in the y-direction
 
    real(rp), allocatable :: lon(:,:)  !< Array to save longitude coordinates of the local mesh nodes
    real(rp), allocatable :: lat(:,:)  !< Array to save latitude coordinates of the local mesh nodes
  end type localmesh2d
 
  public :: localmesh2d_init, localmesh2d_final
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
 
contains
 
  !> Initialize an object to manage a local mesh for 2D domain
!OCL SERIAL
  subroutine localmesh2d_init( this, &
    lcdomID, refElem, myrank )
    
    implicit none
 
    class(localmesh2d), intent(inout) :: this
    integer, intent(in) :: lcdomid
    class(elementbase2d), intent(in), target :: refelem
    integer, intent(in), optional :: myrank
    !-------------------------------------------------
 
    this%refElem2D  => refelem
    call localmeshbase_init(this, lcdomid, refelem, 2, myrank)
    !$acc enter data copyin(this)
    return
  end subroutine localmesh2d_init
 
  !> Finalize an object to manage a local mesh for 2D domain
!OCL SERIAL
  subroutine localmesh2d_final( this, is_generated )
    implicit none
    type(localmesh2d), intent(inout) :: this
    logical, intent(in) :: is_generated
    !-------------------------------------------------
 
    if (is_generated) then
      !$acc exit data delete(this%G_ij, this%GIJ)
      deallocate( this%G_ij, this%GIJ )
      !$acc exit data delete(this%lon, this%lat)
      deallocate( this%lon, this%lat )
    end if
    call localmeshbase_final( this, is_generated )
    !$acc exit data delete(this)
    return
  end subroutine localmesh2d_final
  
end module scale_localmesh_2d