<div dir="ltr"><div>Hi Paul,</div><div><br></div><div>If I have it this way, the hpts.out has pressure data for istep=0 but after that my pressure values are not written. </div><div><br></div><div><div>      call hpts()</div>

<div><br></div><div>      ifpsco(1) = .true.</div><div>      ntot1 = nx1*ny1*nz1*nelv</div><div><br></div><div>      if(mod(istep,iostep).eq.0 .and. istep.gt.0) then</div><div>        call lambda2(t(1,1,1,1,1))</div><div>

        ifto=.true.</div><div>        if(ldimt.gt.1) call comp_vort3(t(1,1,1,1,2),wo1,wo2,vx,vy,vz)</div><div>      endif</div></div><div><br></div><div>It looks like using the temperature field and passive scalar has messed up the o/p in hpts. I require the pressure values but not the passive scalars(however it would be nice to have it). Is there a way to force hpts not to o/p any passive scalar ? </div>

<br clear="all">Regards<br>Shriram<br>
<br></div>