1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Subroutine: ADJUST_ROP C 4 ! Purpose: Remove small negative values of density. C 5 ! C 6 ! C 7 ! Author: M. Syamlal Date: 7-AUG-96 C 8 ! Reviewer: Date: C 9 ! C 10 ! C 11 ! Literature/Document References: C 12 ! C 13 ! Variables referenced: C 14 ! Variables modified: C 15 ! Local variables: C 16 ! C 17 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 18 19 SUBROUTINE ADJUST_ROP(ROP) 20 21 !----------------------------------------------- 22 ! Modules 23 !----------------------------------------------- 24 USE param 25 USE param1 26 USE geometry 27 USE indices 28 USE compar 29 USE functions 30 IMPLICIT NONE 31 !----------------------------------------------- 32 ! Dummy arguments 33 !----------------------------------------------- 34 ! density 35 DOUBLE PRECISION, INTENT(INOUT) :: ROP(DIMENSION_3) 36 !----------------------------------------------- 37 ! Local variables 38 !----------------------------------------------- 39 ! Indices 40 INTEGER :: IJK 41 !----------------------------------------------- 42 43 DO IJK = ijkstart3, ijkend3 44 IF (FLUID_AT(IJK)) ROP(IJK) = DMAX1(ZERO,ROP(IJK)) 45 ENDDO 46 47 RETURN 48 END SUBROUTINE ADJUST_ROP 49 50 51