def __init__(self, initDica=None, t=None, mode='init'):
if t == None and initDica == None:
assert SNConfigDict != None
t = SNConfigDict['t']
else: t = t
self.data={ 'chl': 25.0,
'e_b-v_gal': -1.,
'e_b-v_host': -1.,
'em_high': 6500.0,
'em_low': 1000.0,
'grid': 1.25,
'inthigh': 10000.0,
'intlow': 2500.0,
'itt': 4.0,
'js': 20.0,
'kb': 501.0,
'kr': 250819801106,
'lg_tau': -3.5,
'log_l_low_high': 8.9399999999999995,
'log_lbol': -1.0,
'm-m': -1.,
'mb': 10000.0,
'mu': 1500.0,
'nc': 5.0,
'np5': 4.0,
'options': '1 1 1 1 1',
't': t,
'tb': 10000.0,
'wl': 0.25,
'z': -1.,
'xe1':0.20}
if mode=='init':
self.data['v_ph'] = initialize.time2vph(t)
self.tRise=19.5
self._lockTemp=False
self._lockVph=False
self._lockLum=False
#self.data.update(config.getInitDicaParam(tRise=self.tRise))
if initDica!=None:
self.data.update(initDica)
if SNConfigDict != None:
self.data['e_b-v_gal'] = SNConfigDict['extgal']
self.data['e_b-v_host'] = SNConfigDict['exthost']
self.data['m-m'] = SNConfigDict['mu']
self.data['z'] = SNConfigDict['z']
else:
pass
#pdb.set_trace()
#print "Warning: Dica initialized with non-sensible default values. This might be problematic."
assert self.data['t'] != None
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)
