ReadInitFlowSection.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 flow initialization.
!
! Description: None.
!
! Input: 
!   regions     Data for regions
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************
!
! $Id: ReadInitFlowSection.F90,v 1.12 2008/12/06 08:44:09 mtcampbe Exp $
!
! Copyright: (c) 2002-2006 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE readinitflowsection(regions)

  USE moddatatypes
  USE modglobal, ONLY: t_global
  USE moddatastruct, ONLY: t_region
  USE moderror
  USE modparameters
  USE modmixture, ONLY: t_mixt_input
  
  USE modinterfaces, ONLY: readsection, & 
                           readregionsection  
  
  IMPLICIT NONE

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

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

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

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

  INTEGER :: i,ireg,nvals
#ifdef RFLO
  INTEGER :: brbeg, brend
  INTEGER, PARAMETER :: nvals_max = 1
#endif
#ifdef RFLU
  INTEGER, PARAMETER :: nvals_max = 29
#endif

  CHARACTER(10) :: keys(nvals_max)
  LOGICAL :: defined(nvals_max)
  REAL(RFREAL) :: vals(nvals_max)
  TYPE(t_global), POINTER :: global

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

  global => regions(1)%global

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

! specify keywords and search for them

  nvals = nvals_max

#ifdef RFLO
  keys(1) = 'NDUMMY'
#endif
#ifdef RFLU
  keys( 1) = 'FLAG'
  keys( 2) = 'VELX'
  keys( 3) = 'VELY'
  keys( 4) = 'VELZ'
  keys( 5) = 'PRESS'
  keys( 6) = 'DENS'
  keys( 7) = 'TEMP'
  keys( 8) = 'IVAL1'
  keys( 9) = 'IVAL2'  
  keys(10) = 'IVAL3'
  keys(11) = 'IVAL4'
  keys(12) = 'IVAL5'
  keys(13) = 'IVAL6' 
  keys(14) = 'RVAL1'
  keys(15) = 'RVAL2'  
  keys(16) = 'RVAL3' 
  keys(17) = 'RVAL4'
  keys(18) = 'RVAL5'  
  keys(19) = 'RVAL6'
  keys(20) = 'RVAL7'
  keys(21) = 'RVAL8'         
  keys(22) = 'RVAL9'
  keys(23) = 'RVAL10'
  keys(24) = 'RVAL11' 
  keys(25) = 'RVAL12'
  keys(26) = 'RVAL13'  
  keys(27) = 'RVAL14'
  keys(28) = 'RVAL15'
  keys(29) = 'RVAL16'         
#endif


#ifdef RFLO
  CALL readregionsection( global,if_input,nvals,keys,vals, &
                          brbeg,brend,defined )

  IF (defined(1).eqv..true.) regions(brbeg:brend)%nDumCells = abs(vals(1)+0.5_rfreal)
#endif

