subroutine boundary(type, var)
                                                                 !=begin
  !== Dependency

                                                                 !=end  
  !--- 暗黙の型宣言禁止
  implicit none
                                                                 !=begin  
  !== Input
  character(*), intent(in)   :: type

  !== In/Out
  real(8), intent(inout)     :: var(DimXMin:DimXMax, DimZMin:DimZMax)
                                                                 !=end    
  !--- 内部変数
  integer   :: i

   call BeginSub("boundary",                  fmt="%c",                        c1="Set boundary conditions.")

  if (type == "CC") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do
     
     !--- z 方向: 周期境界     
     do i = 0, Margin
        var(:, RegZMin - i) = var(:, RegZMax - i) 
     end do
     do i = 1, Margin
        var(:, RegZMax + i) = var(:, RegZMin + i) 
     end do

  elseif (type == "CsS") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do

     !--- z 方向: 対称
     do i = 0, Margin
        var(:, RegZMin - i) = var(:, RegZMin + 1 + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax + 1 - i )
     end do
     
  elseif (type == "CfS") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do

     !--- z 方向: 対称
     var(:, RegZMin) = 0.0d0
     var(:, RegZMax) = 0.0d0
     do i = 1, Margin
        var(:, RegZMin - i) = var(:, RegZMin + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax - i )
     end do
     
  elseif (type == "CsA") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do

     !--- z 方向: 反対称     
     do i = 0, Margin
        var(:, RegZMin - i) = - var(:, RegZMin + 1 + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i) = - var(:, RegZMax + 1 - i )
     end do
     
  elseif (type == "CfA") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do
     
     !--- z 方向: 反対称     
     var(:, RegZMin) = 0.0d0
     var(:, RegZMax) = 0.0d0
     do i = 1, Margin
        var(:, RegZMin - i) = - var(:, RegZMin + i )        
    end do
     do i = 1, Margin
        var(:, RegZMax + i) = - var(:, RegZMax - i )
     end do

     
  !--- test     
  elseif (type == "CsAS") then 
     !--- x 方向: 周期境界
     do i = 0, Margin
        var(RegXMin - i , RegZMin:RegZMax)  = var(RegXMax - i , RegZMin:RegZMax) 
     end do
     do i = 1, Margin
        var(RegXMax + i , RegZMin:RegZMax)  = var(RegXMin + i , RegZMin:RegZMax) 
     end do

     !--- z 方向: 対称と反対称     
     do i = 0, Margin
        var(:, RegZMin - i) = - var(:, RegZMin + 1 + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax + 1 - i )
     end do


  elseif (type == "sSsS") then 
     !--- x 方向: 対称
     do i = 0, Margin
        var(RegXMin - i,:) = var(RegXMin + 1 + i,: )
     end do
     do i = 1, Margin
        var(RegXMax + i ,:) = var(RegXMax + 1 - i,: )
     end do

     !--- z 方向: 対称
     do i = 0, Margin
        var(:, RegZMin - i) = var(:, RegZMin + 1 + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax + 1 - i )
     end do
     

  elseif (type == "sSfS") then 
     !--- x 方向: 対称
     do i = 0, Margin
        var(RegXMin - i,:) = var(RegXMin + 1 + i,: )
     end do
     do i = 1, Margin
        var(RegXMax + i ,:) = var(RegXMax + 1 - i,: )
     end do

     !--- z 方向: 対称
     var(:,RegZMin) = var(:,RegZMin + 1)

     do i = 1, Margin
        var(:, RegZMin - i) = var(:, RegZMin + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax - i )
     end do

  elseif (type == "fSsS") then 
     !--- x 方向: 対称
     var(RegXMin,:) = var(RegXMin + 1,:)

     do i = 1, Margin
        var(RegXMin - i,:) = var(RegXMin + i,: )
     end do
     do i = 1, Margin
        var(RegXMax + i ,:) = var(RegXMax - i,: )
     end do

     !--- z 方向: 対称
     do i = 0, Margin
        var(:, RegZMin - i) = var(:, RegZMin + 1 + i )
     end do
     do i = 1, Margin
        var(:, RegZMax + i ) = var(:, RegZMax + 1 - i )
     end do
     

  else
     call MessageNotify("Error", "boundary",                        "Unknown Boundary Type", c1=type)

  end if

  call EndSub("boundary")

end subroutine boundary