#!/bin/sh
echo " "
echo " "
echo " ------------------------- Output Domain ---------------------"
echo " - ESQUIF   [Paris]..........................................1"
echo " - ESCOMPTE [Marseille]......................................2"
echo " "
echo " - France [Nord].............................................3"
echo " - France [Sud]..............................................4"
echo " - Europe ...................................................5"
echo " "
read domain
if [ $domain = "1" ]; then
   latmin=48.5
   latmax=49.2
   longmin=1.7
   longmax=3.
   box="-Bfa0.2/fa0.1WSne"
fi
if [ $domain = "2" ]; then
   latmin=42.7
   latmax=44.2
   longmin=4.1
   longmax=6.2
   box="-Bfa0.2/fa0.1WSne"
    echo "4.2 44.1 15 0 0 0 " > .date3
fi
if [ $domain = "3" ]; then
   latmin=47.5
   latmax=50.5
   longmin=0
   longmax=5
   box="-Bfa0.5/fa0.5WSne"
fi
if [ $domain = "4" ]; then
   latmin=41.5
   latmax=45.5
   longmin=2.5
   longmax=11.5
    box="-Bf1/f.5WSne"
    echo "3.1 45. 15 0 0 0 " > .date3
fi
if [ $domain = "5" ]; then
    latmin=38
    latmax=58
    longmin=-12
    longmax=23.
    box="-Bf2/f2WSne"
fi

dircarto="/home/menut/script"
palette=".palesc"
psca="-D15.5/6./12./1"
pars="-JX15/12.1  -R"$longmin"/"$longmax"/"$latmin"/"$latmax

echo " -- Extension des fichiers a traiter en animation GIF"
echo " "
ls *
echo " -------------------------------------------------------------"
read extens

\rm ps_$extens* $extens*.gif $palette

echo " -- Donner les valeurs MIN et MAX et le pas pour la palette de couleurs"
read vmin vmax pas

echo " "
echo " "
echo " --------------------------- Type de sortie ---------------------"
echo " - Couleur ...................................................1"
echo " - Noir & blanc ..............................................2"
echo " - Uniquement iso-contours ...................................3"
echo " "
read color

if [ $color = "1" ]; then
   makecpt -Crainbow -T${vmin}/${vmax}/${pas} > $palette
fi
if [ $color = "2" ]; then
   makecpt -Cgray -T${vmin}/${vmax}/${pas} -I > $palette
fi
if [ $color = "3" ]; then
   makecpt -Cgray -T${vmin}/${vmax}/${pas} -I > $palette
fi

\rm .dom2
# cadre chimere etendu pour traceur
longminmod=3
longmaxmod=11
latminmod=42.
latmaxmod=45
echo $longminmod $latminmod >> .dom2
echo $longmaxmod $latminmod >> .dom2
echo ">" >> .dom2
echo $longminmod $latmaxmod >> .dom2
echo $longmaxmod $latmaxmod >> .dom2
echo ">" >> .dom2
echo $longminmod $latminmod >> .dom2
echo $longminmod $latmaxmod >> .dom2
echo ">" >> .dom2
echo $longmaxmod $latminmod >> .dom2
echo $longmaxmod $latmaxmod >> .dom2
echo ">" >> .dom2


supr=1
if [ $supr="1" ];then
# suppression du violet dans la palette couleur
gawk '{
       if(FNR==4){
           $2=255
           $3=255
           $4=255
           $6=255
           $7=255
           $8=255
       }
       print $0 >> ".pal2"
       }' $palette
\mv .pal2 $palette
fi

#*****************************************************************************
# Boucle generale sur les fichiers

for date in `ls ${extens}*`
do

\rm .txt .data
gawk '{
     if(FNR==1){print $2 > ".txt"}
     if(FNR>1){print $0 > ".data"}
}' $date

echo "-- Traitement du fichier "$date "   no  " ${count}

