<div dir="ltr">Dear Neks,<div><br></div><div>I am trying to average a box geometry along the xz-plane using the planar_average_s subroutine from turbChannel. But it is not outputting all the y values from my domain.</div><div><br></div><div>I am using the average file from avg_all as a restart file and hope to output the planar averaged velocity and passive scalars along the y direction. Could someone help me fix my code?<br></div><div><br></div><div>Below is what I have so far:</div><div><br></div><div><div>      call planar_average_s(yy,ym1,w1,w2)</div><div><br></div><div>      call planar_average_s(uavg_pl,vx,w1,w2)</div><div>      call planar_average_s(vavg_pl,vy,w1,w2)</div><div>      call planar_average_s(wavg_pl,vz,w1,w2)</div><div>      call planar_average_s(tavg_pl,t(:,1,1,1,1),w1,w2)</div><div>      call planar_average_s(ttavg_pl,t(:,1,1,1,2),w1,w2)</div><div>      call planar_average_s(tttavg_pl,t(:,1,1,1,3),w1,w2)</div></div><div><br></div><div><div>      if(nid.eq.0) then</div><div>      open(unit=57,file='mean_prof.dat')</div><div>      do i=1,ny1*nely/2</div><div>        write(57,3) yy(i)</div><div>     &            , uavg_pl(i)</div><div>     &            , vavg_pl(i)</div><div>     &            , wavg_pl(i)</div><div>     &            , tavg_pl(i)</div><div>     &            , ttavg_pl(i)</div><div>     &            , tttavg_pl(i)</div><div>      enddo</div><div>   3  format(7f17.9)</div><div>      close(57)</div><div>      endif</div></div><div>      call exitt</div><div><br></div><div><br></div><div>Thanks for any help,</div><div>Jon</div><div><br></div></div>