<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Dear Ernesto,</p>
<p>Check the file source/IO/IOMain/IO_writeIntegralQuantities.F90:</p>
<p>Here you see that ENER does only contain the magnetic field
energy if the variable MAGP_VAR is defined. <br>
</p>
<p>#ifdef ENER_VAR<br>
lsum(5) = lsum(5) + solnData(ENER_VAR,i,j,k) * &
<br>
&
solnData(DENS_VAR,i,j,k)*dvol<br>
#ifdef MAGP_VAR<br>
! total plasma energy<br>
!!$ lsum(5) = lsum(5) + (solnData(ENER_VAR,i,j,k) *
& <br>
!!$ & solnData(DENS_VAR,i,j,k) +
solnData(MAGP_VAR,i,j,k))*dvol<br>
<br>
lsum(5) = lsum(5) + solnData(MAGP_VAR,i,j,k)*dvol<br>
#endif<br>
#endif</p>
<p>should likely be replaced by something like</p>
<p>#ifdef ENER_VAR<br>
lsum(5) = lsum(5) + solnData(ENER_VAR,i,j,k) * &<br>
&
solnData(DENS_VAR,i,j,k)*dvol<br>
#ifdef MAGX_VAR<br>
lsum(5) = lsum(5) +
0.5*(solnData(MAGX_VAR,i,j,k)**2+solnData(MAGY_VAR,i,j,k)**2+solnData(MAGZ_VAR,i,j,k)**2)*dvol<br>
#endif<br>
#endif</p>
<p>as there is no need for MAGP_VAR here.<br>
</p>
<p>Best regards,<br>
Dominik<br>
</p>
<br>
<div class="moz-cite-prefix">On 18.10.2017 16:34, ERNESTO ZURBRIGGEN
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAE7HDOJ9LSjFXZaqZ5QHr_asNtiZZcqh4K+3cVEWNOdcOUL-4g@mail.gmail.com">
<div dir="ltr">
<div>Hi everybody!</div>
<div>I'm having a doubt about what data is stored in the UNK
specific total energy ENER_VAR in MHD simulations. </div>
<div><br>
</div>
<div>The flash4.2.2 user's guide equation 14.22 says for the
perfect-gas equation of state that</div>
<div>ENER_VAR= 0.5*v + PRES/(DENS*(gamma-1)) + B**2/(2*DENS)</div>
<div><br>
</div>
<div>But, the outputs seem to show that</div>
<div>ENER_VAR= 0.5*v + PRES/(DENS*(gamma-1))<br>
</div>
<div><br>
</div>
<div>So, I'm confused. Any idea?</div>
<div><br>
</div>
<div><br>
</div>
<div>I'm using flash4.2.2 release in an ideal MHD simulation
with 1T perfect-gas equation of state. </div>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr"><font face="monospace,
monospace"><i><font color="#000000">Ernesto
Zurbriggen</font></i><i> </i></font>
<div><font face="monospace, monospace"><i><br>
</i></font>
<div><font face="monospace, monospace"><i><font
color="#444444">Instituto de
Astronomía Teórica y
Experimental </font><font
color="#274e13">(IATE)</font><font
color="#444444">. </font></i></font></div>
<div><font face="monospace, monospace"><i><font
color="#444444">Observatorio
Astronómico de Córdoba </font><font
color="#274e13">(OAC)</font><font
color="#444444">, </font></i></font><i
style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">Universidad
Nacional de Córdoba </i><i
style="font-family:monospace,monospace;font-size:12.8px"><font
color="#274e13">(UNC)</font></i><i
style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">. </i></div>
</div>
<div><i
style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px">Teléfono:
+54 0351 4331064-5, interno 222. </i></div>
<div><span style="font-size:12.8px"><i><font
color="#444444">Laprida 854, Bº
Observatorio (CP 5000), Córdoba,
Argentina.</font></i></span><font
color="#444444"><br>
</font></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
</body>
</html>