TURB_CheckParamInput.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: Check TURB parameters either specified by user or set in the code.
!
! Description: The checking includes the existency and order of parameters.
!
! Input: regions = input parameters contained in turbInput of all regions.
!
! Output: Error msg.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: TURB_CheckParamInput.F90,v 1.19 2009/08/26 12:28:52 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_checkparaminput( regions )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modmixture, ONLY    : t_mixt_input
  USE modbndpatch, ONLY   : t_patch
#ifdef RFLO
  USE modinterfaces, ONLY : rflo_getdimensdummy
#endif
  USE modturbulence
  USE moderror
  USE modmpi
  USE modparameters
  USE turb_modparameters

  IMPLICIT NONE

! ... parameters
  TYPE(t_region), POINTER :: regions(:)

! ... loop variables
  INTEGER :: ireg, m, ipatch

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

  TYPE(t_mixt_input), POINTER  :: mixtinput  
  TYPE(t_turb_input), POINTER  :: input  
  TYPE(t_patch), POINTER       :: patch1
  LOGICAL :: turbinactive, fixedgrid, movegrid, wransactive, nonswall

#ifdef RFLO
  INTEGER :: idcbeg, idcend, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: ilev, maxrange
#endif

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

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

  global => regions(1)%global
  CALL registerfunction( global,'TURB_CheckParamInput',&
  'TURB_CheckParamInput.F90' )

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_none ) THEN
    WRITE(stdout,'(A,1X,A)') solver_name,'Entering TURB_CheckParamInput...'
  END IF ! global%verbLevel

! check fixed parameters setting ---------------------------------------------

! fixed parameters for general rocturb

  IF ((ycoord - xcoord)/=1 .OR. &
      (zcoord - ycoord)/=1) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'XCOORD,YCOORD,ZCOORD' )
  ENDIF

#ifdef RFLO
  IF (ndir /= 3) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'NDIR /= 3' )
  ENDIF
  IF ((dirk-dirj /= 1).OR.(dirj-diri /= 1).OR.(diri /= 1)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'DIRI,DIRJ,DIRK' )
  ENDIF
#endif

  IF ((cv_turb_vvel - cv_turb_uvel /= 1) .OR. &
      (cv_turb_wvel - cv_turb_vvel /= 1)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'CV_TURB_...' )
  ENDIF

  IF ((cv_turb_xmom - cv_turb_dens /= 1) .OR. &
      (cv_turb_ymom - cv_turb_xmom /= 1) .OR. &
      (cv_turb_zmom - cv_turb_ymom /= 1)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'CV_TURB_...' )
  ENDIF

  IF (((gr_turb_vx - gr_turb_ux)/=1).OR.((gr_turb_wx - gr_turb_vx)/=1).OR. &
      ((gr_turb_uy - gr_turb_wx)/=1).OR.((gr_turb_vy - gr_turb_uy)/=1).OR. &
      ((gr_turb_wy - gr_turb_vy)/=1).OR.((gr_turb_uz - gr_turb_wy)/=1).OR. &
      ((gr_turb_vz - gr_turb_uz)/=1).OR.((gr_turb_wz - gr_turb_vz)/=1)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'GR_TURB_...' )
  ENDIF

#ifdef RFLU
  IF ((cv_turb_dens /= cv_mixt_dens) .OR. &  ! relevant to LesLij, LesMij, LesHij,
      (cv_turb_xmom /= cv_mixt_xmom) .OR. &  ! etc
      (cv_turb_ymom /= cv_mixt_ymom) .OR. &
      (cv_turb_zmom /= cv_mixt_zmom)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'CV_TURB_.../= CV_MIXT_...' )
  ENDIF

!  IF ((CV_TURB_DENS /= 1) .OR. &    ! relevant to RFLU_InterpCells2Faces, etc
!      (E11          /= 1)) THEN
!    CALL ErrorStop( global,ERR_TURB_FIXPARAM,__LINE__,'CV_TURB_DES or E11 /= 1' )
!  ENDIF
#endif

! fixed parameters pertinent to LES

  IF ((filwidth_one  /= 1) .OR. &
      (filwidth_two  /= 2) .OR. &
      (filwidth_four /= 4)) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'FILWIDTH_...' )
  ENDIF

