<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
The sesame table you’re trying to convert is probably missing the 304 tables for the electron EOS. I’ve run into this problem before for another of the sesame tables that only included the 301 tables for the total EOS. Unfortunately in this case there is nothing
 you can do, since FLASH expects the separate ion & electron tables in the tabulated EOS implementation.
<div class=""><br class="">
</div>
<div class="">Adam<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Aug 26, 2021, at 3:07 AM, Ong Jian Fuh <<a href="mailto:jianfuh.ong@eli-np.ro" class="">jianfuh.ong@eli-np.ro</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">
<p class="">Dear all,</p>
<p class="">I am converting sesame data into the one readable by FLASH using opacplot2. It seems that no one is maintaining the GitHub repository. I first extract the sesame data by:</p>
<p class=""><font size="-2" class="">$ sesame-extract xses_asc 3715</font></p>
<p class="">Then</p>
<p class=""><font size="-2" class="">$ opac-convert --Znum 13 xses_3715.ses --outname al-imx-006</font></p>
<p class="">and it returns the following:<br class="">
</p>
<p class=""><font size="-2" class="">Traceback (most recent call last):<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/scripts/opac_convert.py", line 135, in __init__<br class="">
    self.eos_dict = self.handle_dict[args.input]()<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/scripts/opac_convert.py", line 181, in sesame_toEosDict<br class="">
    eos_dict = op.toEosDict(Znum=self.args.Znum,<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/opg_sesame.py", line 290, in toEosDict<br class="">
    opp_ses_data = opp.utils.EosMergeGrids(<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/utils.py", line 233, in __init__<br class="">
    i_grids[key] = eos_data['_'.join((intersect[0],key))]<br class="">
KeyError: 'ele_dens'<br class="">
<br class="">
During handling of the above exception, another exception occurred:<br class="">
<br class="">
Traceback (most recent call last):<br class="">
  File "/home/ong/.anaconda3/bin/opac-convert", line 8, in <module><br class="">
    sys.exit(convert_tables())<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/scripts/opac_convert.py", line 302, in convert_tables<br class="">
    eos_dict = Formats_toEosDict(input_data['args'],<br class="">
  File "/home/ong/.anaconda3/lib/python3.8/site-packages/opacplot2/scripts/opac_convert.py", line 137, in __init__<br class="">
    raise KeyError('Must use valid format name.')<br class="">
KeyError: 'Must use valid format name.'</font></p>
<p class="">Then I try the other element. There was no problem in converting for element 3713 but it does not work for FLASH with</p>
<p class=""> <font size="-2" class="">DRIVER_ABORT: [op_readIonmix4Tables] ERROR: no IONMIX4 energy group grid found</font></p>
<p class="">When I compare the generated IONMIX data (element 3715 & 3713) with the one come with FLASH, there is a lot of zeros</p>
<p class=""><font size="-2" class="">        28        17<br class="">
 atomic #s of gases:         13<br class="">
 relative fractions:   1.00E+00<br class="">
           0<br class="">
numbersnumbersnumbersnumbers<br class="">
numbersnumbersnumbersnumbers<br class="">
</font></p>
<p class=""><font size="-2" class="">.....<br class="">
</font></p>
<p class=""><font size="-2" class="">0.000000E+000.000000E+000.000000E+000.000000E+00<br class="">
0.000000E+000.000000E+000.000000E+000.000000E+00<br class="">
.....</font></p>
<p class="">Does anyone has a successful experience in converting a working EoS data from sesame? Thanks.<br class="">
</p>
<pre class="moz-signature" cols="72">-- 
Regards,

Dr. Ong Jian Fuh
Extreme Light Infrastructure - Nuclear Physics
"Horia Hulubei" National Institute for Physics and Nuclear Engineering
30 Reactorului Street, 077125 Magurele (Ilfov), Romania
<a class="moz-txt-link-abbreviated" href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.eli-2Dnp.ro&d=DwMDaQ&c=kbmfwr1Yojg42sGEpaQh5ofMHBeTl9EI2eaqQZhHbOU&r=1C4d97g-LRBm3GqeRtlZHfnzFQIEWizIeq2b_IZClJA&m=Vx2iTYvpwWtWhJUdr3JGALczzJvBMDf1ZevSJvN6gN4&s=k5XhW15y1fCZwmTzesYoGJqePkkh-hjk0UD0LgKSfns&e=">www.eli-np.ro</a></pre>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>