RFLU_PrintUserInput.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: Write out user input for checking purposes.
!
! Description: None.
!
! Input:
!   regions        Region data
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************
!
! $Id: RFLU_PrintUserInput.F90,v 1.60 2008/12/06 08:44:12 mtcampbe Exp $
!
! Copyright: (c) 2002-2006 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE rflu_printuserinput(regions)

  USE moddatatypes
  USE modglobal, ONLY: t_global
  USE moddatastruct, ONLY: t_region
  USE moderror
  USE modparameters
  USE modmixture, ONLY: t_mixt_input

#ifdef PERI
  USE modinterfacesperiodic, ONLY: peri_printuserinput
#endif

#ifdef PLAG
  USE modinterfaceslagrangian, ONLY: plag_printuserinput
#endif

#ifdef SPEC
  USE modinterfacesspecies, ONLY: spec_printuserinput
#endif

#ifdef TURB
  USE modinterfacesturbulence, ONLY: turb_printuserinput
#endif

#ifdef INRT
  USE modinterfacesinteract, ONLY: inrt_printuserinput, &
                                   inrt_printmaterialinput
#endif

  IMPLICIT NONE

! ******************************************************************************
! Declarations and definitions
! ******************************************************************************

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region), DIMENSION(:), POINTER :: regions

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: ireg
  CHARACTER(CHRLEN) :: rcsidentstring
  TYPE(t_mixt_input), POINTER :: pmixtinput
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  rcsidentstring = '$RCSfile: RFLU_PrintUserInput.F90,v $ $Revision: 1.60 $'

  global => regions(1)%global

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

! ******************************************************************************
! Set pointers and variables
! ******************************************************************************

  pmixtinput => regions(1)%mixtInput

  ireg = 1

! ******************************************************************************
! Echo user input
! ******************************************************************************

  WRITE(stdout,'(A,1X,A)') solver_name
  WRITE(stdout,'(A,1X,A)') solver_name,'Echoing user input...'

! ==============================================================================
! Region-independent variables
! ==============================================================================

! ------------------------------------------------------------------------------
! Formats
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Formats:'
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Grid:',global%gridFormat
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Solution:',global%solutFormat
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Grid source:',global%gridSource

! ------------------------------------------------------------------------------
! Reference values
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')          solver_name,'Reference values:'
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Velocity:',global%refVelocity
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Pressure:',global%refPressure
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Density:',global%refDensity
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name, &
                                    'Specific heat at constant pressure:', &
                                    global%refCp
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Ratio of specific heats:', &
                                    global%refGamma
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Length:',global%refLength
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reynolds number:', &
                                    global%refREnum
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Viscosity:',global%refVisc
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Laminar Prandtl number:', &
                                    global%prLam
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Turbulent Prandtl number:', &
                                    global%prTurb
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Laminar Schmidt number:', &
                                    global%scnLam
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Turbulent Schmidt number:', &
                                    global%scnTurb
                                    
  IF ( pmixtinput%gasModel == gas_model_mixt_gasliq ) THEN
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'BetaPLiq:',global%refBetaPLiq
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'BetaTLiq:',global%refBetaTLiq
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'CvLiq:',global%refCvLiq
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'DensityLiq:', & 
                                      global%refDensityLiq
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'PressLiq:',global%refPressLiq
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'TempLiq:',global%refTempLiq    
  END IF ! pMixtInput%gasModel                                    

! ------------------------------------------------------------------------------
! Forces
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Forces:'
  WRITE(stdout,'(A,5X,A,1X,L1)') solver_name,'Flag:',global%forceFlag
  WRITE(stdout,'(A,5X,A,1X,L1)') solver_name,'Patch coefficient flag:', & 
                                 global%patchCoeffFlag 
                                    
  IF ( global%forceFlag .EQV. .true. ) THEN                                     
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference area:', &
                                      global%forceRefArea
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference length:', &
                                      global%forceRefLength
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference x-coordinate:', &
                                      global%forceRefXCoord
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference y-coordinate:', &
                                      global%forceRefYCoord
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference z-coordinate:', &
                                      global%forceRefZCoord                                 
  END IF ! global%forceFlag

