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

modified: monitor.plt

	modified:   src/post.f90
parent d206f832
No related branches found
No related tags found
No related merge requests found
#set term postscript enhanced eps 10 color linewidth 1 'Times-Roman'
#set term dumb
#set output "post.eps"
#set xrange [0:]
set multiplot layout 3,2 title 'pafiX -- by H. Grosshans'
set multiplot layout 4,4 title 'pafiX -- by H. Grosshans'
f = 'output/monitor.dat'
header=system("head -3 ".f." | tail -1")
set key box opaque width 2 t l
set xlabel offset 0,0.5
set ylabel offset 2
set linetype 1 ps 0.5 lc 3
set linetype 2 ps 0.5 lc 4
set xlabel 's'
set xlabel word(header,3) #header columns are shifted by 1 bc of #
set ylabel 'm/s'
p f u 2:4 t 'u_{cl}', f u 2:5 t 'u_{av}'
set ylabel '-'
p f u 2:8 t 'Re_c', f u 2:9 t 'Re_b'
set ylabel 'm/s'
p f u 2:6 t 'v_{rms}', f u 2:7 t 'w_{rms}'
set ylabel '-'
p f u 2:($6/$11) t 'v_{rms}/u_{tau}', f u 2:($7/$11) t 'w_{rms}/u_{tau}'
#set ylabel '1/s'
#p f u 2:14 t 'gradu'
set ylabel 'N/m**2'
p f u 2:10 t 'tau_w'
set ylabel '-'
p f u 2:12 t 'Re_{tau}'
do for [i=4:16] {
set ylabel word(header,i+1)
p f u 2:i notitle ls 1
}
pause 20
reread
......@@ -119,7 +119,7 @@
write(10,'(a,a10,a14,14(a10))') &
'#','nt','t','dt','ucl','av(u)','av(v)','av(w)','rms(v)','rms(w)', &
'Rec','Reb','tau_w','u_tau','Re_tau', &
'delta_v','avqp'
'delta_v','av(qp)'
endif
write(10,'(x,i10,es14.6e2,20(es10.2e2))') &
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment