PLAG_ReadDisPartnContSection.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 in user input related to discrete particle module.
!
! Description: none.
!
! Input: user input file.
!
! Output: regions = number of constituents and their properties.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: PLAG_ReadDisPartnContSection.F90,v 1.5 2008/12/06 08:44:35 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE plag_readdispartncontsection( regions )

  USE moddatatypes 
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modpartlag, ONLY    : t_plag_input
  USE modinterfaces, ONLY : readlistsection, readprefixedlistsection
  USE modinterfacesinteract, ONLY : inrt_setmaterial
  USE moderror
  USE modparameters
  USE modmaterials, ONLY  : t_material
  USE plag_modparameters
  IMPLICIT NONE

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

! ... loop variables
  INTEGER :: ival, ireg
  
! ... local variables
  CHARACTER(CHRLEN) :: rcsidentstring
  CHARACTER(15)     :: keysncont
  CHARACTER(CHRLEN), POINTER :: strvalsncont(:)

  INTEGER :: brbeg, brend, errorflag, ncols, ncont, nregions, nrows

  LOGICAL :: definedncont
    
  REAL(RFREAL), POINTER :: valsncont(:,:)
  
  TYPE(t_global),   POINTER :: global 
  TYPE(t_material), POINTER :: pmaterial  

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

  rcsidentstring = '$RCSfile: PLAG_ReadDisPartnContSection.F90,v $ $Revision: 1.5 $'

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

! Read section pertinent to constituents --------------------------------------

  definedncont = .false.
  keysncont    = 'NCONT'
!  nCols        = 1
  ncols        = 2
  
  nregions = global%nRegions
   
#ifdef RFLO    
  brbeg = 1
  brend = nregions
#endif
#ifdef RFLU
  brbeg = lbound(regions,1)
  brend = ubound(regions,1)
#endif

!  CALL ReadListSection( global, IF_INPUT,keysnCont,nCols,nRows, &
!                        valsnCont,definednCont )

  CALL readprefixedlistsection( global, if_input,keysncont,ncols,nrows, &
                                valsncont,strvalsncont,definedncont )  

  IF (definedncont) THEN
    regions(brbeg:brend)%plagInput%nCont        = nrows

    DO ireg = brbeg,brend
      ALLOCATE( regions(ireg)%plagInput%molw(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ ,'regions%plagInput%molw' ) 
      END IF ! global%error
    
      ALLOCATE( regions(ireg)%plagInput%dens(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ ,'regions%plagInput%dens' ) 
      END IF ! global%error
    
      ALLOCATE( regions(ireg)%plagInput%spht(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ ,'regions%plagInput%spht' ) 
      END IF ! global%error
    
      ALLOCATE( regions(ireg)%plagInput%surftens(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ ,'regions%plagInput%surftens' ) 
      END IF ! global%error
    
      ALLOCATE( regions(ireg)%plagInput%injcMassFluxRatio(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ , &
                        'regions%plagInput%injcMassFluxRatio' )
      END IF ! global%error
    
      ALLOCATE( regions(ireg)%plagInput%materialIndex(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ , &
                        'regions%plagInput%imaterialIndex' )
      END IF ! global%error

      ALLOCATE( regions(ireg)%plagInput%materialName(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ , &
                       'regions%plagInput%materialName' ) 
      END IF ! global%error
      
      ALLOCATE( regions(ireg)%plagInput%injcTemp(nrows),stat=errorflag )
      global%error = errorflag
      IF (global%error /= err_none) THEN
        CALL errorstop( global, err_allocate,__line__ , &
                        'regions%plagInput%injcTemp' )
      END IF ! global%error
    
    END DO !iReg
    
    DO ireg= brbeg,brend
      ncont = regions(ireg)%plagInput%nCont   
      DO ival=1, ncont  
        CALL inrt_setmaterial(global,pmaterial,strvalsncont(ival))
        
        regions(ireg)%plagInput%materialIndex(ival) = pmaterial%index
        regions(ireg)%plagInput%materialName(ival)  = pmaterial%name

        regions(ireg)%plagInput%molw(ival) = pmaterial%molw
        regions(ireg)%plagInput%dens(ival) = pmaterial%dens
        regions(ireg)%plagInput%spht(ival) = pmaterial%spht
        regions(ireg)%plagInput%surftens(ival) = pmaterial%surftens
        regions(ireg)%plagInput%injcMassFluxRatio(ival) = abs(valsncont(ival,1))
        regions(ireg)%plagInput%injcTemp(ival)          = abs(valsncont(ival,2))
      ENDDO ! ival
    ENDDO ! iReg
    
  ELSE
    CALL errorstop( global, err_missing_value,__line__) 
  ENDIF ! definednCont  

! deallocate pointer arrays

  DEALLOCATE( strvalsncont, stat=errorflag ) 
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop( global, err_deallocate,__line__ ,'strValsnCont' ) 
  END IF ! global%error
  
  DEALLOCATE( valsncont, stat=errorflag ) 
  global%error = errorflag
  IF (global%error /= err_none) THEN
    CALL errorstop( global, err_deallocate,__line__ ,'valsnCont' ) 
  END IF ! global%error
   
! finalize --------------------------------------------------------------------

  CALL deregisterfunction( global )

END SUBROUTINE plag_readdispartncontsection

!******************************************************************************
!
! RCS Revision history:
!
! $Log: PLAG_ReadDisPartnContSection.F90,v $
! Revision 1.5  2008/12/06 08:44:35  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:48  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2006/09/18 20:37:35  fnajjar
! Added injcTemp in particle input
!
! Revision 1.2  2005/12/01 18:39:15  fnajjar
! Added error trap for missing input value
!
! Revision 1.1  2004/12/01 20:58:04  fnajjar
! Initial revision after changing case
!
! Revision 1.7  2004/02/26 21:02:20  haselbac
! Added RFLU support
!
! Revision 1.6  2003/09/26 21:48:20  fnajjar
! Changed interface call for INRT_SetMaterial to ModInterfacesInteract
!
! Revision 1.5  2003/09/13 20:14:22  fnajjar
! Added infrastructure for Breakup model
!
! Revision 1.4  2003/03/24 23:28:41  jferry
! converted SetMaterial to INRT_SetMaterial
!
! Revision 1.3  2003/03/12 21:21:46  fnajjar
! Use Material datastructure for Particle properties
!
! Revision 1.2  2003/01/10 17:03:41  f-najjar
! Bug fix to read properly input data for all regions
!
! Revision 1.1  2002/10/25 14:19:16  f-najjar
! Initial Import of Rocpart
!
!
!******************************************************************************