TURB_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 turbulence
!          for all active regions on current processor.
!
! Description: none.
!
! Input: region = info of current region data
!
! Output: region%turb = info of current turbulence data
!
! Notes: The total size of memory depends on the input selected by user.
!
!******************************************************************************
!
! $Id: TURB_AllocateMemory.F90,v 1.17 2009/08/31 05:17:35 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_allocatememory( region ) ! PUBLIC

  USE moddatatypes
  USE modbndpatch, ONLY   : t_patch
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modturbulence, ONLY : t_turb
  USE moderror
  USE modparameters
  USE turb_modparameters
#ifdef RFLO
  USE modinterfaces, ONLY : rflo_getdimensdummynodes, rflo_getnodeoffset, &
                            rflo_getdimensdummy, rflo_getcelloffset
#include "Indexing.h"
#endif
  IMPLICIT NONE

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

! ... loop variables
  INTEGER :: ipatch, ijbeg, ijend, ibc, iec, ibn, ien
#ifdef RFLO
  INTEGER :: ilev
#endif

! ... local variables
  CHARACTER(CHRLEN)       :: rcsidentstring
  TYPE(t_global), POINTER :: global
  TYPE(t_turb), POINTER   :: turb
  TYPE(t_patch), POINTER  :: patch1, patch

  INTEGER :: turbmodel, modelclass, filtertype, errfl
  INTEGER :: ncv, ndv, nsv, nst, ngrad, nzof, filterwidth(diri:dirk)
#ifdef RFLO
  INTEGER :: idcbeg, idcend, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: idnbeg, idnend, jdnbeg, jdnend, kdnbeg, kdnend
  INTEGER :: n1, n2, ioff, icoff, ijcoff, inoff, ijnoff
  INTEGER :: homdir(diri:dirk)
#endif
#ifdef RFLU
  INTEGER :: npatches, ncellstot, nfaces, nfacestot, nbfaces, nbfacestot
#endif

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

  rcsidentstring = '$RCSfile: TURB_AllocateMemory.F90,v $ $Revision: 1.17 $'

  global => region%global
  CALL registerfunction( global,'TURB_AllocateMemory',&
  'TURB_AllocateMemory.F90' )

! get variables ---------------------------------------------------------------

  turbmodel      = region%mixtInput%turbModel
  modelclass     = region%turbInput%modelClass
  ncv            = region%turbInput%nCv
  ndv            = region%turbInput%nDv
  nsv            = region%turbInput%nSv
  nst            = region%turbInput%nSt
  ngrad          = region%turbInput%nGrad
  nzof           = region%turbInput%nZof
  filterwidth(:) = region%turbInput%filterWidth(:)  ! RFLU only uses 1st comp.

#ifdef RFLO
  filtertype     = region%turbInput%filterType
  homdir(:)      = region%turbInput%homDir(:)

! loop over all grid levels
  
  DO ilev=1,region%nGridLevels

    turb => region%levels(ilev)%turb

! - get cell and node dimensions ----------------------------------------------

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

    CALL rflo_getdimensdummynodes( region,ilev,idnbeg,idnend, &
                                   jdnbeg,jdnend,kdnbeg,kdnend )
    CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
    ibn = indijk(idnbeg,jdnbeg,kdnbeg,inoff,ijnoff)
    ien = indijk(idnend,jdnend,kdnend,inoff,ijnoff)
#endif
#ifdef RFLU
! - get variables and pointers ------------------------------------------------
    turb => region%turb

! - get cell and node dimensions
    ncellstot = region%grid%nCellsTot
    nfaces    = region%grid%nFaces
    nfacestot = region%grid%nFacesTot
    ibc       = 1
    iec       = ncellstot
    ibn       = 1
    ien       = nfaces
    npatches  = region%grid%nPatches

    nbfaces    = 0
    nbfacestot = 0

    DO ipatch = 1,npatches
      patch => region%patches(ipatch)

      nbfaces    = nbfaces    + patch%nBTris    + patch%nBQuads
      nbfacestot = nbfacestot + patch%nBTrisTot + patch%nBQuadsTot
    END DO ! iPatch
#endif

