def __init__(self, filename, ig, fext, intc, tke_diss, data_prefix):
    super(InternalEnergyEquation, self).__init__(ig)
    eht = self.customLoad(filename)
    xzn0 = self.getRAdata(eht, 'xzn0')
    nx = self.getRAdata(eht, 'nx')
    dd = self.getRAdata(eht, 'dd')[intc]
    ux = self.getRAdata(eht, 'ux')[intc]
    pp = self.getRAdata(eht, 'pp')[intc]
    ddux = self.getRAdata(eht, 'ddux')[intc]
    ddei = self.getRAdata(eht, 'ddei')[intc]
    ddeiux = self.getRAdata(eht, 'ddeiux')[intc]
    divu = self.getRAdata(eht, 'divu')[intc]
    ppdivu = self.getRAdata(eht, 'ppdivu')[intc]
    ddenuc1 = self.getRAdata(eht, 'ddenuc1')[intc]
    ddenuc2 = self.getRAdata(eht, 'ddenuc2')[intc]
    t_timec = self.getRAdata(eht, 'timec')
    t_dd = self.getRAdata(eht, 'dd')
    t_ddei = self.getRAdata(eht, 'ddei')
    t_fht_ei = (t_ddei / t_dd)
    fht_ux = (ddux / dd)
    fht_ei = (ddei / dd)
    fei = (ddeiux - ((ddux * ddei) / dd))
    self.minus_dt_dd_fht_ei = (- self.dt((t_dd * t_fht_ei), xzn0, t_timec, intc))
    self.minus_div_dd_fht_ux_fht_ei = (- self.Div(((dd * fht_ux) * fht_ei), xzn0))
    self.minus_div_fei = (- self.Div(fei, xzn0))
    self.minus_div_ftt = (- np.zeros(nx))
    self.minus_pp_div_ux = ((- pp) * self.Div(ux, xzn0))
    self.minus_eht_ppf_df = (- (ppdivu - (pp * divu)))
    self.plus_dd_fht_enuc = (ddenuc1 + ddenuc2)
    self.plus_disstke = (+ tke_diss)
    self.minus_resEiEquation = (- (((((((self.minus_dt_dd_fht_ei + self.minus_div_dd_fht_ux_fht_ei) + self.minus_div_fei) + self.minus_div_ftt) + self.minus_pp_div_ux) + self.minus_eht_ppf_df) + self.plus_dd_fht_enuc) + self.plus_disstke))
    self.data_prefix = data_prefix
    self.xzn0 = xzn0
    self.fht_ei = fht_ei
    self.fext = fext
def __init__(self, filename, ig, fext, intc, tke_diss, data_prefix):
    super(InternalEnergyEquation, self).__init__(ig)
    eht = self.customLoad(filename)
    xzn0 = self.getRAdata(eht, 'xzn0')
    nx = self.getRAdata(eht, 'nx')
    dd = self.getRAdata(eht, 'dd')[intc]
    ux = self.getRAdata(eht, 'ux')[intc]
    pp = self.getRAdata(eht, 'pp')[intc]
    ddux = self.getRAdata(eht, 'ddux')[intc]
    ddei = self.getRAdata(eht, 'ddei')[intc]
    ddeiux = self.getRAdata(eht, 'ddeiux')[intc]
    divu = self.getRAdata(eht, 'divu')[intc]
    ppdivu = self.getRAdata(eht, 'ppdivu')[intc]
    ddenuc1 = self.getRAdata(eht, 'ddenuc1')[intc]
    ddenuc2 = self.getRAdata(eht, 'ddenuc2')[intc]
    t_timec = self.getRAdata(eht, 'timec')
    t_dd = self.getRAdata(eht, 'dd')
    t_ddei = self.getRAdata(eht, 'ddei')
    t_fht_ei = (t_ddei / t_dd)
    fht_ux = (ddux / dd)
    fht_ei = (ddei / dd)
    fei = (ddeiux - ((ddux * ddei) / dd))
    self.minus_dt_dd_fht_ei = (- self.dt((t_dd * t_fht_ei), xzn0, t_timec, intc))
    self.minus_div_dd_fht_ux_fht_ei = (- self.Div(((dd * fht_ux) * fht_ei), xzn0))
    self.minus_div_fei = (- self.Div(fei, xzn0))
    self.minus_div_ftt = (- np.zeros(nx))
    self.minus_pp_div_ux = ((- pp) * self.Div(ux, xzn0))
    self.minus_eht_ppf_df = (- (ppdivu - (pp * divu)))
    self.plus_dd_fht_enuc = (ddenuc1 + ddenuc2)
    self.plus_disstke = (+ tke_diss)
    self.minus_resEiEquation = (- (((((((self.minus_dt_dd_fht_ei + self.minus_div_dd_fht_ux_fht_ei) + self.minus_div_fei) + self.minus_div_ftt) + self.minus_pp_div_ux) + self.minus_eht_ppf_df) + self.plus_dd_fht_enuc) + self.plus_disstke))
    self.data_prefix = data_prefix
    self.xzn0 = xzn0
    self.fht_ei = fht_ei
    self.fext = fext