<div name='mp-default' style='font-family:돋움,Dotum;font-size:10pt;'><p>Hi, Neks</p><p> </p><p>I wonder why planar_average_z in navier5.f shows wrong result.</p><p>here I attached 2 images.</p><p> </p><p>First image shows the z-directional grids points from zm1</p><p>(I wrote a code to extract z coordinate at y = 0 as follows:</p><p>======================================================</p><p>insnum = nid</p><p>insuni = insnum + 1000</p><p>WRITE(insnam,98) 'inst_field_', insnum</p><p>OPEN(unit=insuni,file=insnam</p><p><span style="font-size: 13.3333px;">      insnum = nid</span></p><p><span style="font-size: 13.3333px;">      insuni = insnum + 1000</span></p><p><span style="font-size: 13.3333px;">      WRITE(insnam,98) 'inst_field_', insnum</span></p><p><span style="font-size: 13.3333px;">      Open(unit=insuni,file=insnam)</span></p><p><span style="font-size: 13.3333px;"> </span></p><p><span style="font-size: 13.3333px;">      do i = 1,n</span></p><p><span style="font-size: 13.3333px;">      if (ym1(i,1,1,1).eq.0.0)then</span></p><p><span style="font-size: 13.3333px;">      WRITE(insuni,5) xm1(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">     %                  ,  zm1(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">     %                  ,  vx(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">     %                  ,  vy(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">     %                  ,  vz(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">     %                  ,  phig(i,1,1,1)</span></p><p><span style="font-size: 13.3333px;">      endif</span></p><p><span style="font-size: 13.3333px;"> </span></p><p><span style="font-size: 13.3333px;">      enddo</span></p><p><span style="font-size: 13.3333px;">      CLOSE(insuni)</span></p><p><span style="font-size: 13.3333px;"></span><span style="font-size: 10pt;">======================================================</span>​</p><p><span style="font-size: 13.3333px;"> </span></p><p><span style="font-size: 13.3333px;">Then, I can collect lots of field data from each processor and sort data with ascending order using matlab.</span></p><p>So, I got the result as 1st image. </p><p><span style="font-size: 13.3333px;"> </span></p><p><span style="font-size: 13.3333px;">But when I use plannar_average_z </span></p><p>I got the result as shown in 2nd image.</p><p><span style="font-size: 10pt;">1. I really want to modify the problems, so could you give me your idea on this problem?</span> </p><p>2. Additionally, could you explain 'zgm1' and 'area'?</p><p>I noticed that, planar_average_r or s ... in turbChannel case has also 'area' but, those case are used to calculate on <span style="font-size: 10pt;">different direction.<br>But they used same input parameters. </span></p><p><span style="font-size: 10pt;"> </span></p><p><span style="font-size: 10pt;">What I understand about area is, surface area of element.</span></p><p><span style="font-size: 10pt;">But I thought that different direction needs different indexing </span></p><p><span style="font-size: 10pt;"> </span></p><p>=======================================================</p><p><span style="font-size: 13.3333px;">      subroutine planar_average_z(ua,u,w1,w2)</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">c     Compute r-s planar average of quantity u()</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      include 'SIZE'</span></p><p><span style="font-size: 13.3333px;">      include 'GEOM'</span></p><p><span style="font-size: 13.3333px;">      include 'PARALLEL'</span></p><p><span style="font-size: 13.3333px;">      include 'WZ'</span></p><p><span style="font-size: 13.3333px;">      include 'ZPER'</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      real ua(lz1,nelz),u(lx1*ly1,lz1,nelv),w1(lz1,nelz),w2(lz1,nelz)</span></p><p><span style="font-size: 13.3333px;">      integer e,eg,ez</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      melxy = nelx*nely</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      nz = lz1*nelz</span></p><p><span style="font-size: 13.3333px;">      call rzero(ua,nz)</span></p><p><span style="font-size: 13.3333px;">      call rzero(w1,nz)</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      do e=1,nelt</span></p><p><span style="font-size: 13.3333px;">c </span><span style="font-size: 10pt;">        </span></p><p><span style="font-size: 10pt;">          eg = lglel(e)</span></p><p><span style="font-size: 13.3333px;">c         ez = 1 + (eg-1)/melxy</span></p><p><span style="font-size: 13.3333px;">         call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">         do k=1,lz1</span></p><p><span style="font-size: 13.3333px;">         do i=1,lx1*ly1</span></p><p><span style="font-size: 13.3333px;">            zz = (1.-zgm1(k,3))/2.  ! = 1 for k=1, = 0 for k=lz1</span></p><p><span style="font-size: 13.3333px;">            aa = zz*area(i,1,5,e) + (1-zz)*area(i,1,6,e)  ! wgtd jacobian</span></p><p><span style="font-size: 13.3333px;">            w1(k,ez) = w1(k,ez) + aa</span></p><p><span style="font-size: 13.3333px;">            ua(k,ez) = ua(k,ez) + aa*u(i,k,e)</span></p><p><span style="font-size: 13.3333px;">         enddo</span></p><p><span style="font-size: 13.3333px;">         enddo</span></p><p><span style="font-size: 13.3333px;">      enddo</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      call gop(ua,w2,'+  ',nz)</span></p><p><span style="font-size: 13.3333px;">      call gop(w1,w2,'+  ',nz)</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      do i=1,nz</span></p><p><span style="font-size: 13.3333px;">         ua(i,1) = ua(i,1) / w1(i,1)   ! Normalize</span></p><p><span style="font-size: 13.3333px;">      enddo</span></p><p><span style="font-size: 13.3333px;">c</span></p><p><span style="font-size: 13.3333px;">      return</span></p><p><span style="font-size: 13.3333px;">      end</span></p><p><span style="font-size: 10pt;">======================================================</span></p><p><span style="font-size: 10pt;"> </span></p><p><span style="font-size: 10pt;">Regards,</span></p><p><span style="font-size: 10pt;"> </span></p><p><span style="font-size: 10pt;">Hyunduk.​</span></p>





<p><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">------------------------------</span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">------------------------------------------------------------------------------------</span></p>
<p><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">Hyunduk Seo (Mr.)</span></p>
<p><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">Integrated Ph.D. course, Schoo</span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">l of Mechanical Engineering, </span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">Pusan National University</span><span style="font-size: 10pt;">​</span></p>
<p><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">Experimental Thermo-Flui</span><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">ds Mec</span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">hanics and Energy Systems (ExT</span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">ENsys) Laboratory</span></p>
<p><a href="http://fluid.pusan.ac.kr/" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=ko&q=http://fluid.pusan.ac.kr&source=gmail&ust=1519689597878000&usg=AFQjCNFkiLQ2PCtT4IicKeK4iyZhHEx49Q" style="color: rgb(17, 85, 204); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">http://fluid.pusan.ac.kr</a><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"> </span><br style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">
<span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">T.: +82-51-510-1536 M.: +82-10</span><wbr style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);"><span style="color: rgb(34, 34, 34); font-family: arial, sans-serif; font-size: 14px; background-color: rgb(255, 255, 255);">-2964-0994</span></p>

<p> </p></div><img src="http://webmail.pusan.ac.kr/api/notify.php?home=MjAxODA3MjQxNTI5MDF8bmVrNTAwMC11c2Vyc0BsaXN0cy5tY3MuYW5sLmdvdnxoZHNlbzAxMjNAcHVzYW4uYWMua3I=">