36,39c36,37
<   use x_base_module, only: imin, imax
<   use y_base_module, only: jmin, jmax
<   use z_base_module, only: kmin, kmax
< !  use xyz_base_module, only : imin, imax, jmin, jmax, kmin, kmax
---
>   use xyz_base_module, only : imin, imax, jmin, jmax, kmin, kmax, &
>     &                         xmargin, ymargin, zmargin, im, jm, km 
91c89
<     module procedure BoundaryXCyc_xaa
---
>     module procedure BoundaryXCyc_aaa
95c93
<     module procedure BoundaryXCyc_xaa
---
>     module procedure BoundaryXCyc_aaa
99c97,101
<     module procedure BoundaryXCyc_xaa
---
>     module procedure BoundaryXCyc_aaa
>   end interface
> 
>   interface BoundaryXCyc_pyz
>     module procedure BoundaryXCyc_aaa
127c129
<     module procedure BoundaryYCyc_aya
---
>     module procedure BoundaryYCyc_aaa
131c133
<     module procedure BoundaryYCyc_aya
---
>     module procedure BoundaryYCyc_aaa
135c137,141
<     module procedure BoundaryYCyc_aya
---
>     module procedure BoundaryYCyc_aaa
>   end interface
> 
>   interface BoundaryYCyc_xqz
>     module procedure BoundaryYCyc_aaa
163c169
<     module procedure BoundaryZCyc_aaz
---
>     module procedure BoundaryZCyc_aaa
167c173
<     module procedure BoundaryZCyc_aaz
---
>     module procedure BoundaryZCyc_aaa
171c177,181
<     module procedure BoundaryZCyc_aaz
---
>     module procedure BoundaryZCyc_aaa
>   end interface
> 
>   interface BoundaryZCyc_xyr
>     module procedure BoundaryZCyc_aaa
173a184,185
>   integer :: ix, jy, kz
> 
179d190
<       integer           :: jy, kz
182,185c193,201
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundarySym_x(xaa_Var(:,jy,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do jy = jmin, jmax
> !          call BoundarySym_x(xaa_Var(:,jy,kz))
> !        end do
> !      end do
> 
>       do ix = 1, xmargin
>         xaa_Var(1-ix,:,:)  = xaa_Var(ix,:,:)
>         xaa_Var(im+ix,:,:) = xaa_Var(im+1-ix,:,:)
192d207
<       integer           :: jy, kz
195,198c210,217
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundarySym_p(pyz_Var(:,jy,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do jy = jmin, jmax
> !          call BoundarySym_p(pyz_Var(:,jy,kz))
> !        end do
> !      end do
> 
>       do ix = 1, xmargin-1
>         pyz_Var(-ix,:,:)   = pyz_Var(ix,:,:)      
200a220,224
>       do ix = 1, xmargin
>         pyz_Var(im+ix,:,:) = pyz_Var(im+1-ix,:,:)
>       end do
> 
> 
205d228
<       integer           :: jy, kz
208,211c231,239
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundaryAsym_x(xaa_Var(:,jy,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do jy = jmin, jmax
> !          call BoundaryAsym_x(xaa_Var(:,jy,kz))
> !        end do
> !      end do
> 
>       do ix = 1, xmargin
>         xaa_Var(1-ix,:,:)  = - xaa_Var(ix,:,:)
>         xaa_Var(im+ix,:,:) = - xaa_Var(im+1-ix,:,:)
218d245
<       integer           :: jy, kz
221,224c248,258
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundaryAsym_p(pyz_Var(:,jy,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do jy = jmin, jmax
> !          call BoundaryAsym_p(pyz_Var(:,jy,kz))
> !        end do
> !      end do
> 
>       pyz_Var(0,:,:) = 0.0d0
>       pyz_Var(im,:,:) = 0.0d0
>           
>       do ix = 1, xmargin-1
>         pyz_Var(-ix,:,:) = - pyz_Var(ix,:,:)
227,237c261,262
<     end subroutine BoundaryXAsym_pyz
< !--------------------------------------------------------------------
<     subroutine BoundaryXCyc_xaa(xaa_Var)
<       real(DBKIND),intent(inout) :: xaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
<       integer           :: jy, kz
<       
<       ! x 方向に周期境界条件を適用する
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundaryCyc_x(xaa_Var(:,jy,kz))
<         end do
---
>       do ix = 1, xmargin
>         pyz_Var(im+ix,:,:) = - pyz_Var(im-ix,:,:)
240c265
<     end subroutine BoundaryXCyc_xaa
---
>     end subroutine BoundaryXAsym_pyz
242,244c267,268
<     subroutine BoundaryXCyc_pyz(pyz_Var)
<       real(DBKIND),intent(inout) :: pyz_Var(imin:imax,jmin:jmax,kmin:kmax) 
<       integer           :: jy, kz
---
>     subroutine BoundaryXCyc_aaa(aaa_Var)
>       real(DBKIND),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
247,250c271,279
<       do kz = kmin, kmax
<         do jy = jmin, jmax
<           call BoundaryCyc_p(pyz_Var(:,jy,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do jy = jmin, jmax
> !          call BoundaryCyc_x(xaa_Var(:,jy,kz))
> !        end do
> !      end do
> 
>       do ix = 1, xmargin
>         aaa_Var(1-ix,:,:)  = aaa_Var(im+1-ix,:,:)
>         aaa_Var(im+ix,:,:) = aaa_Var(ix,:,:)
252,253c281,282
< 
<     end subroutine BoundaryXCyc_pyz
---
>  
>    end subroutine BoundaryXCyc_aaa
257d285
<       integer           :: ix, kz
260,263c288,296
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call boundarySym_y(aya_Var(ix,:,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do ix = imin, imax
> !          call boundarySym_y(aya_Var(ix,:,kz))
> !        end do
> !      end do
> 
>       do jy = 1, ymargin
>         aya_Var(:,1-jy,:)  = aya_Var(:,jy,:)
>         aya_Var(:,jm+jy,:) = aya_Var(:,jm+1-jy,:)
270d302
<       integer           :: ix, kz
273,276c305,316
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call boundarySym_q(xqz_Var(ix,:,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do ix = imin, imax
> !          call boundarySym_q(xqz_Var(ix,:,kz))
> !        end do
> !      end do
> 
>       do jy = 1, ymargin-1
>         xqz_Var(:,-jy,:)   = xqz_Var(:,jy,:)      
>       end do
> 
>       do jy = 1, ymargin
>         xqz_Var(:,jm+jy,:) = xqz_Var(:,jm+1-jy,:)
283d322
<       integer           :: ix, kz
286,289c325,333
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call boundaryASym_y(aya_Var(ix,:,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do ix = imin, imax
> !          call boundaryASym_y(aya_Var(ix,:,kz))
> !        end do
> !      end do
> 
>       do jy = 1, ymargin
>         aya_Var(:,1-jy,:)  = - aya_Var(:,jy,:)
>         aya_Var(:,jm+jy,:) = - aya_Var(:,jm+1-jy,:)
299,302c343,353
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call boundaryASym_q(xqz_Var(ix,:,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do ix = imin, imax
> !          call boundaryASym_q(xqz_Var(ix,:,kz))
> !        end do
> !      end do
> 
>       xqz_Var(:,0,:) = 0.0d0
>       xqz_Var(:,jm,:) = 0.0d0
>           
>       do jy = 1, ymargin-1
>         xqz_Var(:,-jy,:) = - xqz_Var(:,jy,:)
305,315c356,357
<     end subroutine BoundaryYAsym_xqz
< !--------------------------------------------------------------------
<     subroutine BoundaryYCyc_aya(aya_Var)
<       real(DBKIND),intent(inout) :: aya_Var(imin:imax,jmin:jmax,kmin:kmax) 
<       integer           :: ix, kz
<       
<       ! y 方向に周期境界条件を適用する
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call BoundaryCyc_y(aya_Var(ix,:,kz))
<         end do
---
>       do jy = 1, ymargin
>         xqz_Var(:,jm+jy,:) = - xqz_Var(:,jm-jy,:)
318c360
<     end subroutine BoundaryYCyc_aya
---
>     end subroutine BoundaryYAsym_xqz
320,321c362,363
<     subroutine BoundaryYCyc_xqz(xqz_Var)
<       real(DBKIND),intent(inout) :: xqz_Var(imin:imax,jmin:jmax,kmin:kmax) 
---
>     subroutine BoundaryYCyc_aaa(aaa_Var)
>       real(DBKIND),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
325,328c367,375
<       do kz = kmin, kmax
<         do ix = imin, imax
<           call BoundaryCyc_q(xqz_Var(ix,:,kz))
<         end do
---
> !      do kz = kmin, kmax
> !        do ix = imin, imax
> !          call BoundaryCyc_y(aya_Var(ix,:,kz))
> !        end do
> !      end do
> 
>       do jy = 1, ymargin
>         aaa_Var(:,1-jy,:)  = aaa_Var(:,jm+1-jy,:)
>         aaa_Var(:,jm+jy,:) = aaa_Var(:,jy,:)
331c378
<     end subroutine BoundaryYCyc_xqz
---
>     end subroutine BoundaryYCyc_aaa
335d381
<       integer           :: ix, jy
338,341c384,392
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundarySym_z(aaz_Var(ix,jy,:))
<         end do
---
> !      do jy = jmin, jmax
> !        do ix = imin, imax
> !          call boundarySym_z(aaz_Var(ix,jy,:))
> !        end do
> !      end do
> 
>       do kz = 1, zmargin
>         aaz_Var(:,:,1-kz)  = aaz_Var(:,:,kz)
>         aaz_Var(:,:,km+kz) = aaz_Var(:,:,km+1-kz)
343a395
> 
348d399
<       integer           :: ix, jy
351,354c402,413
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundarySym_r(xyr_Var(ix,jy,:))
<         end do
---
> !      do jy = jmin, jmax
> !        do ix = imin, imax
> !          call boundarySym_r(xyr_Var(ix,jy,:))
> !        end do
> !      end do
> 
>       do kz = 1, zmargin-1
>         xyr_Var(:,:,-kz)   = xyr_Var(:,:,kz)      
>       end do
> 
>       do kz = 1, zmargin
>         xyr_Var(:,:,km+kz) = xyr_Var(:,:,km+1-kz)
361d419
<       integer           :: ix, jy
364,367c422,430
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundaryAsym_z(aaz_Var(ix,jy,:))
<         end do
---
> !      do jy = jmin, jmax
> !        do ix = imin, imax
> !          call boundaryAsym_z(aaz_Var(ix,jy,:))
> !        end do
> !      end do
> 
>       do kz = 1, zmargin
>         aaz_Var(:,:,1-kz)  = - aaz_Var(:,:,kz)
>         aaz_Var(:,:,km+kz) = - aaz_Var(:,:,km+1-kz)
374d436
<       integer           :: ix, jy
377,380c439,449
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundaryAsym_r(xyr_Var(ix,jy,:))
<         end do
---
> !      do jy = jmin, jmax
> !        do ix = imin, imax
> !          call boundaryAsym_r(xyr_Var(ix,jy,:))
> !        end do
> !      end do
> 
>       xyr_Var(:,:,0) = 0.0d0
>       xyr_Var(:,:,km) = 0.0d0
>           
>       do kz = 1, zmargin-1
>         xyr_Var(:,:,-kz) = - xyr_Var(:,:,kz)
383,393c452,453
<     end subroutine BoundaryZAsym_xyr
< !--------------------------------------------------------------------
<     subroutine BoundaryZCyc_aaz(aaz_Var)
<       real(DBKIND),intent(inout) :: aaz_Var(imin:imax,jmin:jmax,kmin:kmax) 
<       integer           :: ix, jy
<       
<       ! z 方向に周期境界条件を適用する
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundaryCyc_z(aaz_Var(ix,jy,:))
<         end do
---
>       do kz = 1, zmargin
>         xyr_Var(:,:,km+kz) = - xyr_Var(:,:,km-kz)
396c456
<     end subroutine BoundaryZCyc_aaz
---
>     end subroutine BoundaryZAsym_xyr
398,400c458,459
<     subroutine BoundaryZCyc_xyr(xyr_Var)
<       real(DBKIND),intent(inout) :: xyr_Var(imin:imax,jmin:jmax,kmin:kmax) 
<       integer           :: ix, jy
---
>     subroutine BoundaryZCyc_aaa(aaa_Var)
>       real(DBKIND),intent(inout) :: aaa_Var(imin:imax,jmin:jmax,kmin:kmax) 
403,406c462,470
<       do jy = jmin, jmax
<         do ix = imin, imax
<           call boundaryCyc_z(xyr_Var(ix,jy,:))
<         end do
---
> !      do jy = jmin, jmax
> !        do ix = imin, imax
> !          call boundaryCyc_z(aaz_Var(ix,jy,:))
> !        end do
> !      end do
> 
>       do kz = 1, zmargin
>         aaa_Var(:,:,1-kz)  = aaa_Var(:,:,km+1-kz)
>         aaa_Var(:,:,km+kz) = aaa_Var(:,:,kz)
409c473
<     end subroutine BoundaryZCyc_xyr
---
>     end subroutine BoundaryZCyc_aaa