[FLASH-USERS] SESAME reader in FLASH 4.7

[Adam Reyes]

Hi Jian Fuh,

Newer SESAME tables, at least the ones I've had access to, are all double
precision (E22.15) and FLASH expects this.

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.

*************************************************************************
Adam Reyes

Code Group Leader, Flash Center for Computational Science
Research Scientist, Dept. of Physics and Astronomy
University of Rochester
River Campus: Bausch and Lomb Hall, 369
500 Wilson Blvd. PO Box 270171, Rochester, NY 14627
Email adam.reyes at rochester.edu
Web https://flash.rochester.edu
 (he / him / his)

*************************************************************************

On Mon, Jan 9, 2023 at 4:20 PM Jian Fuh Ong <jianfuh.ong at eli-np.ro> wrote:

> Dear Flash Developers,
>
> 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.
>
> We have for material 3717 with the following data structure:
>
>  0  3717   101   160   r    11889    11889   1
>     0
> material. aluminum (z=13.0, a=26.9815) /source. j. barnes, s. lyon /date.
> aug 88
>  /refs. none /comp. al /codes. grizzly  /Classification. Unclassified /
>
>  1  3717   102   552   r    11889    11889   1
>     1
> tfd as calculated by candide in grizzly was used for the electronic
> contribution
> . chartjd nuclear model with jd gamma (igrun=7). gamref = 2.136, dgamma =
> -1.6,
> -2.2, gameps = 0.5, 1.0, debkel = 390.0, tmelt = 933.2. cold curve from
> shock da
> ta and sound velocity. used us-up curve composed of a single straight line
> with
> c0 = 5.35 and s1 = 1.34. also used ecohkc = 69.0, faclj = 0.235, cmat =
> 1.588 an
> d xalpha =0.2. both temperature and density grids enhanced with additional
> point
> s. 301 table has van der waals loops. fits p(rho) hugoniot data well.
>  1  3717   201     5   r    11889    11889   1
>     1
>  1.30000000E+01 2.69815000E+01 2.70000000E+00 0.00000000E+00
> 0.00000000E+0011111
> …….
>
>
> But the data come with FLASH have a longer line:
>
>  0  9999   101   160   r    82803    22704   1
>     0
> This Sesame table was generated using OpacPlot2 to convert an IONMIX EoS
> table t
> o the Sesame format
>
>  1  9999   201   5   r    82803    22704   1
>     1
>  1.300000000000000E+01 2.698153860000000E+01 1.000000000000000E+00
> 1.000000000000000E+00 1.000000000000000E+0011111
> …..
>
> What would you suggest to do so that the shorter data converted from
> sesame can be read correctly. Thanks.
>
>
>
>
>
>
>
> Regards,
>
> Jian Fuh ONG
>
> Research Scientist
> Extreme Light Infrastructure - Nuclear Physics
> Horia Hulubei National Institute for Physics and Nuclear Engineering
> 30 Reactorului Street, ELI-NP,
> 077125 Magurele, jud. Ilfov, Romania
> www.eli-np.ro
>
>
>
>
>
> _______________________________________________
> flash-users mailing list
> flash-users at flash.rochester.edu
>
> For list info, including unsubscribe:
> https://flash.rochester.edu/mailman/listinfo/flash-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230109/d5ac8bb5/attachment.htm>