PLAG_AllocateMemoryPost.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: allocate memory for main variables associated with the Lagrangian 
!          particles (PLAG) for all active regions on current processor.
!
! Description: none.
!
! Input: region = current region
!
! Output: region%plag = plag variables
!
! Notes: none.
!
!******************************************************************************
!
! $Id: PLAG_AllocateMemoryPost.F90,v 1.5 2008/12/06 08:45:07 mtcampbe Exp $
!
! Copyright: (c) 2003 by the University of Illinois
!
!******************************************************************************

SUBROUTINE plag_allocatememorypost( region, iReg )

  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
  IMPLICIT NONE

! ... parameters
  TYPE(t_region) :: region

  INTEGER, INTENT(IN) :: ireg

! ... loop variables
  INTEGER :: icont, ilev

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

  INTEGER :: errorflag, kdnend, naiv, narv, ncont, ncv, ndv, npclsmax, ntv

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

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

  IF ( global%myProcid == masterproc .AND. & 
       global%verbLevel > verbose_none ) THEN 
    WRITE(stdout,'(A,3X,A)') solver_name,'Allocating 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

! - Get particle dimensions ---------------------------------------------------

    npclsmax = region%plagInput%nPclsMax
    ncont    = region%plagInput%nCont

    IF ( ireg == 1 ) &
      print*,' PLAG_AllocateMemoryPost: nPclsMax = ',npclsmax

! - Set pointers --------------------------------------------------------------
        
    pplag%nCv  = cv_plag_last + ncont
    pplag%nDv  = dv_plag_last + ncont
    pplag%nTv  = tv_plag_last
    pplag%nAiv = aiv_plag_last
    pplag%nArv = arv_plag_last
    
    naiv = pplag%nAiv
    narv = pplag%nArv
            
    ncv  = pplag%nCv
    ndv  = pplag%nDv
    ntv  = pplag%nTv

    IF ( global%myProcid == masterproc .AND. & 
         global%verbLevel > verbose_none ) THEN 
      WRITE(stdout,1015) '   nPclsMax',npclsmax
      WRITE(stdout,1010) '   nCont   ',ncont
      WRITE(stdout,1010) '   nAiv    ',naiv
      WRITE(stdout,1010) '   nArv    ',narv
      WRITE(stdout,1010) '   nCv     ',ncv
      WRITE(stdout,1010) '   nDv     ',ndv
      WRITE(stdout,1010) '   nTv     ',ntv
    END IF ! global%verbLevel
     
! - Lagrangian particles variables --------------------------------------------

    ALLOCATE( pplag%aiv(naiv,npclsmax),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%aiv' ) 
    END IF ! global%error

    ALLOCATE( pplag%arv(narv,npclsmax),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%arv' ) 
    END IF ! global%error
    
    ALLOCATE( pplag%cv (ncv,npclsmax),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%cv' ) 
    END IF ! global%error

    ALLOCATE( pplag%dv (ndv,npclsmax),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%dv' ) 
    END IF ! global%error
    
    ALLOCATE( pplag%tv (ntv,npclsmax),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%tv' ) 
    END IF ! global%error
    
! - Lagrangian particles mass and volume indices ------------------------------
          
    ALLOCATE( pplag%cvPlagMass(ncont),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%cvPlagMass' )
    END IF ! global%error
          
    ALLOCATE( pplag%dvPlagVolu(ncont),stat=errorflag )
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop( global, err_allocate,__line__ ,'pPlag%dvPlagVolu' )
    END IF ! global%error
 
    DO icont = 1, ncont
      pplag%cvPlagMass(icont) = cv_plag_last  +icont
      pplag%dvPlagVolu(icont) = dv_plag_last  +icont
    END DO ! iCont

  ENDDO   ! iLev

! finalize

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

  CALL deregisterfunction( global )
  
1010 FORMAT(a,' = ',i2)
1015 FORMAT(a,' = ',i7)

END SUBROUTINE plag_allocatememorypost

!******************************************************************************
!
! RCS Revision history:
!
! $Log: PLAG_AllocateMemoryPost.F90,v $
! Revision 1.5  2008/12/06 08:45:07  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:18:17  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2007/03/06 23:27:44  fnajjar
! Renamed nPclsTot to nPclsMax
!
! Revision 1.2  2005/02/16 14:50:19  fnajjar
! Cosmetic change to include verbosity at end of routine
!
! Revision 1.1  2004/12/01 22:00:46  fnajjar
! Initial revision after changing case
!
! Revision 1.1.1.1  2003/05/06 16:14:38  fnajjar
! Import of postprocessing tool for Rocpart
!
!******************************************************************************