Hi Josh,<br><br>Yes, changing the initial field in the first call of userchk worked, thanks.<br><br>David<br><br><div class="gmail_quote">On Tue, Oct 4, 2011 at 6:21 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div><div style="font-family:Calibri,sans-serif;font-size:11pt">Hi David,<br><br>The best place would probably be in userchk.<br>


<br>You could have a block like<br><br>if (istep.eq.0) then<br><br>C........do stuff<br><br>endif<br><br>If you are doing something else in userchk that would depend on your modification of the field quantities, then you would want  this to come first.<br>


<br>Just in case, the velocity fields are accessed like<br><br>vx(i,j,k,e) ( for u)<br><br>and likewise for vy and vz, where i,j,k go from 1 to nx1 (ny1, nz1) and e goes from 1 to nelv.  Sorry if that was more information than you needed, just wanted to make sure all my bases were covered :-).<br>


<br>Does this help?<br><br>Josh Camp<br></div></div></div></blockquote></div>