/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/reset_fluid_vep.f90

Go to the documentation of this file.

subroutine reset_fluid_vep
  use def_velocity_potential, only : vep
  use grid_parameter, only : nrf, ntf, npf, &
  &                          ntfeq, npfxzp, npfxzm, NS_shape
  implicit none
  integer :: ir, it, ip
!
! Impose symmetry of the shape to velocity potential
! NS_shape = JB : tri-axial configuration
! NS_shape = ML : axi-symmetric configuration
!
  if (NS_shape.eq.'IB') then
! Impose xz and equatorial symmetry
    vep(0:nrf,0,0:npf) = 0.0d0
    vep(0:nrf,0:ntf,npfxzp) = 0.0d0
    vep(0:nrf,0:ntf,npfxzm) = 0.0d0
    vep(0:nrf,0:ntf,npf) = 0.0d0
!                          
    do it = 0, ntfeq
      do ip = 0, npfxzm
        vep(0:nrf,ntf-it,    ip) =   vep(0:nrf,it,ip)
        vep(0:nrf,    it,npf-ip) = - vep(0:nrf,it,ip)
        vep(0:nrf,ntf-it,npf-ip) = - vep(0:nrf,it,ip)
      end do
    end do
  end if
!
end subroutine reset_fluid_vep

---