! fixed parameters pertinent to WLM 

  IF (wlm_vals_xix/=1 .OR. wlm_vals_tauux/=10 .OR. &
     (wlm_vals_etx-wlm_vals_xix)/=1 .OR. (wlm_vals_ztx-wlm_vals_etx)/=1 .OR. &
     (wlm_vals_ety-wlm_vals_xiy)/=1 .OR. (wlm_vals_zty-wlm_vals_ety)/=1 .OR. &
     (wlm_vals_etz-wlm_vals_xiz)/=1 .OR. (wlm_vals_ztz-wlm_vals_etz)/=1 .OR. &
     (wlm_vals_xiy-wlm_vals_ztx)/=1 .OR. (wlm_vals_xiz-wlm_vals_zty)/=1 .OR. &
     (wlm_vals_tauuy-wlm_vals_tauux)/=1.OR.(wlm_vals_tauuz-wlm_vals_tauuy)/=1.OR. &
     (wlm_vals_tauvy-wlm_vals_tauvx)/=1.OR.(wlm_vals_tauvz-wlm_vals_tauvy)/=1.OR. &
     (wlm_vals_tauwy-wlm_vals_tauwx)/=1.OR.(wlm_vals_tauwz-wlm_vals_tauwy)/=1.OR. &
     (wlm_vals_tauvx-wlm_vals_tauuz)/=1.OR.(wlm_vals_tauwx-wlm_vals_tauvz)/=1) THEN
    CALL errorstop( global,err_turb_fixparam,__line__,'WLM_VALS_...' )
  ENDIF

! check turbulence parameters selection

  turbinactive = .false.
  fixedgrid    = .false.
  movegrid     = .false.
  wransactive  = .false.
  nonswall     = .true.

#ifdef RFLO
  DO ireg = 1,global%nRegions
    IF (regions(ireg)%procid==global%myProcid .AND. &   ! region active and
        regions(ireg)%active==active) THEN              ! on my processor
#endif
#ifdef RFLU
  DO ireg = lbound(regions,1),ubound(regions,1)
#endif

      mixtinput => regions(ireg)%mixtInput
      input     => regions(ireg)%turbInput

      IF ((mixtinput%turbModel /= turb_model_none).AND. &
          (mixtinput%turbModel /= turb_model_fixsmag).AND. &
          (mixtinput%turbModel /= turb_model_scalsim).AND. &
          (mixtinput%turbModel /= turb_model_dynsmag).AND. &
          (mixtinput%turbModel /= turb_model_dynmixd).AND. &
          (mixtinput%turbModel /= turb_model_sa)     .AND. &
          (mixtinput%turbModel /= turb_model_dessa)  .AND. &
          (mixtinput%turbModel /= turb_model_hdessa)) THEN
         CALL errorstop( global,err_turb_model,__line__ )
      ENDIF
#ifdef RFLU
      IF (mixtinput%turbModel==turb_model_dynmixd) THEN
        CALL errorstop( global,err_turb_model,__line__, &
          'LES Dynamic Mixed model is not available with unstructured grid' )
      ENDIF
#endif
      IF ((input%modelClass /= model_rans) .AND. (input%nCv > 0)) THEN
        CALL errorstop( global,err_turb_fixparam,__line__, &
             'number of conservative variables > 0 only for RANS/DES' )
      ENDIF

      IF ((input%calcVort /= calcvort_no) .AND. &
          (input%calcVort /= calcvort_fdt) .AND. &
          (input%calcVort /= calcvort_sdt)) THEN 
        CALL errorstop( global,err_turb_input,__line__,'CALCVORTIC: 0, 1 or 2' )
      ENDIF

#ifdef GENX
      IF (input%calcVort == calcvort_no) THEN
        CALL errorstop( global,err_turb_input,__line__, &
        'Vorticities should always be computed in Genx, CALCVORTIC cannot be 0' )
      ENDIF
#endif

      IF ((input%nOutField > 1).AND.(input%calcVort == calcvort_no)) THEN
        CALL errorstop( global,err_turb_input,__line__, &
                        'CALCVORTIC should be > 0 for OUTPUTNUMBER > 1' )
      ENDIF

      IF (input%nZof > zof_nelm) THEN
        CALL errorstop( global,err_turb_fixparam,__line__, &
                        'ZOF_NELM < input%nZof, increase the former' )
      ENDIF

