[Nek5000-users] cfill with temperature

[nek5000-users at lists.mcs.anl.gov]

Neks,

I am attempting to reset the temperature of a certain part of my flow domain
to a constant value using cfill.  I have had success doing this with the x,
y, and z velocities, but temperature doesn't seem to work.  I have attached
the pertinent code from my usrchk, I would appreciate any help!  Thanks,

Josh

      nxyz = nx1*ny1*nz1
      n = nxyz*nelv


      if (istep.eq.0.and.time.lt.1.0) then
      do e=1,nelv
        if (xm1(2,2,2,e).gt.0.9) then
          call cfill(vx(1,1,1,e),1.0,nxyz)
          call cfill(vy(1,1,1,e),4.0,nxyz)
          call cfill(vz(1,1,1,e),0.0,nxyz)
        endif
        if (zm1(2,2,2,e).ge.0) then
          call cfill(t(1,1,1,e,1),980,nxyz)
        else
          call cfill(t(1,1,1,e,1),327,nxyz)
        endif
        if (xm1(2,2,2,e).le.-1.3) then
          call cfill(t(1,1,1,e,1),980,nxyz)
        endif
      enddo

      endif
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101104/89da288b/attachment.html>