def plot_spread_mc_nomext(m, color='green', label='Spread Rate'):
"""Tool for plotting spread rate vs MC for each realization."""
MC = pl.arange(0., 0.51, .01)
y_trace = []
sigList = []
hocList = []
mExtList = []
for modSig in m.modSig.trace():
sigList.append(modSig)
for hoc in m.hoc.trace():
hocList.append(hoc)
for i in range(len(sigList)):
y = roth.flameSpread(sigma=sigList[i],
bulk=data.meanBulk,
hoc=hocList[i],
moist=MC,
mExt=0.3)
pl.plot(MC, y, color='gray', alpha=.75, zorder=-1)
y_trace.append(y)
pl.plot(data.mcs, data.rates, 'bo', label='Experiment')
pl.plot(MC,
pl.mean(y_trace, axis=0),
color=color,
linewidth=5,
label=label)
decorate_plot()
def plot_spread_mc_nomext(m, color='green', label='Spread Rate'):
"""Tool for plotting spread rate vs MC for each realization."""
MC = pl.arange(0., 0.51, .01)
y_trace = []
sigList = []
hocList = []
mExtList = []
for modSig in m.modSig.trace():
sigList.append(modSig)
for hoc in m.hoc.trace():
hocList.append(hoc)
for i in range(len(sigList)):
y = roth.flameSpread(sigma=sigList[i], bulk=data.meanBulk, hoc=hocList[i], moist=MC, mExt=0.3)
pl.plot(MC, y, color='gray', alpha=.75, zorder=-1)
y_trace.append(y)
pl.plot(data.mcs, data.rates, 'bo', label='Experiment')
pl.plot(MC, pl.mean(y_trace, axis=0), color=color, linewidth=5, label=label)
decorate_plot()
def y_mean(modSig=modSig, mBulk=data.meanBulk, hoc=hoc, moistC=data.mcs, moistE=mExt):
return roth.flameSpread(sigma=(modSig+100.), bulk=mBulk, hoc=hoc, moist=moistC, mExt=moistE)
def y_mean(modSig=modSig, mBulk=data.meanBulk, hoc=hoc, moistC=moistC, moistE=0.3):
return roth.flameSpread(sigma=modSig, bulk=mBulk, hoc=hoc, moist=moistC, mExt=moistE)