TURB_floLesGenCoCC.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: Get filter coefficients of cell-cell filtering.
!
! Description: The coefficients are obtained in i, j and k direction.
!              For each direction, distinction is made between the low delta
!              (1 or 2 grid-spacing) and high delta (2 or 4 grid-spacing).
!
! Input: region  = data of current region
!
! Output: turb%ccCofi1,2,4, turb%ccCofj1,2,4, turb%ccCofk1,2,4. 
!
! Notes: This routine is only relevant if non-uniform filter is selected.
!
!******************************************************************************
!
! $Id: TURB_floLesGenCoCC.F90,v 1.5 2008/12/06 08:44:43 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_flolesgencocc( region )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modinterfaces, ONLY : rflo_getdimensdummy,rflo_getdimensdummynodes, &
                            rflo_getcelloffset, rflo_getnodeoffset
  USE turb_modinterfaces, ONLY : turb_flolesgencocclo, turb_flolesgencocchi
  USE modturbulence, ONLY : t_turb
  USE moderror
  USE modparameters
  USE turb_modparameters
  IMPLICIT NONE

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

! ... loop variables

! ... local variables
  CHARACTER(CHRLEN)       :: rcsidentstring
  TYPE(t_global), POINTER :: global

  INTEGER           :: ibeg,iend,jbeg,jend,kbeg,kend
  INTEGER           :: idcbeg,idcend,jdcbeg,jdcend,kdcbeg,kdcend
  INTEGER           :: idnbeg,idnend,jdnbeg,jdnend,kdnbeg,kdnend
  INTEGER           :: ilev,icoff,ijcoff,inoff,ijnoff,minidx,maxidx,segid
  INTEGER           :: homdir(diri:dirk),filterwidth(diri:dirk)
  TYPE(t_turb), POINTER     :: turb
  REAL(RFREAL), POINTER     :: segm(:,:),cccofa(:,:),cccofb(:,:)
  REAL(RFREAL), ALLOCATABLE :: ds(:)

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

  rcsidentstring = '$RCSfile: TURB_floLesGenCoCC.F90,v $'

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

! get indices, parameters and pointers ----------------------------------------

  ilev           =  region%currLevel
  homdir(:)      =  region%turbInput%homDir(:)
  filterwidth(:) =  region%turbInput%filterWidth(:)
  turb           => region%levels(ilev)%turb
  CALL rflo_getdimensdummy( region,ilev,idcbeg,idcend, &
                            jdcbeg,jdcend,kdcbeg,kdcend )
  CALL rflo_getdimensdummynodes( region,ilev,idnbeg,idnend, &
                                 jdnbeg,jdnend,kdnbeg,kdnend )
  CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )

  minidx = min(idnbeg,jdnbeg,kdnbeg)-2   
  maxidx = max(idnend,jdnend,kdnend)+2   
  ALLOCATE( ds(minidx:maxidx) )

  ibeg = idcbeg
  iend = idcend
  jbeg = jdcbeg
  jend = jdcend
  kbeg = kdcbeg
  kend = kdcend

! filter coefficients in I direction

  IF (homdir(diri) == off) THEN
    IF (filterwidth(diri)==filwidth_one) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workK
      cccofa => turb%ccCofi1
      cccofb => turb%ccCofi2
      segid  =  1
      CALL turb_flolesgencocclo( global,diri,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ELSEIF ((filterwidth(diri)==filwidth_two) .OR. &
            (filterwidth(diri)==filwidth_zero)) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workK
      cccofa => turb%ccCofi2
      cccofb => turb%ccCofi4
      segid  =  1
      CALL turb_flolesgencocchi( global,diri,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ENDIF
  ENDIF

! filter coefficients in J direction

  IF (homdir(dirj) == off) THEN
    IF (filterwidth(dirj)==filwidth_one) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workI
      cccofa => turb%ccCofj1
      cccofb => turb%ccCofj2
      segid  =  1
      CALL turb_flolesgencocclo( global,dirj,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ELSEIF ((filterwidth(dirj)==filwidth_two) .OR. &
            (filterwidth(dirj)==filwidth_zero)) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workI
      cccofa => turb%ccCofj2
      cccofb => turb%ccCofj4
      segid  =  1
      CALL turb_flolesgencocchi( global,dirj,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ENDIF
  ENDIF

! filter coefficients in K direction

  IF (homdir(dirk) == off) THEN
    IF (filterwidth(dirk)==filwidth_one) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workJ
      cccofa => turb%ccCofk1
      cccofb => turb%ccCofk2
      segid  =  1
      CALL turb_flolesgencocclo( global,dirk,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ELSEIF ((filterwidth(dirk)==filwidth_two) .OR. &
            (filterwidth(dirk)==filwidth_zero)) THEN

! --- filter coefficients at cell centers
      segm   => region%levels(ilev)%turb%workJ
      cccofa => turb%ccCofk2
      cccofb => turb%ccCofk4
      segid  =  1
      CALL turb_flolesgencocchi( global,dirk,ibeg,iend,jbeg,jend,kbeg,kend, &
                              minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff,ds, &
                              segm,cccofa,cccofb )  

    ENDIF
  ENDIF

! deallocate temporary arrays

  DEALLOCATE( ds )

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

  CALL deregisterfunction( global )

END SUBROUTINE turb_flolesgencocc

!******************************************************************************
!
! RCS Revision history:
!
! $Log: TURB_floLesGenCoCC.F90,v $
! Revision 1.5  2008/12/06 08:44:43  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:55  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2004/08/04 02:50:34  wasistho
! removed turb%avgCoI,J,K replaced with turb%workI,J,K
!
! Revision 1.2  2004/03/12 02:55:36  wasistho
! changed rocturb routine names
!
! Revision 1.1  2004/03/08 23:35:45  wasistho
! changed turb nomenclature
!
! Revision 1.2  2003/08/29 01:41:24  wasistho
! Added TARGET attribute to region variable, since pointers are cached into it
!
! Revision 1.1  2002/10/14 23:55:29  wasistho
! Install Rocturb
!
!
!******************************************************************************