! ------------------------------------------------------------------------------
! Acceleration
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Acceleration:'
  WRITE(stdout,'(A,5X,A,1X,L1)') solver_name,'Flag:',global%accelOn 

  IF ( global%accelOn .EQV. .true. ) THEN 
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'x-component:',global%accelX
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'y-component:',global%accelY
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'z-component:',global%accelZ          
  END IF ! global%accelOn
    
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'gravity:',global%gravity

! ------------------------------------------------------------------------------
! Time stepping
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Time stepping:'
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Solver type:',global%solverType

  IF ( global%flowType == flow_unsteady ) THEN
    WRITE(stdout,'(A,5X,A)')          solver_name,'Flow type: Unsteady'
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Timestep:',global%dtImposed
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Current time:', &
                                      global%timeStamp
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Maximum time:', &
                                      global%maxTime
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Writing interval:', &
                                      global%writeTime
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Printing interval:', &
                                      global%printTime
    WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'Runge-Kutta scheme:', &
                                      global%rkScheme           
  ELSE
    WRITE(stdout,'(A,5X,A)')          solver_name,'Flow type: Steady'
    WRITE(stdout,'(A,5X,A,1X,I8)')    solver_name,'Current iteration:', &
                                      global%currentIter
    WRITE(stdout,'(A,5X,A,1X,I8)')    solver_name,'Maximum iteration:', &
                                      global%maxIter
    WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Residual tolerance:', &
                                      global%resTol
    WRITE(stdout,'(A,5X,A,1X,I8)')    solver_name,'Writing interval:', &
                                      global%writeIter
    WRITE(stdout,'(A,5X,A,1X,I8)')    solver_name,'Printing interval:', &
                                      global%printIter
  END IF ! global%flowType

! ------------------------------------------------------------------------------
! Time zooming
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Time Zooming:'
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Min Plane:',global%tzMinPlane
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Max Plane:',global%tzMaxPlane
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Propellant Density:',global%tzRhos
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Burn A:',global%tzA
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Burn N:',global%tzN
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Throat Radius:',global%tzThroatRad
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Chamber Length:',global%tzLenChamb
  IF(global%tzSubNoz) THEN
     WRITE(stdout,'(A,5X,A)') solver_name, 'Submerged Nozzle: Yes'
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Nozzle Inlet:',global%tzNozInlet
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Nozzle Radius:',global%tzNozRad
  ELSE
     WRITE(stdout,'(A,5X,A)') solver_name, 'Submerged Nozzle: No'
  ENDIF
  IF(global%tzCoordLong == xcoord) THEN
     WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: X'
  ELSE IF (global%tzCoordLong == ycoord) THEN
     WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: Y'
  ELSE
     WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: Z'
  ENDIF

! ------------------------------------------------------------------------------
! Case Constraints
! ------------------------------------------------------------------------------

  IF(global%cnstrCaseRad > 0.0) THEN
     WRITE(stdout,'(A,3X,A)')       solver_name,'Case Constraints:'
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Case Radius:',global%cnstrCaseRad
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Head End:',global%cnstrHeadEnd
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Aft End:',global%cnstrAftEnd
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Ellipse Major:',global%cnstrEllipsL
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Ellipse Minor:',global%cnstrEllipsT
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Sub Nozzle Rad:',global%cnstrNozY
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Tolerance 1:',global%cnstrTol1
     WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Tolerance 2:',global%cnstrTol2
     IF(global%cnstrCoordL == xcoord) THEN
        WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: X'
     ELSE IF (global%cnstrCoordL == ycoord) THEN
        WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: Y'
     ELSE
        WRITE(stdout,'(A,5X,A)')          solver_name,'Rocket Axis: Z'
     ENDIF
  ELSE
     WRITE(stdout,'(A,3X,A)')       solver_name,'No case constraints:'     
  ENDIF

! ------------------------------------------------------------------------------
! Random number generator values
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Random number generator values:'
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'SeedOffset:',&
                                 global%randSeedOffset
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'SeedType:',&
                                 global%randSeedType

