<div dir="ltr">Dear neks,<div>I am trying to user comp_sij to compute dissipation in every point of the domain.</div><div>The data storage is organised according to </div><div><div>      real sij(lx1*ly1*lz1,nij,lelv)</div></div><div>so the index of a component of tensor</div><div>precedes number of element, and </div><div>the is no continuous block of lx1*ly1*lz1,lelv points for an element of tensor.<br></div><div>And so, the global functions like col3 of glsc can not be applied.</div><div>At least I don't understand how.</div><div><br></div><div>To bypass this, I copied comp_sij to my usr file and changed the order of</div><div>last indices in sij</div><div><div>      real sij(lx1*ly1*lz1,lelv,nij)</div></div><div><div><br></div><div><br></div><div>After that I use the followng code to compute viscous dissipation</div><div><br></div><div>       nij = 3</div><div>       if (if3d) nij=6</div><div><br></div><div>      call il_comp_sij(sij,nij,vx,vy,vz,ur,us,ut,vr,vs,vt,wr,ws,wt)</div><div><br></div><div>      call rzero(il_dissip,lv)</div><div><br></div><div>      do j=1,nij</div><div>      call col3(sij2,sij(1,1,1,1,j),sij(1,1,1,1,j),nxyz)</div><div>      il_dissip=il_dissip+sij2</div><div>       if (if3d .and. j.ge.4) then</div><div>      il_dissip=il_dissip+sij2</div><div>       end if</div><div>      if ( .not. if3d .and. j.eq.3 ) then</div><div>      il_dissip=il_dissip+sij2</div><div>      end if</div><div>      il_dissip=2*param(2)*il_dissip</div><div>      enddo</div><div><br></div></div><div>Is is a fair way to compute dissipation?</div><div><br></div><div>Thanks, Ilias</div><div><br></div><div><br></div></div>