RVAV_ReadFileStream2Analyt.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: read File for Stream 2 Data based on Analytical Solutions.
!
! Description: none.
!
! Input: 
!
! Output: Memory location with data for Stream 2
!
! ISSUE: Where do you read the Analytical and Experimental Data on Nodes or Cells
! 
! Notes: none.
!
!******************************************************************************
!
! $Id: RVAV_ReadFileStream2Analyt.F90,v 1.3 2008/12/06 08:45:08 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE rvav_readfilestream2analyt ( global, regionsS2 )

  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 modmpi
  USE modparameters
  USE rvav_modparameters
  USE rvav_modglobal
  USE rvav_moddatastruct
  IMPLICIT NONE

#include "Indexing.h"
  
! ... parameter variables
  TYPE (t_region), POINTER :: regionss2(:)
  TYPE(t_global)  , POINTER  :: global
  
! ... loop variables
  INTEGER :: ireg, ilev

! ... local variables
  CHARACTER(CHRLEN) :: msg, fname
  
  INTEGER :: i, j, k
  INTEGER :: inodes, jnodes, knodes, nvars, ivars
  INTEGER :: ipc, jpc, kpc, ibc, iec, ibn, ien
  INTEGER :: idcbeg, idcend, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: idnbeg, idnend, jdnbeg, jdnend, kdnbeg, kdnend
  INTEGER :: ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
  INTEGER :: ipnbeg, ipnend, jpnbeg, jpnend, kpnbeg, kpnend
  INTEGER :: icoff, ijcoff, inoff, ijnoff
  INTEGER :: nregionss2
  INTEGER :: ijkn0, ijkc0
  INTEGER :: ngridlevels, ivar, errorflag

  REAL(RFREAL), DIMENSION(3)               :: xyzmin, xyzmax
  REAL(RFREAL), DIMENSION(:), ALLOCATABLE  ::  cvmin,  cvmax
  REAL(RFREAL), POINTER                    ::  cv(:,:),xyz(:,:)  
  
  TYPE(t_grid)    , POINTER  :: grid
  TYPE(t_mixt)    , POINTER  :: mixt  
    
!******************************************************************************

  CALL registerfunction( global, 'RVAV_ReadFileStream2Analyt',&
  'RVAV_ReadFileStream2Analyt.F90' )
           
  IF ( global%verbLevel/=verbose_none ) &          
    WRITE(stdout,'(/,A)') 'Reading grid and solution from Stream 2 - ANALYTICAL...' 

  global%casename = trim(globalrvav%casename)//'_s2'

! open file

  IF ( globalrvav%fileTypeS2 == file_analytical ) THEN
    WRITE(fname,'(A,I5.5,A,1PE11.5)') trim(global%casename)//'.anl'
  ENDIF ! fileTypeS2

  OPEN(if_rvav_file_s2,file=fname,form='formatted',status='old',iostat=errorflag)
  global%error = errorflag
  IF (global%error /= 0) &
    CALL errorstop( global, err_file_open,__line__,'File: '//trim(fname) )

! read number of Regions 

  READ(if_rvav_file_s2,*,err=10,end=10) nregionss2
    
  globalrvav%nRegionsS2 = nregionss2 
     
  ALLOCATE( regionss2(globalrvav%nRegionsS2),stat=errorflag )
  global%error = errorflag
  IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )
    
  DO ireg=1,globalrvav%nRegionsS2

    WRITE(stdout,'(A,I5.5)') '  - region ',ireg

    regionss2(ireg)%startLevel = global%startLevel
    regionss2(ireg)%currLevel  = global%startLevel

    ilev =  regionss2(ireg)%currLevel
    regionss2(ireg)%nGridLevels = regionss2(ireg)%startLevel
    ngridlevels = regionss2(ireg)%nGridLevels
      
    IF ( global%verbLevel/=verbose_none ) THEN 
      WRITE(stdout,'(/,A,I5)') 'regionsS2(iReg)%startLevel=',regionss2(ireg)%startLevel
      WRITE(stdout,'(A,I5)')   'regionsS2(iReg)%currLevel =',regionss2(ireg)%currLevel
      WRITE(stdout,'(A,I5)')   'iLev = ',ilev
    END IF ! verbLevel
      
    ALLOCATE( regionss2(ireg)%levels(ngridlevels),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )

  ENDDO ! iReg

  DO ireg=1,globalrvav%nRegionsS2
          
    regionss2(ireg)%startLevel = global%startLevel
    regionss2(ireg)%currLevel  = global%startLevel 

    ilev =  regionss2(ireg)%currLevel
      
    IF ( global%verbLevel/=verbose_none ) THEN 
      WRITE(stdout,'(A,I5.5)') '  - region ',ireg
      WRITE(stdout,'(/,A,I5)') 'regionsS2(iReg)%startLevel=',regionss2(ireg)%startLevel
      WRITE(stdout,'(A,I5)')   'regionsS2(iReg)%currLevel =',regionss2(ireg)%currLevel
      WRITE(stdout,'(A,I5)')   'iLev = ',ilev
    END IF ! verbLevel

    grid => regionss2(ireg)%levels(ilev)%grid
    mixt => regionss2(ireg)%levels(ilev)%mixt
             
