PLAG_InitMemory.F90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
!******************************************************************************
!
! Purpose: initialize memory for all variables associated 
!          with the Lagrangian particles (PLAG).
!
! Description: none.
!
! Input: region = current region
!
! Output: region%plag = plag variables
!
! Notes: none.
!
!******************************************************************************
!
! $Id: PLAG_InitMemory.F90,v 1.5 2008/12/06 08:44:33 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE plag_initmemory( region )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region, t_level
  USE modglobal, ONLY     : t_global
  USE modpartlag, ONLY    : t_plag, t_plag_input
  USE moderror
  USE modparameters
  USE modmpi
  USE plag_modparameters
#ifdef STATS
  USE plag_modeulerian, ONLY: plag_initeulerianfield
  USE plag_modstats, ONLY: plag_initstat
#endif

  IMPLICIT NONE

! ... parameters
  TYPE(t_region) :: region
  
! ... loop variables
  INTEGER :: icont, ilev

! ... local variables
  CHARACTER(CHRLEN) :: rcsidentstring

  INTEGER :: naiv, narv, ncont, ncv, ndv, ntv
  INTEGER :: iaiv, iarv, icv, idv, ipcl, itv, npclsmax

  TYPE(t_level),  POINTER :: plevel
  TYPE(t_plag),   POINTER :: pplag
  TYPE(t_global), POINTER :: global
  
!******************************************************************************

  rcsidentstring = '$RCSfile: PLAG_InitMemory.F90,v $ $Revision: 1.5 $'
  
  global => region%global
  
  CALL registerfunction( global, 'PLAG_InitMemory',&
  'PLAG_InitMemory.F90' )

  IF ( global%myProcid == masterproc .AND. & 
       global%verbLevel > verbose_none ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Initializing memory for PLAG...'
  END IF ! global%verbLevel

! Loop over all grid levels ---------------------------------------------------

  DO ilev=1,region%nGridLevels

    plevel     => region%levels(ilev)
    pplag      => region%levels(ilev)%plag

    ncont     = region%plagInput%nCont
    npclsmax  = region%plagInput%nPclsMax
    naiv      = region%levels(ilev)%plag%nAiv
    narv      = region%levels(ilev)%plag%nArv
    ncv       = region%levels(ilev)%plag%nCv
    ndv       = region%levels(ilev)%plag%nDv
    ntv       = region%levels(ilev)%plag%nTv
    
    DO icont = 1, ncont
      WRITE(stdout,'(A,3(I2,3X))') 'iCont, cvPlagMass, dvPlagVolu ', &
                                    icont, pplag%cvPlagMass(icont),   &
                                    pplag%dvPlagVolu(icont)
    END DO ! iCont   

    pplag%nPcls       = 0
    pplag%nextIdNumber= 0
    pplag%nPclsMax    = npclsmax

    DO ipcl = 1, npclsmax 
      DO iaiv = 1, naiv
!        pPlag%aiv(iAiv,iPcl)    = 0
!        pPlag%aivOld(iAiv,iPcl) = 0
        pplag%aiv(iaiv,ipcl)    = crazy_value_int 
        pplag%aivOld(iaiv,ipcl) = crazy_value_int
      END DO ! iAiv

      DO iarv = 1, narv
!        pPlag%arv(iArv,iPcl)    = 0.0_RFREAL
!        pPlag%arvOld(iArv,iPcl) = 0.0_RFREAL
        pplag%arv(iarv,ipcl)    = -1.0e+30_rfreal
        pplag%arvOld(iarv,ipcl) = -1.0e+30_rfreal
      END DO ! iArv

      DO icv  = 1, ncv
!        pPlag%cv(iCv,iPcl)     = 0.0_RFREAL
!        pPlag%rhs(iCv,iPcl)    = 0.0_RFREAL
!        pPlag%rhsSum(iCv,iPcl) = 0.0_RFREAL 
!        pPlag%cvOld(iCv,iPcl)     = 0.0_RFREAL

        pplag%cv(icv,ipcl)     = -1.0e+30_rfreal
        pplag%rhs(icv,ipcl)    = -1.0e+30_rfreal
        pplag%rhsSum(icv,ipcl) = -1.0e+30_rfreal 
        pplag%cvOld(icv,ipcl)  = -1.0e+30_rfreal
      END DO ! iCv

      DO idv  = 1, ndv
!        pPlag%dv(iDv,iPcl)     = 0.0_RFREAL
        pplag%dv(idv,ipcl)     = -1.0e+30_rfreal
      END DO ! iDv

      DO itv  = 1, ntv
!        pPlag%tv(iTv,iPcl)     = 0.0_RFREAL
        pplag%tv(itv,ipcl)     = -1.0e+30_rfreal
      END DO ! iTv
    ENDDO ! iPcl
    
! - Initialize memory for statistics -------------------------------------------
    
#ifdef STATS  
    IF ( ( global%flowType == flow_unsteady ) .AND. &
         ( global%doStat == active )                ) THEN
      CALL plag_initeulerianfield( region, pplag )
      CALL plag_initstat( region, pplag )
    END IF ! global%flowType        
#endif
 
  ENDDO   ! iLev
  
! finalize --------------------------------------------------------------------

  CALL deregisterfunction( global )

END SUBROUTINE plag_initmemory

!******************************************************************************
!
! RCS Revision history:
!
! $Log: PLAG_InitMemory.F90,v $
! Revision 1.5  2008/12/06 08:44:33  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:46  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2007/03/06 23:13:13  fnajjar
! Renamed nPclsTot to nPclsMax
!
! Revision 1.2  2005/01/08 20:42:13  fnajjar
! Included initialization calls for PLAG statistics
!
! Revision 1.1  2004/12/01 20:57:38  fnajjar
! Initial revision after changing case
!
! Revision 1.5  2004/07/28 18:58:13  fnajjar
! Included new definition for nPclsTot from dynamic memory reallocation
!
! Revision 1.4  2004/04/09 23:11:23  fnajjar
! Made initialization within DO loop construct and added crazy values for proper trapping
!
! Revision 1.3  2003/11/21 22:43:18  fnajjar
! Removed nPclsTot and added nextIdNumber
!
! Revision 1.2  2003/05/27 19:19:57  fnajjar
! Removed distPartBurning and all pertinent LOGICAL datastructure
!
! Revision 1.1  2002/10/25 14:16:31  f-najjar
! Initial Import of Rocpart
!
!
!******************************************************************************