scale_mesh_base2d Module Reference

module FElib / Mesh / Base 2D More...

Data Types
type	meshbase2d
	Derived type to manage a computational mesh (base type for 2D domain) More...
interface	meshbase2d_generate

Functions/Subroutines
subroutine, public	meshbase2d_init (this, refelem, nlocalmeshperprc, nprocs, myrank)
	Initialize an object to manage a 2D computational mesh.
subroutine, public	meshbase2d_final (this)
	Finalize an object to manage a 2D computational mesh.
subroutine, public	meshbase2d_setgeometricinfo (lcmesh, coord_conv, calc_normal)

Variables
integer, public	meshbase2d_dimtype_num = 4
integer, public	meshbase2d_dimtypeid_x = 1
integer, public	meshbase2d_dimtypeid_y = 2
integer, public	meshbase2d_dimtypeid_xy = 3
integer, public	meshbase2d_dimtypeid_xyt = 4

Detailed Description

module FElib / Mesh / Base 2D

Description

Base module to manage 2D meshes for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshbase2d_init()

subroutine, public scale_mesh_base2d::meshbase2d_init	(	class(meshbase2d), intent(inout)	this,
		class(elementbase2d), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in), optional	nprocs,
		integer, intent(in), optional	myrank )

Initialize an object to manage a 2D computational mesh.

