1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: UPDATE_OLD C 4 ! Purpose: Update the stored previous-time-step values of certain C 5 ! field variables C 6 ! *****Remember to modify reset_new also 7 ! C 8 ! Author: M. Syamlal Date: 21-JAN-92 C 9 ! Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti, Date: 29-JAN-92 C 10 ! W. Rogers C 11 ! C 12 ! Revision Number: 1 C 13 ! Purpose: Store old solids velocity values C 14 ! Author: M. Syamlal Date: 17-JUN-92 C 15 ! Reviewer: S. Venkatesan Date: 11-DEC-92 C 16 ! C 17 ! Literature/Document References: C 18 ! C 19 ! Variables referenced: ROP_g, EP_g, ROP_s, IJKMAX2, MMAX, U_s, V_s, C 20 ! W_s C 21 ! C 22 ! Variables modified: ROP_go, ROP_so, IJK, M, U_so, V_so, W_so C 23 ! C 24 ! Local variables: NONE C 25 ! C 26 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 27 ! 28 SUBROUTINE UPDATE_OLD 29 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 30 !...Switches: -xf 31 ! 32 !----------------------------------------------- 33 ! M o d u l e s 34 !----------------------------------------------- 35 USE param 36 USE param1 37 USE parallel 38 USE fldvar 39 USE geometry 40 USE indices 41 USE physprop 42 USE run 43 USE trace 44 USE visc_s 45 USE scalars 46 47 IMPLICIT NONE 48 !----------------------------------------------- 49 ! G l o b a l P a r a m e t e r s 50 !----------------------------------------------- 51 !----------------------------------------------- 52 ! L o c a l P a r a m e t e r s 53 !----------------------------------------------- 54 !----------------------------------------------- 55 ! L o c a l V a r i a b l e s 56 !----------------------------------------------- 57 ! 58 ! Indices 59 ! 60 INTEGER :: M 61 !----------------------------------------------- 62 63 EP_GO(:) = EP_G(:) 64 P_GO(:) = P_G(:) 65 P_STARO(:) = P_STAR(:) 66 RO_GO(:) = RO_G(:) 67 ROP_GO(:) = ROP_G(:) 68 U_GO(:) = U_G(:) 69 V_GO(:) = V_G(:) 70 W_GO(:) = W_G(:) 71 IF (ENERGY_EQ) T_GO(:) = T_G(:) 72 IF (SPECIES_EQ(0)) THEN 73 IF (NMAX(0) > 0) THEN 74 X_GO(:,:NMAX(0)) = X_G(:,:NMAX(0)) 75 ENDIF 76 ENDIF 77 78 IF (NScalar > 0) THEN 79 ScalarO(:,:NScalar) = Scalar(:,:NScalar) 80 ENDIF 81 82 IF (K_Epsilon) THEN 83 K_Turb_GO(:) = K_Turb_G(:) 84 E_Turb_GO(:) = E_Turb_G(:) 85 ENDIF 86 87 !!!!$omp parallel do private(M,IJK,N) 88 DO M = 1, MMAX 89 ROP_SO(:,M) = ROP_S(:,M) 90 ! add by rong 91 IF(Call_DQMOM) D_Po(:,M)=D_P(:,M) 92 ! IF (NScalar>0) ome_o(:,M)=ome(:,M) 93 ! add by rong 94 IF (ENERGY_EQ) T_SO(:,M) = T_S(:,M) 95 IF (GRANULAR_ENERGY) THEN 96 THETA_MO(:,M) = THETA_M(:,M) 97 TRD_S_CO(:,M) = TRD_S_C(:,M) 98 ENDIF 99 U_SO(:,M) = U_S(:,M) 100 V_SO(:,M) = V_S(:,M) 101 W_SO(:,M) = W_S(:,M) 102 IF (SPECIES_EQ(M)) THEN 103 IF (NMAX(M) > 0) THEN 104 X_SO(:,M,:NMAX(M)) = X_S(:,M,:NMAX(M)) 105 RO_SO(:,M) = RO_S(:,M) 106 ENDIF 107 ENDIF 108 END DO 109 110 RETURN 111 END SUBROUTINE UPDATE_OLD 112