[Nek5000-users] Getting planar average and outputting data
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Feb 6 23:42:27 CST 2018
Dear Neks,
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.
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?
Below is what I have so far:
call planar_average_s(yy,ym1,w1,w2)
call planar_average_s(uavg_pl,vx,w1,w2)
call planar_average_s(vavg_pl,vy,w1,w2)
call planar_average_s(wavg_pl,vz,w1,w2)
call planar_average_s(tavg_pl,t(:,1,1,1,1),w1,w2)
call planar_average_s(ttavg_pl,t(:,1,1,1,2),w1,w2)
call planar_average_s(tttavg_pl,t(:,1,1,1,3),w1,w2)
if(nid.eq.0) then
open(unit=57,file='mean_prof.dat')
do i=1,ny1*nely/2
write(57,3) yy(i)
& , uavg_pl(i)
& , vavg_pl(i)
& , wavg_pl(i)
& , tavg_pl(i)
& , ttavg_pl(i)
& , tttavg_pl(i)
enddo
3 format(7f17.9)
close(57)
endif
call exitt
Thanks for any help,
Jon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180206/cd400143/attachment.html>
More information about the Nek5000-users
mailing list