! --- pertinent to LES

      IF (input%modelClass == model_les) THEN
        IF ((mixtinput%turbModel==turb_model_fixsmag).OR. &
            (mixtinput%turbModel==turb_model_scalsim).OR. &
            (mixtinput%turbModel==turb_model_dynsmag).OR. &
            (mixtinput%turbModel==turb_model_dynmixd)) THEN
        ELSE
          CALL errorstop( global,err_turb_lesinput,__line__, &
            'selected turbulence model is not of LES class' )
        ENDIF

        IF ((input%nOutField < 1).OR.(input%nOutField > maxoutfld_les)) THEN
          CALL errorstop( global,err_turb_input,__line__, &
                          'OUTPUTNUMBER for LES out of range' )
        ENDIF

        IF ((input%cSmag < 0._rfreal).OR.(input%cSmag > 0.2_rfreal)) THEN
          CALL errorstop( global,err_turb_lesinput,__line__, &
                          '0<= CSMAGORINSKY <=0.2' )
        ENDIF

#ifdef RFLO
        IF ((input%filterType /= filtype_uniform) .AND. &
            (input%filterType /= filtype_nonunif)) THEN
          CALL errorstop( global,err_turb_lesinput,__line__, &
                          'FILTERTYPE: 0 or 1' )
        ENDIF

        IF ((input%deltaType /= deltype_cbrt) .AND. &
            (input%deltaType /= deltype_sqrt)) THEN
          CALL errorstop( global,err_turb_lesinput,__line__, &
                          'DELTATYPE: 0 or 1' )
        ENDIF

        DO m = diri,dirk
          IF ((input%filterWidth(m) /= filwidth_zero) .AND. &
              (input%filterWidth(m) /= filwidth_one)  .AND. &
              (input%filterWidth(m) /= filwidth_two)) THEN
            CALL errorstop( global,err_turb_lesinput,__line__, &
                            'FILTERWIDTH: 0,1 or 2' )
          ENDIF
          IF ((input%homDir(m) /= off) .AND. &
              (input%homDir(m) /= active)) THEN
            CALL errorstop( global,err_turb_lesinput,__line__, &
                            'HOMOGENDIR: 0 or 1' )
          ENDIF
        ENDDO
#endif
#ifdef RFLU
        IF ((input%filterWidth(diri) /= filwidth_zero) .AND. &
            (input%filterWidth(diri) /= filwidth_one)  .AND. &
            (input%filterWidth(diri) /= filwidth_two)) THEN
          CALL errorstop( global,err_turb_lesinput,__line__, &
                          'FILTERWIDTH: 0,1 or 2' )
        ENDIF
#endif
        IF ((input%engModel /= off) .AND. &
            (input%engModel /= active)) THEN
          CALL errorstop( global,err_turb_lesinput,__line__, &
                          'ENERGYMODEL: 0 or 1' )
        ENDIF

      ENDIF  ! modelClass LES

! --- pertinent to RaNS/DES

      IF (input%modelClass == model_rans) THEN
        IF ((mixtinput%turbModel==turb_model_sa).OR. &
            (mixtinput%turbModel==turb_model_dessa).OR. &
            (mixtinput%turbModel==turb_model_hdessa)) THEN
          IF (input%wDistMethod /= wdist_direct) THEN
            CALL errorstop( global,err_turb_ransinput,__line__, &
              'only direct WALLDISTMETHOD (0) currently valid for SA/DES' )
          ENDIF
!          IF ((mixtInput%moveGrid .eqv. .true.) .AND. (input%wDistMethod==WDIST_DIRECT)) THEN
!            CALL ErrorStop( global,ERR_TURB_RANSINPUT,__LINE__, &
!              'direct WALLDISTMETHOD (0) is not efficient for moving grid' )
!          ENDIF
          IF (input%functV1 /= sa_fv1_pow3 .AND. &
              input%functV1 /= sa_fv1_pow2) THEN
            CALL errorstop( global,err_turb_ransinput,__line__, &
              'selected formula for SA function fv1 is invalid' )
          ENDIF
        ELSE
          CALL errorstop( global,err_turb_ransinput,__line__, &
            'selected turbulence model is not of RANS/DES class' )
        ENDIF

        IF ((input%nOutField < 1).OR.(input%nOutField > maxoutfld_rans)) THEN
          CALL errorstop( global,err_turb_input,__line__, &
                          'OUTPUTNUMBER for RANS out of range' )
        ENDIF

