We figured it out, we think it had something to do with our if statements.  We corrected it by just setting everything to 980.0, since it's a very small portion of the domain that needed 327.0.  Thanks!<br><br>Josh<br>
<br><div class="gmail_quote">On Tue, Nov 9, 2010 at 2:08 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="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Sorry,<br>
<br>
I didn't see this till just now...<br>
<br>
Is it still unresolved?<br>
<br>
Paul<br>
<br>
<br>
I would propose using 980. and 327. instead of integers<br>
as the input value.<br>
<br>
Paul<div><div></div><div class="h5"><br>
<br>
<br>
On Thu, 4 Nov 2010, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Neks,<br>
<br>
I am attempting to reset the temperature of a certain part of my flow domain<br>
to a constant value using cfill.  I have had success doing this with the x,<br>
y, and z velocities, but temperature doesn't seem to work.  I have attached<br>
the pertinent code from my usrchk, I would appreciate any help!  Thanks,<br>
<br>
Josh<br>
<br>
     nxyz = nx1*ny1*nz1<br>
     n = nxyz*nelv<br>
<br>
<br>
     if (istep.eq.0.and.time.lt.1.0) then<br>
     do e=1,nelv<br>
       if (xm1(2,2,2,e).gt.0.9) then<br>
         call cfill(vx(1,1,1,e),1.0,nxyz)<br>
         call cfill(vy(1,1,1,e),4.0,nxyz)<br>
         call cfill(vz(1,1,1,e),0.0,nxyz)<br>
       endif<br>
       if (zm1(2,2,2,e).ge.0) then<br>
         call cfill(t(1,1,1,e,1),980,nxyz)<br>
       else<br>
         call cfill(t(1,1,1,e,1),327,nxyz)<br>
       endif<br>
       if (xm1(2,2,2,e).le.-1.3) then<br>
         call cfill(t(1,1,1,e,1),980,nxyz)<br>
       endif<br>
     enddo<br>
<br>
     endif<br>
<br>
</blockquote></div></div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke (probably)<br>