#ifdef STATS
! ------------------------------------------------------------------------------
! Statistics
! ------------------------------------------------------------------------------

  IF ( global%doStat == active ) THEN
    IF ( global%reStat == active ) THEN
      WRITE(stdout,'(A,3X,A)') solver_name,'Continued statistics:'
    ELSE
      WRITE(stdout,'(A,3X,A)') solver_name,'New statistics:'
    END IF ! global%reStat

    IF ( global%mixtNStat > 0 ) THEN
      WRITE(stdout,'(A,5X,A,1X,I3)') solver_name,'mixture NStat ', &
      global%mixtNStat
      WRITE(stdout,'(A,5X,A,1X,20I3)') solver_name,'mixture StatID', &
      global%mixtStatID(1,:)*10+global%mixtStatID(2,:)
    END IF ! global%mixtNStat

    IF ( global%turbNStat > 0 ) THEN
      WRITE(stdout,'(A,5X,A,1X,I3)') solver_name,'TURB NStat    ', &
      global%turbNStat
      WRITE(stdout,'(A,5X,A,1X,20I3)') solver_name,'TURB StatID   ', &
      global%turbStatID(1,:)*10+global%turbStatID(2,:)
    END IF ! global%turbNStat
  END IF ! global%doStat
#endif

! ------------------------------------------------------------------------------
! Flow specific info
! ------------------------------------------------------------------------------

#ifdef PERI
  IF (regions(1)%periInput%flowKind /= off) THEN
    CALL peri_printuserinput( regions(1) )
  ENDIF  ! periodic flows
#endif

! ------------------------------------------------------------------------------
! Preprocessor
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')       solver_name,'Preprocessor:'
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Partition mode:', &
                                 global%prepPartMode

! ------------------------------------------------------------------------------
! Postprocessor
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)')          solver_name,'Postprocessor:'
  WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'Output format:', &
                                    global%postOutputFormat
  WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'Plot type:', &
                                    global%postPlotType
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Merge flag:', &
                                    global%postMergeFlag
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Write merged files flag:', &
                                    global%postWriteMergeFlag                               
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Plot volume flag:', &
                                    global%postPlotVolFlag
  WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'Interpolation type:', &
                                    global%postInterpType
  WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'Interpolation order:', &
                                    global%postInterpOrder
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Special cell flag:', &
                                    global%postSpecFlag
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Extraction flag:', &
                                    global%postExtractFlag
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Discontinuity flag:', &
                                    global%postDiscFlag
  WRITE(stdout,'(A,5X,A,1X,I3)')    solver_name,'Schlieren type:', &
                                    global%postSchType
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Schlieren exponent:', &
                                    global%postSchExp
  WRITE(stdout,'(A,5X,A,1X,I3)')    solver_name,'Number of fringes:', &
                                    global%postNFringes
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Vorticity flag:', &
                                    global%postVortFlag                                    
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Vortex core flag:', &
                                    global%postVortCoreFlag                                    
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Error flag:', &
                                    global%postCompErrFlag 
  WRITE(stdout,'(A,5X,A,1X,I3)')    solver_name,'Number of servers:', &
                                    global%postNServers
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Patch plotting flag:', &
                                    global%postPlotPatchFlag
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Lagrangian to Eulerian '// &
                                    'flag:',global%postLag2EulFlag   
  WRITE(stdout,'(A,5X,A,1X,L1)')    solver_name,'Gradient flag:', &
                                    global%postGradFlag                                                                        

#ifdef INRT
! ------------------------------------------------------------------------------
! Material input
! ------------------------------------------------------------------------------

  CALL inrt_printmaterialinput(global)
#endif

! ==============================================================================
! Region-dependent variables. NOTE these are not really region-dependent,
! just stored that way for commonality with Rocflo, so only write out for
! first region.
! ==============================================================================

