<div dir="ltr">Hi Jian Fuh,<div><br></div><div>Newer SESAME tables, at least the ones I've had access to, are all double precision (E22.15) and FLASH expects this. </div><div><br></div><div>The simplest solution would be to copy the file "source/physics/Eos/EosMain/Tabulated/eos_tabReadSesameTables.F90" into your simulation directory. Then you can find and replace all occurrences of "E22.15" with "E15.8". Then when you run setup you will use the single-precision version of the code you made in your simulation directory.</div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">*************************************************************************</div><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">Adam Reyes<div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH.jpg"><br></div><div>Code Group Leader, Flash Center for Computational Science <span> </span><br></div><div>Research Scientist, Dept. of Physics and Astronomy</div><div>University of Rochester</div><div>River Campus: Bausch and Lomb Hall, 369 <span> </span><br></div><div>500 Wilson Blvd. PO Box 270171, Rochester, NY 14627</div><div dir="ltr">Email<span> </span><a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a></div><div>Web<span> </span><a href="https://flash.rochester.edu/" target="_blank">https://flash.rochester.edu</a><br></div><div> (he / him / his)</div><div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH-pride-sml.png"><br></div>*************************************************************************</div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 9, 2023 at 4:20 PM Jian Fuh Ong <<a href="mailto:jianfuh.ong@eli-np.ro">jianfuh.ong@eli-np.ro</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><span style="font-style:normal;font-size:14px">Dear Flash Developers,</span><div><span style="font-style:normal;font-size:14px"><br></span></div><div><span style="font-style:normal;font-size:14px">Thank you for the release of FLASH4.7. I have tried out with the new feature for SESAME reader. I face a problem that the sesame data we obtain by using sesame-extract with opacplot2 is not compatible with the fortran format in the code.<br></span><div>
<div dir="auto" style="text-align:start;text-indent:0px"><div dir="auto" style="text-align:start;text-indent:0px"><div dir="auto" style="text-align:start;text-indent:0px"><div dir="auto" style="text-align:start;text-indent:0px"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;text-align:start;text-indent:0px"><span style="font-style:normal;font-size:14px"><br></span></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;text-align:start;text-indent:0px"><span style="font-style:normal;font-size:14px">We have for material 3717 with the following data structure:</span></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;text-align:start;text-indent:0px"><br></div><div dir="auto" style="text-align:start;text-indent:0px"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 0 3717 101 160 r 11889 11889 1 0</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">material. aluminum (z=13.0, a=26.9815) /source. j. barnes, s. lyon /date. aug 88</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> /refs. none /comp. al /codes. grizzly /Classification. Unclassified / </span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 1 3717 102 552 r 11889 11889 1 1</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">tfd as calculated by candide in grizzly was used for the electronic contribution</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">. chartjd nuclear model with jd gamma (igrun=7). gamref = 2.136, dgamma = -1.6, </span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">-2.2, gameps = 0.5, 1.0, debkel = 390.0, tmelt = 933.2. cold curve from shock da</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">ta and sound velocity. used us-up curve composed of a single straight line with </span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">c0 = 5.35 and s1 = 1.34. also used ecohkc = 69.0, faclj = 0.235, cmat = 1.588 an</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">d xalpha =0.2. both temperature and density grids enhanced with additional point</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">s. 301 table has van der waals loops. fits p(rho) hugoniot data well. </span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 1 3717 201 5 r 11889 11889 1 1</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 1.30000000E+01 2.69815000E+01 2.70000000E+00 0.00000000E+00 0.00000000E+0011111</span></font></div><div dir="auto"><font size="1"><font color="#000000">…….</font></font></div><div dir="auto"><font size="1"><font color="#000000"><br></font></font></div><div dir="auto"><br></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><span style="font-style:normal;font-size:14px">But the data come with FLASH have a longer line:</span></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><br></div><div dir="auto"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 0 9999 101 160 r 82803 22704 1 0</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">This Sesame table was generated using OpacPlot2 to convert an IONMIX EoS table t</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal">o the Sesame format </span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 1 9999 201 5 r 82803 22704 1 1</span></font></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font size="1"><span style="font-style:normal"> 1.300000000000000E+01 2.698153860000000E+01 1.000000000000000E+00 1.000000000000000E+00 1.000000000000000E+0011111</span></font></div><div dir="auto"><font size="1"><font color="#000000"><span>…..</span></font></font></div><div dir="auto"><font size="1"><font color="#000000"><span><br></span></font></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px">What would you suggest to do so that the shorter data converted from sesame can be read correctly. Thanks.</span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px"><br></span></font></div><div dir="auto"><font color="#000000"><span style="font-style:normal;font-size:14px">Regards,</span></font></div></div><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><br></div><font color="#000000"><span>Jian Fuh ONG</span></font><br><br><font color="#000000"><span>Research Scientist</span></font><br><font color="#000000"><span>Extreme Light Infrastructure - Nuclear Physics</span></font><br><font color="#000000"><span>Horia Hulubei National Institute for Physics and Nuclear Engineering</span></font><br><font color="#000000"><span>30 Reactorului Street, ELI-NP,</span></font><br><font color="#000000"><span>077125 Magurele, jud. Ilfov, Romania </span></font><br><font color="#000000"><span><a href="http://www.eli-np.ro" target="_blank">www.eli-np.ro</a></span></font><br><br><br><br><br></div></div></div></div></div>
</div>
<br></div></div>_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>