Class

xyz_bc_module

In:

../src/utils/xyz_bc_module.F90

3 次元 (xyz 方向) 不等間隔交互格子有限差分モデル用境界条件モジュール

概要

xyz_bc_module は, 3 次元 (xyz 方向) 不等間隔交互格子を用いた有限差分法に基づく数値モデルのための, 境界条件設定 Fortran 90 副プログラムを提供する.

このモジュールは xyz_module の下位モジュールである. 下請けモジュールとして data_type, xyz_base_module, x_bc_module, y_bc_module, z_bc_module モジュールを用いている.

手続きの命名法

境界条件を設定する副プログラムは

  Boundary[境界条件を示す文字列][...]_(入力配列の次元情報)

のように命名されている. 境界条件を示す文字列の種類は x_bc_module, y_bc_module, z_bc_module の場合と同じである.

Methods

BoundaryXAsym_pyz BoundaryXAsym_xqz BoundaryXAsym_xyr BoundaryXAsym_xyz BoundaryXCyc_pyz BoundaryXCyc_pyz BoundaryXCyc_xqz BoundaryXCyc_xqz BoundaryXCyc_xyr BoundaryXCyc_xyr BoundaryXCyc_xyz BoundaryXCyc_xyz BoundaryXSym_pyz BoundaryXSym_xqz BoundaryXSym_xyr BoundaryXSym_xyz BoundaryYASym_xqz BoundaryYAsym_pyz BoundaryYAsym_xyr BoundaryYAsym_xyz BoundaryYCyc_pyz BoundaryYCyc_xqz BoundaryYCyc_xyr BoundaryYCyc_xyz BoundaryYSym_pyz BoundaryYSym_xqz BoundaryYSym_xyr BoundaryYSym_xyz BoundaryZAsym_pyz BoundaryZAsym_xqz BoundaryZAsym_xyr BoundaryZAsym_xyz BoundaryZCyc_pyz BoundaryZCyc_xqz BoundaryZCyc_xyr BoundaryZCyc_xyz BoundaryZSym_pyz BoundaryZSym_xqz BoundaryZSym_xyr BoundaryZSym_xyz

Included Modules

dc_types gridset mpi_wrapper

Public Instance methods

BoundaryXAsym_pyz( pyz_Var )

