RFLO_AllocateMemory.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 all variables associated with the mixture
!          for all active regions on current processor.
!
! Description: none.
!
! Input: region = current region
!
! Output: region%mixt = mixture variables
!         region%grid = grid variables (also region%gridOld)
!
! Notes: none.
!
!******************************************************************************
!
! $Id: RFLO_AllocateMemory.F90,v 1.24 2008/12/06 08:44:25 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE rflo_allocatememory( region )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region, t_level
  USE modbndpatch, ONLY   : t_patch
  USE modglobal, ONLY     : t_global
  USE modgrid, ONLY       : t_grid
  USE modinterfaces, ONLY : rflo_getdimensdummynodes, rflo_getnodeoffset, &
        rflo_getdimensdummy, rflo_getcelloffset, rflo_getdimensphysnodes
  USE modmixture, ONLY    : t_mixt
  USE moderror
  USE modparameters
  IMPLICIT NONE

#include "Indexing.h"

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

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

! ... local variables
  INTEGER :: idcbeg, idcend, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: idnbeg, idnend, jdnbeg, jdnend, kdnbeg, kdnend
  INTEGER :: ipnbeg, ipnend, jpnbeg, jpnend, kpnbeg, kpnend
  INTEGER :: ibc, iec, ibn, ien, icelloffset, ijcelloffset, &
             inodeoffset, ijnodeoffset, h1, h2, errorflag
  INTEGER :: n1, n2, ijbeg, ijend, ioff

  TYPE(t_level), POINTER  :: level
  TYPE(t_grid) , POINTER  :: grid, gridold, gridold2
  TYPE(t_mixt) , POINTER  :: mixt
  TYPE(t_global), POINTER :: global
  TYPE(t_patch) , POINTER :: patch

!******************************************************************************

  global => region%global

  CALL registerfunction( global,'RFLO_AllocateMemory',&
  'RFLO_AllocateMemory.F90' )

! loop over all grid levels

  DO ilev=1,region%nGridLevels

    level    => region%levels(ilev)
    grid     => region%levels(ilev)%grid
    gridold  => region%levels(ilev)%gridOld
    gridold2 => region%levels(ilev)%gridOld2
    mixt     => region%levels(ilev)%mixt

! - get cell and node dimensions

    CALL rflo_getdimensdummy( region,ilev,idcbeg,idcend, &
                              jdcbeg,jdcend,kdcbeg,kdcend )
    CALL rflo_getcelloffset( region,ilev,icelloffset,ijcelloffset )
    ibc = indijk(idcbeg,jdcbeg,kdcbeg,icelloffset,ijcelloffset)
    iec = indijk(idcend,jdcend,kdcend,icelloffset,ijcelloffset)

    CALL rflo_getdimensdummynodes( region,ilev,idnbeg,idnend, &
                                   jdnbeg,jdnend,kdnbeg,kdnend )
    CALL rflo_getnodeoffset( region,ilev,inodeoffset,ijnodeoffset )
    ibn = indijk(idnbeg,jdnbeg,kdnbeg,inodeoffset,ijnodeoffset)
    ien = indijk(idnend,jdnend,kdnend,inodeoffset,ijnodeoffset)

    CALL rflo_getdimensphysnodes( region,ilev,ipnbeg,ipnend, &
                                  jpnbeg,jpnend,kpnbeg,kpnend )

