&ppn_frame

!       T9 = 0.25D0               ! constant temperature
        T9 = 0.2D0                ! constant temperature
!       T9 = 0.3D0                ! constant temperature
        RHO = 1.0D4               ! /cgs
!       RHO = 1.995D2             ! /cgs

        DT = 5.63D-11              ! initial time step
        T_MAX = 3.16D-03
        DT_MAX = 3.16D-5             ! max number dzeitj can be
        dt_factor = 1.00875d0

!       DT = 1.0D-10 ! 1.0D-8               ! initial time step
!       DT = 1.0D-12 ! 1.0D-8               ! initial time step
!       T_MAX = 2.00000D-03 ! 2.51622D-05 ! 1.62134e-4        ! max age of evolution John
!       T_MAX = 2.0D-1 !  7.08407D-03       ! max age of evolution, 2nd peak (CEMP)
!       T_MAX = 1D-01 ! 2.06633D-02 ! 2.12557D-02 ! 2nd peak, denn3.2d14, 0.2T9
!       T_MAX = 1.28254D-01 ! 2nd peak, denn3.2d13, 0.2T9
!       T_MAX = 1.48254D-01 ! 2nd peak, denn3.16d13, 0.2T9
!       T_MAX = 2.06633D-02 ! 2nd peak, denn3.16d14, 0.2T9
!       T_MAX = 2.06633D-02 ! 2nd peak, denn2.0d15, 0.2T9
!       T_MAX = 2.06633D-02 ! 2nd peak, denn1.0d18, 0.2T9
!       T_MAX = 5D-05
!       T_MAX = 1D-01 ! 2.1256D-02 ! 2nd peak, denn3.16d14, 0.3T9
!       T_MAX = 5.0000D-02 ! 2nd peak, denn3.16d13, 0.3T9
!       T_MAX = 1.00000D+00 ! pure C, Sam and Ondrea project
!       T_MAX = 1.20000D-01 ! pure CNO, Sam and Ondrea project
!       T_MAX = 2.0D+00 ! 1D+00 ! 1.80000D-01 ! pure CNO, Sam and Ondrea project
!       DT_MAX = 1.D-3 ! 1.D-4 ! 1.D-3            ! max number dzeitj can be
!       dt_factor = 1.025d0 ! 1.010d0 ! 1.025d0
        NSOURCE = 1               ! source of outside trajectory
        IABUINI = 100             ! how to initialize
!       ini_filename = "p_ini.dat.CPb"
!       ini_filename = "p_ini.dat.pureCO"
!       ini_filename = "p_ini.dat.pureCNO"
!       ini_filename = "p_ini.dat.pureC"
        ini_filename = "p_ini.dat"
!       ini_filename = "p_ini.dat.2nd_peak_X0"
        IOLEVEL =  0 ! 1          ! controls input/output, 0 prints only the last cycle
        PLIST = 0                 ! new in progress
        CPREFIX = 'iso_massf'     ! prefix of abundance vector output
        iplot_flux_option = 0     ! 0: no print fluxes; 1: print fluxes in flux_*.DAT;

        i_denn_const = 1          ! 1: use constant neutron number density (denn)
        denn = 3.16d+16            ! constant netron number density used when iden_const = 1

        decay_time = 1.0d+09      ! decay time (in years) for abundances used at the end
!       decay_time = 0.0d+00      ! decay time (in years) for abundances used at the end

/

! *** end of namelist

   ! the parameter values are given as fortran95 "namelists"
   ! here's a brief summary of how things work with namelists.
   ! each line in a namelist is of the form   name = value ! comment
   ! each name must be declared as a control in the program,
   ! but you don't need to include all of the controls in the list.
   ! and the order of the list doesn't matter.
   ! blank lines and comment lines can be freely included in the list.
   ! blanks at the start of a line containing a name-value pair are okay too.
   ! all of the controls are given reasonable default values at initialization.
   ! so you only need to set the ones that you actually want to change
   ! if you don't want to change a control, but you do want to keep in the list
   ! either comment out the line, or leave the value field blank.
   ! fortran supports structured values in namelists (arrays and derived types)
   ! but for simplicity, we only use basic types (integers, doubles, logicals, strings).
   ! values are given using the usual fortran syntax.
   ! for example, character string values are enclosed in quotes
   ! ****** REMINDER ******
   ! A '/' character terminates the namelist.
   ! So, do not use '/' outside of comments or strings.
   ! In particular, don't write something like 1d0/3d0 as a value.