Rocburn_PY/test/rocburn_2D.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.                          *
! *********************************************************************
!     ------------------------------------------------------------------------

  PROGRAM burn_initialize

    USE m_rocburn_1d_py
    IMPLICIT NONE

    INTEGER :: nxmax,i
    REAL(DBL) :: to_read,p, patm, tflame,dzero,time
    REAL(DBL) :: toa, fr,rhoc,rhoc_cgs,rb, dt,p_coor(3),fact
    REAL(DBL),ALLOCATABLE :: tn(:),tnp1(:)
    INTEGER :: bflag,ntime_steps
    CHARACTER (LEN=80) pass_dir
    TYPE(parameter_structure),POINTER :: pass_in_out

    
     pass_dir = './'

!----READ INPUT FILE AND GENERATE GRID
     CALL burn_init_0d(pass_in_out,1,trim(pass_dir), nxmax, to_read)

!----ALLOCATION
     ALLOCATE(tn(nxmax))
     ALLOCATE(tnp1(nxmax))

!----USE a burning interface
     bflag =1   
     dzero = zero

!----SET Patm
     patm = 30.0
     p = patm /pa2atm

!----SET rhoc, in genx is passed by Rocfrac
     rhoc_cgs = 1.7026
     rhoc = rhoc*gcc2kgmc

!----INITIALIZE temperature profile
     CALL burn_init_1d(pass_in_out, bflag, &
                    p, to_read, rhoc, p_coor, rb, &
                    toa, fr, tn(:), tflame)
  
     time = 0 
     ntime_steps = 100000
     dt = 1.0e-6

!----RESET P if starting from steady conditions
     if(bflag ==1) then
        fact = 2.0d0
        p = p*fact
     endif

!----MAIN TIME LOOP
     do i = 1,ntime_steps
        time = time+dt
        CALL burn_get_burning_rate1d(pass_in_out, dt, p, &
             to_read, tn(:),   &
             dzero, dzero,     &
             dzero, dzero,     &
             rhoc, dzero,      &
             rb, dzero,        &
             bflag, tnp1(:), tflame)
        tn = tnp1
        print*,i,time,tn(1)
     enddo
 
  stop
END PROGRAM burn_initialize
!*****************************************************************************