def __init__(self, filename, ig, ieos, fext, intc, tke_diss, bconv, tconv, data_prefix):
    super(LuminosityEquation, self).__init__(ig)
    eht = self.customLoad(filename)
    xzn0 = self.getRAdata(eht, 'xzn0')
    yzn0 = self.getRAdata(eht, 'yzn0')
    zzn0 = self.getRAdata(eht, 'zzn0')
    nx = self.getRAdata(eht, 'nx')
    dd = self.getRAdata(eht, 'dd')[intc]
    ux = self.getRAdata(eht, 'ux')[intc]
    pp = self.getRAdata(eht, 'pp')[intc]
    tt = self.getRAdata(eht, 'tt')[intc]
    cp = self.getRAdata(eht, 'cp')[intc]
    gg = self.getRAdata(eht, 'gg')[intc]
    abar = self.getRAdata(eht, 'abar')[intc]
    ddux = self.getRAdata(eht, 'ddux')[intc]
    dduy = self.getRAdata(eht, 'dduy')[intc]
    dduz = self.getRAdata(eht, 'dduz')[intc]
    ddttux = self.getRAdata(eht, 'ddttux')[intc]
    dduxttx = self.getRAdata(eht, 'dduxttx')[intc]
    dduytty = self.getRAdata(eht, 'dduytty')[intc]
    dduzttz = self.getRAdata(eht, 'dduzttz')[intc]
    eiuxddx = self.getRAdata(eht, 'eiuxddx')[intc]
    eiuyddy = self.getRAdata(eht, 'eiuyddy')[intc]
    eiuzddz = self.getRAdata(eht, 'eiuzddz')[intc]
    dduxux = self.getRAdata(eht, 'dduxux')[intc]
    dduyuy = self.getRAdata(eht, 'dduyuy')[intc]
    dduzuz = self.getRAdata(eht, 'dduzuz')[intc]
    dduxuy = self.getRAdata(eht, 'dduxuy')[intc]
    dduxuz = self.getRAdata(eht, 'dduxuz')[intc]
    ddekux = self.getRAdata(eht, 'ddekux')[intc]
    ddek = self.getRAdata(eht, 'ddek')[intc]
    ddei = self.getRAdata(eht, 'ddei')[intc]
    ddeiux = self.getRAdata(eht, 'ddeiux')[intc]
    eiux = self.getRAdata(eht, 'eiux')[intc]
    ddetux = self.getRAdata(eht, 'ddetux')[intc]
    divu = self.getRAdata(eht, 'divu')[intc]
    ppdivu = self.getRAdata(eht, 'ppdivu')[intc]
    dddivu = self.getRAdata(eht, 'dddivu')[intc]
    uxdivu = self.getRAdata(eht, 'uxdivu')[intc]
    ppux = self.getRAdata(eht, 'ppux')[intc]
    ddenuc1 = self.getRAdata(eht, 'ddenuc1')[intc]
    ddenuc2 = self.getRAdata(eht, 'ddenuc2')[intc]
    chim = self.getRAdata(eht, 'chim')[intc]
    chit = self.getRAdata(eht, 'chit')[intc]
    chid = self.getRAdata(eht, 'chid')[intc]
    gamma1 = self.getRAdata(eht, 'gamma1')[intc]
    gascon = 83144629.0
    if (ieos == 1):
        cp = self.getRAdata(eht, 'cp')[intc]
        cv = self.getRAdata(eht, 'cv')[intc]
        gamma1 = (cp / cv)
    t_timec = self.getRAdata(eht, 'timec')
    t_dd = self.getRAdata(eht, 'dd')
    t_tt = self.getRAdata(eht, 'tt')
    t_pp = self.getRAdata(eht, 'pp')
    t_ddei = self.getRAdata(eht, 'ddei')
    t_ddss = self.getRAdata(eht, 'ddss')
    t_ddtt = self.getRAdata(eht, 'ddtt')
    t_ddux = self.getRAdata(eht, 'ddux')
    t_dduy = self.getRAdata(eht, 'dduy')
    t_dduz = self.getRAdata(eht, 'dduz')
    t_dduxux = self.getRAdata(eht, 'dduxux')
    t_dduyuy = self.getRAdata(eht, 'dduyuy')
    t_dduzuz = self.getRAdata(eht, 'dduzuz')
    t_uxux = self.getRAdata(eht, 'uxux')
    t_uyuy = self.getRAdata(eht, 'uyuy')
    t_uzuz = self.getRAdata(eht, 'uzuz')
    t_fht_ek = ((0.5 * ((t_dduxux + t_dduyuy) + t_dduzuz)) / t_dd)
    t_fht_ei = (t_ddei / t_dd)
    t_fht_et = (t_fht_ek + t_fht_ei)
    t_fht_ss = (t_ddss / t_dd)
    t_fht_ux = (t_ddux / t_dd)
    t_fht_uy = (t_dduy / t_dd)
    t_fht_uz = (t_dduz / t_dd)
    t_fht_ui_fht_ui = (((t_fht_ux * t_fht_ux) + (t_fht_uy * t_fht_uy)) + (t_fht_uz * t_fht_uz))
    t_fht_tt = (t_ddtt / t_dd)
    fht_ek = (ddek / dd)
    fht_ux = (ddux / dd)
    fht_uy = (dduy / dd)
    fht_uz = (dduz / dd)
    fht_ei = (ddei / dd)
    fht_et = (fht_ek + fht_ei)
    fht_enuc = ((ddenuc1 + ddenuc2) / dd)
    fht_eiux = (ddeiux / dd)
    fei = (ddeiux - ((ddux * ddei) / dd))
    fekx = (ddekux - (fht_ux * fht_ek))
    fpx = (ppux - (pp * ux))
    fekx = (ddekux - (fht_ux * fht_ek))
    fht_ui_fht_ui = (((fht_ux * fht_ux) + (fht_uy * fht_uy)) + (fht_uz * fht_uz))
    if (self.ig == 1):
        alpha = 1.0
        delta = 1.0
        phi = 1.0
    elif (self.ig == 2):
        alpha = (1.0 / chid)
        delta = ((- chit) / chid)
        phi = (chid / chim)
    fht_rxx = (dduxux - ((ddux * ddux) / dd))
    fdil = (uxdivu - (ux * divu))
    gg = (- gg)
    if (self.ig == 1):
        surface = ((yzn0[(- 1)] - yzn0[0]) * (zzn0[(- 1)] - zzn0[0]))
    elif (self.ig == 2):
        surface = (((+ 4.0) * np.pi) * (xzn0 ** 2.0))
    else:
        print(('ERROR(Properties.py): ' + self.errorGeometry(self.ig)))
        sys.exit()
    self.minus_cp_rho_dTdt = ((- cp) * ((self.dt(t_ddtt, xzn0, t_timec, intc) + self.Div(ddttux, xzn0)) - ((dduxttx + dduytty) + dduzttz)))
    self.plus_delta_dPdt = ((+ delta) * self.dt(t_pp, xzn0, t_timec, intc))
    self.minus_dd_div_eiui = (- (self.Div(ddeiux, xzn0) - ((eiuxddx + eiuyddy) + eiuzddz)))
    self.plus_tke_diss = (+ tke_diss)
    self.minus_resLumExactEquation = (- (((self.minus_cp_rho_dTdt + self.plus_delta_dPdt) + self.minus_dd_div_eiui) + self.plus_tke_diss))
    self.data_prefix = data_prefix
    self.xzn0 = xzn0
    self.fht_et = (fht_ei + fht_ek)
    self.nx = nx
    self.bconv = bconv
    self.tconv = tconv
    self.fext = fext