#ifdef RFLO
        IF (input%spaceDiscr /= rans_discr_cen .AND. &
            input%spaceDiscr /= rans_discr_upw) THEN
          CALL errorstop( global,err_turb_ransinput,__line__, &
            'selected RaNS space discretization is not defined' )
        ENDIF
#endif
#ifdef RFLU
        IF (input%spaceDiscr /= rans_discr_upw) THEN
          CALL errorstop( global,err_turb_ransinput,__line__, &
            'selected RaNS space discretization is not defined in Rocflu' )
        ENDIF
#endif

        IF (input%spaceOrder /= rans_discr_ord1 .AND. &
            input%spaceOrder /= rans_discr_ord2) THEN
          CALL errorstop( global,err_turb_ransinput,__line__, &
            'selected RaNS space discretization order is not defined' )
        ENDIF
      ENDIF

! --- pertinent to WLM
 
#ifdef RFLO
      ilev = 1  ! check input based on finest level
      CALL rflo_getdimensdummy( regions(ireg),ilev,idcbeg,idcend, &
                                jdcbeg,jdcend,kdcbeg,kdcend )

      DO ipatch=1,regions(ireg)%nPatches
        patch1 => regions(ireg)%levels(ilev)%patches(ipatch)
#endif
#ifdef RFLU
      DO ipatch=1,regions(ireg)%grid%nPatches
        patch1 => regions(ireg)%patches(ipatch)
#endif
        IF (patch1%bcType>=bc_noslipwall .AND. &
            patch1%bcType<=bc_noslipwall+bc_range) THEN
          
          IF (patch1%valBola%nSwitches <= 0) THEN
            CALL errorstop( global,err_turb_fixparam,__line__, &
                           'nSwitches in bcvalues type valBola should be > 0' )
          ENDIF
          
          IF (patch1%valBola%switches(wlm_input_model) < 0 .OR. &
              patch1%valBola%switches(wlm_input_model) > wlm_model_extern) THEN
            CALL errorstop( global,err_turb_wlminput,__line__,'MODEL: 0-3' )
          ENDIF
          
          IF (patch1%valBola%switches(wlm_input_model) == wlm_model_extern) THEN
            CALL errorstop( global,err_turb_wlminput,__line__, &
            'Ext. tau_wall not ready yet' )
          ENDIF

          IF (patch1%valBola%switches(wlm_input_model)/=wlm_model_nomodel) THEN

            IF (minval(patch1%valBola%vals(:,wlm_vals_rough)) < 0._rfreal) THEN
              CALL errorstop( global,err_turb_wlminput,__line__, &
                           'ROUGHNESS: > 0.0' )
            ENDIF
#ifdef RFLO
            IF (patch1%valBola%switches(wlm_input_refpoint) < 1) THEN
              CALL errorstop( global,err_turb_wlminput,__line__, &
                             'WLM reference point should be >= 1' )
            ENDIF
            maxrange = idcend-idcbeg+ 1
            IF ((patch1%lbound == 1 .OR. patch1%lbound == 2) .AND. &
                 patch1%valBola%switches(wlm_input_refpoint) > maxrange) THEN
              CALL errorstop( global,err_turb_wlminput,__line__, &
                             'WLM reference point exceeds region max. range' )
            ENDIF
            maxrange = jdcend-jdcbeg+ 1
            IF ((patch1%lbound == 3 .OR. patch1%lbound == 4) .AND. &
                 patch1%valBola%switches(wlm_input_refpoint) > maxrange) THEN
              CALL errorstop( global,err_turb_wlminput,__line__, &
                             'WLM reference point exceeds region max. range' )
            ENDIF
            maxrange = kdcend-kdcbeg+ 1
            IF ((patch1%lbound == 5 .OR. patch1%lbound == 6) .AND. &
                 patch1%valBola%switches(wlm_input_refpoint) > maxrange) THEN
              CALL errorstop( global,err_turb_wlminput,__line__, &
                             'WLM reference point exceeds region max. range' )
            ENDIF
#endif
#ifdef RFLU
            CALL errorstop( global,err_turb_wlminput,__line__, &
                           'Wall Layer Model is not available yet in Rocflu' )
