set_param.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Subroutinee: SET_PARAMETERS                                         !
!  Purpose: Set parameters used in array allocations.                  !
!                                                                      !
!  Author: J.Musser                                  Date: 17-APR-14   !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE set_parameters

! Global Varialbes:
!---------------------------------------------------------------------//
! Number of TFM solids
      use physprop, only: mmax
! Number of DES solids
      use discretelement, only: des_mmax
! Number of species
      use physprop, only: nmax
! Number of scalar equations
      use scalars, only: nscalar

! Domain indices.
      use geometry, only: imax3, jmax3, kmax3, ijkmax3

! Rank specific decompositions.
      use compar, only: ijksize3_all
      USE compar, only: istart3, iend3, istart4, iend4
      USE compar, only: jstart3, jend3, jstart4, jend4
      USE compar, only: kstart3, kend3, kstart4, kend4

! Global Parameters:
!---------------------------------------------------------------------//
! Total number of solids phases
      use param, only: dimension_m
! Maximum number of species.
      use param, only: dimension_n_g ! Gas
      use param, only: dimension_n_s ! Solids
! Maximum number of user-defined scalars
      use param, only: dimension_scalar
      use param, only: dim_scalar2
! Dimension for the upper triangle of an MxM matrix
      use param1, only: dimension_lm
! Total number of species
      use param1, only: dimension_n_all
! Parameter constants.
      use param1, only: undefined_i
! Axis decomposition
      USE param, only: dimension_i, dimension_j, dimension_k
      USE param, only: dimension_3, dimension_4
      USE param, only: dimension_3g, dimension_3l, dimension_3p

! MPI-Domain decompoint and rank flags.
      use compar, only: mype
! Flag for POST_MFIX
      use cdist, only: bdoing_postmfix


      IMPLICIT NONE

! Local Variables:
!---------------------------------------------------------------------//
! Loop index
      INTEGER :: M, MMAX_TOT
!......................................................................!

! The total number of solids.
      dimension_m = max(1, merge(mmax,mmax+des_mmax,bdoing_postmfix))

! Number of gas phase species.
      dimension_n_g = merge(nmax(0), 1, nmax(0) /= undefined_i)

! Number of solids phase species. (Max over all phases)
      dimension_n_s = 1

      mmax_tot = merge(mmax, mmax+des_mmax, bdoing_postmfix)
      DO m = 1, mmax_tot
         IF(nmax(m) /= undefined_i) &
            dimension_n_s = max(dimension_n_s, nmax(m))
      ENDDO

! Max number of species over all phases
      dimension_n_all = max(dimension_n_g, dimension_n_s)

! Size of MxM upper triangular matrix.
      dimension_lm = (dimension_m * (dimension_m-1)/2)+1

! Number of scalar equations.
      dimension_scalar = nscalar
      dim_scalar2 = 2*nscalar



! Set DIMENSION_x variables.
      dimension_i = imax3
      dimension_j = jmax3
      dimension_k = kmax3

      dimension_3 = (kend3-kstart3+1)*(jend3-jstart3+1)*(iend3-istart3+1)
      dimension_4 = (kend4-kstart4+1)*(jend4-jstart4+1)*(iend4-istart4+1)

      dimension_3g = ijkmax3            ! Global IJK array
      dimension_3l = ijksize3_all(mype) ! Local IJK array
      dimension_3p = merge(dimension_3, 1, .NOT.bdoing_postmfix)


      RETURN
      END SUBROUTINE set_parameters