RFLU_InterpolateWrapper.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: Compute variables other than conserved variables at vertices
!
! Description: None.
!
! Input:
!   pRegion        Pointer to region data
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************
!
! $Id: RFLU_InterpolateWrapper.F90,v 1.15 2008/12/06 08:45:05 mtcampbe Exp $
!
! Copyright: (c) 2003-2006 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE rflu_interpolatewrapper(pRegion)

  USE moddatatypes
  USE modglobal, ONLY: t_global
  USE modparameters
  USE moderror
  USE moddatastruct, ONLY: t_region
  USE modgrid, ONLY: t_grid

  USE rflu_modinterpolation, ONLY: rflu_interpcells2verts, &
                                   rflu_interpsimplecells2verts

  USE modinterfaces, ONLY: rflu_computevertexvariables

  IMPLICIT NONE

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

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

  TYPE(t_region), POINTER :: pregion

! ==============================================================================
! Local variables
! ==============================================================================

  CHARACTER(CHRLEN) :: rcsidentstring
  TYPE(t_global), POINTER :: global
  TYPE(t_grid), POINTER :: pgrid

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

  rcsidentstring = '$RCSfile: RFLU_InterpolateWrapper.F90,v $ $Revision: 1.15 $'

  global => pregion%global

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

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

  pgrid => pregion%grid

! ******************************************************************************
! Mixture
! ******************************************************************************

! ==============================================================================
! Conserved variables
! ==============================================================================

  IF ( global%postInterpType == interp_type_proper ) THEN
    CALL rflu_interpcells2verts(pregion,global%postInterpOrder, &
                                pregion%mixtInput%nCv,pregion%mixt%cv, &
                                pregion%mixt%cvVert)
  ELSE IF ( global%postInterpType == interp_type_simple ) THEN
    CALL rflu_interpsimplecells2verts(pregion,pregion%mixtInput%nCv, & 
                                      pregion%mixt%cv,pregion%mixt%cvVert)
  ELSE
    CALL errorstop(global,err_reached_default,__line__)
  END IF ! global%postInterpType

! ==============================================================================
! Gas variables
! ==============================================================================

  SELECT CASE ( pregion%mixtInput%fluidModel )
    CASE ( fluid_model_comp ) 
      IF ( pregion%mixtInput%nGvAct == 0 ) THEN
        pregion%mixt%gvVert(gv_mixt_cp ,:) = pregion%mixt%gv(gv_mixt_cp ,:)
        pregion%mixt%gvVert(gv_mixt_mol,:) = pregion%mixt%gv(gv_mixt_mol,:)
      ELSE 
        IF ( global%postInterpType == interp_type_proper ) THEN
          CALL rflu_interpcells2verts(pregion,global%postInterpOrder, &
                                      pregion%mixtInput%nGv,pregion%mixt%gv, &
                                      pregion%mixt%gvVert)
        ELSE IF ( global%postInterpType == interp_type_simple ) THEN
          CALL rflu_interpsimplecells2verts(pregion,pregion%mixtInput%nGv, & 
                                            pregion%mixt%gv,pregion%mixt%gvVert)
        ELSE
          CALL errorstop(global,err_reached_default,__line__)
        END IF ! global%postInterpType
      END IF ! pRegion%mixtInput%nGvAct
    CASE default
      CALL errorstop(global,err_reached_default,__line__)
  END SELECT ! pRegion%mixtInput%fluidModel 

! ******************************************************************************
! Physical modules
! ******************************************************************************

#ifdef SPEC
  IF ( global%specUsed .EQV. .true. ) THEN

! ==============================================================================
!   Conserved variables
! ==============================================================================

    IF ( global%postInterpType == interp_type_proper ) THEN
      CALL rflu_interpcells2verts(pregion,global%postInterpOrder, &
                                  pregion%specInput%nSpecies, &
                                  pregion%spec%cv,pregion%spec%cvVert)
    ELSE IF ( global%postInterpType == interp_type_simple ) THEN
      CALL rflu_interpsimplecells2verts(pregion,pregion%specInput%nSpecies, &
                                        pregion%spec%cv,pregion%spec%cvVert)
    ELSE
      CALL errorstop(global,err_reached_default,__line__)
    END IF ! global%postInterpType
  END IF ! global%specUsed
#endif

! ******************************************************************************
! Compute remaining vertex variables
! ******************************************************************************

  CALL rflu_computevertexvariables(pregion)

! ******************************************************************************
! Interpolate plotting variables from cell centroids to vertices
! ******************************************************************************

  IF ( pregion%mixtInput%fluidModel == fluid_model_comp ) THEN 
    IF ( global%postInterpType == interp_type_proper ) THEN
      CALL rflu_interpcells2verts(pregion,global%postInterpOrder, &
                                  pregion%plot%nPv,pregion%plot%pv, &
                                  pregion%plot%pvVert)
    ELSE IF ( global%postInterpType == interp_type_simple ) THEN
      CALL rflu_interpsimplecells2verts(pregion,pregion%plot%nPv, & 
                                        pregion%plot%pv,pregion%plot%pvVert)
    ELSE
      CALL errorstop(global,err_reached_default,__line__)
    END IF ! global%postInterpType
  END IF ! pRegion%mixtInput%fluidModel

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

  CALL deregisterfunction(global)

END SUBROUTINE rflu_interpolatewrapper

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: RFLU_InterpolateWrapper.F90,v $
! Revision 1.15  2008/12/06 08:45:05  mtcampbe
! Updated license.
!
! Revision 1.14  2008/11/19 22:18:16  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.13  2007/03/19 21:41:44  haselbac
! Adapted to changes related to plotting variables
!
! Revision 1.12  2006/03/26 20:22:33  haselbac
! Added support for GL model, changed order
!
! Revision 1.11  2005/11/14 17:05:01  haselbac
! Generalized to support pseudo-gas model
!
! Revision 1.10  2005/11/10 02:47:52  haselbac
! Added support for gas models
!
! Revision 1.9  2005/10/31 21:09:39  haselbac
! Changed specModel and SPEC_MODEL_NONE
!
! Revision 1.8  2005/05/28 18:02:33  haselbac
! Replaced CV_MIXT_ENER by pRegion%mixtInput%nCv
!
! Revision 1.7  2005/05/01 14:22:19  haselbac
! Added interp of plotting vars
!
! Revision 1.6  2004/11/14 19:57:14  haselbac
! Added IF for fluid model
!
! Revision 1.5  2004/07/28 15:29:21  jferry
! created global variable for spec use
!
! Revision 1.4  2004/07/21 14:56:59  haselbac
! Added capability of using simple interpolation
!
! Revision 1.3  2004/03/18 03:34:11  haselbac
! Changed call to renamed interpolation function
!
! Revision 1.2  2003/12/04 03:29:32  haselbac
! Adapted to changes in RFLU_ModInterpolation
!
! Revision 1.1  2003/11/25 21:04:05  haselbac
! Initial revision
!
! ******************************************************************************