APN_input_0d.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.                          *
! *********************************************************************
!                                                                             
! ---------------------------------------------------------------------------  
!                                                                              
!   SUBROUTINE  : APN_input_0d                                                  
!                                                                              
!   This subroutine read inputs for burning rate model Rocburn_1D_APN
!
!   Authors          : 
!
!   Creation Date    :  Sep. 10, 2002
!   
!   Modifications    :
!
!    No.     Date         Authors       Description
!
!                                                                              
! ---------------------------------------------------------------------------  
!                                                                              
!                                                                              
!   arguments   :                                                              
!                                                                              
!      G_APN     : Global variables for Rocburn_1D_APN
!      Indir     : directory for input data file
!                                                                              
! ---------------------------------------------------------------------------  
!
  SUBROUTINE apn_input_0d(G_APN, Indir)

    USE m_rocburn_apn_global_data

    IMPLICIT NONE
    include 'mpif.h'

!
!   Global data for Rocburn_1D_APN passed as a pointer
!

    TYPE (g_burn_1d), POINTER :: g_apn
!
!   arguments
!
    CHARACTER(*), INTENT(IN)  :: indir
    CHARACTER(*), PARAMETER :: controlfile = "RocburnAPNControl.txt"

!
!   
!
!
! ----------------------------------------------------------------
!   local variables
    CHARACTER(LEN=80) :: infile

    INTEGER   :: ir, ioerr
    INTEGER   :: mat

!
!       read propellant thermophysical properties
!

    ir = 10

    if (indir(len_trim(indir):len_trim(indir)) == '/') then
       infile= trim(indir) // controlfile
    else
       infile= trim(indir) // '/' // controlfile
    endif

    OPEN (unit=ir,file=infile,status='old')

!RAF
!RAF Extend for multiple propellants, but keep it backward compatible
!RAF Specify same value for nxmax and To for all materials.
!RAF

    DO mat = 1, matmax
      READ(ir,*,iostat=ioerr)   g_apn%a_p(mat)
      IF (ioerr /= 0) THEN
        g_apn%nmat = mat - 1
        EXIT
      ENDIF
      g_apn%nmat = mat
      READ(ir,*)   g_apn%n_p(mat)
      READ(ir,*)   g_apn%nxmax
      READ(ir,*)   g_apn%Tf_adiabatic(mat)
      READ(ir,*)   g_apn%To
      READ(ir,*,iostat=ioerr) g_apn%xmax(mat)
      IF (ioerr /= 0) THEN
        g_apn%xmax(mat) = 1.0e+15
      ENDIF
      IF (ioerr /= 0) EXIT
    END DO

    READ(ir,*,iostat=ioerr) g_apn%verbosity
    IF(ioerr /= 0) THEN
      g_apn%verbosity = 1
    ENDIF

    CLOSE(ir)

    IF(g_apn%rank .eq. 0 .AND. g_apn%verbosity .gt. 0) THEN
       WRITE(6,'(A)') 'RocburnAPN: *********** Using APN Burn Model **************'
    ENDIF
    IF(g_apn%rank .eq. 0 .AND. g_apn%verbosity .gt. 1) THEN
       WRITE(6,'(A,i3,A)') 'RocburnAPN:  Found a total of ',g_apn%nmat,' materials'
       DO mat = 1, g_apn%nmat
          WRITE(6,'(A,i3,A,f12.4)') 'RocburnAPN:  a_p(',mat,')  =',g_apn%a_p(mat)
          WRITE(6,'(A,i3,A,f12.4)') 'RocburnAPN:  n_p(',mat,')  =',g_apn%n_p(mat)
          WRITE(6,'(A,i3)') 'RocburnAPN:  nxmax= ',g_apn%nxmax
          WRITE(6,'(A,i3,A,f12.4)') 'RocburnAPN:  Tf_adiabatic(',mat,') = ',&
                                     g_apn%Tf_adiabatic(mat)
          WRITE(6,'(A,f12.4)') 'RocburnAPN:  To= ',g_apn%To
          WRITE(6,'(A,i3,A,f12.4)') 'RocburnAPN:  xmax(',mat,')  =',&
                                    g_apn%xmax(mat)
       END DO
    END IF

    IF (g_apn%verbosity.gt.2) print *,'RocburnAPN: rank=',g_apn%rank,       &
            ' done input ', controlfile

    RETURN 

  END SUBROUTINE apn_input_0d