<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Hi Ilias,
<div><br>
</div>
<div>opcopy is for vector to vector copies - so that would not be correct usage.</div>
<div><br>
</div>
<div>I would just look carefully at the hpts constructs if you want to interpolate your </div>
<div>dissipation.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF542981" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Monday, November 09, 2015 9:34 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] dissipation<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>An afterquesion;</div>
Am I right that  I use the obtained field of dissipation 
<div>for hpts output through</div>
<div>
<div>      call opcopy(vx,vy,pr,<span style="color:rgb(0,0,0); font-family:Courier; font-size:13.3333330154419px">sij2,vy,vz</span>) </div>
</div>
<div>      hpts()</div>
<div>?</div>
<div><br>
</div>
<div>Thank you, Ilias</div>
<div><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Mon, Nov 9, 2015 at 5:03 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">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">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div><br>
</div>
Hi Ilias,
<div><br>
</div>
<div><br>
</div>
<div>How about:</div>
<div><br>
</div>
<div>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      subroutine dissipation(sij2)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      include 'SIZE'</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      include 'TOTAL'</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      real sij2(lx1*ly1*lz1,lelt)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      common /scrns/         sij (lx1*ly1*lz1,6,lelv)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      parameter (lr=lx1*ly1*lz1)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      common /scruz/         ur(lr),us(lr),ut(lr)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
     $                     , vr(lr),vs(lr),vt(lr)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
     $                     , wr(lr),ws(lr),wt(lr)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      integer e</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      nij = 3+3*(ndim-2)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      call comp_sij(sij,nij,vx,vy,vz,ur,us,ut,vr,vs,vt,wr,ws,wt)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      nxyz = nx1*ny1*nz1</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      n    = nx1*ny1*nz1*nelv</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      call rzero(sij2,n)</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      do e=1,nelv</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      do j=1,nij</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      do i=1,nxyz</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
         sij2(i,e)=sij2(i,e)+sij(i,j,e)**2</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      enddo</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      enddo</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      enddo</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier; min-height:16px">
<br>
</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      return</p>
<p style="margin-right:0px; margin-left:0px; line-height:normal; font-family:Courier">
      end</p>
<div><br>
</div>
<div>?</div>
<div><br>
</div>
<div>This should function correctly.</div>
<div><br>
</div>
<div>Best, Paul</div>
<div><br>
</div>
<div><br>
</div>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr>
<div style="direction:ltr"><font face="Tahoma" size="2" color="#000000"><b>From:</b>
<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of
<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>]<br>
<b>Sent:</b> Sunday, November 08, 2015 4:37 PM<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
<b>Subject:</b> [Nek5000-users] dissipation<br>
</font><br>
</div>
<div></div>
<div>
<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>
</div>
</div>
</div>
</div>
</div>
<br>
_______________________________________________<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" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</body>
</html>