*Subroutine* :
`pyz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXAsym_pyz(pyz_Var)
    real(DP),intent(inout) :: pyz_Var(imin:imax,jmin:jmax,kmin:kmax) 
    
    ! x 方向に反対称境界条件を適用する
    pyz_Var(0,:,:) = 0.0d0
    pyz_Var(nx,:,:) = 0.0d0
          
    do ix = 1, xmg-1
      pyz_Var(-ix,:,:) = - pyz_Var(ix,:,:)
    end do
    
    do ix = 1, xmg
      pyz_Var(nx+ix,:,:) = - pyz_Var(nx-ix,:,:)
    end do
    
  end subroutine BoundaryXAsym_pyz

BoundaryXAsym_xqz( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXAsym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に反対称境界条件を適用する

    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = - xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = - xaa_Var(nx+1-ix,:,:)
    end do
        
  end subroutine BoundaryXAsym_xaa

BoundaryXAsym_xyr( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXAsym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に反対称境界条件を適用する

    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = - xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = - xaa_Var(nx+1-ix,:,:)
    end do
        
  end subroutine BoundaryXAsym_xaa

BoundaryXAsym_xyz( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXAsym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に反対称境界条件を適用する

    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = - xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = - xaa_Var(nx+1-ix,:,:)
    end do
        
  end subroutine BoundaryXAsym_xaa

BoundaryXCyc_pyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax, jmin:jmax, kmin:kmax)` :	`real(8), intent(inout)`

  subroutine MPIBoundaryXCyc_aaa(aaa_Var)
    
    implicit none
    
    !変数定義
    real(8), intent(inout)  :: aaa_Var(imin:imax, jmin:jmax, kmin:kmax)
    integer            :: idest_a, idep_a, idest_b, idep_b
    integer, parameter :: nvars = 1
    real(8)            :: sbuf_a( Xmg, ny, nz )
    real(8)            :: rbuf_a( Xmg, ny, nz )
    real(8)            :: sbuf_b( Xmg, ny, nz )
    real(8)            :: rbuf_b( Xmg, ny, nz )
    integer            :: ireqs_a, ireqr_a, ireqs_b, ireqr_b
    integer            :: i, j, k
    
    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_a)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_a( i, j, k ) = aaa_var( nx + 1 - i , j, k ) 
        end do
      end do
    end do

    idest_a = mod(( myrank+1 )       , nprocs)   !送信先
    idep_a  = mod(( myrank-1 )+nprocs, nprocs)   !受信元

    call MPIWrapperISend( idest_a, xmg, ny, nz, sbuf_a, ireqs_a ) !送信
    call MPIWrapperIRecv( idep_a , xmg, ny, nz, rbuf_a, ireqr_a ) !受信
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_b)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_b( i, j, k ) = aaa_var( i, j, k ) 
        end do
      end do
    end do
    
    idest_b = mod(( myrank-1 )+nprocs, nprocs)    !送信先
    idep_b  = mod(( myrank+1 )       , nprocs)    !受信元

    call MPIWrapperISend( idest_b, xmg, ny, nz, sbuf_b, ireqs_b ) !送信
    call MPIWrapperIRecv( idep_b , xmg, ny, nz, rbuf_b, ireqr_b ) !受信

    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_a)を代入する.
    call MPIWrapperWait( ireqs_a )
    call MPIWrapperWait( ireqr_a )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var( 1-i , j, k ) = rbuf_a( i, j, k ) 
        end do
      end do
    end do
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_b)を代入する.
    call MPIWrapperWait( ireqs_b )
    call MPIWrapperWait( ireqr_b )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var(nx + i , j, k ) = rbuf_b( i, j, k ) 
        end do
      end do
    end do
    
  end subroutine MPIBoundaryXCyc_aaa

BoundaryXCyc_pyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

   subroutine BoundaryXCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! x 方向に周期境界条件を適用する

      do ix = 1, xmg
        aaa_Var(1-ix,:,:)  = aaa_Var(nx+1-ix,:,:)
        aaa_Var(nx+ix,:,:) = aaa_Var(ix,:,:)
      end do
 
   end subroutine BoundaryXCyc_aaa