def __init__(self, filename, ig, ieos, fext, intc, tke_diss, bconv, tconv, data_prefix):
    super(LuminosityEquation, self).__init__(ig)
    eht = self.customLoad(filename)
    xzn0 = self.getRAdata(eht, 'xzn0')
    yzn0 = self.getRAdata(eht, 'yzn0')
    zzn0 = self.getRAdata(eht, 'zzn0')
    nx = self.getRAdata(eht, 'nx')
    dd = self.getRAdata(eht, 'dd')[intc]
    ux = self.getRAdata(eht, 'ux')[intc]
    pp = self.getRAdata(eht, 'pp')[intc]
    tt = self.getRAdata(eht, 'tt')[intc]
    cp = self.getRAdata(eht, 'cp')[intc]
    gg = self.getRAdata(eht, 'gg')[intc]
    abar = self.getRAdata(eht, 'abar')[intc]
    ddux = self.getRAdata(eht, 'ddux')[intc]
    dduy = self.getRAdata(eht, 'dduy')[intc]
    dduz = self.getRAdata(eht, 'dduz')[intc]
    ddttux = self.getRAdata(eht, 'ddttux')[intc]
    dduxttx = self.getRAdata(eht, 'dduxttx')[intc]
    dduytty = self.getRAdata(eht, 'dduytty')[intc]
    dduzttz = self.getRAdata(eht, 'dduzttz')[intc]
    eiuxddx = self.getRAdata(eht, 'eiuxddx')[intc]
    eiuyddy = self.getRAdata(eht, 'eiuyddy')[intc]
    eiuzddz = self.getRAdata(eht, 'eiuzddz')[intc]
    dduxux = self.getRAdata(eht, 'dduxux')[intc]
    dduyuy = self.getRAdata(eht, 'dduyuy')[intc]
    dduzuz = self.getRAdata(eht, 'dduzuz')[intc]
    dduxuy = self.getRAdata(eht, 'dduxuy')[intc]
    dduxuz = self.getRAdata(eht, 'dduxuz')[intc]
    ddekux = self.getRAdata(eht, 'ddekux')[intc]
    ddek = self.getRAdata(eht, 'ddek')[intc]
    ddei = self.getRAdata(eht, 'ddei')[intc]
    ddeiux = self.getRAdata(eht, 'ddeiux')[intc]
    eiux = self.getRAdata(eht, 'eiux')[intc]
    ddetux = self.getRAdata(eht, 'ddetux')[intc]
    divu = self.getRAdata(eht, 'divu')[intc]
    ppdivu = self.getRAdata(eht, 'ppdivu')[intc]
    dddivu = self.getRAdata(eht, 'dddivu')[intc]
    uxdivu = self.getRAdata(eht, 'uxdivu')[intc]
    ppux = self.getRAdata(eht, 'ppux')[intc]
    ddenuc1 = self.getRAdata(eht, 'ddenuc1')[intc]
    ddenuc2 = self.getRAdata(eht, 'ddenuc2')[intc]
    chim = self.getRAdata(eht, 'chim')[intc]
    chit = self.getRAdata(eht, 'chit')[intc]
    chid = self.getRAdata(eht, 'chid')[intc]
    gamma1 = self.getRAdata(eht, 'gamma1')[intc]
    gascon = 83144629.0
    if (ieos == 1):
        cp = self.getRAdata(eht, 'cp')[intc]
        cv = self.getRAdata(eht, 'cv')[intc]
        gamma1 = (cp / cv)
    t_timec = self.getRAdata(eht, 'timec')
    t_dd = self.getRAdata(eht, 'dd')
    t_tt = self.getRAdata(eht, 'tt')
    t_pp = self.getRAdata(eht, 'pp')
    t_ddei = self.getRAdata(eht, 'ddei')
    t_ddss = self.getRAdata(eht, 'ddss')
    t_ddtt = self.getRAdata(eht, 'ddtt')
    t_ddux = self.getRAdata(eht, 'ddux')
    t_dduy = self.getRAdata(eht, 'dduy')
    t_dduz = self.getRAdata(eht, 'dduz')
    t_dduxux = self.getRAdata(eht, 'dduxux')
    t_dduyuy = self.getRAdata(eht, 'dduyuy')
    t_dduzuz = self.getRAdata(eht, 'dduzuz')
    t_uxux = self.getRAdata(eht, 'uxux')
    t_uyuy = self.getRAdata(eht, 'uyuy')
    t_uzuz = self.getRAdata(eht, 'uzuz')
    t_fht_ek = ((0.5 * ((t_dduxux + t_dduyuy) + t_dduzuz)) / t_dd)
    t_fht_ei = (t_ddei / t_dd)
    t_fht_et = (t_fht_ek + t_fht_ei)
    t_fht_ss = (t_ddss / t_dd)
    t_fht_ux = (t_ddux / t_dd)
    t_fht_uy = (t_dduy / t_dd)
    t_fht_uz = (t_dduz / t_dd)
    t_fht_ui_fht_ui = (((t_fht_ux * t_fht_ux) + (t_fht_uy * t_fht_uy)) + (t_fht_uz * t_fht_uz))
    t_fht_tt = (t_ddtt / t_dd)
    fht_ek = (ddek / dd)
    fht_ux = (ddux / dd)
    fht_uy = (dduy / dd)
    fht_uz = (dduz / dd)
    fht_ei = (ddei / dd)
    fht_et = (fht_ek + fht_ei)
    fht_enuc = ((ddenuc1 + ddenuc2) / dd)
    fht_eiux = (ddeiux / dd)
    fei = (ddeiux - ((ddux * ddei) / dd))
    fekx = (ddekux - (fht_ux * fht_ek))
    fpx = (ppux - (pp * ux))
    fekx = (ddekux - (fht_ux * fht_ek))
    fht_ui_fht_ui = (((fht_ux * fht_ux) + (fht_uy * fht_uy)) + (fht_uz * fht_uz))
    if (self.ig == 1):
        alpha = 1.0
        delta = 1.0
        phi = 1.0
    elif (self.ig == 2):
        alpha = (1.0 / chid)
        delta = ((- chit) / chid)
        phi = (chid / chim)
    fht_rxx = (dduxux - ((ddux * ddux) / dd))
    fdil = (uxdivu - (ux * divu))
    gg = (- gg)
    if (self.ig == 1):
        surface = ((yzn0[(- 1)] - yzn0[0]) * (zzn0[(- 1)] - zzn0[0]))
    elif (self.ig == 2):
        surface = (((+ 4.0) * np.pi) * (xzn0 ** 2.0))
    else:
        print(('ERROR(Properties.py): ' + self.errorGeometry(self.ig)))
        sys.exit()
    self.minus_cp_rho_dTdt = ((- cp) * ((self.dt(t_ddtt, xzn0, t_timec, intc) + self.Div(ddttux, xzn0)) - ((dduxttx + dduytty) + dduzttz)))
    self.plus_delta_dPdt = ((+ delta) * self.dt(t_pp, xzn0, t_timec, intc))
    self.minus_dd_div_eiui = (- (self.Div(ddeiux, xzn0) - ((eiuxddx + eiuyddy) + eiuzddz)))
    self.plus_tke_diss = (+ tke_diss)
    self.minus_resLumExactEquation = (- (((self.minus_cp_rho_dTdt + self.plus_delta_dPdt) + self.minus_dd_div_eiui) + self.plus_tke_diss))
    self.data_prefix = data_prefix
    self.xzn0 = xzn0
    self.fht_et = (fht_ei + fht_ek)
    self.nx = nx
    self.bconv = bconv
    self.tconv = tconv
    self.fext = fext