#endif
          ENDIF   ! switch
        ENDIF     ! bcType
      ENDDO       ! iPatch    

! --- check turbulence statistics input

      IF (input%nSt > st_turb_nvar) THEN
        CALL errorstop( global,err_turb_fixparam,__line__, &
                       'increase ST_TURB_NVAR in parameters module file' )
      ENDIF

#ifdef STATS
      IF (global%turbNStat > (input%nFixSt+input%nSv+input%nSt)) THEN
        CALL errorstop( global,err_turb_statsinput,__line__, &
            'TURBNSTAT larger than allocated; nSv (stress), nSt (stats) may be 0' )
      ENDIF

      IF ((global%doStat==1) .AND. &
          (mixtinput%turbModel/=turb_model_dynsmag) .AND. &
          (mixtinput%turbModel/=turb_model_dynmixd)) THEN
        DO m = 1,global%turbNStat
          IF (global%turbStatId(1,m)==0 .AND. global%turbStatId(2,m)>2 ) THEN
            CALL errorstop( global,err_turb_statsinput,__line__, &
            'TURBSTATID > 02 only for dyn.LES; chg setting if otherwise desired' )
          ENDIF
        ENDDO
      ENDIF
#endif

#ifdef RFLO
    ENDIF ! region active
#endif
  ENDDO   ! iReg

! global checking --------------------------------------------------------

#ifdef RFLO
  DO ireg=1,global%nRegions    ! serial process
#endif
#ifdef RFLU
  DO ireg = lbound(regions,1),ubound(regions,1)
#endif
    mixtinput => regions(ireg)%mixtInput
    input     => regions(ireg)%turbInput

! - general: turbInactive and moveGrid vs fixedGrid

    IF (mixtinput%turbModel == turb_model_none) THEN
      turbinactive = .true.
    ENDIF
    IF       (mixtinput%moveGrid .eqv. .true.)  movegrid = .true.
    IF (.NOT. (mixtinput%moveGrid .eqv. .true.)) fixedgrid = .true.

! - RaNS/DES: wRansActive and noNsWall

    IF (input%modelClass == model_rans) THEN
      IF ((mixtinput%turbModel==turb_model_sa).OR. &
          (mixtinput%turbModel==turb_model_dessa).OR. &
          (mixtinput%turbModel==turb_model_hdessa)) THEN
        wransactive = .true.          
      ENDIF
    ENDIF

# ifdef RFLO
    DO ipatch=1,regions(ireg)%nPatches
      patch1 => regions(ireg)%levels(ilev)%patches(ipatch)
#endif
# ifdef RFLU
    DO ipatch=1,regions(ireg)%grid%nPatches
      patch1 => regions(ireg)%patches(ipatch)
#endif
      IF (patch1%bcType>=bc_noslipwall .AND. &
          patch1%bcType<=bc_noslipwall+bc_range) THEN
        nonswall = .false.
      ENDIF
    ENDDO

  ENDDO ! iReg

#ifdef GENX
  IF (turbinactive .eqv. .true.) THEN
    CALL errorstop( global,err_turb_region,__line__, &
      'For Genx, turbulence should be active in all regions or none at all' )
  ENDIF
#endif
  IF ((fixedgrid .eqv. .true.) .AND. (movegrid .eqv. .true.)) THEN
    CALL errorstop( global,err_turb_input,__line__, &
      'QUESTION: Is it safe to only update metrics in regions that move?' )
  ENDIF

  IF ((wransactive .eqv. .true.) .AND. (nonswall .eqv. .true.)) THEN
    CALL errorstop( global,err_turb_input,__line__, &
      'RaNS model selected needs wall distance but no ns-wall presents' )
  ENDIF

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

  IF ( global%myProcid == masterproc .AND. &
       global%verbLevel > verbose_none ) THEN
    WRITE(stdout,'(A,1X,A)') solver_name,'Leaving TURB_CheckParamInput.'
  END IF ! global%verbLevel

  CALL deregisterfunction( global )

END SUBROUTINE turb_checkparaminput

