search

ClusterMonteCarlo / CMC-COSMIC

14 stars 19 forks source link

BSE setting #9

Closed scottcoughlin2014 closed 3 years ago

scottcoughlin2014 commented 4 years ago 
pts1=0.001

can lead NaN radius

scottcoughlin2014 commented 4 years ago

The condition that gives a Nan

kstar_1,kstar_2,mass_1,mass_2,porb,ecc,metallicity,tphysf,mass0_1,mass0_2,rad_1,rad_2,lum_1,lum_2,massc_1,massc_2,radc_1,radc_2,menv_1,menv_2,renv_1,renv_2,omega_spin_1,omega_spin_2,B_1,B_2,bacc_1,bacc_2,tacc_1,tacc_2,epoch_1,epoch_2,tms_1,tms_2,bhspin_1,bhspin_2,tphys
7,1,1.1761,5.69906,0.780364,0,0.001,73.3077,1.1761,5.69906,0.239908,3.45774,454.242,1936.42,0,0,0,0,1.1761,1e-10,0.239908,1e-10,10.4184,2925.68,0,0,0,0,0,0,58.2163,4.37994,14.3918,70.7224,0,0,59.669
import pandas
from cosmic.sample.initialbinarytable import InitialBinaryTable

from cosmic.evolve import Evolve

single_binary =  pandas.read_csv("initc.csv")
single_binary["binfrac"] = 0
bpp, bcm, initC, kick_info = Evolve.evolve(initialbinarytable=single_binary, params="Params.ini")
print(bpp)
carlrodriguez commented 3 years ago

Wait is this a BSE issue or a CMC issue? Would love to move this over to the COSMIC page...

carlrodriguez commented 3 years ago

Fixed by changing the default back to pts1=0.05 in CMC