def initLumVphGrid(step=1,gridSize=10,doPickle=True): t=config.getTimeFromExplosion() vph=initialize.time2vph(t) mainConfigDir=config.getMainConfigDir() modelW7=initialize.readW7Data(os.path.join(mainConfigDir,'w7.combined.dat')) initDica=config.getVanillaDica() initDica['t']=t #Preparing the normalization initComp=config.getVanillaComp() initComp.update(initialize.getW7Comp(modelW7,t)) initComp=abund.setNiDecay(initComp,t) initComp=abund.setCONe(initComp) initComp=abund.normAbundances(initComp) initRunDir='init_lumvph_run/' initialVph=[0.6*vph,1.4*vph] return lumVphGrid(initialLum,initialVph,1,initRunDir,initComp=initComp,initDica=initDica,gridSize=gridSize,doPickle=doPickle)
def getInitComp(model="w7"): if model == "w7": w7Data = initialize.readW7Data() initialize.getW7Comp(w7Data)