RFLU_ModDeallocateMemory.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: Suite of routines to deallocate memory.
!
! Description: None.
!
! Notes: None.
!
!*******************************************************************************
!
! $Id: RFLU_ModDeallocateMemory.F90,v 1.11 2008/12/06 08:44:21 mtcampbe Exp $
!
! Copyright: (c) 2004 by the University of Illinois
!
!*******************************************************************************

MODULE rflu_moddeallocatememory

  USE moddatatypes
  USE moderror
  USE modglobal, ONLY: t_global
  USE modparameters
  USE modgrid, ONLY: t_grid
  USE modbndpatch, ONLY: t_patch
  USE moddatastruct, ONLY: t_region
  USE modmixture, ONLY: t_mixt_input
  USE modmpi

  USE modinterfaces, ONLY: rflu_decideneedbgradface

  IMPLICIT NONE

  PRIVATE
  PUBLIC :: rflu_deallocatememorygspeeds, &
            rflu_deallocatememorysol, &
            rflu_deallocatememorysolcv, &
            rflu_deallocatememorysoldv, &
            rflu_deallocatememorysolgv, &
            rflu_deallocatememorysoltv, &
            rflu_deallocatememorytstep, & 
            rflu_deallocatememorytstep_c, & 
            rflu_deallocatememorytstep_i

! ******************************************************************************
! Declarations and definitions
! ******************************************************************************

  CHARACTER(CHRLEN) :: RCSIdentString = &
    '$RCSfile: RFLU_ModDeallocateMemory.F90,v $ $Revision: 1.11 $'

! ******************************************************************************
! Routines
! ******************************************************************************

  CONTAINS





