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
  ! 
  type, extends(localmeshbase), public :: localmesh2d
 
    type(elementbase2d), pointer :: refelem2d
    real(rp) :: xmin, xmax
    real(rp) :: ymin, ymax
    integer :: nex, ney
 
    real(rp), allocatable :: lon(:,:)     
    real(rp), allocatable :: lat(:,:)     
  end type localmesh2d
 
  public :: localmesh2d_init, localmesh2d_final
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
 
contains
!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)
 
    return
  end subroutine localmesh2d_init
 
!OCL SERIAL
  subroutine localmesh2d_final( this, is_generated )
    implicit none
    type(localmesh2d), intent(inout) :: this
    logical, intent(in) :: is_generated
    !-------------------------------------------------
 
    if (is_generated) then
      deallocate( this%lon, this%lat )
    end if
    call localmeshbase_final( this, is_generated )
 
    return
  end subroutine localmesh2d_final
  
end module scale_localmesh_2d