!******************************************************************************
!
! RCS Revision history:
!
! $Log: TURB_CheckParamInput.F90,v $
! Revision 1.19  2009/08/26 12:28:52  mtcampbe
! Ported to Hera.   Fixed logical expression syntax errors.  Replaced all
! IF (logical_variable)  with IF (logical_variable .eqv. .true.) as
! consistent with the specification.  Also changed: IF( ASSOCIATED(expr) )
! to IF ( ASSOCIATED(expr) .eqv. .true. ).   Intel compilers produce code
! which silently fails for some mal-formed expressions, so these changes
! are a net which should ensure that they are evaluated as intended.
!
! Revision 1.18  2008/12/06 08:44:41  mtcampbe
! Updated license.
!
! Revision 1.17  2008/11/19 22:17:53  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.16  2006/02/04 05:00:18  wasistho
! added enter and leave statements
!
! Revision 1.15  2006/01/12 09:50:07  wasistho
! enabled tripping fixed Smagorinsky
!
! Revision 1.14  2005/12/29 19:47:14  wasistho
! bug fixed start/end index of ireg for rflu
!
! Revision 1.13  2005/03/09 06:34:45  wasistho
! incorporated HDESSA
!
! Revision 1.12  2005/03/08 01:54:56  wasistho
! fixed bug, added if dostat cheking in turb-statistics
!
! Revision 1.11  2004/12/10 01:03:13  wasistho
! trap turbStatId error in connection with turbmodel selected
!
! Revision 1.10  2004/08/12 20:58:39  wasistho
! check roughness value only when WLM is active
!
! Revision 1.9  2004/08/07 01:08:03  wasistho
! error msg if Dyn.Mixed selected in Rocflu
!
! Revision 1.8  2004/07/03 01:44:00  wasistho
! removed restriction to only fixed Smag in Rocflu
!
! Revision 1.7  2004/05/28 01:59:43  wasistho
! update unstructured grid LES
!
! Revision 1.6  2004/04/08 04:00:27  wasistho
! allow DES with moving grid
!
! Revision 1.5  2004/03/27 02:16:42  wasistho
! compiled with Rocflu
!
! Revision 1.4  2004/03/24 03:37:02  wasistho
! prepared for RFLU
!
! Revision 1.3  2004/03/20 03:28:29  wasistho
! prepared for RFLU
!
! Revision 1.2  2004/03/19 02:45:12  wasistho
! prepared for RFLU
!
! Revision 1.1  2004/03/05 21:08:02  wasistho
! changed nomenclature
!
! Revision 1.18  2004/02/26 21:25:50  wasistho
! delete energy sgs warning
!
! Revision 1.17  2004/02/24 21:03:23  wasistho
! modified the warning previously checked in
!
! Revision 1.16  2004/02/24 02:52:18  wasistho
! added warning for unbalance load MPI if LES energy model is used
!
! Revision 1.15  2004/02/19 04:03:42  wasistho
! added new rans/SA parameter VISCFUNCTION
!
! Revision 1.14  2004/02/14 03:43:07  wasistho
! added new WLM parameter: reference point
!
! Revision 1.13  2004/02/11 03:24:40  wasistho
! added feature: variable number of turbulence output fields
!
! Revision 1.12  2003/12/24 00:04:36  wasistho
! modify WLM checking to satisfy Turing
!
! Revision 1.11  2003/10/26 00:56:14  wasistho
! bug fixed
!
! Revision 1.10  2003/10/26 00:09:26  wasistho
! added multiple discr.types and order
!
! Revision 1.9  2003/10/07 02:05:04  wasistho
! initial installation of RaNS-SA and DES
!
! Revision 1.8  2003/09/11 03:07:30  wasistho
! removed ASSOCIATED(patch1%valBola%switches)
!
! Revision 1.7  2003/08/29 20:30:35  wasistho
! Added check for vorticity computation in Genx
!
! Revision 1.6  2003/08/06 15:56:05  wasistho
! added vorticities computation
!
! Revision 1.5  2003/07/18 20:12:56  wasistho
! put ifdef STATS
!
! Revision 1.4  2003/07/17 01:23:10  wasistho
! prepared rocturb for Genx
!
! Revision 1.3  2003/05/31 01:46:21  wasistho
! installed turb. wall layer model
!
! Revision 1.2  2003/05/24 02:11:05  wasistho
! turbulence statistics expanded
!
! Revision 1.1  2002/10/14 23:55:29  wasistho
! Install Rocturb
!
!
!******************************************************************************