pscoast ${pars} ${box} -S205/255/255WS -A10 -Ia/1/0/0/255 -W1 -Df -Na -K -P  > ps_$date

\rm wind.data
gawk '{
       ff=sqrt($5*$5+$4*$4)
       dir=atan2($5,$4)
       if(dir >= 0.)dir2=dir/3.14159*180.
       if(dir < 0.)dir2=360.+dir/3.14159*180.
       if(FNR%30){
          printf("%-f %-f %f %-f \n",$1,$2,dir2,ff/5.) >> "wind.data"
       }
       }' .data


if [ $color = "1" ]; then
  pscontour .data  ${pars} ${box} -C$palette -I -O -P -K -Af1  >> ps_$date
# suppression de la colorbar pour les traceurs
#  psscale -C$palette ${psca} -L -O -K  -P >> ps_$date
fi
if [ $color = "2" ]; then
  pscontour .data  ${pars} ${box} -C$palette -I -O -P -K -Af1  >> ps_$date
  psscale -C$palette ${psca} -L -O -K  >> ps_$date
fi
if [ $color = "3" ]; then
  pscontour .data  ${pars} ${box} -C$palette -K -W2/0 -O -P -Af10  >> ps_$date
fi
#psxy wind.data ${pars}  -Sv0.02/0.1/0.05 -G200 -M -O -P -K >> ps_$date


if [ $domain = "1" ]; then
   psxy $dircarto/ESQ_Relief_Roads ${pars}  -K -M -O -P -W1/20  >> ps_$date
   psxy $dircarto/ESQ_Paris.pge ${pars}  -K -M -O -P -W10/0 >> ps_$date
fi
if [ $domain = "2" ]; then
   pstext $dircarto/names_poiesc ${pars}  -K -O -P  >> ps_$date
fi

# coordonnees des traceurs
\rm .dot
echo "4.4  43.85" >> ".dot"
echo "4.4 43.85 8 0 0 0  Nimes"  > .tit1     
pstext ".tit1"  ${pars} -K  -O -P  >> ps_$date


echo "4.18   43.56" >> ".dot"
echo "4.18 43.56 8 0 0 0  AiguesMortes"  > .tit1     
pstext ".tit1"  ${pars} -K  -O -P  >> ps_$date

echo "5.15  43.45" >> ".dot"
echo "5.15  43.45 8 0 0 0  PntdeBerre"  > .tit1     
pstext ".tit1"  ${pars} -K  -O -P  >> ps_$date

echo "5.38  43.33" >> ".dot"
echo "5.38  43.33 8 0 0 0  Marseille"  > .tit1     
pstext ".tit1"  ${pars} -K  -O -P  >> ps_$date

psxy ".dot" ${pars}  -Sc0.15 -G255/0/0 -M -O -P -K >> ps_$date

psxy ".dom2" ${pars}  -M  -O -P -K -W4/255/0/0 >> ps_$date

pscoast ${pars} ${box}  -A10 -Ia/1/0/0/255 -W1 -Df -Na -K -P -O >> ps_$date

# le titre de la figure
paste .date3 .txt > .date2
pstext .date2 ${pars}  -W240 -O -N -P -K  >>  ps_$date

pstext .txt ${pars} -W255/0/0  -O -P  >>  ps_$date

# Bug dans le set-up des post-scripts
gawk '{
       if($0 !~ "/ImagingBBox")print $0 >> "test.ps"
       }' ps_$date
\mv test.ps ps_$date


echo "              Conversion en format jpg"
convert -crop 0x0  ps_$date " " $date".gif"

done
#*****************************************************************************

echo " ------------------------ MONTAGE --------------------"
# Montage sur unique post-script
#\rm $extens.ps
#montage +frame +shadow +label  -crop 0x0 -geometry 100x50+1+1 -tile 4x6 #ps_$extens* $extens.ps
#ghostview  -font 6x13 -geometry 800x800 $extens.ps &




