
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">

<p>I guess I will answer my own question. It turns out that by loading the file the way I was doing it (calling load_fld() in userchk()), the program tries to calculate the convection terms using previous timesteps, which were not specified in my case. This

 produces the NaN during the calculation of the RHS. </p>

<p><br>

</p>

<p>I solved my problem by setting my *.rea file to:</p>

<p><br>

</p>

<p></p>

<div><span>    </span>1 PRESOLVE/RESTART OPTIONS  *****</div>

<div>new0.f00001</div>

<div><br>

</div>

<div>Then in userchk() I copy the initial field to auxiliary common arrays (ub and vb, similar to those used in the "blasius" example) during istep=0 which remain unchanged throughout the simulation and are used to set the appropriate boundary conditions. </div>

<div><br>

</div>

<div>Juan Diego</div>

<div><br>

</div>

<p></p>

</div>

</body>

</html>