PEUL_ReadInputFile.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 for PEUL information (done on all processors).
!
! Description: none.
!
! Input: user input file.
!
! Output: regions = PEUL information
!
! Notes: none.
!
!******************************************************************************
!
! $Id: PEUL_ReadInputFile.F90,v 1.3 2008/12/06 08:44:39 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE peul_readinputfile( regions )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region
  USE modglobal,     ONLY : t_global
  USE moderror
  USE modparameters
  USE peul_modparameters

  USE peul_modinterfaces, ONLY : peul_readconpartsection, &
                                 peul_readconpartptypesection
  IMPLICIT NONE

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

! ... loop variables
  INTEGER :: iread, ireg

! ... local variables
  CHARACTER(CHRLEN)   :: rcsidentstring
  CHARACTER(CHRLEN+4) :: fname
  CHARACTER(256)      :: line

  LOGICAL :: usedsomewhere, unusedsomewhere

  INTEGER :: errorflag, iptype, nptypes, nsecs, brbeg, brend, readstatus

  TYPE(t_global), POINTER :: global

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

  rcsidentstring = '$RCSfile: PEUL_ReadInputFile.F90,v $ $Revision: 1.3 $'

  global => regions(1)%global

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

! begin -----------------------------------------------------------------------

! search for CONPART and CONPART_PTYPE sections

  fname = trim(global%inDir)//trim(global%casename)//'.inp'

  DO iread = 1,2

! - open file

    OPEN(if_input,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 file looking for keywords

    SELECT CASE (iread)

    CASE (1) ! on first pass, find # of sections and Ptypes

      nsecs = 0 ! initialize counter for multiply referenced sections

      DO
        READ(if_input,'(A256)',err=10,end=86) line

        SELECT CASE(trim(line))

        CASE ('# CONPART')
          IF (nsecs == 1) THEN
            nptypes = iptype
          ELSE IF (nsecs > 1) THEN
            IF (iptype /= nptypes) & ! nPtypes same in each sec?
              CALL errorstop( global,err_peul_nptypes,__line__ )
          END IF ! nSecs
          nsecs = nsecs + 1
          iptype = 0

        CASE ('# CONPART_PTYPE')
          IF (nsecs == 0) &
            CALL errorstop( global,err_peul_ptype,__line__ )
          iptype = iptype + 1

        END SELECT ! line
      END DO

86    CONTINUE

      IF (nsecs == 1) THEN
        nptypes = iptype
      ELSE IF (nsecs > 1) THEN
        IF (iptype /= nptypes) & ! nPtypes same in each sec?
          CALL errorstop( global,err_peul_nptypes,__line__ )
      END IF ! nSecs

    CASE (2) ! on second pass, search for real- and string-valued keys

      DO
        READ(if_input,'(A256)',err=10,end=87) line

        SELECT CASE(trim(line))

        CASE ('# CONPART')
          CALL peul_readconpartsection( regions,nptypes,brbeg,brend )
          iptype = 0

        CASE ('# CONPART_PTYPE')
          iptype = iptype + 1
          CALL peul_readconpartptypesection( regions,brbeg,brend,iptype )

        END SELECT ! line
      END DO

87    CONTINUE

    END SELECT ! iRead

! close file ------------------------------------------------------------------

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

    IF (nsecs == 0) EXIT

  END DO ! iRead

! set global%peulUsed -------------------------------------------------------

  usedsomewhere   = .false.
  unusedsomewhere = .false.

  DO ireg = lbound(regions,1),ubound(regions,1)
    readstatus = regions(ireg)%peulInput%readStatus
    usedsomewhere   = usedsomewhere  .OR.(readstatus == 1)
    unusedsomewhere = unusedsomewhere.OR.(readstatus /= 1) ! == 0 or -1
  END DO ! iReg

  IF (usedsomewhere.AND.unusedsomewhere) THEN
    CALL errorstop( global,err_mp_allornone,__line__ )
  END IF ! usedSomewhere.AND.unusedSomewhere

  global%peulUsed = usedsomewhere

! finalization & error handling -----------------------------------------------

  CALL deregisterfunction( global )
  goto 999

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

999 CONTINUE

END SUBROUTINE peul_readinputfile

!******************************************************************************
!
! RCS Revision history:
!
! $Log: PEUL_ReadInputFile.F90,v $
! Revision 1.3  2008/12/06 08:44:39  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:17:51  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 21:09:47  haselbac
! Initial revision after changing case
!
! Revision 1.7  2004/03/05 22:09:04  jferry
! created global variables for peul, plag, and inrt use
!
! Revision 1.6  2003/04/14 16:32:24  jferry
! minor edits
!
! Revision 1.5  2003/03/17 18:42:45  jblazek
! Added inDir to path of the input file.
!
! Revision 1.4  2003/03/11 16:04:57  jferry
! Created data type for material properties
!
! Revision 1.3  2003/02/26 23:38:30  jferry
! eliminated end=999 convention to ensure that files get closed
!
! Revision 1.2  2003/02/12 23:34:48  jferry
! Replaced [io]stat=global%error with local errorFlag
!
! Revision 1.1  2003/02/11 22:52:51  jferry
! Initial import of Rocsmoke
!
!******************************************************************************