<div dir="ltr">Hello,<div><br></div><div>I am having an issue with hpts data being overwritten while using the comp_vort3() and the hpts() routine. However the data dumps (blah0.f0000) are uncorrupted. Eddy viscosity is computed and stored in temp array, lambda2 and vorticity components are stored in 4 passive scalar arrays. It seems like the last two PS values of previous history points are over-written on vx and vy. I have attached an image of the hpts.out file and also pasted below is a piece of the usr file :</div>

<div><br></div><div><div>      common /gaaa/    wo1(lx1,ly1,lz1,lelv)</div><div>     &              ,  wo2(lx1,ly1,lz1,lelv)</div><div>     &              ,  wo3(lx1,ly1,lz1,lelv)</div></div><div><br></div><div><div>

      if(mod(istep,iostep).eq.0) then</div><div>        ! compute lambda2 vortex</div><div>        call lambda2(t(1,1,1,1,2))</div><div>        ifpsco(1)=.true.</div><div>        ! compute vorticity</div><div>        if(ldimt.ge.5) call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)</div>

<div>        call hpts()</div><div>       ifpsco(2)=.true.</div><div>       ifpsco(3)=.true.</div><div>       ifpsco(4)=.true.</div><div>      endif</div></div><div><br></div><div>t(1,1,1,1,1) has eddy viscosity and ldimt is set to 5. Also, I am wondering what ldimt1 and ldimt3 means in the SIZE file ? Thanks.</div>

<div><br></div><div>Regards<br>Shriram<br>
</div></div>