<div dir="ltr">Hi, <div><br></div><div>Thanks so much for your reply. </div><div><br></div><div>The visit file obtained after running the usr file attached last email is called pipe0.f00001. Visualization of pipe0.f00001 shows the peak value of streamwise velocity is 0.0003 which should be around 1.3 in non-dimension. So I guess the averaging process is wrong in usr file. Would you please have a look at the averaging calculation in usr file is correct. </div><div><br></div><div>The 2D mesh was generated firstly using pretex and extruded to 3D via n2to3. </div><div><br></div><div>If the averaging process is correct, may I know how I can use Paraview or VisIt to get the line data. By the way, is there any way to get profile like channel flow as shown below,</div><div><br></div><div><div>      if (nid.eq.0 .and. istep.gt.0 .and. </div><div>     &   mod(istep,iostep_avg).eq.0) then</div><div>         write(6,*) 'Dumping statistics ....'</div><div>         open(unit=56,file='reystresses.dat')</div><div>         write(56,'(A,1pe14.7)') '#time = ', time</div><div>         write(56,'(A)')</div><div>     &   '#  y     y+     R_uu     R_vv     R_ww     R_uv'</div><div>         open(unit=57,file='means.dat')</div><div>         write(57,'(A,1pe14.7)') '#time = ', time</div><div>         write(57,'(A)')</div><div>     &   '#  y     y+    Umean'</div></div><div><br></div><div><br></div><div>Kind regards,</div><div><br></div><div>ZJ</div></div>