<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=">