! - RANS/DES ----------------------------------------------------------------- 
! - RANS transport equations variables: cv, cvOld, rhs, rhsSum, diss, lens,
!                                       srad, epsIrs (only RFLO)

    IF ((modelclass==model_rans) .AND. (ncv > 0)) THEN
      ALLOCATE( turb%cv(    ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%cvOld( ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%rhs(   ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%rhsSum(ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%diss(  ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%dsterm(ncv,ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      ALLOCATE( turb%lens(      ibc:iec),stat=errfl )
      IF (errfl>0) goto 88
      IF (global%solverType==solv_implicit) THEN  ! Dual Tst
         ALLOCATE( turb%cvn(  ncv,ibc:iec),stat=errfl )
         IF (errfl>0) goto 88
         ALLOCATE( turb%cvn1( ncv,ibc:iec),stat=errfl )
         IF (errfl>0) goto 88
         ALLOCATE( turb%cvn2( ncv,ibc:iec),stat=errfl )
         IF (errfl>0) goto 88
         ALLOCATE( turb%sDual(ncv,ibc:iec),stat=errfl ); 
         IF (errfl>0) goto 88
      ELSE
         nullify( turb%cvn,turb%cvn1,turb%cvn2,turb%sDual )
!         NULLIFY( turb%cv,turb%cvOld,turb%rhs,turb%rhsSum,turb%diss)
      ENDIF  ! implicit
   ELSE
      nullify( turb%cvn,turb%cvn1,turb%cvn2,turb%sDual )
      nullify( turb%cv,turb%cvOld,turb%rhs,turb%rhsSum,turb%diss)
   ENDIF  ! rans
   
#ifdef RFLO
   IF ((modelclass==model_rans) .AND. (ncv > 0)) THEN
       ALLOCATE( turb%srad(diri:dirk,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
       IF ((global%flowType==flow_steady) .AND. &
            (region%turbInput%smoocf > 0._rfreal)) THEN
          ALLOCATE( turb%epsIrs(diri:dirk,ibc:iec), stat=errfl )
          IF (errfl>0) goto 88
       ELSE
          nullify( turb%epsIrs ) 
       ENDIF
    ELSE
       nullify( turb%srad,turb%epsIrs )
    ENDIF
#endif

! - RANS gradient variables
! - note: put within MODEL_RANS for not allocated twice by LES as they may be
!         allocated somewhere as LES work variables to save memory

    IF (modelclass==model_rans) THEN 
       IF (ngrad > 0) THEN
#ifdef RFLO
          ALLOCATE( turb%gradi(ngrad,ibn:ien),stat=errfl )
          IF (errfl>0) goto 88  
          ALLOCATE( turb%gradj(ngrad,ibn:ien),stat=errfl )
          IF (errfl>0) goto 88  
          ALLOCATE( turb%gradk(ngrad,ibn:ien),stat=errfl )
          IF (errfl>0) goto 88
       ELSE
          nullify( turb%gradi,turb%gradj,turb%gradk )
#endif
#ifdef RFLU
          ALLOCATE( turb%gradi( xcoord:zcoord,ngrad, nfaces),stat=errfl )
          IF (errfl>0) goto 88
          ALLOCATE( turb%bGradi(xcoord:zcoord,ngrad,nbfaces),stat=errfl )
          IF (errfl>0) goto 88
       ELSE
          nullify( turb%gradi, turb%bGradi )
#endif
       ENDIF
    ENDIF
    
! - total number of no-slip faces for wall dist., only on finest grid
    
#ifdef RFLO
    IF (ilev == 1) THEN
       DO ipatch=1,region%nPatches
          patch  => region%levels(ilev)%patches(ipatch)
          
          IF ((patch%bcType>=bc_noslipwall .AND. &
               patch%bcType<=bc_noslipwall+bc_range) .OR. &
               (patch%bcType>=bc_injection .AND. &
               patch%bcType<=bc_injection+bc_range)) 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)
             global%turbWallDim = global%turbWallDim + ijend-ijbeg+1
          ENDIF    ! bcType 
       ENDDO      ! iPatch
    ENDIF        ! iLev
#endif
#ifdef RFLU
    DO ipatch=1,region%grid%nPatches
       patch  => region%patches(ipatch)
       
       IF ((patch%bcType>=bc_noslipwall .AND. &
            patch%bcType<=bc_noslipwall+bc_range) .OR. &
            (patch%bcType>=bc_injection .AND. &
            patch%bcType<=bc_injection+bc_range)) THEN
          
          global%turbWallDim = global%turbWallDim + patch%nBFaces
       ENDIF    ! bcType 
    ENDDO      ! iPatch
#endif

! - GENERAL -------------------------------------------------------------------
! - TURB derived variables
    
    IF (ndv > 0) THEN
       ALLOCATE( turb%dv(ndv,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%dv )
    ENDIF
    
! - turbulent stress components
    
    IF (nsv > 0) THEN
       ALLOCATE( turb%sv(nsv,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%sv )
    ENDIF
    
! - vorticity components
    
    IF ((region%turbInput%calcVort /= calcvort_no) .OR. &
         (modelclass == model_rans)) THEN
       ALLOCATE( turb%vort(xcoord:zcoord,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%vort )
    ENDIF
    
! - zero-one switch field
    
    IF (nzof > 0) THEN
#ifdef RFLO
       ALLOCATE( turb%zofi(xcoord:zcoord,nzof,ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
       ALLOCATE( turb%zofj(xcoord:zcoord,nzof,ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
       ALLOCATE( turb%zofk(xcoord:zcoord,nzof,ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
#endif
#ifdef RFLU
       ALLOCATE( turb%zofi(xcoord:zcoord,nzof,nfaces),stat=errfl )
       IF (errfl>0) goto 88
       ALLOCATE( turb%bZofi(xcoord:zcoord,nzof,nbfaces),stat=errfl )
       IF (errfl>0) goto 88
#endif
    ENDIF
    
! - statistics of TURB pertinent variables (s.u. dynamic model coefficient)
    
#ifdef STATS
    IF ((global%flowType == flow_unsteady) .AND. &
         (global%doStat == active) .AND. &
         (global%turbNStat > 0)) THEN
       ALLOCATE( turb%tav(global%turbNStat,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%tav )
    ENDIF
    
    IF (nst > 0) THEN
       ALLOCATE( turb%st(nst,ibc:iec),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%st )
    ENDIF
#endif
    
! - LES -----------------------------------------------------------------------
! - LES filter coefficients, fix arrays
    
#ifdef RFLO
    IF ((turbmodel==turb_model_fixsmag) .OR. &  ! involve eddy vis.type
         (turbmodel==turb_model_dynsmag) .OR. &
         (turbmodel==turb_model_dynmixd)) THEN
       ALLOCATE( turb%fvolI(ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
       ALLOCATE( turb%fvolJ(ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
       ALLOCATE( turb%fvolK(ibn:ien),stat=errfl )
       IF (errfl>0) goto 88
    ELSE
       nullify( turb%fvolI, turb%fvolJ, turb%fvolK )
    ENDIF
    
    IF (((turbmodel==turb_model_scalsim)  .OR. & ! involve filtering
         (turbmodel==turb_model_dynsmag)  .OR. &
         (turbmodel==turb_model_dynmixd)) .AND. &
         (filtertype == filtype_nonunif)) THEN
       
       IF (homdir(diri) == off) THEN
          IF (filterwidth(diri) == filwidth_one) THEN
             ALLOCATE( turb%ccCofi1( 3,ibc:iec),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofi2( 3,ibc:iec),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi1I(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi1J(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi1K(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2I(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2J(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2K(3,ibn:ien),stat=errfl )
             IF (errfl>0) goto 88
          ELSEIF ((filterwidth(diri) == filwidth_two) .OR. &
               (filterwidth(diri) == filwidth_zero)) THEN
             ALLOCATE( turb%ccCofi2( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofi4( 5,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi2K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi4I(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi4J(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofi4K(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
          ENDIF
       ENDIF
       IF (homdir(dirj) == off) THEN
          IF (filterwidth(dirj) == filwidth_one) THEN
             ALLOCATE( turb%ccCofj1( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofj2( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj1I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj1J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj1K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
          ELSEIF ((filterwidth(dirj) == filwidth_two) .OR. &
               (filterwidth(dirj) == filwidth_zero)) THEN
             ALLOCATE( turb%ccCofj2( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofj4( 5,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj2K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj4I(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj4J(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofj4K(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
          ENDIF
       ENDIF
       IF (homdir(dirk) == off) THEN
          IF (filterwidth(dirk) == filwidth_one) THEN
             ALLOCATE( turb%ccCofk1( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofk2( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk1I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk1J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk1K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
          ELSEIF ((filterwidth(dirk) == filwidth_two) .OR. &
               (filterwidth(dirk) == filwidth_zero)) THEN
             ALLOCATE( turb%ccCofk2( 3,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ccCofk4( 5,ibc:iec),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2I(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2J(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk2K(3,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk4I(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk4J(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
             ALLOCATE( turb%ffCofk4K(5,ibn:ien),stat=errfl ); IF (errfl>0) goto 88
          ENDIF
       ENDIF
    ELSE
       nullify( turb%ccCofi1 ,turb%ccCofi2 ,turb%ccCofi4  )
       nullify( turb%ccCofj1 ,turb%ccCofj2 ,turb%ccCofj4  )
       nullify( turb%ccCofk1 ,turb%ccCofk2 ,turb%ccCofk4  )
       nullify( turb%ffCofi1I,turb%ffCofi2I,turb%ffCofi4I )
       nullify( turb%ffCofi1J,turb%ffCofi2J,turb%ffCofi4J )
       nullify( turb%ffCofi1K,turb%ffCofi2K,turb%ffCofi4K )
       nullify( turb%ffCofj1I,turb%ffCofj2I,turb%ffCofj4I )
       nullify( turb%ffCofj1J,turb%ffCofj2J,turb%ffCofj4J )
       nullify( turb%ffCofj1K,turb%ffCofj2K,turb%ffCofj4K )
       nullify( turb%ffCofk1I,turb%ffCofk2I,turb%ffCofk4I )
       nullify( turb%ffCofk1J,turb%ffCofk2J,turb%ffCofk4J )
       nullify( turb%ffCofk1K,turb%ffCofk2K,turb%ffCofk4K )
    ENDIF
#endif
#ifdef RFLU
    IF ((turbmodel==turb_model_fixsmag) .OR. &  ! involve eddy vis.type
         (turbmodel==turb_model_dynsmag) .OR. &
         (turbmodel==turb_model_dynmixd)) THEN
       ALLOCATE( turb%fvolI(  nfaces),stat=errfl ); IF (errfl>0) goto 88
       ALLOCATE( turb%bfVolI(nbfaces),stat=errfl ); IF (errfl>0) goto 88
    ELSE
       nullify( turb%fvolI, turb%bfVolI )
    ENDIF
    
    IF (modelclass == model_les) THEN
       ALLOCATE( turb%avgCoI( 2, nfaces),stat=errfl ); IF (errfl>0) goto 88
       ALLOCATE( turb%bAvgCoI(2,nbfaces),stat=errfl ); IF (errfl>0) goto 88
    ELSE
       nullify( turb%avgCoI, turb%bAvgCoI )
    ENDIF
#endif
    
! - nullify unused arrays
    
    IF (modelclass == model_rans) THEN
       nullify( turb%trace, turb%lij, turb%mij )
       nullify( turb%fVar, turb%ffVar, turb%ccVar, turb%coef, turb%mueT )
#ifdef RFLO
       nullify( turb%fISij, turb%fJSij, turb%fKSij )
#endif
#ifdef RFLU
       nullify( turb%bLij, turb%bMij )
       nullify( turb%bfVar, turb%bffVar, turb%bCoef, turb%bMueT )
       nullify( turb%fISij, turb%bfISij )
#endif
    ENDIF
    
! - WLM (wall layer model)-----------------------------------------------------
! - allocate arrays pertinent to WLM on coarser grid level (RFLO only)
    
#ifdef RFLO
    IF (ilev > 1) THEN
       DO ipatch=1,region%nPatches
          patch1 => region%levels(1)%patches(ipatch)
          patch  => region%levels(ilev)%patches(ipatch)
          
          IF (patch%bcType>=bc_noslipwall .AND. &
               patch%bcType<=bc_noslipwall+bc_range .AND. &
               patch1%valBola%switches(wlm_input_model) /= wlm_model_nomodel) THEN
             patch%valBola%nData     = patch1%valBola%nData
             patch%valBola%nSwitches = patch1%valBola%nSwitches
             patch%valBola%distrib   = patch1%valBola%distrib
             
             ALLOCATE( patch%valBola%switches(patch%valBola%nSwitches), &
                  stat=errfl ); IF (errfl>0) goto 88
             patch%valBola%switches  = patch1%valBola%switches
             
             IF (patch%valBola%nData > 0) THEN  ! vals at level 1 are defined
                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)
                
                ALLOCATE( patch%valBola%vals(ijbeg:ijend,patch%valBola%nData), &
                     stat=errfl ); IF (errfl>0) goto 88
                
             ENDIF  ! nData
          ENDIF    ! bcType 
       ENDDO      ! iPatch
    ENDIF        ! iLev
    
 ENDDO  !iLev
#endif
 
 goto 999
 
! finalize ----------------------------------------------------------
 
88 CONTINUE
 
 global%error = errfl
 CALL errorstop( global,err_allocate,__line__ )
 
999 CONTINUE
 
 CALL deregisterfunction( global )
 
END SUBROUTINE turb_allocatememory

!******************************************************************************
!
! RCS Revision history:
!
! $Log: TURB_AllocateMemory.F90,v $
! Revision 1.17  2009/08/31 05:17:35  mtcampbe
! Fix DES allocation bug
!
! Revision 1.16  2009/04/07 18:51:50  mtcampbe
! Fixed intel compiler complaint about line continuations
!
! Revision 1.15  2008/12/06 08:44:41  mtcampbe
! Updated license.
!
! Revision 1.14  2008/11/19 22:17:53  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.13  2006/01/12 09:49:44  wasistho
! enabled tripping fixed Smagorinsky
!
! Revision 1.12  2004/12/09 22:18:05  wasistho
! allocated added data structures
!
! Revision 1.11  2004/10/25 02:11:45  wasistho
! shift allocation of turb%st within ifdef STATS
!
! Revision 1.10  2004/09/23 22:22:58  wasistho
! added nullify turb%dsterm
!
! Revision 1.9  2004/08/04 02:49:06  wasistho
! removed turb%avgCoI,J,K as it is defined as grid%c2fCoI,J,K
!
! Revision 1.8  2004/06/03 02:10:25  wasistho
! enabled non-uniform fix-Smagorinsky
!
! Revision 1.7  2004/05/28 01:57:21  wasistho
! update unstructured grid LES
!
! Revision 1.6  2004/03/27 02:16:42  wasistho
! compiled with Rocflu
!
! Revision 1.5  2004/03/25 04:40:41  wasistho
! prepared for RFLU
!
! Revision 1.4  2004/03/20 03:28:29  wasistho
! prepared for RFLU
!
! Revision 1.3  2004/03/19 02:42:19  wasistho
! prepared for RFLU
!
! Revision 1.2  2004/03/13 03:10:45  wasistho
! prepared for RFLU
!
! Revision 1.1  2004/03/05 04:36:59  wasistho
! changed nomenclature
!
! Revision 1.15  2004/02/04 22:30:53  wasistho
! move MPIsum of global%turbWallDim from allocateMemory to calcMetrics
!
! Revision 1.14  2004/01/21 03:51:08  wasistho
! add injection wall condition for counting number of wall faces
!
! Revision 1.13  2003/10/21 20:31:23  wasistho
! added dt relaxation in steady flow due to RANS source term
!
! Revision 1.12  2003/10/07 20:31:24  wasistho
! allocate vorticities also for RaNS/DES
!
! Revision 1.10  2003/08/29 01:40:25  wasistho
! Added TARGET attribute to region variable, since pointers are cached into it
!
! Revision 1.9  2003/08/06 15:55:58  wasistho
! added vorticities computation
!
! Revision 1.8  2003/06/09 23:19:09  wasistho
! extra condition for wlm data allocation
!
! Revision 1.7  2003/05/31 01:46:14  wasistho
! installed turb. wall layer model
!
! Revision 1.6  2003/05/24 02:08:00  wasistho
! turbulence statistics expanded
!
! Revision 1.5  2003/05/15 02:57:05  jblazek
! Inlined index function.
!
! Revision 1.4  2002/11/02 02:05:37  wasistho
! Added TURB statistics
!
! Revision 1.3  2002/10/16 07:48:53  wasistho
! Enable Fix Smagorinsky
!
! Revision 1.2  2002/10/16 01:58:44  wasistho
! Changed global%error flag
!
! Revision 1.1  2002/10/14 23:55:29  wasistho
! Install Rocturb
!
!
!******************************************************************************