<span style="font-family:SimSun;"> 
<p>
        <span style="font-family:SimSun;">Hi all!</span> 
</p>
<p>
        <span style="font-family:SimSun;"> </span> 
</p>
<p>
        <span style="font-family:SimSun;">I met some problems when trying to add the resistivity η to the simulation
of <strong>magnetoHD/NohCylindrical(user manual 32.2.4)</strong>,</span> 
</p>
<p>
        <span style="font-family:SimSun;"> </span> 
</p>
<p>
        <span style="font-family:SimSun;">First, the units of resistivity troubles me. In the CGS UnitSystem,
whether setting </span><i><span style="font-family:SimSun;"> <strong>resistivity = 1 </strong></span></i><span style="font-family:SimSun;">in </span><i><span style="font-family:SimSun;">“flash.par” </span></i><span style="font-family:SimSun;">means set magnetic resistivity η=1
cm</span><sup><span style="font-family:SimSun;">2</span></sup><span style="font-family:SimSun;">/s=10</span><sup><span style="font-family:SimSun;">-11 </span></sup><span style="font-family:SimSun;">Ω·m and magnetic viscosity
υ</span><sub><span style="font-family:SimSun;">m</span></sub><span style="font-family:SimSun;">=η·c·c/(4·pi), where c=3·10</span><sup><span style="font-family:SimSun;">10 </span></sup><span style="font-family:SimSun;">cm/s?</span> 
</p>
<p>
        <span style="font-family:SimSun;"> </span> 
</p>
<p>
        <span style="font-family:SimSun;">Second, the calculation about resistivity troubles me. Looking through
the kernels and helpers of Diffuse Unit, I can’t find any calculation about resistivity.
Is it in the kernels and helpers of <strong>Hydro Unit</strong>? I can’t find the calculation about <strong>magnetic diffusion</strong> in the Hydro Unit. Is it needed to </span><em><span style="font-family:Times New Roman;"><strong>REQUIRE Diffuse</strong></span></em><span style="font-family:SimSun;">?</span> 
</p>
<p>
        <span style="font-family:SimSun;"> </span> 
</p>
<p>
        <span style="font-family:SimSun;">Finally, When I set these parameter in </span><i><span style="font-family:SimSun;">“flash.par”</span></i><span style="font-family:SimSun;"> as below and setup with </span><i><span style="font-family:Times New Roman;"><strong>+usm</strong></span></i> 
</p>
<blockquote style="padding:0px;margin-top:0px;margin-right:0px;margin-bottom:0px;" dir="ltr">
        <p>
                <i><span style="font-family:Times New Roman;"><strong>#     Units
(none, SI or CGS)</strong></span></i> 
        </p>
        <p>
                <i><span style="font-family:Times New Roman;"><strong>UnitSystem        = "CGS"</strong></span></i> 
        </p>
        <p>
                <i><span style="font-family:Times New Roman;"><strong>##   
Resistive MHD parameters     ##</strong></span></i> 
        </p>
        <p>
                <i><span style="font-family:Times New Roman;"><strong>useDiffuse    =
.true.</strong></span></i> 
        </p>
        <p>
                <i><span style="font-family:Times New Roman;"><strong> dt_diff_factor 
= 0.8</strong></span></i> 
        </p>
        <p>
                <i><strong><span style="font-family:Times New Roman;">useMagneticResistivity     = .true.</span></strong></i> 
        </p>
        <p>
                <i><strong><span style="font-family:Times New Roman;"><span class="ke-content-forecolor" style="color:#E53333;"> resistivity          = 0.0</span> # magnetic Resistivity</span></strong></i> 
        </p>
</blockquote>
<p>
        mpirun stopped with the same message, even if setting a large range of <i>resistivity</i>, from 10<sup>20</sup> to 10<sup>-40</sup>
</p>
<p>
        <b><i>      DRIVER_ABORT:
[Hydro]: Computed dt is not positive! Aborting!</i></b> 
</p>
<p>
        <b><i>      dtCheck=   0.0000000000000000</i></b><b></b> 
</p>
<p>
        <br>
</p>
<p>
        How to solve these problem? Really need someone to help!
</p>
<p>
         
</p>
<p>
        Zhang
</p>
</span><br>