[Nek5000-users] Hpts output & lambda2

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

Are wo1 and wo2 properly declared in the calling routine ?

Paul


On Thu, 9 Sep 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hi,
>
> I am trying to compute lambda2 vortex just like how it is done in the
> turbchannel example. But it seems like that it messes up my output of hpts.
> The values of pressure are not written and I suppose instead of pressure the
> t(1,1,1,1,1) or t(1,1,1,1,2) is written. I tried setting the ifto and ifpsco
> value to false but it doesn't seem to work. Any thoughts on where I am going
> wrong ?
>
>     if(mod(istep,iostep).eq.0 .and. istep.gt.0) then
>        call lambda2(t(1,1,1,1,1))
>        ifto=.true.
>        if(ldimt.gt.1) call comp_vort3(t(1,1,1,1,2),wo1,wo2,vx,vy,vz)
>        ifpsco(1)=.true.
>      endif
>
>      ifto=.false.
>      ifpo=.true.
>      ifpsco(1)=.false.
>      if(istep.gt.0) call hpts()
>
> Thanks
> Shriram
>