<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- p { margin-top: 0px; margin-bottom: 0px; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dear all,<br>
</p>
<p><br>
</p>
<p>I've been trying to integrate ionization into my simulation, but am running into an error early in the run.<br>
</p>
<p><br>
</p>
<p>I looked at the NeiTest problem, which works fine, and added what I thought was necessary to the various setup files in my simulation.<br>
</p>
<p>I added the following to the Config:<br>
</p>
<p><br>
</p>
<div>REQUIRES physics/Eos/EosMain/Multigamma</div>
<div>REQUIRES Simulation/SimulationComposition/Ionize<br>
</div>
<div>REQUIRES physics/sourceTerms/Ionize<br>
</div>
<div><br>
</div>
<div>
<div>SPECIES he TO 2</div>
<div>SPECIES c TO 6</div>
<div>SPECIES n TO 7<br>
</div>
<div>#SPECIES o TO 8</div>
<div>#SPECIES ne TO 10</div>
<div>#SPECIES mg TO 12</div>
<div>#SPECIES si TO 14</div>
<div>#SPECIES s TO 16</div>
<div>#SPECIES ar TO 18</div>
<div>#SPECIES ca TO 20</div>
<div>#SPECIES fe TO 26</div>
<div>#SPECIES ni TO 28<br>
</div>
<br>
</div>
<p>Some species were commented out because I received an error regarding file size and memory, and that seemed to be fixed by removing some species and reducing the checkpoint file size.<br>
</p>
<p><br>
</p>
<p>The setup was as follows, and is fairly standard:<br>
</p>
<p>setup.py MySimulation -2d -auto +spherical -objdir balaramasim +parallelio -maxblocks=10000 +pm4dev<br>
</p>
<p><br>
</p>
<p>The logfile is attached, and the error there is:<br>
</p>
<div> [ 04-26-2017 17:16:36.447 ] step: n=198 t=6.445519E+03 dt=2.535561E+01</div>
<div> [ 04-26-2017 17:16:41.217 ] [DRIVER_ABORT]: Driver_abort() called by PE 0</div>
<div> [ 04-26-2017 17:16:41.230 ] abort_message: error in routine neimn: wrong sum of pop. fraction<br>
</div>
<div><br>
</div>
<div>In the actual simulation, there was mention of elsum:<br>
</div>
<div>
<div>Element I.D.: 1 elsum = 0.60024578323989231 </div>
<div> error in routine neimn: wrong sum of population fraction<br>
</div>
<div><br>
</div>
<div>In the initBlock, there was an existing section that deals with fractions, and it divided fractions evenly.<br>
</div>
<div>Meaning the various ionization variables (elec,h,he,he2, and so on) were all equal to 0.05 at the start of the simulation.<br>
</div>
<div><br>
</div>
<div>Does anyone know what is causing the problem, and how to solve it?<br>
</div>
<div><br>
</div>
<div>Kind Regards,<br>
</div>
<div>Michiel Bustraan<br>
</div>
<div><br>
<br>
</div>
<br>
</div>
<p><br>
</p>
</body>
</html>