RVAV_ComputeAnalyticalSolution.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: Generate analytical data stream (S2) for various flow fields.
!
! Description: currently supported formats are:
!              - ASCII 
!              - Binary
!
! Input: case name from the list of arguments
!
! Output: File for Analytical solution
!
! Notes: 
!
!******************************************************************************
!
! $Id: RVAV_ComputeAnalyticalSolution.F90,v 1.3 2008/12/06 08:45:08 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE rvav_computeanalyticalsolution( similarityType,regionsS1 )

  USE moddatatypes
  USE moderror
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modgrid, ONLY       : t_grid
  USE modmixture, ONLY    : t_mixt
  USE modinterfaces, ONLY : rflo_getdimensdummy, rflo_getcelloffset, &
        rflo_getdimensphys, rflo_copygeometrydummy, rflo_calccellcentroids
  USE modmpi
  USE modparameters
  USE rvav_modglobal
  USE rvav_modparameters
  USE rvav_modinterfaces, ONLY : rvav_blasiussolution, rvav_gammbumpsolution, &
                                 rvav_proudmanculicksolution
  IMPLICIT NONE

#include "Indexing.h"
  
! ... parameters
  INTEGER :: similaritytype
  TYPE (t_region), POINTER :: regionss1(:)
  
! ... loop variables
  INTEGER :: i, j, k 
  
! ... local variables
  CHARACTER(CHRLEN)    :: msg
  CHARACTER(CHRLEN+23) :: fname

  TYPE(t_grid)   , POINTER :: grid
  TYPE(t_mixt)   , POINTER :: mixt
  TYPE(t_global) , POINTER :: global

  INTEGER :: ireg, ilev  
  INTEGER :: ipc, jpc, kpc, ibc, iec
  INTEGER :: idcbeg, idcend, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: icoff, ijcoff
  INTEGER :: icelloffset, ijcelloffset
  INTEGER :: ibpc,iepc
  INTEGER :: ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
  INTEGER :: ivars, inodes, jnodes, knodes, nvari, nvars, errorflag

  REAL(RFREAL), POINTER :: cofg(:,:)
  REAL(RFREAL), DIMENSION(:,:), ALLOCATABLE :: cvsa
  
!******************************************************************************

  global => regionss1(1)%global
  
  CALL registerfunction( global, 'RVAV_ComputeAnalyticalSolution',&
  'RVAV_ComputeAnalyticalSolution.F90' )
   
