<div dir="ltr">Hi Ketan,<div><br></div><div>Yes there are 288 elements indeed. </div><div><br></div><div>As far as I can see it outpost the rms and TKE in the code. May I know where it output the mean velocity?</div><div><br></div><div><div>      lt = lx1*ly1*lz1*nelt</div><div>      ifld  = 1</div><div><br></div><div>ccc   avg all the avgp10.f0* files</div><div>      call auto_averager(fname1)</div><div><br></div><div>ccc   copy uk=vx**2, vk=vy**2, wk=vz**2</div><div>      call col3(uk,vx,vx,lt)</div><div>      call col3(vk,vy,vy,lt)</div><div>      call col3(wk,vz,vz,lt)</div><div><br></div><div>ccc   avg all the rmsp10.f0* files</div><div>      call auto_averager(fname2)</div><div><br></div><div>ccc   u'^2,v'^2,w'^2 = rms - avg**2</div><div>      call sub2(vx,uk,lt)</div><div>      call sub2(vy,vk,lt)</div><div>      call sub2(vz,wk,lt)</div><div><br></div><div>ccc   Average in the z direction</div><div>      call my_z_avg(uk,vx,gs_avg_hndl,nelxy,ifld)</div><div>      call my_z_avg(vk,vy,gs_avg_hndl,nelxy,ifld)</div><div>      call my_z_avg(wk,vz,gs_avg_hndl,nelxy,ifld)</div><div><span style="white-space:pre">   </span>  </div><div>ccc   TKE = 0.5*(u'^2+v'^2+w'^2)</div><div>      call add4(pk,uk,vk,wk,lt)</div><div>      call cmult(pk,0.5,lt)</div><div><span style="white-space:pre">       </span>  </div><div>ccc   sqrt uk,vk,wk to get urms, vrms and wrms</div><div>      call vsqrt(uk,lt)</div><div>      call vsqrt(vk,lt)</div><div>      call vsqrt(wk,lt)</div><div><span style="white-space:pre">        </span>  </div><div>ccc   output urms**2,vrms**2 and wrms**2 along with the TKE </div><div>ccc   TKE is in pressure field</div><div>      call outpost(uk,vk,wk,pk,t,'   ')</div></div><div><br></div><div>Kind regards,</div><div><br></div><div>ZJ</div><div><br></div><div><br></div></div>