def plot_luminosity_equation_exact(self, LAXIS, xbl, xbr, ybu, ybd, ilg):
    'Plot luminosity equation in the model'
    grd1 = self.xzn0
    rhs0 = self.minus_cp_rho_dTdt
    rhs1 = self.plus_delta_dPdt
    rhs2 = self.minus_dd_div_eiui
    rhs3 = self.plus_tke_diss
    res = self.minus_resLumExactEquation
    plt.figure(figsize=(7, 6))
    plt.gca().yaxis.get_major_formatter().set_powerlimits((0, 0))
    to_plot = [rhs0, rhs1, res]
    self.set_plt_axis(LAXIS, xbl, xbr, ybu, ybd, to_plot)
    self.bconv = 400000000.0
    self.tconv = 1200000000.0
    xlimitrange = np.where(((grd1 > self.bconv) & (grd1 < self.tconv)))
    xlimitbottom = np.where((grd1 < self.bconv))
    xlimittop = np.where((grd1 > self.tconv))
    plt.title('standard luminosity equation exact')
    if (self.ig == 1):
        plt.plot(grd1[xlimitrange], rhs0[xlimitrange], color='#FF8C00', label='$-c_P \\overline{\\rho \\partial_t T}$')
        plt.plot(grd1[xlimitrange], rhs1[xlimitrange], color='y', label='$+\\delta \\overline{\\partial_t P}$')
        plt.plot(grd1[xlimitrange], rhs2[xlimitrange], color='r', label='$-\\overline{\\rho \\nabla \\cdot \\epsilon_I {\\bf u}}$')
        plt.plot(grd1[xlimitrange], rhs3[xlimitrange], color='g', label='$+\\varepsilon_K$')
        plt.plot(grd1, res, color='k', linestyle='--', label='res $\\sim N$')
        zeros = np.zeros(self.nx)
        plt.plot(grd1, zeros, color='k', linewidth=0.6, label='zero')
    elif (self.ig == 2):
        plt.plot(grd1[xlimitrange], rhs0[xlimitrange], color='#FF8C00', label='$-c_P \\rho dT/dt$')
        plt.plot(grd1[xlimitrange], rhs1[xlimitrange], color='y', label='$+\\delta dP/dt$')
        plt.plot(grd1, res, color='k', linestyle='--', label='res $\\sim N$')
        zeros = np.zeros(self.nx)
        plt.plot(grd1, zeros, color='k', linewidth=0.6, label='zero')
    plt.axvline(self.bconv, linestyle='--', linewidth=0.7, color='k')
    plt.axvline(self.tconv, linestyle='--', linewidth=0.7, color='k')
    if (self.ig == 1):
        setxlabel = 'x (cm)'
        setylabel = 'erg g$^{-1}$ s$^{-1}$'
        plt.xlabel(setxlabel)
        plt.ylabel(setylabel)
    elif (self.ig == 2):
        setxlabel = 'r (cm)'
        setylabel = 'erg g$^{-1}$ s$^{-1}$'
        plt.xlabel(setxlabel)
        plt.ylabel(setylabel)
    plt.legend(loc=ilg, prop={'size': 10}, ncol=2)
    plt.show(block=False)
    if (self.fext == 'png'):
        plt.savefig((('RESULTS/' + self.data_prefix) + 'standard_luminosity_exact_eq.png'))
    elif (self.fext == 'eps'):
        plt.savefig((('RESULTS/' + self.data_prefix) + 'standard_luminosity_exact_eq.eps'))
Plot luminosity equation in the model
def plot_luminosity_equation_exact(self, LAXIS, xbl, xbr, ybu, ybd, ilg):
    
    grd1 = self.xzn0
    rhs0 = self.minus_cp_rho_dTdt
    rhs1 = self.plus_delta_dPdt
    rhs2 = self.minus_dd_div_eiui
    rhs3 = self.plus_tke_diss
    res = self.minus_resLumExactEquation
    plt.figure(figsize=(7, 6))
    plt.gca().yaxis.get_major_formatter().set_powerlimits((0, 0))
    to_plot = [rhs0, rhs1, res]
    self.set_plt_axis(LAXIS, xbl, xbr, ybu, ybd, to_plot)
    self.bconv = 400000000.0
    self.tconv = 1200000000.0
    xlimitrange = np.where(((grd1 > self.bconv) & (grd1 < self.tconv)))
    xlimitbottom = np.where((grd1 < self.bconv))
    xlimittop = np.where((grd1 > self.tconv))
    plt.title('standard luminosity equation exact')
    if (self.ig == 1):
        plt.plot(grd1[xlimitrange], rhs0[xlimitrange], color='#FF8C00', label='$-c_P \\overline{\\rho \\partial_t T}$')
        plt.plot(grd1[xlimitrange], rhs1[xlimitrange], color='y', label='$+\\delta \\overline{\\partial_t P}$')
        plt.plot(grd1[xlimitrange], rhs2[xlimitrange], color='r', label='$-\\overline{\\rho \\nabla \\cdot \\epsilon_I {\\bf u}}$')
        plt.plot(grd1[xlimitrange], rhs3[xlimitrange], color='g', label='$+\\varepsilon_K$')
        plt.plot(grd1, res, color='k', linestyle='--', label='res $\\sim N$')
        zeros = np.zeros(self.nx)
        plt.plot(grd1, zeros, color='k', linewidth=0.6, label='zero')
    elif (self.ig == 2):
        plt.plot(grd1[xlimitrange], rhs0[xlimitrange], color='#FF8C00', label='$-c_P \\rho dT/dt$')
        plt.plot(grd1[xlimitrange], rhs1[xlimitrange], color='y', label='$+\\delta dP/dt$')
        plt.plot(grd1, res, color='k', linestyle='--', label='res $\\sim N$')
        zeros = np.zeros(self.nx)
        plt.plot(grd1, zeros, color='k', linewidth=0.6, label='zero')
    plt.axvline(self.bconv, linestyle='--', linewidth=0.7, color='k')
    plt.axvline(self.tconv, linestyle='--', linewidth=0.7, color='k')
    if (self.ig == 1):
        setxlabel = 'x (cm)'
        setylabel = 'erg g$^{-1}$ s$^{-1}$'
        plt.xlabel(setxlabel)
        plt.ylabel(setylabel)
    elif (self.ig == 2):
        setxlabel = 'r (cm)'
        setylabel = 'erg g$^{-1}$ s$^{-1}$'
        plt.xlabel(setxlabel)
        plt.ylabel(setylabel)
    plt.legend(loc=ilg, prop={'size': 10}, ncol=2)
    plt.show(block=False)
    if (self.fext == 'png'):
        plt.savefig((('RESULTS/' + self.data_prefix) + 'standard_luminosity_exact_eq.png'))
    elif (self.fext == 'eps'):
        plt.savefig((('RESULTS/' + self.data_prefix) + 'standard_luminosity_exact_eq.eps'))