Hi Nek's,<div><br></div><div>I am using Nek  in perturbation mode. At a certain time step i need to replace the velocity field with a different  one. </div><div><br></div><div>In usrchk I replace (vxp,vyp,vzp) with my field (vxb,vyb,vzb). I also  reset vxlagp,vylagp,vzlagp and the extrapolation terms exxp1,exyp1,exzp1 and exxp2,exyp2,exzp2.</div>
<div><br></div><div><div>               call opcopy(vxp,vyp,vzp,vxb,vyb,vzb)</div><div>               call opzero(exx1p,exy1p,exz1p)</div><div>               call opzero(exx2p,exy2p,exz2p)</div><div>               do i=1,lorder-1</div>
<div>                  call opzero(vxlagp(1,i,1),vylagp(1,i,1), vzlagp(1,i,1))</div><div>               enddo</div></div><div><br></div><div>It works well if I set Torder=1, but for Torder=2 it seems that there are data contamination from previous time steps, so I think I'm missing some variables which i need to modify.</div>
<div><br></div><div>Andrea</div>