! open file in which stream2 is stored
          
  WRITE(fname,'(A,A)') trim(globalrvav%casename)//'_s2.anl'
  OPEN (if_rvav_file_s2,file=fname,form='formatted',status='unknown',&
        iostat=errorflag)

  global%error = errorflag
  IF (global%error /= 0) &
  CALL errorstop( global, err_file_open,__line__,'File: '//trim(fname) )  

! write number of regions 
  
  WRITE(if_rvav_file_s2,*,err=10) global%nRegions 
  
! calculate cell centroids ----------------------------------------------------
      
  IF (global%verbLevel /= verbose_none) &
    WRITE(stdout,'(/,A)') 'Compute Cell Centroids...'
        
  DO ireg=1, global%nRegions

    ilev =  regionss1(ireg)%currLevel
    grid => regionss1(ireg)%levels(ilev)%grid
    mixt => regionss1(ireg)%levels(ilev)%mixt

! - get cell and node dimensions

    CALL rflo_getdimensdummy( regionss1(ireg),ilev,idcbeg,idcend, &
                              jdcbeg,jdcend,kdcbeg,kdcend )
    CALL rflo_getcelloffset( regionss1(ireg),ilev,icelloffset,ijcelloffset )
    ibc = indijk(idcbeg,jdcbeg,kdcbeg,icelloffset,ijcelloffset)
    iec = indijk(idcend,jdcend,kdcend,icelloffset,ijcelloffset)
      
    CALL rflo_getdimensphys( regionss1(ireg),ilev,ipcbeg,ipcend, &
                             jpcbeg,jpcend,kpcbeg,kpcend )
    CALL rflo_getcelloffset( regionss1(ireg),ilev,icoff,ijcoff )
    
    ibpc = indijk(ipcbeg,jpcbeg,kpcbeg,icoff,ijcoff)
    iepc = indijk(ipcend,jdcend,kpcend,icoff,ijcoff)
      
    ALLOCATE( grid%cofg(3,ibc:iec),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )
        
    IF (global%verbLevel /= verbose_none) &
       WRITE(stdout,'(/,A,2X,I5)') ' nGridLevels = ',regionss1(ireg)%nGridLevels

    CALL rflo_copygeometrydummy( regionss1(ireg) )
    CALL rflo_calccellcentroids( regionss1(ireg) )
      
  ENDDO   ! iReg

! start case selection for analytical solutions -------------------------------

  SELECT CASE (similaritytype)
                                
! generates Input DataStream for C0
      
    CASE (rvav_culick)
                
      IF (global%verbLevel /= verbose_none) &
        WRITE(stdout,'(/,A)') 'Generates Input Deck for C0 ...'
        
      CALL rvav_proudmanculicksolution( fname,regionss1 )
         
! generates Input DataStream for GAMM Bump
           
    CASE (rvav_gammbump)
        
      IF (global%verbLevel /= verbose_none) &
        WRITE(stdout,'(/,A)') 'Generates Input Deck for GAMM Bump ...'
  
      CALL rvav_gammbumpsolution( fname,regionss1 )
          
! generates Input DataStream for Laminar Blasius
           
    CASE (rvav_blasius) 
        
      IF (global%verbLevel /= verbose_none) &
        WRITE(stdout,'(/,A)') 'Generates Input Deck for Laminar Blasius ...'

      CALL rvav_blasiussolution( fname,regionss1 )
          
  END SELECT ! similarityType  

! close file ------------------------------------------------------------------
      
  CLOSE(if_rvav_file_s2,iostat=errorflag)
  global%error = errorflag
  IF (global%error /= 0) &
  CALL errorstop( global, err_file_close,__line__,'File: '//trim(fname) )
  
  goto 999

10 CONTINUE
   CALL errorstop( global, err_file_read,__line__,'File: '//trim(fname) ) 
      
! finalize --------------------------------------------------------------------

999 CONTINUE

  CALL deregisterfunction( global )
  
1000 FORMAT('Region ',i5,', grid level= ',i2,'.')

END SUBROUTINE rvav_computeanalyticalsolution

!******************************************************************************
!
! RCS Revision history:
!
! $Log: RVAV_ComputeAnalyticalSolution.F90,v $
! Revision 1.3  2008/12/06 08:45:08  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:18:18  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 22:46:35  fnajjar
! Initial revision after changing case
!
! Revision 1.15  2003/11/20 16:40:41  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.11  2003/05/15 02:57:07  jblazek
! Inlined index function.
!
! Revision 1.10  2002/10/12 03:20:51  jblazek
! Replaced [io]stat=global%error with local errorFlag for Rocflo.
!
! Revision 1.9  2002/09/10 00:01:45  f-najjar
! Variable global moved into regions()
!
! Revision 1.8  2002/08/15 19:48:06  jblazek
! Implemented grid deformation capability.
!
! Revision 1.7  2002/07/31 02:34:57  f-najjar
! Split Analytical Solutions into individual routines
!
! Revision 1.6  2002/06/24 16:05:05  f-najjar
! Cleanup for file writing
!
! Revision 1.5  2002/06/24 15:48:09  f-najjar
! Included Blasius Profile Computation
!
! Revision 1.4  2002/06/21 14:11:15  f-najjar
! Clean up for Consistency
!
! Revision 1.3  2002/06/19 20:29:02  f-najjar
! Included GAMM Bump Computation
!
! Revision 1.2  2002/06/19 14:40:14  f-najjar
! Included verbLevel calls for cleanup
!
! Revision 1.1  2002/06/18 03:19:15  f-najjar
! Initial Import
!
!******************************************************************************