geom_mod.F90

Go to the documentation of this file.

module geom
  use kinds, only : dp
  implicit none

  real(DP), allocatable, target, dimension(:,:,:,:) :: xm1, ym1, zm1
  real(DP), pointer,             dimension(:,:,:,:) :: xm2, ym2, zm2

  real(DP), allocatable :: jacmi (:,:,:,:)

  real(DP), allocatable, target, dimension(:,:,:,:) :: &
     rxm1, sxm1, txm1, rym1, sym1, tym1, rzm1, szm1, tzm1, jacm1 
  real(DP), pointer,             dimension(:,:,:,:) :: &
     rxm2, sxm2, txm2, rym2, sym2, tym2, rzm2, szm2, tzm2, jacm2

  real(DP), allocatable :: rx(:,:,:) !>

  real(DP), allocatable, dimension(:,:,:,:) :: &
     g1m1, g2m1, g3m1, g4m1, g5m1, g6m1 

  real(DP), allocatable, dimension(:,:,:,:) :: &
     unx, uny, unz, area 

  real(DP) :: dlam

  real(DP), allocatable, dimension(:,:,:,:) :: &
     vnx, vny, vnz, v1x, v1y, v1z, v2x, v2y, v2z 

  real(DP), allocatable, target :: bm1(:,:,:,:) 
  real(DP), pointer             :: bm2(:,:,:,:) 

  real(DP), allocatable, dimension(:,:,:,:) :: &
    binvm1, bintm1, bm2inv, baxm1, yinvm1

  real(DP), allocatable :: bm1lag(:,:,:,:,:)

  real(DP) :: volvm1,volvm2,voltm1,voltm2

  logical :: ifgeom,ifgmsh3,ifvcor,ifsurt,ifmelt,ifwcno, ifbcor 
  logical, allocatable :: ifrzer(:), ifqinp(:,:), ifeppm(:,:)
  logical, allocatable :: iflmsf(:), iflmse(:), iflmsc(:)
  logical, allocatable :: ifmsfc(:,:,:), ifmseg(:,:,:),ifmscr(:,:,:)
  logical, allocatable :: ifnskp(:,:) 

  logical :: bm1_compress !>

  contains

  subroutine init_geom()
    use size_m
    use mesh, only : if_ortho
    implicit none

    allocate( xm1(lx1,ly1,lz1,lelt), ym1(lx1,ly1,lz1,lelt), zm1(lx1,ly1,lz1,lelt) )


    if (if_ortho) then
      allocate(rx(lxd*lyd*lzd,ldim,lelv)) !verified
      allocate(sxm1(lx1,ly1,lz1,1), txm1(lx1,ly1,lz1,1) &
             , rym1(lx1,ly1,lz1,1), tym1(lx1,ly1,lz1,1) &
             , rzm1(lx1,ly1,lz1,1), szm1(lx1,ly1,lz1,1) )
    else
      allocate(rx(lxd*lyd*lzd,ldim*ldim,lelv)) !verified
      allocate(sxm1(lx1,ly1,lz1,lelt), txm1(lx1,ly1,lz1,lelt) &
             , rym1(lx1,ly1,lz1,lelt), tym1(lx1,ly1,lz1,lelt) &
             , rzm1(lx1,ly1,lz1,lelt), szm1(lx1,ly1,lz1,lelt) )
    endif
    allocate( rxm1(lx1,ly1,lz1,lelt), sym1(lx1,ly1,lz1,lelt), tzm1(lx1,ly1,lz1,lelt) ) 
    allocate( jacm1(lx1,ly1,lz1,lelt), jacmi(lx1,ly1,lz1,lelt) )

    allocate( g1m1(lx1,ly1,lz1,lelt), g2m1(lx1,ly1,lz1,lelt), g3m1(lx1,ly1,lz1,lelt) &
            , g4m1(lx1,ly1,lz1,lelt), g5m1(lx1,ly1,lz1,lelt), g6m1(lx1,ly1,lz1,lelt) )

    allocate( unx(lx1,lz1,6,lelt), uny(lx1,lz1,6,lelt), unz(lx1,lz1,6,lelt) )
    allocate( area(lx1,lz1,6,lelt) )

    allocate( &
     ifrzer(lelt),ifqinp(6,lelv),ifeppm(6,lelv) &
    ,iflmsf(0:1),iflmse(0:1),iflmsc(0:1) &
    ,ifmsfc(6,lelt,0:1) &
    ,ifmseg(12,lelt,0:1) &
    ,ifmscr(8,lelt,0:1) &
    ,ifnskp(8,lelt) &
    )

    allocate( bm1(lx1,ly1,lz1,lelt), binvm1(lx1,ly1,lz1,lelv) )
    iflmsf = .false.

  end subroutine init_geom

#if 0
  subroutine compress_geom()
    use size_m, only : lx1,ly1,lz1,nelt, nid
    implicit none

    integer, parameter :: lxyz = lx1*ly1*lz1
    real(DP) :: thresh
    real(DP), external :: dnrm2
    integer :: ie

    bm1_compress = .true.
    thresh = 1.d-4 * dnrm2(lxyz*nelt, bm1, 1) / nelt
    
    do ie = 2, nelt
      if (dnrm2(lxyz,bm1(:,:,:,ie) - bm1(:,:,:,ie-1),1) > thresh) then
        bm1_compress = .false.
        if (nid == 0) then
        write(*,*) bm1(:,1,1,ie)
        write(*,*) bm1(:,1,1,ie-1)
        endif
      endif
    enddo
    write(*,*) "MAX: bm1_compress =", bm1_compress

  end subroutine compress_geom
#endif

end module geom