<div dir="ltr">Hello Andrea,<div><br></div><div><br></div><div>  Are you calling EOS with the proper mode? Here is how I call my EOS with mode_dens_ei :</div><div><br></div><div><div>!!!!!</div><div><div><div><br></div><div>
#include "constants.h"</div><div><br></div><div><div>  real :: rhoZone, presZone, eintZone</div><div><br></div><div>! Variables for Eos<br></div><div>  real, dimension(EOS_NUM) :: eosData</div></div><div><br></div>
<div> ...</div><div><br></div><div><br></div><div>    rhoZone = solnData(DENS_VAR,i,j,k)</div></div><div><div>    eintZone = solnData(EINT_VAR,i,j,k)</div><div><br></div></div><div>    eosData(EOS_DENS) = rhoZone</div><div>
    eosData(EOS_EINT) = eintZone</div><div><br></div><div>    call Eos(MODE_DENS_EI, 1, eosData)</div></div><div><br></div><div>    presZone = eosData(EOS_PRES)</div><div>    solnData(PRES_VAR,i,j,k) = presZone</div><div>
<br></div><div>!!!!!</div></div><div><br></div><div><br></div><div><div><div><br></div></div></div><div>Where EOS_NUM and MODE_DENS_EI are defined in constants.h . eosData is then filled according to MODE_DENS_EI. </div><div>
<br></div><div><br></div><div>Looking at the code, what you're looking for should be:</div><div><br></div><div><br></div><div>!!!!!</div><div><br></div><div><div><div><div>#include "constants.h"</div><div><br>
</div><div><div>  real :: rhoZone, presZone</div><div><br></div><div>! Variables for Eos<br></div><div>  real, dimension(EOS_NUM) :: eosData</div></div><div><br></div><div>...</div><div><br></div><div>    ! assuming density and pressure are known</div>
<div>    rhoZone = solnData(DENS_VAR,i,j,k)</div></div><div><div>    presZone = solnData(PRES_VAR,i,j,k)</div><div><br></div></div><div>    eosData(EOS_DENS) = rhoZone</div><div>    eosData(EOS_PRES) = presZone</div><div>
<br></div><div><div>    call Eos(MODE_DENS_PRES, 1, eosData)</div></div><div><br></div><div>    eintZone = eosData(EOS_EINT)</div><div>    solnData(EINT_VAR,i,j,k) = eintZone</div><div><br></div><div><br></div></div><div>
!!!!!</div><div><br></div></div><div><br></div><div>Hope that helps,</div><div>   -Mark</div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Jan 17, 2014 at 10:42 AM, andrea costa <span dir="ltr"><<a href="mailto:acosta@oac.uncor.edu" target="_blank">acosta@oac.uncor.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all, I would like to know how to set up the mode_dens_pres.<div>I define in Simulation_initBlock and Simulation_adjustEvolution the density and the pressure.</div>
<div>However, printing the mode value it turns that the mode used is mode_dens_ei.</div>
<div><br></div><div>thanks </div><span class="HOEnZb"><font color="#888888"><div>andrea</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Mark Richardson, <br><a href="mailto:Mark.L.Richardson@asu.edu" target="_blank">Mark.L.Richardson@asu.edu</a><div>Ph.D. Candidate: Astrophysics<br>ISTB4 623<br>
School of Earth and Space Exploration<br>Arizona State University<br>480 318-4449<br><a href="http://www.public.asu.edu/~mlricha4" target="_blank">www.public.asu.edu/~mlricha4</a><br></div>
</div>