#!/bin/sh
# Script de lecture des resultats meteorologiques et chimiques
# Abscisse = Heure, Ordonnee = Maillage vertical en metres
# Version pour GMT 3.3.1 - Juin 1999

xmin=3.
xmax=11.
zmin=0.
zmax=3.
box="-Bf0.5/f0.2/WeSn"
pars="-JX15/10  -R"$xmin"/"$xmax"/"$zmin"/"$zmax
psca="-D15.5/5./10./.7"
palette="GMT_palette_Temp.cpt"
echo ${xmin} ${zmax} " 15 0 0 0 " > .date3

\rm $palette wind.data box_wind wind2.data

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 " ----------------------- Type de sortie -------------------"
echo " - Couleur (rainbow)......................................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
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
makecpt -Cgray -T-999/0/999  > .palblack

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

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

\rm .txt .data .topo

gawk '{
     if(FNR==1){print $2 > ".txt"}
     if(FNR>1 && $3!=-999){print $0 >> ".data"}
     if(FNR>1 && $3==-999 && $2!=0.){print $0 >> ".topo"}
}' $date
echo ${xmax} "   0  -999." >> .topo
echo ${xmin} "   0  -999." >> .topo


echo "-- Traitement du fichier "$date

if [ $color = "1" ]; then
   pscontour .data  ${pars} -C$palette -I ${box} -P -K  > ps_$date
#   pscontour .topo  ${pars} -C".palblack" -I ${box} -P -K -O >> ps_$date
    psxy .topo ${pars}  -G0 -P -K -O >> ps_$date
#   psscale -C$palette ${psca} -I -L -O -K >> ps_$date
fi
if [ $color = "2" ]; then
   pscontour .data  ${pars} -C$palette -I ${box} -P -K  >> ps_$date
   psscale -C$palette ${psca} -L -O -K >> ps_$date
fi
if [ $color = "3" ]; then
   pscontour .data  ${pars} -C$palette  ${box} -P -K -W2/0 >> ps_$date
fi

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

echo " --  Correction of GMT 3.3.1 bug"
# 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 gif"
convert -crop 0x0  ps_$date " " $date".gif"

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