#ifdef RFLU
  CALL readsection(global,if_input,nvals,keys,vals,defined )

  IF ( defined(1) .EQV. .false. ) THEN
    CALL errorstop(global,err_val_undefined,__line__,'INITFLOW-FLAG')
  ELSE
    IF ( nint(vals(1)) == initflow_fromscratch ) THEN 
      global%initFlowFlag = initflow_fromscratch
    
      DO ireg = lbound(regions,1),ubound(regions,1) 
        SELECT CASE ( regions(ireg)%mixtInput%fluidModel ) 
          CASE ( fluid_model_incomp ) 
            DO i = 2,5
              IF ( defined(i) .EQV. .false. ) THEN 
                CALL errorstop(global,err_val_undefined,__line__,keys(i))
              END IF ! defined
            END DO ! i          
          
            regions(ireg)%mixtInput%iniVelX  = vals(2)
            regions(ireg)%mixtInput%iniVelY  = vals(3)   
            regions(ireg)%mixtInput%iniVelZ  = vals(4)  
            regions(ireg)%mixtInput%iniPress = vals(5)            
          CASE ( fluid_model_comp ) 
            SELECT CASE ( regions(ireg)%mixtInput%gasModel ) 
              CASE ( gas_model_tcperf, & 
                     gas_model_mixt_tcperf, & 
                     gas_model_mixt_tperf, & 
                     gas_model_mixt_pseudo ) 
                DO i = 2,6
                  IF ( defined(i) .EQV. .false. ) THEN 
                    CALL errorstop(global,err_val_undefined,__line__,keys(i))
                  END IF ! defined
                END DO ! i          

                regions(ireg)%mixtInput%iniVelX  = vals(2)
                regions(ireg)%mixtInput%iniVelY  = vals(3)   
                regions(ireg)%mixtInput%iniVelZ  = vals(4)  
                regions(ireg)%mixtInput%iniPress = vals(5)
                regions(ireg)%mixtInput%iniDens  = vals(6) 
              CASE ( gas_model_mixt_gasliq )                         
                DO i = 2,4
                  IF ( defined(i) .EQV. .false. ) THEN 
                    CALL errorstop(global,err_val_undefined,__line__,keys(i))
                  END IF ! defined
                END DO ! i

                DO i = 6,7
                  IF ( defined(i) .EQV. .false. ) THEN 
                    CALL errorstop(global,err_val_undefined,__line__,keys(i))
                  END IF ! defined
                END DO ! i

                regions(ireg)%mixtInput%iniVelX  = vals(2)
                regions(ireg)%mixtInput%iniVelY  = vals(3)   
                regions(ireg)%mixtInput%iniVelZ  = vals(4)  
                regions(ireg)%mixtInput%iniPress = vals(5)
                regions(ireg)%mixtInput%iniTemp  = vals(7)
              CASE default
                CALL errorstop(global,err_reached_default,__line__)
            END SELECT ! regions(iReg)%mixtInput%gasModel
          CASE default
            CALL errorstop(global,err_reached_default,__line__)
        END SELECT ! regions(iReg)%mixtInput%fluidModel       
      END DO ! iReg
    ELSE IF ( nint(vals(1)) == initflow_fromfile ) THEN 
      global%initFlowFlag = initflow_fromfile
    ELSE IF ( nint(vals(1)) == initflow_fromhardcode ) THEN  
      global%initFlowFlag = initflow_fromhardcode      
    ELSE IF ( nint(vals(1)) == initflow_fromcombo_serial ) THEN
      global%initFlowFlag = initflow_fromcombo_serial
    ELSE IF ( nint(vals(1)) == initflow_fromcombo_parallel ) THEN
      global%initFlowFlag = initflow_fromcombo_parallel     
    ELSE 
      CALL errorstop(global,err_reached_default,__line__)
    END IF ! NINT(vals(1))     
  END IF ! defined

  DO ireg = lbound(regions,1),ubound(regions,1)  
    IF ( defined(8) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal1 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal1 = vals(8)
    END IF ! defined  

    IF ( defined(9) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal2 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal2 = vals(9)
    END IF ! defined  

    IF ( defined(10) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal3 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal3 = vals(10)
    END IF ! defined
    
    IF ( defined(11) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal4 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal4 = vals(11)
    END IF ! defined  
    
    IF ( defined(12) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal5 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal5 = vals(12)
    END IF ! defined  
    
    IF ( defined(13) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepIntVal6 = crazy_value_int
    ELSE
      regions(ireg)%mixtInput%prepIntVal6 = vals(13)
    END IF ! defined                

    IF ( defined(14) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal1 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal1 = vals(14)
    END IF ! defined  

    IF ( defined(15) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal2 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal2 = vals(15)
    END IF ! defined  

    IF ( defined(16) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal3 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal3 = vals(16)
    END IF ! defined 

    IF ( defined(17) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal4 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal4 = vals(17)
    END IF ! defined 

    IF ( defined(18) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal5 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal5 = vals(18)
    END IF ! defined  

    IF ( defined(19) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal6 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal6 = vals(19)
    END IF ! defined
    
    IF ( defined(20) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal7 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal7 = vals(20)
    END IF ! defined
    
    IF ( defined(21) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal8 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal8 = vals(21)
    END IF ! defined        

    IF ( defined(22) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal9 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal9 = vals(22)
    END IF ! defined        

    IF ( defined(23) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal10 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal10 = vals(23)
    END IF ! defined        

    IF ( defined(24) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal11 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal11 = vals(24)
    END IF ! defined        

    IF ( defined(25) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal12 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal12 = vals(25)
    END IF ! defined        

    IF ( defined(26) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal13 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal13 = vals(26)
    END IF ! defined        

    IF ( defined(27) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal14 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal14 = vals(27)
    END IF ! defined        

    IF ( defined(28) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal15 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal15 = vals(28)
    END IF ! defined        

    IF ( defined(29) .EQV. .false. ) THEN 
      regions(ireg)%mixtInput%prepRealVal16 = REAL(crazy_value_int,kind=rfreal)
    ELSE
      regions(ireg)%mixtInput%prepRealVal16 = vals(29)
    END IF ! defined        
  END DO ! iReg           
#endif

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

  CALL deregisterfunction(global)

END SUBROUTINE readinitflowsection

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: ReadInitFlowSection.F90,v $
! Revision 1.12  2008/12/06 08:44:09  mtcampbe
! Updated license.
!
! Revision 1.11  2008/11/19 22:17:23  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.10  2008/10/23 18:20:55  mtcampbe
! Crazy number of changes to track and fix initialization and
! restart bugs.  Many improperly formed logical expressions
! were fixed, and bug in allocation for data associated with
! the BC_INFLOWVELTEMP boundary condition squashed in
! RFLO_ReadBcInflowVelSection.F90.
!
! Revision 1.9  2007/04/05 00:56:57  haselbac
! Added additional RVALxy params for 2p shocktube problems
!
! Revision 1.8  2006/04/07 15:19:15  haselbac
! Removed tabs
!
! Revision 1.7  2006/03/26 20:21:20  haselbac
! Added TEMP input argument
!
! Revision 1.6  2005/11/17 14:37:24  haselbac
! Added more RVAL variables
!
! Revision 1.5  2005/09/13 21:36:45  haselbac
! Adapted to new init option
!
! Revision 1.4  2005/04/20 14:38:36  haselbac
! Added more int and real vals
!
! Revision 1.3  2005/03/29 22:28:31  haselbac
! Added setting of initFlowFlag for combo option
!
! Revision 1.2  2005/03/22 03:32:39  haselbac
! Added initialization of integer and real helper variables
!
! Revision 1.1  2004/12/01 16:50:26  haselbac
! Initial revision after changing case
!
! Revision 1.13  2004/11/14 19:35:42  haselbac
! Added initialization for incompressible fluid model
!
! Revision 1.12  2004/07/28 16:41:32  haselbac
! Bug fix: Initial values not assigned to all regions
!
! Revision 1.11  2004/04/08 03:15:24  wasistho
! nDummyCells in Rocflo read from INITFLOW section
!
! Revision 1.10  2003/11/20 16:40:35  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.7  2003/09/15 00:36:19  haselbac
! Added hard-code option as input
!
! Revision 1.6  2003/05/16 02:27:43  haselbac
! Removed KIND=RFREAL from NINT statements
!
! Revision 1.5  2003/05/15 02:57:02  jblazek
! Inlined index function.
!
! Revision 1.4  2003/03/25 19:15:17  haselbac
! Fixed bug in RCSIdentString
!
! Revision 1.3  2003/03/15 16:27:39  haselbac
! Added KIND qualifyer
!
! Revision 1.2  2003/02/13 22:30:54  jferry
! removed RFLU_ prefix in RegisterFunction
!
! Revision 1.1  2003/01/28 16:12:56  haselbac
! Moved here from rfluprep
!
! Revision 1.4  2002/10/27 19:23:46  haselbac
! Removed tabs
!
! Revision 1.3  2002/10/07 14:11:29  haselbac
! Removed tabs
!
! Revision 1.2  2002/09/09 16:40:13  haselbac
! global and mixtInput now under regions
!
! Revision 1.1  2002/04/11 19:13:22  haselbac
! Initial revision
!
! ******************************************************************************