<div dir="ltr">Dear flash users,<div>I would really appreciate it if somebody helps me with the following issue.</div><div>I am experimenting with the equation of state (EOS). Don't ask me how it is possible but my pressure becomes negative sometimes. To make sure that sound speed is real, after calculating Gamma=Cs^2*rho/P (Gamma is GAMC, P is pressure, rho is density, Cs is sound speed). I am checking if  Gamma*P is positive (density is always positive). And 

Gamma*P   is always positive in subroutine EOS. However  

I receive an error  message from subroutine hy_uhd_shockDetec because Gamma*P/rho becomes negative so sound speed becomes imaginary. What I know is that it happens in guard cells only. However I do not understand how it can happen. Somewhere on the way from subroutine EOS to subroutine hy_uhd_shockDetect 

Gamma*P in guard cells becomes negative yet as far as I understand both Gamma and P are calculated only in EOS.</div><div>Many thanks!</div><div>MAksim</div></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br> <table style="border-top:1px solid #d3d4de">
        <tr>
      <td style="width:55px;padding-top:18px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
                <td style="width:470px;padding-top:17px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Без вирусов. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color:#4453ea">www.avast.com</a>                </td>
        </tr>
</table>
<a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div>