! - read dimensions

    READ(if_rvav_file_s2,*,err=10,end=10) inodes,jnodes,knodes,nvars 
            
    IF ( global%verbLevel/=verbose_none ) THEN
      WRITE(stdout,'(A,3(I5,2X))')   'nRegionsS2 = ',nregionss2
      WRITE(stdout,'(A,3(I5,2X))')   'iNodes,jNodes,kNodes = ',inodes,jnodes,knodes
      WRITE(stdout,'(A,1(I5,2X))')   'nVars = ',nvars
    END IF ! verbLevel

! - find array extent
      
    ipnbeg = 1
    jpnbeg = 1
    kpnbeg = 1
      
    ipnend = inodes
    jpnend = jnodes
    kpnend = knodes
 
    inoff  = ipnend - ipnbeg + 1
    ijnoff = inoff*(jpnend-jpnbeg+1)
      
    ibn = indijk(ipnbeg,jpnbeg,kpnbeg,inoff,ijnoff)
    ien = indijk(ipnend,jpnend,kpnend,inoff,ijnoff) 
                              
    globalrvav%iCOffS2  = inoff
    globalrvav%ijCOffS2 = ijnoff
      
! - allocate memory

    ALLOCATE( grid%xyz(3,ibn:ien),stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )

    ALLOCATE( mixt%cv(nvars-3,ibn:ien) ,stat=errorflag )
    global%error = errorflag
    IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )
    
! - initialize variables 

    grid%xyz = 0.0_rfreal
    mixt%cv  = 0.0_rfreal

! - create new pointers for reading data 

    xyz => grid%xyz
     cv => mixt%cv
        
! - read stream2 data 

    READ(if_rvav_file_s2,*,err=10,end=10) &
    (((xyz(xcoord,indijk(i,j,k,inoff,ijnoff)),    i=ipnbeg,ipnend), &
                                                  j=jpnbeg,jpnend), &
                                                  k=kpnbeg,kpnend), &
    (((xyz(ycoord,indijk(i,j,k,inoff,ijnoff)),    i=ipnbeg,ipnend), &
                                                  j=jpnbeg,jpnend), &
                                                  k=kpnbeg,kpnend), &
    (((xyz(zcoord,indijk(i,j,k,inoff,ijnoff)),    i=ipnbeg,ipnend), &
                                                  j=jpnbeg,jpnend), &
                                                  k=kpnbeg,kpnend), &
    ((((cv(ivars, indijk(i,j,k,inoff,ijnoff)),    i=ipnbeg,ipnend), &
                                                  j=jpnbeg,jpnend), &
                                                  k=kpnbeg,kpnend), &
                                              ivars=1,nvars-zcoord)

