ReadTimeZoomingSection.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 timezooming.
!
! Description: none.
!
! Input: user input file.
!
! Output: global = timezooming parameters
!
! Notes: none.
!
!******************************************************************************


SUBROUTINE readtimezoomingsection( global )

  USE moddatatypes
  USE modglobal, ONLY     : t_global
  USE modinterfaces, ONLY : readsection
  USE moderror
  USE modparameters
  USE modmpi

  IMPLICIT NONE

! ... parameters
  TYPE(t_global), POINTER :: global

! ... local variables
  INTEGER :: nvals
  INTEGER, PARAMETER :: nvals_max = 9

  CHARACTER(10) :: keys(nvals_max)

  LOGICAL :: defined(nvals_max)

  REAL(RFREAL) :: vals(nvals_max)

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

  CALL registerfunction( global,'ReadTimeZoomingSection',&
  'ReadTimeZoomingSection.F90' )
! specify keywords and search for them

  nvals = nvals_max

  keys(1) = 'MINPLANE'
  keys(2) = 'MAXPLANE'
  keys(3) = 'NOZINLET'
  keys(4) = 'NOZRAD'
  keys(5) = 'AXIS'
  keys(6) = 'THROATRAD'
  keys(7) = 'N'
  keys(8) = 'RHO'
  keys(9) = 'A'

  
  CALL readsection( global,if_input,nvals,keys,vals,defined )

  IF (defined( 1)) THEN
     global%tzMinPlane        = vals( 1)
  ELSE
     global%tzMinPlane = -9d10
  ENDIF

  IF (defined( 2)) THEN
     global%tzMaxPlane        = vals( 2)
  ELSE
     global%tzMaxPlane = 0.85_rfreal
  ENDIF

  IF (defined( 8)) THEN
     global%tzRhos        = vals( 8)
  ELSE
     global%tzRhos = 1702.68_rfreal
  ENDIF

  IF (defined( 9)) THEN
     global%tzA           = vals( 9)
  ELSE
     global%tzA = .000003789917_rfreal
  ENDIF

  IF (defined( 7)) THEN
     global%tzN           = vals( 7)
  ELSE
     global%tzN = 0.461_rfreal
  ENDIF

  IF (defined( 6)) THEN
     global%tzThroatRad   = vals( 6)
  ELSE
     global%tzThroatRad = 0.013208_rfreal
  ENDIF

  IF (defined( 5)) THEN
     IF(vals(5) == 1.0) THEN
        global%tzCoordLong   = xcoord
        global%tzCoordTrans1 = ycoord
        global%tzCoordTrans2 = zcoord
     ELSE IF (vals(5) == 2.0) THEN
        global%tzCoordLong   = ycoord
        global%tzCoordTrans1 = zcoord
        global%tzCoordTrans2 = xcoord
     ELSE
        global%tzCoordLong   = zcoord
        global%tzCoordTrans1 = xcoord
        global%tzCoordTrans2 = ycoord
     ENDIF
  ELSE
     global%tzCoordLong   = xcoord
     global%tzCoordTrans1 = ycoord
     global%tzCoordTrans2 = zcoord
  ENDIF

  global%tzSubNoz = .false.
  IF (defined(3)) THEN
     global%tzNozInlet   = vals( 3)
     global%tzSubNoz = .true.
  ELSE
     global%tzNozInlet = global%tzMaxPlane
  ENDIF

  IF (defined(4)) THEN
     global%tzNozRad   = vals(4)
     global%tzSubNoz = .true.
  ELSE
     global%tzNozRad = global%tzThroatRad
  ENDIF

! finalize

  CALL deregisterfunction( global )

END SUBROUTINE readtimezoomingsection