[FLASH-USERS] Reading in grid coordinate data in the EOS unit

[Max Katz]

A more direct way of doing this may be to simply disallow changes to the
internal energy in the hydro update. If for example you're using the split
PPM solver, you can comment out

              solnData(EINT_VAR,i,j,k) = einternal

in source/physics/Hydro/HydroMain/split/PPM/hy_ppm_updateSoln.F90.

Max Katz
Ph.D. Candidate, Department of Physics and Astronomy
Stony Brook University
http://astro.sunysb.edu/mkatz/

On Wed, Jun 10, 2015 at 1:20 PM, Andrea Derdzinski <
aderdzinski at astro.columbia.edu> wrote:

> Hi FLASH users,
>
>  I am running a simulation in which I want to hold the internal energy
> fixed. The internal energy has a radial dependence, so I need to read in
> the grid coordinate information at every time step and set the internal
> energy to its initial profile eint(r).   My idea was to modify the EOS unit
> (using the ideal gamma EOS, so changing eos_idealGamma.F90) such that it
> reads in the grid coordinate information and resets the internal energy to
> it's initial state. However, this module does not receive the block
> information as an input, so I am unsure of the best way to retrieve the
> coordinate information.
> Does anyone have any suggestions for how to do this?
>
>
> Thanks,
> -Andrea
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20150610/01e551b9/attachment.htm>