<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"Times New Roman \(Body CS\)";
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:Menlo;
panose-1:2 11 6 9 3 8 4 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Arial",sans-serif;
color:windowtext;
font-weight:normal;
font-style:normal;}
p.p1, li.p1, div.p1
{mso-style-name:p1;
margin:0in;
background:black;
font-size:8.5pt;
font-family:Menlo;
color:#E6E6E6;}
p.p2, li.p2, div.p2
{mso-style-name:p2;
margin:0in;
background:black;
font-size:8.5pt;
font-family:Menlo;
color:#E6E6E6;}
span.apple-tab-span
{mso-style-name:apple-tab-span;}
span.s1
{mso-style-name:s1;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:"Arial",sans-serif">Hi everyone,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:"Arial",sans-serif"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal">I’m setting up a new simulation that incorporates the MGD, multispecies, and diffusion (for electron temperature conduction) units – see attached. This test problem is a simple 2/3D test that has two materials – one that is optically thick
and another that is optically thin. The material is initially cold with a hot electron and radiation temperature at x \approx 0.0, resulting in radiation flowing through the optically thin material. It is a simple variation of the “crooked pipe” radiation
transport test problem.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">This setup compiles successfully, but I am getting a seg fault at runtime related to the multispecies opacity. Here is the backtrace.<o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-converted-space"> </span><span class="s1">MaterialProperties initialized</span><o:p></o:p></p>
<p class="p2" style="margin-left:.5in"><o:p> </o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">Program received signal SIGSEGV: Segmentation fault - invalid memory reference.</span><o:p></o:p></p>
<p class="p2" style="margin-left:.5in"><o:p> </o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">Backtrace for this error:</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#0</span><span class="apple-converted-space">
</span><span class="s1">0x7f3b095353ff in ???</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#1</span><span class="apple-converted-space">
</span><span class="s1">0x44b551 in radtrans_mgdgetbound_</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/RadTrans_mgdGetBound.F90:29</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#2</span><span class="apple-converted-space">
</span><span class="s1">0x4ddd04 in op_setenergygroupboundaries_</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/op_setEnergyGroupBoundaries.F90:46</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#3</span><span class="apple-converted-space">
</span><span class="s1">0x445df7 in opacity_init_</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/Opacity_init.F90:141</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#4</span><span class="apple-converted-space">
</span><span class="s1">0x4126f6 in driver_initflash_</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/Driver_initFlash.F90:143</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#5</span><span class="apple-converted-space">
</span><span class="s1">0x407521 in flash</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/Flash.F90:49</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">#6</span><span class="apple-converted-space">
</span><span class="s1">0x407521 in main</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="apple-tab-span"> </span>
<span class="s1">at /nfs/depot/nerc_projects/dtra_rom/FLASH4.6.2/MGD_Tunnel_2D/Flash.F90:43</span><o:p></o:p></p>
<p class="p1" style="margin-left:.5in"><span class="s1">Segmentation fault (core dumped)</span><o:p></o:p></p>
<p class="MsoNormal" style="margin-left:.5in"><o:p></o:p></p>
<p class="MsoNormal">After marching through this backtrace, it seems as though the opacity initialization occurs
<i>before </i>RadTrans_init.F90 and results in the array of reals <span style="font-size:10.0pt;font-family:Menlo">
rt_mgdBounds</span> being empty when it is accessed on line 29 of RadTrans_mgdGetBound.F90…
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Does anyone have any suggestions on how to adjust the setup to fix this?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks in advance,<o:p></o:p></p>
<p class="MsoNormal">Tony <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:"Arial",sans-serif">-----------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black;background:white">Tony Alberti</span><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black"><br>
<span style="background:white">Postdoctoral Scholar</span><br>
<span style="background:white">School of Nuclear Science and Engineering</span><br>
</span><span style="font-size:10.5pt;font-family:"Arial",sans-serif"><a href="https://ne.oregonstate.edu/tony-alberti" target="_blank"><span style="color:#1155CC;background:white">https://ne.oregonstate.edu/tony-alberti</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Arial",sans-serif"><o:p> </o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>