! - grid coordinates

    ALLOCATE( grid%xyz(3,ibn:ien)  ,stat=errorflag )
    ALLOCATE( grid%ijkDgen(ibn:ien),stat=errorflag )
    IF (region%mixtInput%moveGrid) THEN
      ALLOCATE( gridold%xyzOld(3,ibn:ien),stat=errorflag )
      ALLOCATE( gridold%xyz(3,ibn:ien)   ,stat=errorflag )
      IF (ilev == 1) THEN     ! only on the finest grid
        ALLOCATE( grid%arcLen12(jpnbeg:jpnend,kpnbeg:kpnend),stat=errorflag )
        ALLOCATE( grid%arcLen34(kpnbeg:kpnend,ipnbeg:ipnend),stat=errorflag )
        ALLOCATE( grid%arcLen56(ipnbeg:ipnend,jpnbeg:jpnend),stat=errorflag )
      ENDIF
      IF (global%moveGridScheme/=movegrid_blocks .AND. ilev==1) THEN
        ALLOCATE( grid%xyzOld(3,ibn:ien),stat=errorflag )
      ENDIF
      IF ((global%moveGridScheme==movegrid_frame .OR. &
           global%moveGridScheme==movegrid_foms  .OR. &
           global%moveGridScheme==movegrid_elframe) .AND. ilev==1) THEN
        ALLOCATE( grid%nCorns(global%nRegions), stat=errorflag )
      ENDIF
      IF (global%moveGridScheme==movegrid_foms .AND. ilev==1) THEN
        ALLOCATE( grid%xyzOrth(3,ibn:ien), stat=errorflag )
        ALLOCATE( grid%xyzTemp(3,ibn:ien), stat=errorflag )
      ENDIF
      IF ((global%moveGridScheme==movegrid_elglobal .OR. &
           global%moveGridScheme==movegrid_elframe) .AND. ilev==1) THEN
        ALLOCATE( grid%stu(      3,ibn:ien),stat=errorflag )
        ALLOCATE( grid%stuOld(   3,ibn:ien),stat=errorflag )
        ALLOCATE( grid%stui(     3,ibn:ien),stat=errorflag )
        ALLOCATE( grid%stuj(     3,ibn:ien),stat=errorflag )
        ALLOCATE( grid%stuk(     3,ibn:ien),stat=errorflag )
        ALLOCATE( grid%pmat(   3,6,ibn:ien),stat=errorflag )

        ALLOCATE( grid%aijk(       ibn:ien),stat=errorflag )
        ALLOCATE( grid%aimjk(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aipjk(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijmk(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijpk(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijkm(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijkp(      ibn:ien),stat=errorflag )
        ALLOCATE( grid%aimjmk(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aipjmk(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aimjpk(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aipjpk(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aimjkm(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aipjkm(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aimjkp(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aipjkp(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijmkm(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijpkm(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijmkp(     ibn:ien),stat=errorflag )
        ALLOCATE( grid%aijpkp(     ibn:ien),stat=errorflag )

        DO ipatch=1,region%nPatches
          patch  => region%levels(ilev)%patches(ipatch)

          h1 = abs(patch%l1end-patch%l1beg) + 2
          h2 = abs(patch%l2end-patch%l2beg) + 2

! ------- elliptic PDE grid motion

          ALLOCATE( patch%st(     3,h1,h2),stat=errorflag )
          ALLOCATE( patch%stOld(  3,h1,h2),stat=errorflag )
          ALLOCATE( patch%sti(    3,h1,h2),stat=errorflag )
          ALLOCATE( patch%stj(    3,h1,h2),stat=errorflag )
          ALLOCATE( patch%stii(   3,h1,h2),stat=errorflag )
          ALLOCATE( patch%stjj(   3,h1,h2),stat=errorflag )
          ALLOCATE( patch%stij(   3,h1,h2),stat=errorflag )
          ALLOCATE( patch%pfun( 2,3,h1,h2),stat=errorflag )

          ALLOCATE( patch%aimjm( h1,h2),stat=errorflag )
          ALLOCATE( patch%aijm(  h1,h2),stat=errorflag )
          ALLOCATE( patch%aipjm( h1,h2),stat=errorflag )
          ALLOCATE( patch%aimj(  h1,h2),stat=errorflag )
          ALLOCATE( patch%aij(   h1,h2),stat=errorflag )
          ALLOCATE( patch%aipj(  h1,h2),stat=errorflag )
          ALLOCATE( patch%aimjp( h1,h2),stat=errorflag )
          ALLOCATE( patch%aijp(  h1,h2),stat=errorflag )
          ALLOCATE( patch%aipjp( h1,h2),stat=errorflag )

        ENDDO      ! iPatch
      ENDIF        ! moveGridScheme
    ENDIF          ! moveGrid
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

! - averaging coefficients, face vectors, volumes, grid speeds

    ALLOCATE( grid%c2fCoI(2,ibn:ien),stat=errorflag )
    ALLOCATE( grid%c2fCoJ(2,ibn:ien),stat=errorflag )
    ALLOCATE( grid%c2fCoK(2,ibn:ien),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    ALLOCATE( grid%c2eCoI(4,ibn:ien),stat=errorflag )
    ALLOCATE( grid%c2eCoJ(4,ibn:ien),stat=errorflag )
    ALLOCATE( grid%c2eCoK(4,ibn:ien),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    ALLOCATE( grid%si(4,ibn:ien),stat=errorflag )
    ALLOCATE( grid%sj(4,ibn:ien),stat=errorflag )
    ALLOCATE( grid%sk(4,ibn:ien),stat=errorflag )
    ALLOCATE( grid%vol(ibc:iec) ,stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    IF ((region%mixtInput%spaceOrder == discr_order_4) .OR. &
        global%calcCellCtr) THEN
      ALLOCATE( grid%cofg(3,ibc:iec),stat=errorflag )
      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )
    ELSE
      nullify( grid%cofg )
    ENDIF

    IF (global%calcFaceCtr) THEN
      ALLOCATE( grid%cfcI(3,ibn:ien),stat=errorflag )
      ALLOCATE( grid%cfcJ(3,ibn:ien),stat=errorflag )
      ALLOCATE( grid%cfcK(3,ibn:ien),stat=errorflag )
      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )
    ELSE
      nullify( grid%cfcI, grid%cfcJ, grid%cfcK )
    ENDIF

    IF (region%mixtInput%moveGrid) THEN
      ALLOCATE( grid%siVel(ibn:ien)  ,stat=errorflag )
      ALLOCATE( grid%sjVel(ibn:ien)  ,stat=errorflag )
      ALLOCATE( grid%skVel(ibn:ien)  ,stat=errorflag )
      ALLOCATE( gridold%si(4,ibn:ien),stat=errorflag )
      ALLOCATE( gridold%sj(4,ibn:ien),stat=errorflag )
      ALLOCATE( gridold%sk(4,ibn:ien),stat=errorflag )
      ALLOCATE( gridold%vol(ibc:iec) ,stat=errorflag )
      IF (global%solverType==solv_implicit .AND. &
          global%tstepOrder==3 .AND. ilev == 1) THEN
        ALLOCATE( gridold2%vol(ibc:iec),stat=errorflag )
      ENDIF
    ELSE
      ALLOCATE( grid%siVel(0:1),stat=errorflag )
      ALLOCATE( grid%sjVel(0:1),stat=errorflag )
      ALLOCATE( grid%skVel(0:1),stat=errorflag )
    ENDIF
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    grid%siVel(:) = 0._rfreal   ! zero out grid speeds
    grid%sjVel(:) = 0._rfreal
    grid%skVel(:) = 0._rfreal

! - time step & spectral radii

    ALLOCATE( level%dt(ibc:iec)   ,stat=errorflag )
    ALLOCATE( mixt%srad(3,ibc:iec),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

! - mixture variables

    ALLOCATE( mixt%cv   (cv_mixt_neqs,ibc:iec),stat=errorflag )
    ALLOCATE( mixt%cvOld(cv_mixt_neqs,ibc:iec),stat=errorflag )
    ALLOCATE( mixt%dv   (mixt%nDv    ,ibc:iec),stat=errorflag )
    ALLOCATE( mixt%rhs  (cv_mixt_neqs,ibc:iec),stat=errorflag )
    ALLOCATE( mixt%diss (cv_mixt_neqs,ibc:iec),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    IF (mixt%nTv > 0) THEN
      ALLOCATE( mixt%tv(mixt%nTv,ibc:iec),stat=errorflag )
    ELSE
      nullify( mixt%tv )
    ENDIF

    IF (mixt%nGrad > 0) THEN
      ALLOCATE( mixt%gradi(mixt%nGrad,ibn:ien),stat=errorflag )
      ALLOCATE( mixt%gradj(mixt%nGrad,ibn:ien),stat=errorflag )
      ALLOCATE( mixt%gradk(mixt%nGrad,ibn:ien),stat=errorflag )
    ELSE
      nullify( mixt%gradi )
      nullify( mixt%gradj )
      nullify( mixt%gradk )
    ENDIF

    IF (region%mixtInput%gasModel == gas_model_tcperf) THEN
      ALLOCATE( mixt%gv(mixt%nGv,0:1),stat=errorflag )
    ELSE
      ALLOCATE( mixt%gv(mixt%nGv,ibc:iec),stat=errorflag )
    ENDIF

#ifdef STATS
    IF ((global%flowType == flow_unsteady) .AND. &
        (global%doStat == active) .AND. &
        (global%mixtNStat > 0)) THEN
      ALLOCATE( mixt%tav(global%mixtNStat,ibc:iec),stat=errorflag )
    ELSE
      nullify( mixt%tav )
    ENDIF
#endif

    IF (ilev>1 .AND. global%cycleType/=mgcycle_no) THEN
      ALLOCATE( mixt%fterm(cv_mixt_neqs,ibc:iec),stat=errorflag )
    ELSE
      nullify( mixt%fterm )
    ENDIF

    IF (global%flowType == flow_unsteady) THEN
      ALLOCATE( mixt%rhsSum(cv_mixt_neqs,ibc:iec),stat=errorflag )
    ELSE
      nullify( mixt%rhsSum )
    ENDIF

    IF (global%flowType==flow_unsteady .AND. &
        global%solverType==solv_implicit .AND. &
        ilev == 1) THEN
      ALLOCATE( mixt%cvn  (cv_mixt_neqs,ibc:iec),stat=errorflag )
      ALLOCATE( mixt%cvn1 (cv_mixt_neqs,ibc:iec),stat=errorflag )
      ALLOCATE( mixt%cvn2 (cv_mixt_neqs,ibc:iec),stat=errorflag )
      ALLOCATE( mixt%sDual(cv_mixt_neqs,ibc:iec),stat=errorflag )
    ELSE
      nullify( mixt%cvn   )
      nullify( mixt%cvn1  )
      nullify( mixt%cvn2  )
      nullify( mixt%sDual )
    ENDIF

    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

! - patch data allocation -----------------------------------------------------

    DO ipatch=1,region%nPatches
      patch  => region%levels(ilev)%patches(ipatch)

! --- aerodynamic coefficients

      IF (global%aeroCoeffs == active) THEN
        n1    = abs(patch%l1end-patch%l1beg)
        n2    = abs(patch%l2end-patch%l2beg)
        ioff  = n1 + 1
        ijbeg = indij( 0, 0,ioff)
        ijend = indij(n1,n2,ioff)
      ELSE
        ijbeg = 0
        ijend = 1
      ENDIF
      ALLOCATE( patch%cp(ijbeg:ijend),   stat=errorflag )
      ALLOCATE( patch%ch(ijbeg:ijend),   stat=errorflag )
      ALLOCATE( patch%cf(3,ijbeg:ijend), stat=errorflag )
      ALLOCATE( patch%forceCoeffs(3,forces_press:forces_visc), stat=errorflag )
      ALLOCATE( patch%momentCoeffs(3,forces_press:forces_visc),stat=errorflag )

      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

! --- patch arclengths

      h1 = abs(patch%l1end-patch%l1beg) + 2
      h2 = abs(patch%l2end-patch%l2beg) + 2
      ALLOCATE( patch%arclen1( h2 ), stat=errorflag )
      ALLOCATE( patch%arclen2( h1 ), stat=errorflag )

      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )

    ENDDO ! iPatch

  ENDDO   ! iLev

! data allocation independent of grid level -----------------------------------

! coefficients of implicit residual smoothing

  IF ((global%flowType==flow_steady .OR. &
      (global%flowType==flow_unsteady .AND. global%solverType==solv_implicit)) &
      .AND. region%mixtInput%smoocf>0.) THEN
    ALLOCATE( mixt%epsIrs(3,ibc:iec),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global,err_allocate,__line__ )
  ELSE
    nullify( mixt%epsIrs )
  ENDIF

! finalize --------------------------------------------------------------------

  CALL deregisterfunction( global )

END SUBROUTINE rflo_allocatememory

!******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLO_AllocateMemory.F90,v $
! Revision 1.24  2008/12/06 08:44:25  mtcampbe
! Updated license.
!
! Revision 1.23  2008/11/19 22:17:36  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.22  2006/03/25 23:25:27  wasistho
! adjusted size of all si,j,k
!
! Revision 1.21  2006/03/24 05:03:27  wasistho
! increased size of si, sj, sk
!
! Revision 1.20  2006/03/18 08:20:58  wasistho
! allocate patch%arclen1,2
!
! Revision 1.19  2006/03/16 22:02:08  wasistho
! bug fixed aero coeffs allocation
!
! Revision 1.18  2006/03/13 03:40:40  wasistho
! allocate aero coeffs
!
! Revision 1.17  2006/03/08 23:33:23  wasistho
! made gridOld%xyzOld available for all type gm
!
! Revision 1.16  2006/03/02 01:27:17  wasistho
! split movegrid_epde to elglobal and elframe
!
! Revision 1.15  2006/02/09 00:26:23  wasistho
! allocate nCorns for EPDE too
!
! Revision 1.14  2005/12/07 08:47:53  wasistho
! allocate arrays for surface mesh motion EPDE
!
! Revision 1.13  2005/12/03 09:33:00  wasistho
! allocate arrays for movegrid EPDE
!
! Revision 1.12  2005/11/28 20:04:14  wasistho
! allocate gridOld%xyzOld
!
! Revision 1.11  2005/11/11 07:18:18  wasistho
! allocate ijkDgen
!
! Revision 1.10  2005/10/31 21:09:36  haselbac
! Changed specModel and SPEC_MODEL_NONE
!
! Revision 1.9  2005/10/28 07:37:55  wasistho
! corrected nullify cfcI,J,K
!
! Revision 1.8  2005/10/27 05:12:01  wasistho
! allocate xyzOrth
!
! Revision 1.7  2005/10/20 06:49:47  wasistho
! allocate cfcI,J,K
!
! Revision 1.6  2005/06/10 19:32:40  wasistho
! move memory allocation of movegridFrame variables to modflo/RFLO_ModMoveGridFrame
!
! Revision 1.5  2005/06/02 03:21:32  wasistho
! shuffle MoveGridVms with MoveGridFrame
!
! Revision 1.4  2005/05/29 20:45:01  wasistho
! added regCornOrig
!
! Revision 1.3  2005/05/28 08:07:40  wasistho
! allocate data for moveGridFrame
!
! Revision 1.2  2005/05/21 05:08:29  wasistho
! changed == MOVEGRID_GLOBAL to /=MOVEGRID_BLOCKS
!
! Revision 1.1  2004/11/29 20:51:38  wasistho
! lower to upper case
!
! Revision 1.33  2004/08/03 00:33:03  wasistho
! added allocation of grid%c2eConI,J,K
!
! Revision 1.32  2004/08/02 19:32:24  wasistho
! changed grid%avgCo to grid%c2fCo
!
! Revision 1.31  2004/07/30 17:25:18  wasistho
! provide cell2face averaging coefficients
!
! Revision 1.30  2003/11/20 16:40:36  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.26  2003/10/03 20:18:07  wasistho
! added wall dist. criterion for cellcentroids allocation
!
! Revision 1.25  2003/10/01 23:52:10  jblazek
! Corrected bug in moving noslip wall BC and grid speeds.
!
! Revision 1.24  2003/08/11 21:51:18  jblazek
! Added basic global grid smoothing scheme.
!
! Revision 1.23  2003/07/03 21:48:45  jblazek
! Implemented dual-time stepping.
!
! Revision 1.22  2003/05/15 02:57:03  jblazek
! Inlined index function.
!
! Revision 1.21  2003/04/23 22:51:28  olawlor
! Added TARGET attribute to region variable, since
! pointers are cached into it.
!
! Revision 1.20  2003/03/14 22:05:11  jblazek
! Improved mesh motion algorithm - node movement exchaged between blocks.
!
! Revision 1.19  2002/11/02 01:56:04  wasistho
! Added TURB statistics
!
! Revision 1.18  2002/10/12 03:20:50  jblazek
! Replaced [io]stat=global%error with local errorFlag for Rocflo.
!
! Revision 1.17  2002/09/05 17:40:20  jblazek
! Variable global moved into regions().
!
! Revision 1.16  2002/08/30 01:47:58  jblazek
! Added support for moving grids.
!
! Revision 1.15  2002/08/15 19:48:05  jblazek
! Implemented grid deformation capability.
!
! Revision 1.14  2002/06/14 21:38:45  wasistho
! Added time avg statistics
!
! Revision 1.13  2002/03/29 23:15:22  jblazek
! Corrected bug in MPI send.
!
! Revision 1.12  2002/03/22 19:26:23  jblazek
! Unused pointers are nullified.
!
! Revision 1.11  2002/03/18 23:11:32  jblazek
! Finished multiblock and MPI.
!
! Revision 1.10  2002/02/27 18:38:20  jblazek
! Changed extrapol. to dummy cells at injection boundaries and slip walls.
!
! Revision 1.9  2002/02/21 23:25:05  jblazek
! Blocks renamed as regions.
!
! Revision 1.8  2002/02/09 01:47:01  jblazek
! Added multi-probe option, residual smoothing, physical time step.
!
! Revision 1.7  2002/02/06 00:15:39  jblazek
! Improved injection BC. Added pointers to gradients.
!
! Revision 1.6  2002/01/28 23:55:22  jblazek
! Added flux computation (central scheme).
!
! Revision 1.5  2002/01/23 03:51:25  jblazek
! Added low-level time-stepping routines.
!
! Revision 1.4  2002/01/10 00:02:07  jblazek
! Added calculation of mixture properties.
!
! Revision 1.3  2002/01/08 22:09:17  jblazek
! Added calculation of face vectors and volumes.
!
! Revision 1.2  2001/12/22 00:09:39  jblazek
! Added routines to store grid and solution.
!
! Revision 1.1  2001/12/11 21:59:29  jblazek
! memory allocation added.
!
!******************************************************************************