Skip to content
Snippets Groups Projects
Commit 2eadffe3 authored by holger's avatar holger
Browse files

modified: src/post.f90

parent 46880916
No related branches found
No related tags found
No related merge requests found
......@@ -31,7 +31,7 @@
!> @author Holger Grosshans
subroutine monitor
use var
real(kind=pr) :: ucl,avu,tke,gradu,graduy,graduz, &
real(kind=pr) :: ucl,avu,avv,avw,tke,gradu,graduy,graduz, &
Rec,Reb,avvp,avwp,avyp,avzp,avvf,avwf,avxp,avqp, &
syncAv,syncSum
integer :: i,j,l,m
......@@ -44,7 +44,11 @@
ucl= sum(u(imin:imax,j,l))/(imax-imin+1)
avu= sum(abs(u(imin:imax,jmin:jmax,lmin:lmax)) &
avu= sum(u(imin:imax,jmin:jmax,lmin:lmax) &
*vol(imin:imax,jmin:jmax,lmin:lmax)) / volTot
avv= sum(v(imin:imax,jmin:jmax,lmin:lmax) &
*vol(imin:imax,jmin:jmax,lmin:lmax)) / volTot
avw= sum(w(imin:imax,jmin:jmax,lmin:lmax) &
*vol(imin:imax,jmin:jmax,lmin:lmax)) / volTot
rmsv = sum((v(imin:imax,jmin:jmax,lmin:lmax))**2 &
......@@ -87,6 +91,8 @@
ucl=syncAv(ucl)
avu=syncAv(avu)
avv=syncAv(avv)
avw=syncAv(avw)
ubulk=avu
rmsv=syncAv(rmsv)
rmsw=syncAv(rmsw)
......@@ -111,14 +117,14 @@
write(10,'(a,a)') '# ',version
write(10,'(a,i10,i14,15(i10))') '#',(m,m=1,16)
write(10,'(a,a10,a14,15(a10))') &
'#','nt','t','dt','ucl','av(u)','rms(v)','rms(w)', &
'#','nt','t','dt','ucl','av(u)','av(v)','av(w)', &
'Rec','Reb','tau_w','u_tau','Re_tau', &
'delta_v','gradu','dpdx','av(qp)'
'delta_v','gradu','av(qp)'
endif
write(10,'(x,i10,es14.6e2,21(es10.2e2))') &
nt,t,dt,ucl,avu,rmsv,rmsw,Rec,Reb, &
tau_w,u_tau,Re_tau,delta_v,gradu,dpdx,avqp
nt,t,dt,ucl,avu,avv,avw,Rec,Reb, &
tau_w,u_tau,Re_tau,delta_v,gradu,avqp
close(10)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment