#!/bin/sh
export LANG=C

echo "efface tout, ne reste que run.sh	pc1_fourier.gfs	et v.gfv"  
rm -rf TURE/*
rmdir TURE
rm *.eps *.mpg
rm vals.data cut.dat
rm pc1_fourierBi.gfs




mkdir TURE

file="cut.dat"
echo "----" $file
 awk 'BEGIN{
 for (a = 0.; a <= 2; a += 0.025) print a " " " 0.0 0.0 "; }' > $file


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

echo "calcul T imposee standard"
gerris2D pc1_fourier.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set key right bottom
set term postscript eps color 14
set output "ture_erf.eps"
p[0:5] "< awk '{if(\$1<0.1)print \$0}' vals.data"  u ((1-\$2)/(\$1**.5)):9 t'gerris'w p ,(erf((x)/2)) t'exact'w l
FdF

convert ture_erf.eps ../img/ture_erf.png

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

echo "calcul T imposee grand Biot=1000"
sed -e 's/{ BcDirichlet T 0 }/{ GfsBcNavier T 0 .00100 }/g' ./pc1_fourier.gfs  >  ./pc1_fourierBi.gfs
gerris2D pc1_fourierBi.gfs | gfsview2D v.gfv


cat <<FdF | gnuplot
set key right bottom
set key right bottom
set term postscript eps color 14
set output "ture_erf_GB.eps"
p[0:5] "< awk '{if(\$1<0.1)print \$0}' vals.data"  u ((1-\$2)/(\$1**.5)):9 t'gerris'w p ,(erf((x)/2)) t'erf Bi=infinity'w l
FdF

convert ture_erf_GB.eps ../img/ture_erf_GB.png

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

echo "calcul dT/dy=-Bi T, Bi=10"
sed -e 's/{ BcDirichlet T 0 }/{ GfsBcNavier T 0 .100 }/g' ./pc1_fourier.gfs  >  ./pc1_fourierBi.gfs
gerris2D pc1_fourierBi.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set key right bottom
set ylabel "T(x,t) Bi=10"
Bi=1/0.1
set term postscript eps color 14
set output "ture_erf_B10.eps"
p "< awk 'BEGIN{x=0;}{if(x<\$2){if(\$1<0.2)print \$0}else{print \" \"};x=\$2}' vals.data" u ((1-\$2)):9t'gerris' w lp,\
''u 2:(erf((\$2/\$1**.5)/2))+exp(Bi*(\$2+\$1*Bi))*erfc(\$2/(\$1**.5)/2 +Bi*\$1**.5)  t'exact Bi=0.1'w l
FdF

convert ture_erf_B10.eps ../img/ture_erf_B10.png

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

echo "calcul dT/dy=-Bi T, Bi=0.1"
sed -e 's/{ BcDirichlet T 0 }/{ GfsBcNavier T 0 10.000 }/g' ./pc1_fourier.gfs  >  ./pc1_fourierBi.gfs
gerris2D pc1_fourierBi.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set key right bottom
set ylabel "T(x,t) Bi=0.1"
Bi=1/10.
set term postscript eps color 14
set output "ture_erf_Bp1.eps"
p "< awk 'BEGIN{x=0;}{if(x<\$2){if(\$1<0.15)print \$0}else{print \" \"};x=\$2}' vals.data" u ((1-\$2)):9t'gerris' w lp,\
''u 2:(erf((\$2/\$1**.5)/2))+exp(Bi*(\$2+\$1*Bi))*erfc(\$2/(\$1**.5)/2 +Bi*\$1**.5)  t'exact Bi=10'w l
FdF

convert ture_erf_Bp1.eps ../img/ture_erf_Bp1.png

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

echo "calcul dT/dy=-Bi T, Bi=1."
sed -e 's/{ BcDirichlet T 0 }/{ GfsBcNavier T 0 1.000 }/g' ./pc1_fourier.gfs  >  ./pc1_fourierBi.gfs
gerris2D pc1_fourierBi.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set key right bottom
set xlabel "x"
set ylabel "T(x,t) Bi=1"
Bi=1.
set term postscript eps color 14
set output "ture_erf_Bi1.eps"
p "< awk 'BEGIN{x=0;}{if(x<\$2){if(\$1<0.15)print \$0}else{print \" \"};x=\$2}' vals.data" u ((1-\$2)):9t'gerris' w lp,\
''u 2:(erf((\$2/\$1**.5)/2))+exp(Bi*(\$2+\$1*Bi))*erfc(\$2/(\$1**.5)/2 +Bi*\$1**.5)  t'exact Bi=10'w l
FdF

convert ture_erf_Bi1.eps ../img/ture_erf_Bi1.png
 


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

echo "calcul dT/dy=-Bi T, Bi=1"
#sed -e 's/{ BcDirichlet T 0 }/{ GfsBcNavier T 0 1.0000 }/g' ./pc1_fourier.gfs  >  ./pc1_fourierBi.gfs
#gerris2D pc1_fourierBi.gfs | gfsview2D v.gfv

cat <<FdF | gnuplot
set key right bottom
set term postscript eps color 14
set output "ture_erf_Bi1tl.eps"
set xlabel "x"
set ylabel "T(x,t) Bi=1"
TBi1(t,x)= 1.1192*exp(-(0.8603)**2 * t)*cos(0.8603*x)\
        -0.1517*exp(-(3.4256)**2 * t)*cos(3.4256*x)\
        +0.0466*exp(-(6.4376)**2 * t)*cos(6.4376*x)\
        -0.0217*exp(-(9.5293)**2 * t)*cos(9.5293*x); 
p[][0:1] "< awk 'BEGIN{x=0;}{if(x<\$2){if(\$1>0.05)print \$0}else{print \" \"};x=\$2}' vals.data" u 2:9t'gerris' w lp,\
'' u 2:(TBi1(\$1,\$2)) t'4 modes' w l
FdF

convert ture_erf_Bi1tl.eps ../img/ture_erf_Bi1tl.png
        
cat <<FdF | gnuplot
set key right bottom
set term postscript eps color 14
set output "ture_erf_Bi1tlc.eps"
set xlabel "t"
set ylabel "T(x=0,t) Bi=1"
TBi1(t,x)= 1.1192*exp(-(0.8603)**2 * t)*cos(0.8603*x)\
        -0.1517*exp(-(3.4256)**2 * t)*cos(3.4256*x)\
        +0.0466*exp(-(6.4376)**2 * t)*cos(6.4376*x)\
        -0.0217*exp(-(9.5293)**2 * t)*cos(9.5293*x); 
p "< awk '{if(\$2==0.0)print \$0}' vals.data" u 1:9t'gerris' w lp,\
'' u 1:(TBi1(\$1,0)) t'4 modes' w l
FdF

convert ture_erf_Bi1tlc.eps ../img/ture_erf_Bi1tlc.png