TURB_WlmInitia.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: Set initial values of variables needed in wlm computations.
!
! Description: Friction velocity u_tau is set to 0.05 u_edge.
!
! Input: region  = data of current region.
!        patch   = current patch.
!
! Output: initial u_tau.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: TURB_WlmInitia.F90,v 1.5 2008/12/06 08:44:42 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_wlminitia( region,patch )

  USE moddatatypes
  USE modbndpatch, ONLY   : t_patch
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
#ifdef RFLO
  USE modinterfaces, ONLY : rflo_getpatchindices, rflo_getcelloffset

#include "Indexing.h"
#endif
  USE moderror
  USE modparameters
  USE turb_modparameters
  IMPLICIT NONE

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

! ... loop variables
  INTEGER :: i, j, k

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

  INTEGER :: ibeg, iend, jbeg, jend, kbeg, kend, ijkc, ijkval
  REAL(RFREAL)          :: u, v, w
  REAL(RFREAL), POINTER :: vals(:,:)

#ifdef RFLO
  INTEGER :: ilev, lbound, icoff, ijcoff
  REAL(RFREAL), POINTER :: dv(:,:)
#endif
#ifdef RFLU
  REAL(RFREAL)          :: rdens
  REAL(RFREAL), POINTER :: cv(:,:)
#endif

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

  rcsidentstring = '$RCSfile: TURB_WlmInitia.F90,v $ $Revision: 1.5 $'

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

! get common pointers ---------------------------------------------------------

  vals => patch%valBola%vals 

#ifdef RFLO
! get dimensions and parameters -----------------------------------------------

  ilev   =  region%currLevel
  lbound =  patch%lbound

  CALL rflo_getpatchindices( region,patch,ilev,ibeg,iend,jbeg,jend,kbeg,kend )
  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )

! get pointers

  dv   => region%levels(ilev)%mixt%dv

! estimate initial friction velocity utau

  DO k=kbeg,kend
    DO j=jbeg,jend
      DO i=ibeg,iend

        ijkc = indijk(i      ,j      ,k      ,icoff,ijcoff)  

        IF (lbound==1 .OR. lbound==2) THEN
          ijkval = indij(j-jbeg ,k-kbeg ,jend-jbeg+1)
        ELSEIF (lbound==3 .OR. lbound==4) THEN
          ijkval = indij(k-kbeg ,i-ibeg ,kend-kbeg+1)
        ELSEIF (lbound==5 .OR. lbound==6) THEN
          ijkval = indij(i-ibeg ,j-jbeg ,iend-ibeg+1)
        ENDIF

        u = dv(dv_mixt_uvel,ijkc)     
        v = dv(dv_mixt_vvel,ijkc)     
        w = dv(dv_mixt_wvel,ijkc)
#endif
#ifdef RFLU
! specific Rocflu, check the state of cv first --------------------------------
  IF (region%mixt%cvState /= cv_mixt_state_cons) &
                            CALL errorstop(global,err_reached_default,__line__)

! get dimensions and parameters -----------------------------------------------

  ibeg =  1
  iend =  patch%nBFaces

! get pointers
  
  cv   => region%mixt%cv 

! estimate initial friction velocity utau

  DO i=ibeg,iend
        ijkc   = patch%bf2c(i) 
        ijkval = i

        rdens  = 1._rfreal/cv(cv_mixt_dens,ijkc) 
        u = cv(cv_mixt_xmom,ijkc)*rdens
        v = cv(cv_mixt_ymom,ijkc)*rdens
        w = cv(cv_mixt_zmom,ijkc)*rdens
#endif
        vals(ijkval,wlm_vals_utau) = 0.05*sqrt(u*u+v*v+w*w)

#ifdef RFLO
      ENDDO   ! i
    ENDDO     ! j
  ENDDO       ! k
#endif
#ifdef RFLU
  ENDDO       ! i
#endif

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

  CALL deregisterfunction( global )

END SUBROUTINE turb_wlminitia

!******************************************************************************
!
! RCS Revision history:
! 
! $Log: TURB_WlmInitia.F90,v $
! Revision 1.5  2008/12/06 08:44:42  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:55  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2004/03/27 02:16:42  wasistho
! compiled with Rocflu
!
! Revision 1.2  2004/03/24 03:37:03  wasistho
! prepared for RFLU
!
! Revision 1.1  2004/03/05 04:37:01  wasistho
! changed nomenclature
!
! Revision 1.2  2004/03/02 03:49:45  wasistho
! forgot colon after Id and Log
!
!
!******************************************************************************