Esempio n. 1
0
          'c30p5': 'c',
          'c30p10': 'k'}

symbols = {'c30p0': 's',
           'c20p0': 'v',
           'c40p0': '^',
           'c30p5': 'p',
           'c30p10': 'H'}
#############################################################
# First plot: Just Shen, frac diffs
#############################################################
filters = ('arealR<100',)
baseline = cstSequence(shen_cold, uniformMaxRotSlice, filters).getSeqPlot([xVar], [yVar], filters)

for script in colors.keys():
    thisSet = cstDataset(script, "HShenEOS", ye, sourceDb)
    thisSeq = cstSequence(thisSet, uniformMaxRotSlice, filters)
    thisPlot = thisSeq.getSeqPlot([xVar], [yVar], filters)

    result = reduceTwoSeqPlots(baseline, thisPlot, lambda a, b: 100 * b / a)

    plt.loglog(*result, c=colors[script],  label=script)

plt.loglog([1.e14, 3.e15], [100., 100.], ls=':', linewidth=1, c='k')
plt.xlabel(xLabel)
plt.axes().yaxis.set_minor_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
plt.axes().yaxis.set_major_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
plt.ylabel("$M_g$/$M_g(\mathrm{cold}),\, \%$", labelpad=5)
removeExponentialNotationOnAxis('y')
plt.legend()
plt.show()
Esempio n. 2
0
c0p0.getSecInstabilitySeq(0.0, 0.7, 1e15, 3e15)

exit()

sliceDict = {'a': 1.0, 'rpoe': 'min'}
filters = ()  # ('ToverW<.12', 'RedMax=0.0')#('ToverW>0.1',)

thisSeq40 = cstSequence(c40p0, sliceDict, filters)

thisSeq0 = cstSequence(c0p0, sliceDict, filters)

#thisSeq.getColumnData(['J', 'gravMass'], lambda j, m: j / (m * m), ('ToverW<.25',))
#
# guy40 = thisSeq40.getSeqPlot(['rpoe'], ['J', 'gravMass'], filters,
#                              ycolFunc=lambda j, m: j / (m * m))
# guy0 = thisSeq0.getSeqPlot(['rpoe'], ['J', 'gravMass'], filters,
#                            ycolFunc=lambda j, m: j / (m * m))
guy40 = thisSeq40.getSeqPlot(['edMax'], ['gravMass'], filters,
                             ycolFunc=lambda j: j)
guy0 = thisSeq0.getSeqPlot(['edMax'], ['gravMass'], filters,
                           ycolFunc=lambda j: j)
result = reduceTwoSeqPlots(guy0, guy40, lambda a, b: b / a - 1.0)

print result

plt.plot(*result)
plt.show()


#thisSeq = cstSequence(c30p0, {'rpoe': 1.0, 'edMax': 323076923000000.})