primary stops evolving after becoming naked helium HG star.

[poojanagrawal]

I am trying to create some BH-star systems using the cosmic's develop branch. However, I encountered many systems where the primary star stops evolving halfway after becoming a naked helium Hertzsprung-Gap star.

For example,

	tphys	mass_1	mass_2	kstar_1	kstar_2	sep	evol_type
0	0.000000	30.460000	15.000000	1	1	6.728360e+02	1
0	5.884438	27.937086	14.961129	2	1	7.130208e+02	2
0	5.889382	27.918930	14.961595	2	1	4.992910e+02	3
0	5.889382	27.928196	14.961595	2	1	6.265277e-13	5
0	5.889382	27.928196	14.961595	2	1	6.265277e-13	7
0	5.889382	0.000000	14.961595	8	1	0.000000e+00	8
0	12.677995	0.000000	14.596350	8	2	0.000000e+00	2
0	12.704624	0.000000	14.587939	8	4	0.000000e+00	2
0	14.144233	0.000000	10.426111	8	5	0.000000e+00	2
0	14.199501	0.000000	10.001305	8	5	0.000000e+00	16
0	14.199501	0.000000	1.620586	8	13	0.000000e+00	2
0	13700.000000	0.000000	1.620586	8	13	0.000000e+00	10

Here are the IC and BSEDict for reproducing the above system:

single_binary = InitialBinaryTable.InitialBinaries(m1=30.46, m2=15.0, porb=300, ecc=0.3, tphysf=13700.0, kstar1=1, kstar2=1, metallicity=0.02)

BSEDict = {'xi': 1.0, 'bhflag': 0, 'neta': 0.5, 'windflag': 3, 'wdflag': 1, 'alpha1': [1.0,1.0], \ 'pts1': 0.001, 'pts3': 0.02, 'pts2': 0.01, 'epsnov': 0.001, 'hewind': 0.5, 'ck': 1000,\ 'bwind': 0.0, 'lambdaf': 0.0, 'mxns': 3.0, 'beta': -1.0, 'tflag': 1, 'acc2': 1.5, \ 'grflag' : 1, 'remnantflag': 4, 'ceflag': 0, 'eddfac': 1.0, 'ifflag': 0, 'bconst': 3000,\ 'sigma': 265.0, 'gamma': -2.0, 'pisn': 45.0, 'natal_kick_array' : [[-100.0,-100.0,-100.0,-100.0,0.0], [-100.0,-100.0,-100.0,-100.0,0.0]],\ 'bhsigmafrac' : 1.0, 'polar_kick_angle' : 90, 'qcrit_array' : [0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0], \ 'cekickflag' : 2, 'cehestarflag' : 0, 'cemergeflag' : 0, 'ecsn' : 2.25, 'ecsn_mlow' : 1.6, 'aic' : 1, 'ussn' : 0, 'sigmadiv' :-20.0, \ 'qcflag' : 1, 'eddlimflag' : 0, 'fprimc_array' : [2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0], \ 'bhspinflag' : 0, 'bhspinmag' : 0.0, 'rejuv_fac' : 1.0, 'rejuvflag' : 0, 'htpmb' : 1,\ 'ST_cr' : 1, 'ST_tide' : 1, 'bdecayfac' : 1, 'rembar_massloss' : 0.5, 'kickflag' : -1, \ 'zsun' : 0.02, 'bhms_coll_flag' : 0, 'don_lim' : -1, 'acc_lim' : [-1,-1],'rtmsflag' : 0,'wd_mass_lim' : 1}

SSEDict = {'stellar_engine': 'sse' }