cfupdateold.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
!
!  Subroutine: CFUPDATEOLD
!  Purpose: Update old arrays to store current time step information
!           before the new arrays are updated
!
!  Author: Jay Boyalakuntla                           Date: 12-Jun-04
!  Reviewer:                                          Date:
!
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

      SUBROUTINE cfupdateold

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE discretelement, only: max_pip, des_pos_old, des_pos_new, des_vel_old, des_vel_new, omega_old, omega_new
      USE functions, ONLY: is_nonexistent, is_ghost
      USE run
      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
! Loop counters (no. particles)
      INTEGER LL
!-----------------------------------------------

!!$omp parallel do if(max_pip .ge. 10000) default(shared)         &
!!$omp private(ll)                    &
!!$omp schedule (guided,50)
      DO ll = 1, max_pip
         IF(is_nonexistent(ll) .or. is_ghost(ll)) cycle

         des_pos_old(ll,:)  = des_pos_new(ll,:)
         des_vel_old(ll,:)  = des_vel_new(ll,:)
         omega_old(ll,:)    = omega_new(ll,:)

      ENDDO
!!$omp end parallel do

      RETURN
      END SUBROUTINE cfupdateold