! - Write Min-Max Values

    IF ( global%verbLevel/=verbose_none ) THEN
      
      WRITE(stdout,'(/,A,2(I5,2X))') 'ibn-ien ',ibn,ien
      WRITE(stdout,'(A,3(I5,2X))')   'ipnbeg,jpnbeg,kpnbeg',ipnbeg,jpnbeg,kpnbeg
      WRITE(stdout,'(A,3(I5,2X))')   'ipnend,jpnend,kpnend',ipnend,jpnend,kpnend
      WRITE(stdout,'(A,2(I5,2X))')   'iNOff,ijNOff',inoff,ijnoff
      WRITE(stdout,'(A,2(I5,2X))')   'iCOffS2,ijCOffS2',globalrvav%iCOffS2,globalrvav%ijCOffS2
        
      xyzmin = +1.0e+30_rfreal
      xyzmax = -1.0e+30_rfreal

      DO i = ibn, ien
        xyzmin(1) = min(grid%xyz(1,i), xyzmin(1))
        xyzmin(2) = min(grid%xyz(2,i), xyzmin(2))
        xyzmin(3) = min(grid%xyz(3,i), xyzmin(3))

        xyzmax(1) = max(grid%xyz(1,i), xyzmax(1))
        xyzmax(2) = max(grid%xyz(2,i), xyzmax(2))
        xyzmax(3) = max(grid%xyz(3,i), xyzmax(3))
      END DO ! i

      DO k = kpnbeg, kpnbeg
      DO j = jpnbeg, jpnbeg
      DO i = ipnbeg, ipnend
        ijkn0  = indijk(i  ,j,k,inoff,ijnoff)
        WRITE(stdout,'(I5,3(3X,E12.5))') i,grid%xyz(1,ijkn0),grid%xyz(2,ijkn0),grid%xyz(3,ijkn0)
      END DO ! i 
      END DO ! j
      END DO ! k

      WRITE(stdout,'(/,A,2(E12.5,2X))') 'Min-Max of X     ', xyzmin(1),xyzmax(1)
      WRITE(stdout,'(A,2(E12.5,2X))')   'Min-Max of Y     ', xyzmin(2),xyzmax(2)
      WRITE(stdout,'(A,2(E12.5,2X))')   'Min-Max of Z     ', xyzmin(3),xyzmax(3)

      ALLOCATE(cvmin(nvars-zcoord),stat=errorflag )
      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )
      ALLOCATE(cvmax(nvars-zcoord),stat=errorflag )
      global%error = errorflag
      IF (global%error /= 0) CALL errorstop( global, err_allocate,__line__ )

      cvmin = +1.0e+30_rfreal
      cvmax = -1.0e+30_rfreal

      DO i = ibn, ien
        DO ivar = 1, nvars-zcoord
          cvmin(ivar) = min(mixt%cv(ivar,i), cvmin(ivar))
          cvmax(ivar) = max(mixt%cv(ivar,i), cvmax(ivar))
        END DO ! ivar 
      END DO ! i

        WRITE(stdout,'(/,A)') 'Distribution of Stream2 in I-Direction '
        DO k = kpnbeg, kpnbeg
        DO j = jpnbeg, jpnbeg
        DO i = ipnbeg, ipnend
          ijkn0  = indijk(i  ,j,k,inoff,ijnoff)
          WRITE(stdout,'(I5,4(3X,1PE15.7))') i,(mixt%cv(ivar,ijkn0),ivar=1,nvars-zcoord)
        END DO ! i
        END DO ! j 
        END DO ! k

        WRITE(stdout,'(/,A)') 'Distribution of Stream2 in J-Direction '
        DO k = kpnbeg, kpnbeg
        DO j = jpnbeg, jpnend
        DO i = ipnbeg, ipnbeg
          ijkn0  = indijk(i  ,j,k,inoff,ijnoff)
          WRITE(stdout,'(I5,4(3X,1PE15.7))') j,(mixt%cv(ivar,ijkn0),ivar=1,nvars-zcoord)
        END DO ! i
        END DO ! j
        END DO ! k

        DO ivar = 1, nvars-zcoord
          WRITE(stdout,'(/,A,I5)') 'ivar = ',ivar
          WRITE(stdout,'(/,A,2(E12.5,2X))')'Min-Max of CV ', cvmin(ivar),cvmax(ivar)
        END DO ! ivar

        DEALLOCATE(cvmin,stat=errorflag )
        global%error = errorflag
        IF (global%error /= 0) CALL errorstop( global,err_deallocate,__line__ )
        DEALLOCATE(cvmax,stat=errorflag )
        global%error = errorflag
        IF (global%error /= 0) CALL errorstop( global,err_deallocate,__line__ )
        
      END IF ! verbLevel
      
    END DO ! iReg 
  
  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_readfilestream2analyt

!******************************************************************************
!
! RCS Revision history:
!
! $Log: RVAV_ReadFileStream2Analyt.F90,v $
! Revision 1.3  2008/12/06 08:45:08  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:18:19  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 22:46:35  fnajjar
! Initial revision after changing case
!
! Revision 1.4  2003/05/15 02:57:08  jblazek
! Inlined index function.
!
! Revision 1.3  2002/10/12 03:20:51  jblazek
! Replaced [io]stat=global%error with local errorFlag for Rocflo.
!
! Revision 1.2  2002/09/10 00:01:45  f-najjar
! Variable global moved into regions()
!
! Revision 1.1  2002/07/16 22:34:16  f-najjar
! Initial Import
!
!
!******************************************************************************