calc_dtpic.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!  Module name: MAKE_ARRAYS_DES                                        !
!  Author: Jay Boyalakuntla                           Date: 12-Jun-04  !
!                                                                      !
!  Purpose: DES - allocating DES arrays
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE calc_dtpic

      USE compar
      USE cutcell
      USE des_rxns
      USE stl_functions_des
      USE des_thermo
      USE discretelement
      USE functions
      USE funits
      USE geometry
      USE param1
      USE run
      use desmpi
      use error_manager
      use mfix_pic, only: cfl_pic, dtpic_cfl, dtpic_taup
      use mfix_pic, only: dtpic_max
      use mpi_utility

      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
      INTEGER :: L, PC


! MPPIC related quantities
      DOUBLE PRECISION :: DTPIC_TMPX, DTPIC_TMPY, DTPIC_TMPZ
!-----------------------------------------------
! Include statement functions
!-----------------------------------------------

      dtpic_cfl = large_number

      pc = 1
      DO l = 1, max_pip
      IF(pc.GT.pip) EXIT
         IF(is_nonexistent(l)) cycle
         pc = pc+1
         IF(is_ghost(l) .or. is_entering_ghost(l) .or. is_exiting_ghost(l)) cycle

         dtpic_tmpx = (cfl_pic*dx(pijk(l,1)))/&
            (abs(des_vel_new(l,1))+small_number)
         dtpic_tmpy = (cfl_pic*dy(pijk(l,2)))/&
            (abs(des_vel_new(l,2))+small_number)
         dtpic_tmpz = large_number
         IF(do_k) dtpic_tmpz = (cfl_pic*dz(pijk(l,3)))/&
            (abs(des_vel_new(l,3))+small_number)

         dtpic_cfl = min(dtpic_tmpx, dtpic_tmpy, dtpic_tmpz)
      ENDDO

      CALL global_all_max(dtpic_cfl)

      dtpic_max = min(dtpic_cfl, dtpic_taup)
      dtsolid = dtpic_max

      WRITE(err_msg,2000) dtpic_cfl, dtpic_taup, dtsolid
      CALL flush_err_msg(header=.false., footer=.false.)

 2000 FORMAT('DTPIC BASED ON CFL AND TAUP:', 2x, 2(2x,g11.4),          &
         /'DTSOLID set to ', g11.4)

      RETURN
      END SUBROUTINE calc_dtpic