<div dir="ltr">Hi Paul,<br><br>Yes. ldimt is set to 5.<br><br clear="all">Regards<br>Shriram Jagannathan<br>
<br><br><div class="gmail_quote">On 13 December 2010 02:00,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<br>
<br>
How big is ldimt in SIZE ?<br>
<br>
For the given usage, it looks like it should be 5?<br>
<br>
Paul<div><div></div><div class="h5"><br>
<br>
<br>
<br>
 Sun, 12 Dec 2010, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I am having an issue with hpts data being overwritten while using the<br>
comp_vort3() and the hpts() routine. However the data dumps (blah0.f0000)<br>
are uncorrupted. Eddy viscosity is computed and stored in temp array,<br>
lambda2 and vorticity components are stored in 4 passive scalar arrays. It<br>
seems like the last two PS values of previous history points are<br>
over-written on vx and vy. I have attached an image of the hpts.out file and<br>
also pasted below is a piece of the usr file :<br>
<br>
     common /gaaa/    wo1(lx1,ly1,lz1,lelv)<br>
    &              ,  wo2(lx1,ly1,lz1,lelv)<br>
    &              ,  wo3(lx1,ly1,lz1,lelv)<br>
<br>
     if(mod(istep,iostep).eq.0) then<br>
       ! compute lambda2 vortex<br>
       call lambda2(t(1,1,1,1,2))<br>
       ifpsco(1)=.true.<br>
       ! compute vorticity<br>
       if(ldimt.ge.5) call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)<br>
       call hpts()<br>
      ifpsco(2)=.true.<br>
      ifpsco(3)=.true.<br>
      ifpsco(4)=.true.<br>
     endif<br>
<br>
t(1,1,1,1,1) has eddy viscosity and ldimt is set to 5. Also, I am wondering<br>
what ldimt1 and ldimt3 means in the SIZE file ? Thanks.<br>
<br>
Regards<br>
Shriram<br>
<br>
</blockquote></div></div>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br></div>