opacity_conversion_exemple

import opacplot2 as opp

# create a dictionary op with opacity data from a TOPS file 
# check that units are consistent

# write to Ionmix format
opp.writeIonmixFile(fname,
                            op['Znum'], op['Xnum'],
                            numDens=op['idens'][:], # electron density
                            temps=op['temp'][:],
                            ngroups=Ng,
                            opac_bounds=groups_array,
                            planck_absorb=res['opp_mg'][:],
                            rosseland=res['opr_mg'][:],
                            planck_emiss=res['emp_mg'][:])