Definition at line 76 of file scale_mesh_base2d.F90.

    
    implicit none
    class(MeshBase2D), intent(inout) :: this
    class(ElementBase2D), intent(in), target :: refElem
    integer, intent(in) :: NLocalMeshPerPrc
    integer, intent(in), optional :: nprocs
    integer, intent(in), optional :: myrank
 
    integer :: n
    !-----------------------------------------------------------------------------
    
    this%refElem2D => refelem
    call meshbase_init( this, &
      meshbase2d_dimtype_num, refelem, &
      nlocalmeshperprc, 4,             &
      nprocs                           )
 
    allocate( this%lcmesh_list(this%LOCAL_MESH_NUM) )
    do n=1, this%LOCAL_MESH_NUM
      call localmesh2d_init( this%lcmesh_list(n), n, refelem, myrank )
    end do
 
    call this%SetDimInfo( meshbase2d_dimtypeid_x, "x", "m", "X-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_y, "y", "m", "Y-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_xy, "xy", "m", "XY-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_xyt, "xyt", "m", "XY-coordinate" )
 
    return

References scale_localmesh_2d::localmesh2d_init(), meshbase2d_dimtype_num, meshbase2d_dimtypeid_x, meshbase2d_dimtypeid_xy, meshbase2d_dimtypeid_xyt, meshbase2d_dimtypeid_y, and scale_mesh_base::meshbase_init().

Referenced by scale_mesh_base2d::meshbase2d_generate::meshbase2d_generate(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init(), and scale_mesh_rectdom2d::meshrectdom2d_init().

meshbase2d_final()

subroutine, public scale_mesh_base2d::meshbase2d_final

(

class(meshbase2d), intent(inout)

this

)

Finalize an object to manage a 2D computational mesh.

Definition at line 111 of file scale_mesh_base2d.F90.

    implicit none
    class(MeshBase2D), intent(inout) :: this
    integer :: n
    !-----------------------------------------------------------------------------
  
    if ( allocated ( this%lcmesh_list ) ) then 
      do n=1, this%LOCAL_MESH_NUM
        call localmesh2d_final( this%lcmesh_list(n), this%isGenerated )
      end do
  
      deallocate( this%lcmesh_list )
    end if
    
    call meshbase_final(this)
 
    return

References scale_localmesh_2d::localmesh2d_final(), and scale_mesh_base::meshbase_final().

Referenced by scale_mesh_base2d::meshbase2d_generate::meshbase2d_generate(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init(), and scale_mesh_rectdom2d::meshrectdom2d_init().

meshbase2d_setgeometricinfo()

subroutine, public scale_mesh_base2d::meshbase2d_setgeometricinfo	(	type(localmesh2d), intent(inout)	lcmesh,
		external subroutine(real(rp), dimension(elem%np), intent(out) x, real(rp), dimension(elem%np), intent(out) y, real(rp), dimension(elem%np), intent(out) xr, real(rp), dimension(elem%np), intent(out) xs, real(rp), dimension(elem%np), intent(out) yr, real(rp), dimension(elem%np), intent(out) ys, real(rp), dimension(elem%nv), intent(in) vx, real(rp), dimension(elem%nv), intent(in) vy, type(elementbase2d), intent(in) elem)	coord_conv,
		external subroutine(real(rp), dimension(elem%nfptot,2), intent(out) normal_fn, real(rp), dimension(elem%nfptot,2,2), intent(in) escale_f, integer, dimension(elem%nfp,elem%nfaces), intent(in) fid, type(elementbase2d), intent(in) elem)	calc_normal )

Definition at line 144 of file scale_mesh_base2d.F90.

    use scale_mesh_base, only: meshbase_setgeometricinfo
    implicit none
    type(LocalMesh2D), intent(inout) :: lcmesh
    interface
      subroutine coord_conv( x, y, xr, xs, yr, ys, &
        vx, vy, elem )
        import elementbase2d
        import rp
        type(ElementBase2D), intent(in) :: elem
        real(RP), intent(out) :: x(elem%Np), y(elem%Np)
        real(RP), intent(out) :: xr(elem%Np), xs(elem%Np), yr(elem%Np), ys(elem%Np)
        real(RP), intent(in) :: vx(elem%Nv), vy(elem%Nv)        
      end subroutine coord_conv
      subroutine calc_normal( normal_fn, &
        Escale_f, fid, elem )
        import elementbase2d
        import rp
        type(ElementBase2D), intent(in) :: elem
        real(RP), intent(out) :: normal_fn(elem%NfpTot,2)
        integer, intent(in) :: fid(elem%Nfp,elem%Nfaces)
        real(RP), intent(in) :: Escale_f(elem%NfpTot,2,2)
      end subroutine calc_normal      
    end interface
 
    class(ElementBase2D), pointer :: refElem
    integer :: ke
    integer :: f
    integer :: i, j
    integer :: d
    integer :: node_ids(lcmesh%refElem%Nv)
    real(RP) :: vx(lcmesh%refElem%Nv), vy(lcmesh%refElem%Nv)
    real(RP) :: xr(lcmesh%refElem%Np), xs(lcmesh%refElem%Np)
    real(RP) :: yr(lcmesh%refElem%Np), ys(lcmesh%refElem%Np)
    integer :: fmask(lcmesh%refElem%NfpTot)
    integer :: fid(lcmesh%refElem2D%Nfp,lcmesh%refElem2D%Nfaces)
    real(RP) :: Escale_f(lcmesh%refElem%NfpTot,2,2)
 
    !-----------------------------------------------------------------------------
 
    refelem => lcmesh%refElem2D
 
    call meshbase_setgeometricinfo(lcmesh, 2)
 
    allocate( lcmesh%G_ij(refelem%Np,lcmesh%Ne,2,2) )
    allocate( lcmesh%GIJ (refelem%Np,lcmesh%Ne,2,2) )
    allocate( lcmesh%lon(refelem%Np,lcmesh%Ne) )
    allocate( lcmesh%lat(refelem%Np,lcmesh%Ne) )
    !$acc enter data create(lcmesh%G_ij, lcmesh%GIJ, lcmesh%lon, lcmesh%lat)
 
    fmask(:) = reshape(refelem%Fmask, shape(fmask))
    do f=1, refelem%Nfaces
    do i=1, refelem%Nfp
       fid(i,f) = i + (f-1)*refelem%Nfp
    end do
    end do
 
    do ke=1, lcmesh%Ne
      node_ids(:) = lcmesh%EToV(ke,:)
      vx(:) = lcmesh%pos_ev(node_ids(:),1)
      vy(:) = lcmesh%pos_ev(node_ids(:),2)
      call coord_conv( &
        lcmesh%pos_en(:,ke,1), lcmesh%pos_en(:,ke,2), xr, xs, yr, ys, & ! (out)
        vx, vy, refelem )                                             ! (in)
      
      lcmesh%J(:,ke) = - xs*yr + xr*ys
 
      lcmesh%Escale(:,ke,1,1) =   ys/lcmesh%J(:,ke)
      lcmesh%Escale(:,ke,1,2) = - xs/lcmesh%J(:,ke)
      lcmesh%Escale(:,ke,2,1) = - yr/lcmesh%J(:,ke)
      lcmesh%Escale(:,ke,2,2) =   xr/lcmesh%J(:,ke)
      
      !* Face
      
      !
      !mesh%fx(:,n) = mesh%x(fmask(:),n)
      !mesh%fy(:,n) = mesh%y(fmask(:),n)
 
      ! Calculate normal vectors
      do j=1, 2
      do i=1, 2
        escale_f(:,i,j) = lcmesh%Escale(fmask(:),ke,i,j)
      end do
      end do
      
      call calc_normal( lcmesh%normal_fn(:,ke,:), & ! (out)
        escale_f, fid, refelem )                   ! (in)
      
      lcmesh%sJ(:,ke) = sqrt( lcmesh%normal_fn(:,ke,1)**2 + lcmesh%normal_fn(:,ke,2)**2 )
      do d=1, 2
        lcmesh%normal_fn(:,ke,d) = lcmesh%normal_fn(:,ke,d)/lcmesh%sJ(:,ke)
      end do
      lcmesh%sJ(:,ke) = lcmesh%sJ(:,ke)*lcmesh%J(fmask(:),ke)
 
      lcmesh%Fscale(:,ke) = lcmesh%sJ(:,ke)/lcmesh%J(fmask(:),ke)       
    end do
    !$acc update device(lcmesh%pos_en, lcmesh%normal_fn, lcmesh%sJ, lcmesh%J, lcmesh%Escale, lcmesh%Fscale)
 
    !$omp parallel
    !$omp workshare
    lcmesh%Gsqrt (:,:)     = 1.0_rp
    lcmesh%GIJ   (:,:,1,1) = 1.0_rp
    lcmesh%GIJ   (:,:,2,1) = 0.0_rp
    lcmesh%GIJ   (:,:,1,2) = 0.0_rp
    lcmesh%GIJ   (:,:,2,2) = 1.0_rp
    lcmesh%G_ij  (:,:,1,1) = 1.0_rp
    lcmesh%G_ij  (:,:,2,1) = 0.0_rp
    lcmesh%G_ij  (:,:,1,2) = 0.0_rp
    lcmesh%G_ij  (:,:,2,2) = 1.0_rp
    !$omp end workshare
    !$omp end parallel
    !$acc update device(lcmesh%Gsqrt, lcmesh%GIJ, lcmesh%G_ij)
 
    return

References scale_mesh_base::meshbase_setgeometricinfo().

Referenced by scale_mesh_base2d::meshbase2d_generate::meshbase2d_generate(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_setuplocaldom(), and scale_mesh_rectdom2d::meshrectdom2d_setuplocaldom().

Variable Documentation

meshbase2d_dimtype_num

integer, public scale_mesh_base2d::meshbase2d_dimtype_num = 4

Definition at line 57 of file scale_mesh_base2d.F90.

  integer, public :: MESHBASE2D_DIMTYPE_NUM   = 4

Referenced by scale_file_history_meshfield::file_history_meshfield_in::file_history_meshfield_in3d(), and meshbase2d_init().

meshbase2d_dimtypeid_x

integer, public scale_mesh_base2d::meshbase2d_dimtypeid_x = 1

Definition at line 58 of file scale_mesh_base2d.F90.

  integer, public :: MESHBASE2D_DIMTYPEID_X   = 1

Referenced by scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d(), scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d_cubedsphere(), scale_file_history_meshfield::file_history_meshfield_in::file_history_meshfield_in3d(), meshbase2d_init(), and scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init().

meshbase2d_dimtypeid_y

integer, public scale_mesh_base2d::meshbase2d_dimtypeid_y = 2

Definition at line 59 of file scale_mesh_base2d.F90.

  integer, public :: MESHBASE2D_DIMTYPEID_Y   = 2

Referenced by scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d(), scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d_cubedsphere(), scale_file_history_meshfield::file_history_meshfield_in::file_history_meshfield_in3d(), meshbase2d_init(), and scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init().

meshbase2d_dimtypeid_xy

integer, public scale_mesh_base2d::meshbase2d_dimtypeid_xy = 3

Definition at line 60 of file scale_mesh_base2d.F90.

  integer, public :: MESHBASE2D_DIMTYPEID_XY  = 3

Referenced by mod_atmos_mesh_gm::atmosmeshgm_init(), mod_atmos_mesh_rm::atmosmeshrm_init(), scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d(), scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d_cubedsphere(), meshbase2d_init(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init(), and mod_mktopo::mktopo_write().

meshbase2d_dimtypeid_xyt

integer, public scale_mesh_base2d::meshbase2d_dimtypeid_xyt = 4

Definition at line 61 of file scale_mesh_base2d.F90.

  integer, public :: MESHBASE2D_DIMTYPEID_XYT = 4

Referenced by scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d(), scale_file_common_meshfield::file_common_meshfield_get_dims::file_common_meshfield_get_dims2d_cubedsphere(), meshbase2d_init(), and scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init().