[FLASH-USERS] Mode Dens-Pres

[Mark L Richardson]

Hello Andrea,


  Are you calling EOS with the proper mode? Here is how I call my EOS with
mode_dens_ei :

!!!!!

#include "constants.h"

  real :: rhoZone, presZone, eintZone

! Variables for Eos
  real, dimension(EOS_NUM) :: eosData

 ...


    rhoZone = solnData(DENS_VAR,i,j,k)
    eintZone = solnData(EINT_VAR,i,j,k)

    eosData(EOS_DENS) = rhoZone
    eosData(EOS_EINT) = eintZone

    call Eos(MODE_DENS_EI, 1, eosData)

    presZone = eosData(EOS_PRES)
    solnData(PRES_VAR,i,j,k) = presZone

!!!!!



Where EOS_NUM and MODE_DENS_EI are defined in constants.h . eosData is then
filled according to MODE_DENS_EI.


Looking at the code, what you're looking for should be:


!!!!!

#include "constants.h"

  real :: rhoZone, presZone

! Variables for Eos
  real, dimension(EOS_NUM) :: eosData

...

    ! assuming density and pressure are known
    rhoZone = solnData(DENS_VAR,i,j,k)
    presZone = solnData(PRES_VAR,i,j,k)

    eosData(EOS_DENS) = rhoZone
    eosData(EOS_PRES) = presZone

    call Eos(MODE_DENS_PRES, 1, eosData)

    eintZone = eosData(EOS_EINT)
    solnData(EINT_VAR,i,j,k) = eintZone


!!!!!


Hope that helps,
   -Mark






On Fri, Jan 17, 2014 at 10:42 AM, andrea costa <acosta at oac.uncor.edu> wrote:

> Hi all, I would like to know how to set up the mode_dens_pres.
> I define in Simulation_initBlock and Simulation_adjustEvolution the
> density and the pressure.
> However, printing the mode value it turns that the mode used is
> mode_dens_ei.
>
> thanks
> andrea
>



-- 
Mark Richardson,
Mark.L.Richardson at asu.edu
Ph.D. Candidate: Astrophysics
ISTB4 623
School of Earth and Space Exploration
Arizona State University
480 318-4449
www.public.asu.edu/~mlricha4
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20140117/e6aa17ac/attachment-0001.htm>