BoundaryXCyc_xqz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax, jmin:jmax, kmin:kmax)` :	`real(8), intent(inout)`

  subroutine MPIBoundaryXCyc_aaa(aaa_Var)
    
    implicit none
    
    !変数定義
    real(8), intent(inout)  :: aaa_Var(imin:imax, jmin:jmax, kmin:kmax)
    integer            :: idest_a, idep_a, idest_b, idep_b
    integer, parameter :: nvars = 1
    real(8)            :: sbuf_a( Xmg, ny, nz )
    real(8)            :: rbuf_a( Xmg, ny, nz )
    real(8)            :: sbuf_b( Xmg, ny, nz )
    real(8)            :: rbuf_b( Xmg, ny, nz )
    integer            :: ireqs_a, ireqr_a, ireqs_b, ireqr_b
    integer            :: i, j, k
    
    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_a)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_a( i, j, k ) = aaa_var( nx + 1 - i , j, k ) 
        end do
      end do
    end do

    idest_a = mod(( myrank+1 )       , nprocs)   !送信先
    idep_a  = mod(( myrank-1 )+nprocs, nprocs)   !受信元

    call MPIWrapperISend( idest_a, xmg, ny, nz, sbuf_a, ireqs_a ) !送信
    call MPIWrapperIRecv( idep_a , xmg, ny, nz, rbuf_a, ireqr_a ) !受信
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_b)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_b( i, j, k ) = aaa_var( i, j, k ) 
        end do
      end do
    end do
    
    idest_b = mod(( myrank-1 )+nprocs, nprocs)    !送信先
    idep_b  = mod(( myrank+1 )       , nprocs)    !受信元

    call MPIWrapperISend( idest_b, xmg, ny, nz, sbuf_b, ireqs_b ) !送信
    call MPIWrapperIRecv( idep_b , xmg, ny, nz, rbuf_b, ireqr_b ) !受信

    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_a)を代入する.
    call MPIWrapperWait( ireqs_a )
    call MPIWrapperWait( ireqr_a )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var( 1-i , j, k ) = rbuf_a( i, j, k ) 
        end do
      end do
    end do
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_b)を代入する.
    call MPIWrapperWait( ireqs_b )
    call MPIWrapperWait( ireqr_b )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var(nx + i , j, k ) = rbuf_b( i, j, k ) 
        end do
      end do
    end do
    
  end subroutine MPIBoundaryXCyc_aaa

BoundaryXCyc_xqz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

   subroutine BoundaryXCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! x 方向に周期境界条件を適用する

      do ix = 1, xmg
        aaa_Var(1-ix,:,:)  = aaa_Var(nx+1-ix,:,:)
        aaa_Var(nx+ix,:,:) = aaa_Var(ix,:,:)
      end do
 
   end subroutine BoundaryXCyc_aaa

BoundaryXCyc_xyr( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax, jmin:jmax, kmin:kmax)` :	`real(8), intent(inout)`

  subroutine MPIBoundaryXCyc_aaa(aaa_Var)
    
    implicit none
    
    !変数定義
    real(8), intent(inout)  :: aaa_Var(imin:imax, jmin:jmax, kmin:kmax)
    integer            :: idest_a, idep_a, idest_b, idep_b
    integer, parameter :: nvars = 1
    real(8)            :: sbuf_a( Xmg, ny, nz )
    real(8)            :: rbuf_a( Xmg, ny, nz )
    real(8)            :: sbuf_b( Xmg, ny, nz )
    real(8)            :: rbuf_b( Xmg, ny, nz )
    integer            :: ireqs_a, ireqr_a, ireqs_b, ireqr_b
    integer            :: i, j, k
    
    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_a)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_a( i, j, k ) = aaa_var( nx + 1 - i , j, k ) 
        end do
      end do
    end do

    idest_a = mod(( myrank+1 )       , nprocs)   !送信先
    idep_a  = mod(( myrank-1 )+nprocs, nprocs)   !受信元

    call MPIWrapperISend( idest_a, xmg, ny, nz, sbuf_a, ireqs_a ) !送信
    call MPIWrapperIRecv( idep_a , xmg, ny, nz, rbuf_a, ireqr_a ) !受信
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_b)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_b( i, j, k ) = aaa_var( i, j, k ) 
        end do
      end do
    end do
    
    idest_b = mod(( myrank-1 )+nprocs, nprocs)    !送信先
    idep_b  = mod(( myrank+1 )       , nprocs)    !受信元

    call MPIWrapperISend( idest_b, xmg, ny, nz, sbuf_b, ireqs_b ) !送信
    call MPIWrapperIRecv( idep_b , xmg, ny, nz, rbuf_b, ireqr_b ) !受信

    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_a)を代入する.
    call MPIWrapperWait( ireqs_a )
    call MPIWrapperWait( ireqr_a )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var( 1-i , j, k ) = rbuf_a( i, j, k ) 
        end do
      end do
    end do
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_b)を代入する.
    call MPIWrapperWait( ireqs_b )
    call MPIWrapperWait( ireqr_b )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var(nx + i , j, k ) = rbuf_b( i, j, k ) 
        end do
      end do
    end do
    
  end subroutine MPIBoundaryXCyc_aaa

