reset_new.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: RESET_NEW                                              C
!  Purpose: Reset the new variables with the stored previous-time-step C
!           values of field variables.                                 C
!    *****Remember to modify update_old also
!                                                                      C
!  Author: M. Syamlal                                 Date: FEB-6-97   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 reset_new

!...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 fldvar
      USE physprop, only: mmax, nmax
      USE run, only: energy_eq, call_dqmom, granular_energy, k_epsilon, species_eq
      USE scalars, only: nscalar
      USE trace, only: trd_s_c, trd_s_co

      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
!
!                    error index
      INTEGER :: IER

!-----------------------------------------------

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

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

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

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

          ro_s(:,m) = ro_so(:,m)
        ENDIF
      END DO

!     Recalculate all coefficients
      CALL calc_coeff_all (0, ier)

      RETURN
      END SUBROUTINE reset_new

!// Comments on the modifications for DMP version implementation
!// 120 Replaced the index for initialization: (:IJKMAX2) to just (:)