update_old.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: UPDATE_OLD                                             C
!  Purpose: Update the stored previous-time-step values of certain     C
!           field variables                                            C
!    *****Remember to modify reset_new also
!                                                                      C
!  Author: M. Syamlal                                 Date: 21-JAN-92  C
!  Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti,  Date: 29-JAN-92  C
!           W. Rogers                                                  C
!                                                                      C
!  Revision Number: 1                                                  C
!  Purpose: Store old solids velocity values                           C
!  Author: M. Syamlal                                 Date: 17-JUN-92  C
!  Reviewer: S. Venkatesan                            Date: 11-DEC-92  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced: ROP_g, EP_g, ROP_s, IJKMAX2, MMAX, U_s, V_s,  C
!                        W_s                                           C
!                                                                      C
!  Variables modified: ROP_go, ROP_so, IJK, M, U_so, V_so, W_so C
!                                                                      C
!  Local variables: NONE                                               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE update_old
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE parallel
      USE fldvar
      USE geometry
      USE indices
      USE physprop
      USE run
      USE trace
      USE visc_s
      USE scalars

      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   V a r i a b l e s
!-----------------------------------------------
!
!                    Indices
!
      INTEGER ::  M
!-----------------------------------------------

      ep_go(:) = ep_g(:)
      p_go(:) = p_g(:)
      p_staro(:) = p_star(:)
      ro_go(:) = ro_g(:)
      rop_go(:) = rop_g(:)
      u_go(:) = u_g(:)
      v_go(:) = v_g(:)
      w_go(:) = w_g(:)
      IF (energy_eq) t_go(:) = t_g(:)
        IF (species_eq(0)) THEN
          IF (nmax(0) > 0) THEN
             x_go(:,:nmax(0)) = x_g(:,:nmax(0))
          ENDIF
      ENDIF

      IF (nscalar > 0) THEN
        scalaro(:,:nscalar) = scalar(:,:nscalar)
      ENDIF

      IF (k_epsilon) THEN
        k_turb_go(:) = k_turb_g(:)
        e_turb_go(:) = e_turb_g(:)
      ENDIF

!!!!$omp parallel do private(M,IJK,N)
      DO m = 1, mmax
            rop_so(:,m) = rop_s(:,m)
! add by rong
       IF(call_dqmom) d_po(:,m)=d_p(:,m)
!      IF (NScalar>0) ome_o(:,M)=ome(:,M)
! add by rong
            IF (energy_eq) t_so(:,m) = t_s(:,m)
            IF (granular_energy) THEN
               theta_mo(:,m) = theta_m(:,m)
               trd_s_co(:,m) = trd_s_c(:,m)
            ENDIF
            u_so(:,m) = u_s(:,m)
            v_so(:,m) = v_s(:,m)
            w_so(:,m) = w_s(:,m)
            IF (species_eq(m)) THEN
               IF (nmax(m) > 0) THEN
                  x_so(:,m,:nmax(m)) = x_s(:,m,:nmax(m))
                  ro_so(:,m) = ro_s(:,m)
               ENDIF
            ENDIF
      END DO

      RETURN
      END SUBROUTINE update_old