PlasmaCalcs.mhd.mhd_er_tables.eos_file_tables

PlasmaCalcs.mhd.mhd_er_tables.eos_file_tables(eos_file, *, epermcoords, rcoords, u=None)

return erTableManager dict (of {var: Equation Of State lookup table}), for vars in file.

file should have quantities:

‘P’ (pressure),

‘T’ (temperature),

‘ne’ (electron number density),

‘kappaR’ (Rosseland opacity).

file values of P, ne, and kappaR should be in ln-space (e.g. ln(P) instead of P).

file should have values in cgs

(despite note in helita implying ne would be in si… I’m 99% sure it is in cgs here.

e.g. ne close to n (== r / (amu * self.M_AMU_FLUID)) in the corona,

where plasma should be mostly ionized.

As opposed to ne 1e-6 times smaller than n, if ne was truly in si in table.)

eos_file: str

path to eos table file.

epermcoords, rcoords: 1D array

coords for eperm (energy per unit mass) and r (mass density).

Units should be ‘raw’ units.

u: None or UnitsManager

used to calculate scale_quants to convert to ‘raw’ units.

None –> all scale quants will be 1.