#!/bin/sh
export LANG=C

echo "efface tout, ne reste que run.sh	turb.gfs et v.gfv"  
rm -rf SIM/*
rmdir SIM
rm *.eps *.mpg
rm vals.data cuty.dat 




mkdir SIM

file="cuty.dat"
echo "----" $file
 awk 'BEGIN{q=2
 for (a = 0.000; a <= 1.01; a=a+q*0.0005) print "0 " a " 0.0 "; }' > $file


####################################

echo "calcul U longueur de Prandtl"
gerris2D -m turb.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set term postscript eps color 14
set output "mupvp.eps"
set xlabel "y"
p'sim.data' u (\$2):9 t"-<u\'v\'> gerris",1-x t '1-y'
FdF

convert mupvp.eps -resize 1600x1200  mupvp.png

cat <<FdF | gnuplot
set key right bottom
set logscale x
ulo(x)=log(0.82*x + pow(1 + 4*pow(0.41,2)*pow(x,2),0.5))*pow(0.41,-1) + (pow(0.41,-2)*pow(x,-1)*(1 - pow(1 + 4*pow(0.41,2)*pow(x,2),0.5)))/2.
pow(x,n)=x**n
set xlabel "y"
set term postscript eps color 14
set output "u_log_turb.eps"
e=2**-8.
p[e:][0:15]'sim.data' u (\$2*500):6 t 'gerris',\
  x t 'u+=y+',ulo(x) t'exact in y+',log(x)/.4-1.23245 t'u+=log(y+)/0.41-1.23245' 
FdF
convert u_log_turb.eps -resize 1600x1200  u_log_turb.png
 
#   p"< awk '{if(fabs($1)<0.01){print $0}}' sim.data" u 2:6
#awk '{if(($1<0.01)&&($1>0)){print $2 " " $6}}' sim.data > profil500.txt
# p'vals.data'u ($3*500):7



tf=`tail -n 1 vals.data  | awk '{print $1}'`   
awk -v ttf=$tf '{if($1==ttf) {print $3 " " $7}}' vals.data> U.data
moyenne=` awk 'BEGIN{y=0;Intg=0}{Intg=Intg+($1-y)*$2;y=$1}END{print Intg}' U.data`