! ******************************************************************************
!
! Purpose: Deallocate memory for grid speeds.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorygspeeds(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag,ipatch
  TYPE(t_grid), POINTER :: pgrid
  TYPE(t_mixt_input), POINTER :: pmixtinput
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemoryGSpeeds',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pgrid      => pregion%grid
  pmixtinput => pregion%mixtInput

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

! ==============================================================================
! Grid motion active
! ==============================================================================

  IF ( pmixtinput%moveGrid .EQV. .true. ) THEN

! ------------------------------------------------------------------------------
!   Interior faces
! ------------------------------------------------------------------------------

    DEALLOCATE(pgrid%gs,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pGrid%gs')
    END IF ! global%error

! ------------------------------------------------------------------------------
!   Patch faces
! ------------------------------------------------------------------------------

    IF ( pgrid%nPatches > 0 ) THEN
      DO ipatch = 1,pgrid%nPatches
        ppatch => pregion%patches(ipatch)

        DEALLOCATE(ppatch%gs,stat=errorflag)
        global%error = errorflag
        IF ( global%error /= err_none ) THEN
          CALL errorstop(global,err_deallocate,__line__,'pPatch%gs')
        END IF ! global%error
      END DO ! iPatch
    END IF ! pGrid%nPatches

! ==============================================================================
! Grid motion not active
! ==============================================================================

  ELSE

! ------------------------------------------------------------------------------
!   Interior faces
! ------------------------------------------------------------------------------

    DEALLOCATE(pgrid%gs,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pGrid%gs')
    END IF ! global%error

! ------------------------------------------------------------------------------
!   Patch faces
! ------------------------------------------------------------------------------

    IF ( pgrid%nPatches > 0 ) THEN
      DO ipatch = 1,pgrid%nPatches
        ppatch => pregion%patches(ipatch)

        DEALLOCATE(ppatch%gs,stat=errorflag)
        global%error = errorflag
        IF ( global%error /= err_none ) THEN
          CALL errorstop(global,err_deallocate,__line__,'pPatch%gs')
        END IF ! global%error
      END DO ! iPatch
    END IF ! pGrid%nPatches
  END IF ! pMixtInput%moveGrid

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorygspeeds






! ******************************************************************************
!
! Purpose: Deallocate memory for mixture solution.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorysol(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  CHARACTER(CHRLEN) :: rcsidentstring
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemorySol',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Allocate memory
! ******************************************************************************

  CALL rflu_deallocatememorysolcv(pregion)
  CALL rflu_deallocatememorysoldv(pregion)
  CALL rflu_deallocatememorysolgv(pregion)
  CALL rflu_deallocatememorysoltv(pregion)

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorysol






! ******************************************************************************
!
! Purpose: Deallocate memory for conserved variables of mixture.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorysolcv(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  CHARACTER(CHRLEN) :: rcsidentstring
  INTEGER :: errorflag
  TYPE(t_global), POINTER :: global
  TYPE(t_grid), POINTER :: pgrid

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemorySolCv',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pgrid => pregion%grid

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

  DEALLOCATE(pregion%mixt%cv,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%cv')
  END IF ! global%error

  DEALLOCATE(pregion%mixt%cvInfo,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%cvInfo')
  END IF ! global%error

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorysolcv






! ******************************************************************************
!
! Purpose: Deallocate memory for dependent variables of mixture.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorysoldv(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  CHARACTER(CHRLEN) :: rcsidentstring
  INTEGER :: errorflag
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemorySolDv',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

  DEALLOCATE(pregion%mixt%dv,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%dv')
  END IF ! global%error

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorysoldv






! ******************************************************************************
!
! Purpose: Deallocate memory for gas variables of mixture.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorysolgv(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemorySolGv',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

  DEALLOCATE(pregion%mixt%gv,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%gv')
  END IF ! global%error

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorysolgv





! ******************************************************************************
!
! Purpose: Deallocate memory for transport variables of mixture.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorysoltv(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag
  TYPE(t_mixt_input), POINTER :: pmixtinput
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemorySolTv',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pmixtinput => pregion%mixtInput

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

  IF ( pmixtinput%computeTv .EQV. .true. ) THEN
    DEALLOCATE(pregion%mixt%tv,stat=errorflag)
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%tv')
    END IF ! global%error
  END IF ! pMixtInput%computeTv

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorysoltv






! ******************************************************************************
!
! Purpose: Deallocate memory for mixture time-stepping.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorytstep(pRegion)

  USE rflu_modoles

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag,ipatch,nbfaces,nbfacestot
  TYPE(t_grid), POINTER :: pgrid,pgridold,pgridold2
  TYPE(t_mixt_input), POINTER :: pmixtinput
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemoryTStep',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pgrid      => pregion%grid
  pgridold   => pregion%gridOld
  pgridold2  => pregion%gridOld2
  pmixtinput => pregion%mixtInput

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

! ==============================================================================
! Old solutions
! ==============================================================================

  DEALLOCATE(pregion%mixt%cvOld,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%cvOld')
  END IF ! global%error

  IF ( global%solverType == solv_implicit_nk ) THEN 
    DEALLOCATE(pregion%mixt%cvOld1,stat=errorflag)
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%cvOld1')
    END IF ! global%error

    DEALLOCATE(pregion%mixt%cvOld2,stat=errorflag)
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%cvOld2')
    END IF ! global%error
  END IF ! global%solverType

! ==============================================================================
! Time step
! ==============================================================================

  DEALLOCATE(pregion%dt,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%dt')
  END IF ! global%error

! ==============================================================================
! Residuals
! ==============================================================================

  DEALLOCATE(pregion%mixt%rhs,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%rhs')
  END IF ! global%error

  DEALLOCATE(pregion%mixt%diss,stat=errorflag)
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%diss')
  END IF ! global%error

  IF ( global%flowType == flow_unsteady ) THEN
    DEALLOCATE(pregion%mixt%rhsSum,stat=errorflag)
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%rhsSum')
    END IF ! global%error
  END IF ! global%flowType

! ==============================================================================
! Gradients
! ==============================================================================

! ------------------------------------------------------------------------------
! Cell gradients
! ------------------------------------------------------------------------------

  IF ( (pmixtinput%spaceDiscr == discr_upw_roe     ) .OR. &
       (pmixtinput%spaceDiscr == discr_upw_hllc    ) .OR. & 
       (pmixtinput%spaceDiscr == discr_upw_ausmplus) ) THEN
    IF ( pmixtinput%spaceOrder > 1 ) THEN
      DEALLOCATE(pregion%mixt%gradCell,stat=errorflag)
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%gradCell')
      END IF ! global%error
    END IF ! pMixtInput%spaceOrder
  ELSE IF ( pmixtinput%spaceDiscr == discr_opt_les ) THEN
    DEALLOCATE(pregion%mixt%gradCell,stat=errorflag)
    global%error = errorflag
    IF (global%error /= err_none) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%gradCell')
    END IF ! global%error
  END IF ! pMixtInput%spaceDiscr

! ------------------------------------------------------------------------------
! Face gradients
! ------------------------------------------------------------------------------

  IF ( pmixtinput%flowModel == flow_navst ) THEN
    DEALLOCATE(pregion%mixt%gradFace,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%gradFace')
    END IF ! global%error
  END IF ! pMixtInput%flowModel

  IF ( pgrid%nPatches > 0 ) THEN
    DO ipatch = 1,pregion%grid%nPatches
      ppatch => pregion%patches(ipatch)
  
      IF ( rflu_decideneedbgradface(pregion,ppatch) .EQV. .true. ) THEN
        DEALLOCATE(ppatch%mixt%gradFace,stat=errorflag)
        global%error = errorflag
        IF ( global%error /= err_none ) THEN
          CALL errorstop(global,err_deallocate,__line__,'pPatch%mixt%gradFace')
        END IF ! global%error
      END IF ! RFLU_DecideNeedBGradFace
    END DO ! iPatch
  END IF ! pGrid%nPatches

! ==============================================================================
! Grid motion. NOTE grid speeds are allocated separately because they are
! written into grid file, and hence they need to be allocated in pre- and
! postprocessors also.
! ==============================================================================

  IF ( pmixtinput%moveGrid .EQV. .true. ) THEN

! ------------------------------------------------------------------------------
!   Residual
! ------------------------------------------------------------------------------

    DEALLOCATE(pgrid%rhs,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pGrid%rhs')
    END IF ! global%error

! ------------------------------------------------------------------------------
!   Displacement
! ------------------------------------------------------------------------------

    IF ( pmixtinput%moveGridType /= movegrid_type_xyz ) THEN
      DEALLOCATE(pgrid%disp,stat=errorflag)
      global%error = errorflag
      IF ( global%error /= err_none ) THEN
        CALL errorstop(global,err_deallocate,__line__,'pGrid%disp')
      END IF ! global%error
    END IF ! pMixtInput%moveGridType

! ------------------------------------------------------------------------------
!   Old coordinates
! ------------------------------------------------------------------------------

    DEALLOCATE(pgridold%xyz,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%gridOld%xyz')
    END IF ! global%error

    IF ( global%solverType == solv_implicit_nk ) THEN 
      DEALLOCATE(pgridold2%xyz,stat=errorflag)
      global%error = errorflag
      IF ( global%error /= err_none ) THEN
        CALL errorstop(global,err_deallocate,__line__,'pRegion%gridOld2%xyz')
      END IF ! global%error
    END IF ! global%solverType

! ------------------------------------------------------------------------------
!   Old volume
! ------------------------------------------------------------------------------

    DEALLOCATE(pgridold%vol,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%gridOld%vol')
    END IF ! global%error

    IF ( global%solverType == solv_implicit_nk ) THEN 
      DEALLOCATE(pgridold2%vol,stat=errorflag)
      global%error = errorflag
      IF ( global%error /= err_none ) THEN
        CALL errorstop(global,err_deallocate,__line__,'pRegion%gridOld2%vol')
      END IF ! global%error
    END IF ! global%solverType

#ifndef GENX
! ------------------------------------------------------------------------------
!   Patch displacements. NOTE allocate here only if not running inside GENX,
!   because when running inside GENX allocate displacements also for virtual
!   vertices.
! ------------------------------------------------------------------------------

    IF ( pgrid%nPatches > 0 ) THEN
      DO ipatch = 1,pgrid%nPatches
        ppatch => pregion%patches(ipatch)

        DEALLOCATE(ppatch%dXyz,stat=errorflag)
        global%error = errorflag
        IF ( global%error /= err_none ) THEN
          CALL errorstop(global,err_deallocate,__line__,'pPatch%dXyz')
        END IF ! global%error

      END DO ! iPatch
    END IF ! pGrid%nPatches
#endif
  END IF ! pMixtInput%moveGrid

#ifdef STATS
! ==============================================================================
! Time averaged statistics
! ==============================================================================

  IF ( (global%flowType == flow_unsteady) .AND. &
       (global%doStat == active) ) THEN
    DEALLOCATE(pregion%mixt%tav,stat=errorflag)
    global%error = errorflag
    IF ( global%error /= err_none ) THEN
      CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%tav')
    END IF ! global%error
  END IF ! global%flowType
#endif

! ==============================================================================
! Optimal LES
! ==============================================================================

  IF ( pmixtinput%spaceDiscr == discr_opt_les ) THEN
    CALL rflu_destroystencilsweightsoles(pregion)
  END IF ! pMixtInput%spaceDiscr

! ==============================================================================
! Substantial derivative
! ==============================================================================

  DEALLOCATE(pregion%mixt%sd,stat=errorflag)
  global%error = errorflag
  IF ( global%error /= err_none ) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%sd')
  END IF ! global%error

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorytstep






! ******************************************************************************
!
! Purpose: Deallocate memory for mixture time-stepping for compressible fluid.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorytstep_c(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag,ipatch
  TYPE(t_grid), POINTER :: pgrid
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemoryTStep_C',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pgrid => pregion%grid

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

! ==============================================================================
! Mass fluxes
! ==============================================================================

  DEALLOCATE(pregion%mixt%mfMixt,stat=errorflag)
  global%error = errorflag
  IF ( global%error /= err_none ) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%mfMixt')
  END IF ! global%error

  IF ( pregion%grid%nPatches > 0 ) THEN
    DO ipatch = 1,pregion%grid%nPatches
      ppatch => pregion%patches(ipatch)

      DEALLOCATE(ppatch%mfMixt,stat=errorflag)
      global%error = errorflag
      IF ( global%error /= err_none ) THEN
        CALL errorstop(global,err_deallocate,__line__,'pPatch%mfMixt')
      END IF ! global%error
    END DO ! iPatch
  END IF ! pRegion%grid%nPatches

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorytstep_c






! ******************************************************************************
!
! Purpose: Deallocate memory for mixture time-stepping for incompressible fluid.
!
! Description: None.
!
! Input:
!   pRegion        Region pointer
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************

SUBROUTINE rflu_deallocatememorytstep_i(pRegion)

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: errorflag,ipatch
  TYPE(t_grid), POINTER :: pgrid
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => pregion%global

  CALL registerfunction(global,'RFLU_DeallocateMemoryTStep_I',&
  'RFLU_ModDeallocateMemory.F90')

! ******************************************************************************
! Set pointers
! ******************************************************************************

  pgrid => pregion%grid

! ******************************************************************************
! Deallocate memory
! ******************************************************************************

! ==============================================================================
! Face velocities
! ==============================================================================

  DEALLOCATE(pregion%mixt%vfMixt,stat=errorflag)
  global%error = errorflag
  IF ( global%error /= err_none ) THEN
    CALL errorstop(global,err_deallocate,__line__,'pRegion%mixt%vfMixt')
  END IF ! global%error

  IF ( pregion%grid%nPatches > 0 ) THEN
    DO ipatch = 1,pregion%grid%nPatches
      ppatch => pregion%patches(ipatch)

      DEALLOCATE(ppatch%vfMixt,stat=errorflag)
      global%error = errorflag
      IF ( global%error /= err_none ) THEN
        CALL errorstop(global,err_deallocate,__line__,'pPatch%vfMixt')
      END IF ! global%error
    END DO ! iPatch
  END IF ! pRegion%grid%nPatches

! ******************************************************************************
! End
! ******************************************************************************

  CALL deregisterfunction(global)

END SUBROUTINE rflu_deallocatememorytstep_i





! ******************************************************************************
! End
! ******************************************************************************

END MODULE rflu_moddeallocatememory


! ******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLU_ModDeallocateMemory.F90,v $
! Revision 1.11  2008/12/06 08:44:21  mtcampbe
! Updated license.
!
! Revision 1.10  2008/11/19 22:17:32  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.9  2006/11/01 15:50:00  haselbac
! Changed so implicit-solver arrays dealt with properly
!
! Revision 1.8  2006/08/19 15:39:03  mparmar
! Moved region%mixt%bGradFace to patch%mixt%gradFace
!
! Revision 1.7  2006/02/08 21:03:47  hdewey2
! Added old2 quantities
!
! Revision 1.6  2005/09/22 17:11:04  hdewey2
! Added deallocation of cvOld1 and cvOld2 for transient implicit solver.
!
! Revision 1.5  2005/07/14 21:43:26  haselbac
! Added AUSM flux function to memory deallocation IF statement
!
! Revision 1.4  2004/12/19 15:47:08  haselbac
! Added memory deallocation for incompressible solver
!
! Revision 1.3  2004/10/19 19:38:48  haselbac
! Updated for GEN3
!
! Revision 1.2  2004/07/30 22:47:36  jferry
! Implemented Equilibrium Eulerian method for Rocflu
!
! Revision 1.1  2004/03/19 21:15:19  haselbac
! Initial revision
!
! ******************************************************************************