#!/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


zmin=0.
zmax=3.
xmin=4.3
xmax=11.


echo " -------------------------------------------------------------"
echo " -- Destruction des fichiers temporaires "
\rm $palette wind.data box_wind wind2.data

echo " "
echo " -- Name of data file"
echo " "
ls c2f_cs*
echo " -------------------------------------------------------------"
read date



box="-Bf1/f0.2/WeSn"
pars="-JX15/10  -R"$xmin"/"$xmax"/"$zmin"/"$zmax
psca="-D15.5/5./10./.7"
palette="GMT_palette_Temp.cpt"

echo " -- Nettoyage du precedent fichier "$date".ps"
\rm $date.ps .nbline

gawk '{
   if(FNR==1)vmin=10000
   if(FNR==1)vmax=0
   if($3 <= vmin && $3 != -999. && $3 != 999.){vmin=$3}
   if($3 >= vmax && $3 != -999. && $3 != 999.){vmax=$3}
   nline=FNR
   if(nline % 50 == 0){
   print vmin "  " vmax   "  "(vmax-vmin)/20 > ".nbline"
   }
}' $date


echo "-- Choix conseille pour la palette de couleurs:"       
tail -1 .nbline
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


echo " -- Creating scale " $palette
if [ $color = "1" ]; then
   pscontour $date  ${pars} -C$palette -I ${box} -P -K  >> $date.ps
   psscale -C$palette ${psca} -L -O >> $date.ps
fi
if [ $color = "2" ]; then
   pscontour $date  ${pars} -C$palette -I ${box} -P -K  >> $date.ps
   psscale -C$palette ${psca} -L -O -K >> $date.ps
fi
if [ $color = "3" ]; then
   pscontour $date  ${pars} -C$palette  ${box} -P -K -W2/0 >> $date.ps
fi

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

\rm GMT_palette_Temp.cpt wind*.data box_wind

echo " -------------------------------------------------------------"
echo " Fichier cree = " $date".ps   & Affichage Ghostview"
echo " -------------------------------------------------------------"
ghostview  -geometry 800x800 $date.ps &


