def energy_profile(self,ixaxis):
mass = self.get('<STR_LIT>')<EOL>radius = self.get('<STR_LIT>') * ast.rsun_cm<EOL>eps_nuc = self.get('<STR_LIT>')<EOL>eps_neu = self.get('<STR_LIT>')<EOL>if ixaxis == '<STR_LIT>':<EOL><INDENT>xaxis = mass<EOL>xlab = '<STR_LIT>'<EOL><DEDENT>else:<EOL><INDENT>xaxis = old_div(radius, <NUM_LIT>) <EOL>xlab = '<STR_LIT>'<EOL><DEDENT>pl.plot(xaxis, np.log10(eps_nuc),<EOL>'<STR_LIT>',<EOL>label='<STR_LIT>')<EOL>pl.plot(xaxis, np.log10(-eps_nuc),<EOL>'<STR_LIT>',<EOL>label='<STR_LIT>')<EOL>pl.plot(xaxis, np.log10(eps_neu),<EOL>'<STR_LIT>',<EOL>label='<STR_LIT>')<EOL>pl.xlabel(xlab)<EOL>pl.ylabel('<STR_LIT>')<EOL>pl.legend(loc='<STR_LIT>').draw_frame(False)<EOL>
Plot radial profile of key energy generations eps_nuc,
eps_neu etc.

Parameters
----------
ixaxis : 'mass' or 'radius'