BoundaryXCyc_xyr( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

   subroutine BoundaryXCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! x 方向に周期境界条件を適用する

      do ix = 1, xmg
        aaa_Var(1-ix,:,:)  = aaa_Var(nx+1-ix,:,:)
        aaa_Var(nx+ix,:,:) = aaa_Var(ix,:,:)
      end do
 
   end subroutine BoundaryXCyc_aaa

BoundaryXCyc_xyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax, jmin:jmax, kmin:kmax)` :	`real(8), intent(inout)`

  subroutine MPIBoundaryXCyc_aaa(aaa_Var)
    
    implicit none
    
    !変数定義
    real(8), intent(inout)  :: aaa_Var(imin:imax, jmin:jmax, kmin:kmax)
    integer            :: idest_a, idep_a, idest_b, idep_b
    integer, parameter :: nvars = 1
    real(8)            :: sbuf_a( Xmg, ny, nz )
    real(8)            :: rbuf_a( Xmg, ny, nz )
    real(8)            :: sbuf_b( Xmg, ny, nz )
    real(8)            :: rbuf_b( Xmg, ny, nz )
    integer            :: ireqs_a, ireqr_a, ireqs_b, ireqr_b
    integer            :: i, j, k
    
    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_a)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_a( i, j, k ) = aaa_var( nx + 1 - i , j, k ) 
        end do
      end do
    end do

    idest_a = mod(( myrank+1 )       , nprocs)   !送信先
    idep_a  = mod(( myrank-1 )+nprocs, nprocs)   !受信元

    call MPIWrapperISend( idest_a, xmg, ny, nz, sbuf_a, ireqs_a ) !送信
    call MPIWrapperIRecv( idep_a , xmg, ny, nz, rbuf_a, ireqr_a ) !受信
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 送信する部分配列(sbuf_b)を用意する
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          sbuf_b( i, j, k ) = aaa_var( i, j, k ) 
        end do
      end do
    end do
    
    idest_b = mod(( myrank-1 )+nprocs, nprocs)    !送信先
    idep_b  = mod(( myrank+1 )       , nprocs)    !受信元

    call MPIWrapperISend( idest_b, xmg, ny, nz, sbuf_b, ireqs_b ) !送信
    call MPIWrapperIRecv( idep_b , xmg, ny, nz, rbuf_b, ireqr_b ) !受信

    !-------------------------------
    ! 配列の右側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_a)を代入する.
    call MPIWrapperWait( ireqs_a )
    call MPIWrapperWait( ireqr_a )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var( 1-i , j, k ) = rbuf_a( i, j, k ) 
        end do
      end do
    end do
    
    !-------------------------------
    ! 配列の左側を, ノード間で通信する.
    ! 受信した部分配列(rbuf_b)を代入する.
    call MPIWrapperWait( ireqs_b )
    call MPIWrapperWait( ireqr_b )
    
    do k = 1, nz
      do j = 1, ny
        do i = 1, Xmg
          aaa_var(nx + i , j, k ) = rbuf_b( i, j, k ) 
        end do
      end do
    end do
    
  end subroutine MPIBoundaryXCyc_aaa

BoundaryXCyc_xyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

   subroutine BoundaryXCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! x 方向に周期境界条件を適用する

      do ix = 1, xmg
        aaa_Var(1-ix,:,:)  = aaa_Var(nx+1-ix,:,:)
        aaa_Var(nx+ix,:,:) = aaa_Var(ix,:,:)
      end do
 
   end subroutine BoundaryXCyc_aaa

BoundaryXSym_pyz( pyz_Var )

*Subroutine* :
`pyz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXSym_pyz(pyz_Var)
    real(DP),intent(inout) :: pyz_Var(imin:imax,jmin:jmax,kmin:kmax) 
    
    ! x 方向に対称境界条件を適用する

    do ix = 1, xmg-1
      pyz_Var(-ix,:,:)   = pyz_Var(ix,:,:)      
    end do
    
    do ix = 1, xmg
      pyz_Var(nx+ix,:,:) = pyz_Var(nx+1-ix,:,:)
    end do
    
    
  end subroutine BoundaryXSym_pyz

BoundaryXSym_xqz( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXSym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に周期境界条件を適用する
    
    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = xaa_Var(nx+1-ix,:,:)
    end do
    
  end subroutine BoundaryXSym_xaa

BoundaryXSym_xyr( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXSym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に周期境界条件を適用する
    
    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = xaa_Var(nx+1-ix,:,:)
    end do
    
  end subroutine BoundaryXSym_xaa

BoundaryXSym_xyz( xaa_Var )

*Subroutine* :
`xaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

  subroutine BoundaryXSym_xaa(xaa_Var)
    real(DP),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 

    ! x 方向に周期境界条件を適用する
    
    do ix = 1, xmg
      xaa_Var(1-ix,:,:)  = xaa_Var(ix,:,:)
      xaa_Var(nx+ix,:,:) = xaa_Var(nx+1-ix,:,:)
    end do
    
  end subroutine BoundaryXSym_xaa

BoundaryYASym_xqz( xqz_Var )

*Subroutine* :
`xqz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYASym_xqz(xqz_Var)
      real(DP),intent(inout) :: xqz_Var(imin:imax,jmin:jmax,kmin:kmax) 

      ! y 方向に反対称境界条件を適用する

      if (ymg == 0) return 

      xqz_Var(:,0,:) = 0.0d0
      xqz_Var(:,ny,:) = 0.0d0
          
      do jy = 1, ymg-1
        xqz_Var(:,-jy,:) = - xqz_Var(:,jy,:)
      end do

      do jy = 1, ymg
        xqz_Var(:,ny+jy,:) = - xqz_Var(:,ny-jy,:)
      end do

    end subroutine BoundaryYAsym_xqz

BoundaryYAsym_pyz( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYASym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に反対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = - aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = - aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYAsym_aya

BoundaryYAsym_xyr( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYASym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に反対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = - aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = - aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYAsym_aya

BoundaryYAsym_xyz( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYASym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に反対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = - aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = - aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYAsym_aya

BoundaryYCyc_pyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に周期境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aaa_Var(:,1-jy,:)  = aaa_Var(:,ny+1-jy,:)
        aaa_Var(:,ny+jy,:) = aaa_Var(:,jy,:)
      end do

    end subroutine BoundaryYCyc_aaa

BoundaryYCyc_xqz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に周期境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aaa_Var(:,1-jy,:)  = aaa_Var(:,ny+1-jy,:)
        aaa_Var(:,ny+jy,:) = aaa_Var(:,jy,:)
      end do

    end subroutine BoundaryYCyc_aaa

BoundaryYCyc_xyr( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に周期境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aaa_Var(:,1-jy,:)  = aaa_Var(:,ny+1-jy,:)
        aaa_Var(:,ny+jy,:) = aaa_Var(:,jy,:)
      end do

    end subroutine BoundaryYCyc_aaa

BoundaryYCyc_xyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に周期境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aaa_Var(:,1-jy,:)  = aaa_Var(:,ny+1-jy,:)
        aaa_Var(:,ny+jy,:) = aaa_Var(:,jy,:)
      end do

    end subroutine BoundaryYCyc_aaa

BoundaryYSym_pyz( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYSym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYsym_aya

BoundaryYSym_xqz( xqz_Var )

*Subroutine* :
`xqz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYSym_xqz(xqz_Var)
      real(DP),intent(inout) :: xqz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg-1
        xqz_Var(:,-jy,:)   = xqz_Var(:,jy,:)      
      end do

      do jy = 1, ymg
        xqz_Var(:,ny+jy,:) = xqz_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYsym_xqz

BoundaryYSym_xyr( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYSym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYsym_aya

BoundaryYSym_xyz( aya_Var )

*Subroutine* :
`aya_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryYSym_aya(aya_Var)
      real(DP),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! y 方向に対称境界条件を適用する

      if (ymg == 0) return 

      do jy = 1, ymg
        aya_Var(:,1-jy,:)  = aya_Var(:,jy,:)
        aya_Var(:,ny+jy,:) = aya_Var(:,ny+1-jy,:)
      end do

    end subroutine BoundaryYsym_aya

BoundaryZAsym_pyz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZAsym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に反対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = - aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = - aaz_Var(:,:,nz+1-kz)
      end do

    end subroutine BoundaryZAsym_aaz

BoundaryZAsym_xqz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZAsym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に反対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = - aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = - aaz_Var(:,:,nz+1-kz)
      end do

    end subroutine BoundaryZAsym_aaz

BoundaryZAsym_xyr( xyr_Var )

*Subroutine* :
`xyr_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZAsym_xyr(xyr_Var)
      real(DP),intent(inout) :: xyr_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に反対称境界条件を適用する

      xyr_Var(:,:,0) = 0.0d0
      xyr_Var(:,:,nz) = 0.0d0
          
      do kz = 1, zmg-1
        xyr_Var(:,:,-kz) = - xyr_Var(:,:,kz)
      end do

      do kz = 1, zmg
        xyr_Var(:,:,nz+kz) = - xyr_Var(:,:,nz-kz)
      end do

    end subroutine BoundaryZAsym_xyr

BoundaryZAsym_xyz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZAsym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に反対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = - aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = - aaz_Var(:,:,nz+1-kz)
      end do

    end subroutine BoundaryZAsym_aaz

BoundaryZCyc_pyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に周期境界条件を適用する

      do kz = 1, zmg
        aaa_Var(:,:,1-kz)  = aaa_Var(:,:,nz+1-kz)
        aaa_Var(:,:,nz+kz) = aaa_Var(:,:,kz)
      end do

    end subroutine BoundaryZCyc_aaa

BoundaryZCyc_xqz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に周期境界条件を適用する

      do kz = 1, zmg
        aaa_Var(:,:,1-kz)  = aaa_Var(:,:,nz+1-kz)
        aaa_Var(:,:,nz+kz) = aaa_Var(:,:,kz)
      end do

    end subroutine BoundaryZCyc_aaa

BoundaryZCyc_xyr( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に周期境界条件を適用する

      do kz = 1, zmg
        aaa_Var(:,:,1-kz)  = aaa_Var(:,:,nz+1-kz)
        aaa_Var(:,:,nz+kz) = aaa_Var(:,:,kz)
      end do

    end subroutine BoundaryZCyc_aaa

BoundaryZCyc_xyz( aaa_Var )

*Subroutine* :
`aaa_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZCyc_aaa(aaa_Var)
      real(DP),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に周期境界条件を適用する

      do kz = 1, zmg
        aaa_Var(:,:,1-kz)  = aaa_Var(:,:,nz+1-kz)
        aaa_Var(:,:,nz+kz) = aaa_Var(:,:,kz)
      end do

    end subroutine BoundaryZCyc_aaa

BoundaryZSym_pyz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZSym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = aaz_Var(:,:,nz+1-kz)
      end do


    end subroutine BoundaryZSym_aaz

BoundaryZSym_xqz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZSym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = aaz_Var(:,:,nz+1-kz)
      end do


    end subroutine BoundaryZSym_aaz

BoundaryZSym_xyr( xyr_Var )

*Subroutine* :
`xyr_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZSym_xyr(xyr_Var)
      real(DP),intent(inout) :: xyr_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に対称境界条件を適用する

      do kz = 1, zmg-1
        xyr_Var(:,:,-kz)   = xyr_Var(:,:,kz)      
      end do

      do kz = 1, zmg
        xyr_Var(:,:,nz+kz) = xyr_Var(:,:,nz+1-kz)
      end do

    end subroutine BoundaryZSym_xyr

BoundaryZSym_xyz( aaz_Var )

*Subroutine* :
`aaz_Var(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP),intent(inout)`

    subroutine BoundaryZSym_aaz(aaz_Var)
      real(DP),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
      
      ! z 方向に対称境界条件を適用する

      do kz = 1, zmg
        aaz_Var(:,:,1-kz)  = aaz_Var(:,:,kz)
        aaz_Var(:,:,nz+kz) = aaz_Var(:,:,nz+1-kz)
      end do


    end subroutine BoundaryZSym_aaz