<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><br>Hi,<br>

    <br>

    I am trying to setup a simulation using FLASH code. It runs well in
    a test problem, Sod Shock-Tube, using the original initial
    parameters. But when I run it using physical parameters in units of
    "cgs", it raises a warning:<br>

    "after gc filling: min.
    unk(DENS_VAR)=0.1019999999999999981868E-24     PE=0    
    block=21                               type=1      ",<br>

    without any output file except the "0000" one. If I change the
    "plotfileIntervalStep" to 1, it produces some 'plt' files, but the
    result is obviously wrong. Following is what the screen output looks
    like. The "flash.par" file and a result density image are in the
    attachments.<br>

    <br>

    
<blockquote><small>16  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23<br>
         15  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
         14  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
         13  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
         12  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
         11  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
         10  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          9  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          8  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          7  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          6  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          5  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          4  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          3  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          2  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
          1  0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23 0.10E-23
        0.10E-23 0.10E-23 0.10E-23<br>
        WARNING after gc filling: min.
        unk(DENS_VAR)=0.1019999999999999981868E-24     PE=0    
        block=21                               type=1              
                        </small><br>
      ...<br>
        <small>1  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
        WARNING after gc filling: min.
        unk(DENS_VAR)=0.9999999999999999071852E-10     PE=0    
        block=37                              
        type=1                               <br>
         16  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         15  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         14  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         13  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         12  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         11  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
         10  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          9  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          8  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          7  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          6  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          5  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          4  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 ^C0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          3  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          2  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
          1  0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09 0.10E-09
        0.10E-09 0.10E-09 0.10E-09<br>
        WARNING after gc filling: min.
        unk(DENS_VAR)=0.9999999999999999071852E-10     PE=0    
        block=41                              
        type=1                               <br>
      </small><br>
    </blockquote>

    From the warning message, I think it's because of the value of
    density. Is it too small? Or any thoughts on how to make it better?<br>

    <br>

    Thanks in advance.<br>

    <br>

    Cheers,<br>

    Gao-Yuan
<pre class="moz-signature" cols="72">-- 

ZHANG Gao-Yuan

Department of Astronomy, Nanjing University<span style="font-family: Arial; line-height: 1.7;"> </span></pre></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"><span title="neteasefooter"><span id="netease_mail_footer"><a href="#" target="_blank"></a></span></span>
</span></span>