[Nek5000-users] cfill with temperature

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

Sorry,

I didn't see this till just now...

Is it still unresolved?

Paul


I would propose using 980. and 327. instead of integers
as the input value.

Paul


On Thu, 4 Nov 2010, nek5000-users at lists.mcs.anl.gov wrote:

> 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
>