! Fiducial model.
 
&lc_modelpar
afile=lc_model0.dat
! afile=lc_model0_fine.dat
alpha=1.d0
beta=1.d0
irad=10          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44, 10 is with Si-32.
itrap=1         ! 0 for increasing partial escape;  1 for complete trapping.
t_del=1.d0
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
! t_max=40.d0
/

&lc_modelpar1
afile=lc_model1.dat
alpha=1.d0
beta=1.d0
irad=2          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44.
itrap=0         ! 0 for increasing partial escape;  1 for complete trapping.
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_del=1.d0
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
/


&lc_modelpar1
afile=lc_model2.dat
alpha=1.d0
beta=1.d0
irad=4          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44.
itrap=0         ! 0 for increasing partial escape;  1 for complete trapping.
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_del=1.d0
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
/

&lc_modelpar1
afile=lc_model3.dat
alpha=1.d0
beta=1.d0
irad=6          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44.
itrap=0         ! 0 for increasing partial escape;  1 for complete trapping.
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_del=1.d0
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
/


&lc_modelpar1
afile=lc_model4.dat
alpha=1.d0
beta=1.d0
irad=10          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44.
itrap=0         ! 0 for increasing partial escape;  1 for complete trapping.
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_del=1.d0
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
/

&lc_modelpar1
afile=lc_model5.dat
alpha=1.d0
beta=1.d0
irad=10          ! 2 is for just Ni-56, 4 is with Co-57, 6 is with Fe55, 8 is with Ti-44, 10 for Si-32, 12 for Ni-59, 14 for Fe-60, 16 for Mn-53.
itrap=0         ! 0 for increasing partial escape;  1 for complete trapping.
! t_del=.1d0    ! Fine, but maybe too fine for the SM algorithm for 4000. days
t_del=1.d0
t_lcmax=19.5d0  ! conley2006 mean Bmax rise time.
t_max=6000.d0
/




! Sort of Fitted model.  Which probably no use.

&lc_modelpar1
alpha=3.d0
beta=3.d0
t_del=.1d0
t_lcmax=15.d0  ! conley2006 mean Bmax rise time.
t_max=500.d0
/