! ------------------------------------------------------------------------------
! Flow model and grid motion
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)') solver_name,'Flow type:'

  IF ( pmixtinput%flowModel == flow_euler ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name,'Flow model: Euler'
  ELSE
    WRITE(stdout,'(A,5X,A)') solver_name,'Flow model: Navier-Stokes'
  END IF ! pMixtInput%flowModel

  IF ( pmixtinput%moveGrid .EQV. .true. ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name,'Grid motion: Active'
  ELSE
    WRITE(stdout,'(A,5X,A)') solver_name,'Grid motion: Not active'
  END IF ! pMixtInput%moveGrid

! ------------------------------------------------------------------------------
! Mixture
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)') solver_name,'Mixture:'
  WRITE(stdout,'(A,5X,A,1X,L1)') solver_name,'Frozen flag:', &
                                 pmixtinput%frozenFlag
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Gas model:',pmixtinput%gasModel
  
! ------------------------------------------------------------------------------
! Viscosity model
! ------------------------------------------------------------------------------

  IF ( pmixtinput%computeTv ) THEN
    WRITE(stdout,'(A,3X,A)') solver_name,'Viscosity model:'

    IF ( pmixtinput%viscModel == visc_suthr ) THEN
      WRITE(stdout,'(A,5X,A)') solver_name,'Viscosity model: Sutherland'
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference viscosity:', &
                                        pmixtinput%refVisc
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference temperature:', &
                                        pmixtinput%refTemp
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Sutherland coefficient:', &
                                        pmixtinput%suthCoef
    ELSE IF ( pmixtinput%viscModel == visc_fixed ) THEN
      WRITE(stdout,'(A,5X,A)') solver_name,'Viscosity model: Fixed'
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Viscosity value:', &
                                        pmixtinput%refVisc
    ELSE IF ( pmixtinput%viscModel == visc_antib ) THEN
      WRITE(stdout,'(A,5X,A)') solver_name,'Viscosity model: Antibes'
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Reference viscosity:', &
                                        pmixtinput%refVisc
    END IF ! pMixtInput%viscModel
  END IF ! pMixtInput%computeTv

! ------------------------------------------------------------------------------
! Numerics
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)') solver_name,'Numerics:'

  IF ( pmixtinput%timeScheme == tst_hyb5rk ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name, &
          'Time discretization: Explicit multistage'
  ELSE IF ( pmixtinput%timeScheme == tst_std4rk ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name, &
          'Time discretization: Explicit classical Runge-Kutta'
  END IF ! pMixtInput%timeScheme

  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'CFL number:',pmixtinput%cfl

  IF ( pmixtinput%spaceDiscr == discr_upw_roe ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name, &
          'Inviscid flux function: Roe'
  ELSE IF ( pmixtinput%spaceDiscr == discr_upw_hllc ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name, &
          'Inviscid flux function: HLLC'
  ELSE IF ( pmixtinput%spaceDiscr == discr_upw_ausmplus ) THEN
    WRITE(stdout,'(A,5X,A)') solver_name, &
          'Inviscid flux function: AUSM+'          
  END IF ! pMixtInput%spaceDiscr

  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Spatial order of accuracy:', &
                                 pmixtinput%spaceOrder
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Reconstruction:', &
                                 pmixtinput%reconst
  WRITE(stdout,'(A,5X,A)')       solver_name,'Constrained reconstruction:'
  WRITE(stdout,'(A,7X,A,1X,I2,1X,E12.5)') solver_name,'Cells:', & 
                                          pmixtinput%cReconstCells, &
                                          pmixtinput%cReconstCellsWeight
  WRITE(stdout,'(A,7X,A,1X,I2,1X,E12.5)') solver_name,'Faces:', & 
                                          pmixtinput%cReconstFaces, &
                                          pmixtinput%cReconstFacesWeight
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Entropy fix coefficient:', &
                                    pmixtinput%epsentr
  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'Dissipation factor:', &
                                    pmixtinput%dissFact                                    
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Dimensionality:', &
                                 pmixtinput%dimens 
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Cell-stencil dimensionality:', &
                                 pmixtinput%stencilDimensCells  
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name,'Face-stencil dimensionality:', &
                                 pmixtinput%stencilDimensFaces
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name, &
                                 'Boundary-face-stencil dimensionality:', &
                                 pmixtinput%stencilDimensBFaces                 
  WRITE(stdout,'(A,5X,A,1X,I2)') solver_name, &
                                 'Boundary-face-stencil space Order', &
                                 pmixtinput%spaceOrderBFaces 

  WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'In-cell test tolerance:', &
                                    pmixtinput%tolerICT 

! ------------------------------------------------------------------------------
! Grid motion
! ------------------------------------------------------------------------------

  IF ( pmixtinput%moveGrid .EQV. .true. ) THEN
    WRITE(stdout,'(A,3X,A)') solver_name,'Grid motion:'
    WRITE(stdout,'(A,5X,A,2X,I2)')    solver_name,'Type:', &
                                      pmixtinput%moveGridType
                                      
    IF ( pmixtinput%moveGridType /= movegrid_type_genx ) THEN                             
      WRITE(stdout,'(A,5X,A,1X,I2)')    solver_name,'NIter:', &
                                        pmixtinput%moveGridNIter
      WRITE(stdout,'(A,5X,A,1X,E12.5)') solver_name,'SFact:', &
                                        pmixtinput%moveGridSFact
    END IF ! pMixtInput%moveGridType                                    
  END IF ! pMixtInput%moveGrid

! ------------------------------------------------------------------------------
! Initialization
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)') solver_name,'Initialization:'
  WRITE(stdout,'(A,5X,A,2X,I2)')    solver_name,'Flag:     ', &
                                    global%initFlowFlag
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 1:', &
                                    pmixtinput%prepIntVal1 
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 2:', &
                                    pmixtinput%prepIntVal2 
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 3:', &
                                    pmixtinput%prepIntVal3   
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 4:', &
                                    pmixtinput%prepIntVal4 
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 5:', &
                                    pmixtinput%prepIntVal5 
  WRITE(stdout,'(A,5X,A,2X,I12)')   solver_name,'Integer 6:', &
                                    pmixtinput%prepIntVal6                                                                                                        
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  1:  ', &
                                    pmixtinput%prepRealVal1 
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  2:  ', &
                                    pmixtinput%prepRealVal2 
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  3:  ', &
                                    pmixtinput%prepRealVal3  
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  4:  ', &
                                    pmixtinput%prepRealVal4
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  5:  ', &
                                    pmixtinput%prepRealVal5 
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  6:  ', &
                                    pmixtinput%prepRealVal6 
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  7:  ', &
                                    pmixtinput%prepRealVal7 
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  8:  ', &
                                    pmixtinput%prepRealVal8
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real  9:  ', &
                                    pmixtinput%prepRealVal9
  WRITE(stdout,'(A,5X,A,2X,E12.5)') solver_name,'Real 10:  ', &
                                    pmixtinput%prepRealVal10

! ------------------------------------------------------------------------------
! Multi-physics modules
! ------------------------------------------------------------------------------

  WRITE(stdout,'(A,3X,A)') solver_name,'Multi-physics modules:'

! Turbulence -------------------------------------------------------------------

  IF ( pmixtinput%flowModel == flow_navst ) THEN
    IF ( pmixtinput%turbModel /= turb_model_none ) THEN
#ifdef TURB
      WRITE(stdout,'(A,5X,A)') solver_name,'Turbulence module: Active'
      CALL turb_printuserinput(regions(ireg))
#endif
    ELSE
      WRITE(stdout,'(A,5X,A)') solver_name,'Turbulence module: Not active'
    END IF ! MixtInput%turbModel
  END IF ! pMixtInput%flowModel

! Species ----------------------------------------------------------------------

  IF ( global%specUsed .EQV. .true. ) THEN
#ifdef SPEC
    WRITE(stdout,'(A,5X,A)') solver_name,'Species module: Active'
    CALL spec_printuserinput(regions(ireg))
#endif
  ELSE
    WRITE(stdout,'(A,5X,A)') solver_name,'Species module: Not active'
  END IF ! global%specUsed

! Particles --------------------------------------------------------------------

  IF ( global%plagUsed .EQV. .true. ) THEN
#ifdef PLAG
    WRITE(stdout,'(A,5X,A)') solver_name,'Particle module: Active'
    CALL plag_printuserinput(regions(ireg))
#endif
  ELSE
    WRITE(stdout,'(A,5X,A)') solver_name,'Particle module: Not active'
  END IF ! global%plagUsed

! Radiation --------------------------------------------------------------------

  IF ( pmixtinput%radiUsed .EQV. .true. ) THEN
#ifdef RADI
    WRITE(stdout,'(A,5X,A)') solver_name,'Radiation module: Active'
    CALL radi_printuserinput(regions(ireg))
#endif
  ELSE
    WRITE(stdout,'(A,5X,A)') solver_name,'Radiation module: Not active'
  END IF ! pMixtInput%radiUsed

#ifdef INRT
! Interactions -----------------------------------------------------------------

  IF ( global%inrtUsed .EQV. .true. ) THEN
    CALL inrt_printuserinput( regions(ireg) )
  ENDIF

#endif

! ******************************************************************************
! End
! ******************************************************************************

  WRITE(stdout,'(A,1X,A)') solver_name,'Echoing user input done.'
  WRITE(stdout,'(A,1X,A)') solver_name

  CALL deregisterfunction(global)

END SUBROUTINE rflu_printuserinput

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLU_PrintUserInput.F90,v $
! Revision 1.60  2008/12/06 08:44:12  mtcampbe
! Updated license.
!
! Revision 1.59  2008/11/19 22:17:25  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.58  2007/04/23 14:47:51  mtcampbe
! Added notification for NO case constraints
!
! Revision 1.57  2007/04/14 23:27:46  mtcampbe
! Updated for generalized geometry and support for submerged nozzles
!
! Revision 1.56  2007/04/14 21:13:35  mtcampbe
! Updated for TZ and Constraints
!
! Revision 1.55  2006/10/20 21:26:52  mparmar
! Added printing of gravity
!
! Revision 1.54  2006/08/19 15:38:46  mparmar
! Added printing of mixtInput%spaceOrderBFaces
!
! Revision 1.53  2006/04/07 15:19:16  haselbac
! Removed tabs
!
! Revision 1.52  2006/04/07 14:42:45  haselbac
! Added printing of stencilDimens params
!
! Revision 1.51  2006/03/26 20:21:39  haselbac
! Added printing of GL ref params, cosmetics
!
! Revision 1.50  2006/01/06 22:07:13  haselbac
! Added printing of postGradFlag and stencilDimens
!
! Revision 1.49  2005/12/25 15:27:56  haselbac
! Added printing for constrained reconstruction
!
! Revision 1.48  2005/12/24 21:26:26  haselbac
! Added printing of ICT tolerance
!
! Revision 1.47  2005/12/22 19:49:12  gzheng
! fixed a wrong format using L to print an integer, changed to 'I3'
!
! Revision 1.46  2005/12/10 23:28:26  haselbac
! Removed geom post variables
!
! Revision 1.45  2005/12/10 16:58:24  haselbac
! Added printing of postLag2EulFlag
!
! Revision 1.44  2005/12/01 17:09:53  fnajjar
! Added appropriate initialization, checking and printing of random seed type
!
! Revision 1.43  2005/11/17 14:37:53  haselbac
! Added printing of new prepRealVal vars
!
! Revision 1.42  2005/11/10 22:22:05  fnajjar
! ACH: Added frozenFlag
!
! Revision 1.41  2005/11/10 02:02:29  haselbac
! Added printing of acceleration components
!
! Revision 1.40  2005/10/31 19:25:48  haselbac
! Added printing of gasModel
!
! Revision 1.39  2005/10/28 19:17:29  haselbac
! Added printing of postPlotPatchFlag
!
! Revision 1.38  2005/10/27 18:56:51  haselbac
! Added printing of constr variable
!
! Revision 1.37  2005/10/05 20:03:09  haselbac
! Added printing of post output format and nservers
!
! Revision 1.36  2005/08/24 01:35:58  haselbac
! Fixed bug in writing solverType
!
! Revision 1.35  2005/08/10 00:32:08  haselbac
! Added printing of postVortCoreFlag
!
! Revision 1.34  2005/08/09 00:55:01  haselbac
! Added printing of patchCoeffFlag, postWriteMergeFlag, postCompErrFlag
!
! Revision 1.33  2005/08/03 18:19:24  hdewey2
! Added printing solverType
!
! Revision 1.32  2005/07/25 12:21:00  haselbac
! Added postVortFlag, changed format of postSchExp
!
! Revision 1.31  2005/07/14 21:58:45  haselbac
! Added output for AUSM flux function
!
! Revision 1.30  2005/07/11 19:24:13  mparmar
! Added printing of reconst option
!
! Revision 1.29  2005/07/05 19:26:58  haselbac
! Added printing of postSchType and postSchExp
!
! Revision 1.28  2005/05/01 14:19:04  haselbac
! Added printing of postDiscFlag and postNFringes
!
! Revision 1.27  2005/04/20 14:39:39  haselbac
! Added printing of additional int and read vals
!
! Revision 1.26  2005/03/22 03:33:14  haselbac
! Added printing of prep init helper variables
!
! Revision 1.25  2005/03/09 14:53:46  haselbac
! Added printing of dimensionality
!
! Revision 1.24  2004/11/17 16:28:11  haselbac
! Added printing of rkScheme
!
! Revision 1.23  2004/10/26 15:16:35  haselbac
! Added printing of postExtractFlag
!
! Revision 1.22  2004/10/19 19:25:14  haselbac
! Removed printing of surfDiverFlag
!
! Revision 1.21  2004/07/28 15:29:19  jferry
! created global variable for spec use
!
! Revision 1.20  2004/07/23 22:43:15  jferry
! Integrated rocspecies into rocinteract
!
! Revision 1.19  2004/07/21 14:54:49  haselbac
! Added printing of postInterpType
!
! Revision 1.18  2004/07/08 02:17:43  haselbac
! Added printing of dissFact
!
! Revision 1.17  2004/06/17 23:04:09  wasistho
! added PERI_PrintUserInput
!
! Revision 1.16  2004/06/16 20:00:22  haselbac
! Added force variables, clean-up, cosmetics
!
! Revision 1.15  2004/03/05 22:09:01  jferry
! created global variables for peul, plag, and inrt use
!
! Revision 1.14  2004/03/02 21:49:21  jferry
! Added inrtUsed flag to mixture data structure
!
! Revision 1.13  2004/02/02 22:49:40  haselbac
! Added printing of input for materials and particles
!
! Revision 1.12  2004/01/29 22:56:36  haselbac
! Cosmetic changes only
!
! Revision 1.11  2003/12/04 03:23:59  haselbac
! Cosmetic changes only
!
! Revision 1.10  2003/11/25 21:02:50  haselbac
! Added support for rocspecies, some clean-up
!
! Revision 1.9  2003/08/07 15:31:20  haselbac
! Added vars and changed some var names
!
! Revision 1.8  2003/07/22 01:57:51  haselbac
! Added writing of global%postInterpOrder
!
! Revision 1.7  2003/05/16 22:06:15  haselbac
! Fixed bug: LOGICALs should be written out using L format
!
! Revision 1.6  2003/05/05 18:40:15  haselbac
! Added printing of plotType, changed pltVolFlag to plotVolFlag
!
! Revision 1.5  2003/04/29 21:48:32  haselbac
! Added printing of surfDiverFlag
!
! Revision 1.4  2003/04/28 22:40:46  haselbac
! Added output of post and prep flags
!
! Revision 1.3  2003/04/10 23:29:59  fnajjar
! Added printouts for viscosity models
!
! Revision 1.2  2003/03/31 16:12:39  haselbac
! Cosmetics, added printing of grid-motion type
!
! Revision 1.1  2003/01/28 15:53:32  haselbac
! Initial revision, moved from rocflu
!
! Revision 1.7  2002/11/04 22:12:16  haselbac
! Bug fix: Added ifdef STATS
!
! Revision 1.6  2002/11/02 02:04:34  wasistho
! Added TURB statistics
!
! Revision 1.5  2002/09/17 22:51:23  jferry
! Removed Fast Eulerian particle type
!
! Revision 1.4  2002/09/09 15:51:56  haselbac
! global now under region and deleted debug output
!
! Revision 1.3  2002/06/17 13:34:12  haselbac
! Prefixed SOLVER_NAME to all screen output
!
! Revision 1.2  2002/05/04 17:11:25  haselbac
! Cosmetic changes
!
! Revision 1.1  2002/03/26 19:25:09  haselbac
! Initial revision
!
! ******************************************************************************