ReadRocketSection.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 rockets
!
! Description: none.
!
! Input: user input file.
!
! Output: global = constraint params
!
! Notes: none.
!
!******************************************************************************


SUBROUTINE readrocketsection( 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 = 15

  CHARACTER(10) :: keys(nvals_max)
  LOGICAL :: hasacase
  LOGICAL :: defined(nvals_max)
  REAL(RFREAL) :: vals(nvals_max)

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

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

  nvals = nvals_max

  keys(1)  =  'CASERAD'
  keys(2)  =  'HEADEND'
  keys(3)  =  'AFTEND'
  keys(4)  =  'COORDL'
  keys(5)  =  'TOL1'
  keys(6)  =  'TOL2'
  keys(7)  =  'ELLIPSL'
  keys(8)  =  'ELLIPST'
  keys(9)  =  'NOZY'
  keys(10) =  'LMIN'
  keys(11) =  'LMAX'
  keys(12) =  'T1MIN'
  keys(13) =  'T1MAX'
  keys(14) =  'T2MIN'
  keys(15) =  'T2MAX'
  
  CALL readsection( global,if_input,nvals,keys,vals,defined )

  IF (defined( 1)) THEN
     hasacase = .true.
     global%cnstrCaseRad        = vals( 1)
  ELSE
     global%cnstrCaseRad = -999999_rfreal
  ENDIF

  IF (defined( 2)) THEN
     global%cnstrHeadEnd = vals( 2)
  ELSE
     global%cnstrHeadEnd = -9999999_rfreal
  ENDIF

  IF (defined( 3)) THEN
     global%cnstrAftEnd = vals( 3)
  ELSE
     global%cnstrAftEnd = 9999999_rfreal
  ENDIF

  IF (defined( 4)) THEN
     IF(vals( 4) == 1.0) THEN
        global%cnstrCoordL  = xcoord
        global%cnstrCoordT1 = ycoord
        global%cnstrCoordT2 = zcoord
     ELSE IF (vals(4) == 2.0) THEN
        global%cnstrCoordL  = ycoord
        global%cnstrCoordT1 = zcoord
        global%cnstrCoordT2 = xcoord
     ELSE
        global%cnstrCoordL  = zcoord
        global%cnstrCoordT1 = xcoord
        global%cnstrCoordT2 = ycoord
     ENDIF
  ELSE
     global%cnstrCoordL  = xcoord
     global%cnstrCoordT1 = ycoord
     global%cnstrCoordT2 = zcoord
  ENDIF

  IF (defined( 5)) THEN
     global%cnstrTol1 = vals( 5)
  ELSE
     global%cnstrTol1 = 0.000001_rfreal
  ENDIF
  IF (defined( 6)) THEN
     global%cnstrTol2 = vals(6)
  ELSE
     global%cnstrTol2 = 0.00000001_rfreal
  ENDIF

  IF (defined( 7)) THEN
     global%cnstrEllipsL   = vals(7)
  ELSE
     global%cnstrEllipsL = 1.0_rfreal
  ENDIF

  IF (defined( 8)) THEN
     global%cnstrEllipsT   = vals(8)
  ELSE
     global%cnstrEllipsT = global%cnstrCaseRad
  ENDIF

  IF (defined( 9)) THEN
     global%cnstrNozY   = vals(9)
  ELSE
     global%cnstrNozY = 1.0_rfreal
  ENDIF
  IF (defined( 10)) THEN
     global%cnstrLMinPlane   = vals(10)
  ELSE
     global%cnstrLMinPlane = -9999999_rfreal
  ENDIF
  IF (defined( 11)) THEN
     global%cnstrLMaxPlane   = vals(11)
  ELSE
     global%cnstrLMaxPlane = 9999999_rfreal
  ENDIF
  IF (defined( 12)) THEN
     global%cnstrT1MinPlane   = vals(12)
  ELSE
     global%cnstrT1MinPlane = -9999999_rfreal
  ENDIF
  IF (defined( 13)) THEN
     global%cnstrT1MaxPlane   = vals(13)
  ELSE
     global%cnstrT1MaxPlane = 9999999_rfreal
  ENDIF
  IF (defined( 14)) THEN
     global%cnstrT2MinPlane   = vals(14)
  ELSE
     global%cnstrT2MinPlane = -9999999_rfreal
  ENDIF
  IF (defined( 15)) THEN
     global%cnstrT2MaxPlane   = vals(15)
  ELSE
     global%cnstrT2MaxPlane = 9999999_rfreal
  ENDIF

     
! finalize

  CALL deregisterfunction( global )

END SUBROUTINE readrocketsection