#!/bin/ksh
####################################################################
#   Procedure cep_readsond.sh
#       This shell-script run the following commands on a workstation :
#         Retrieve observations from MARS, and unpack them.                     
#                                                                               
# modele de programme pour la recuperation des donnees d'OBSERVATION
#      avec BUFREX au lieu de GRIBEX
#     REGION:  A CHOISIR
#  Guy Cautenet  Anne-Marie Lanquette     Y.P.    7/95  7/96   
####################################################################
#  setting the parameters for LoadLeveler
# ---------------------------
#@ shell           = /usr/bin/ksh
#@ class = normal
#@ queue
# 
# 
#-------------------------------
# commands to be executed
#-------------------------------
set -x
PATH=$PATH:.                    # allows you to run any of your programs or
                                # scripts held in the current directory
export PATH
# modifier en fonction du nom {tpp} utilisateur
SCRATCH="/ws/scratch/ms/fr/auc/scratchdir"
cd $SCRATCH  
#
# Tidy up at the end of the script or if the scripts end abnormally.
trap 'eval cd $CURRENT_DIR;eval /bin/rm -r $SCRATCH/input_datafile ' 0 1 2 3 6 15
#
####################################################################
#                          MARS request
####################################################################
#    1 - Type de donnees  (parmi d'autres)
#
# type=obs      donnees d'observation
# type=fb       donnees des fichiers "feedback"
####################################################################
#    2 - Type d'observations  (parmi d'autres)
#
# obstype=t               donnees de radiosondage
# obstype=s,ident=03955   donnees de surface ?
# obstype=vsns,           donnees de vertical sounding non satellite
####################################################################
#    3 - Selection de la zone
#
#  area=60/-40.5/30/20.25, =latitude Nord/longitude Ouest/latitude 
#                                                         Sud/longitude Est
# 
#  area=60/-40.5/30/20.25
#  area=56./-15./35./25.   # aire Puy de Dome
#  area=50./-80./-60./30.  # aire NARE
#
#  (latitude negative vers le sud, positive vers le nord)
#  (longitude negative vers l'ouest, positive vers l'est)
####################################################################
#    4 - Date(s) et Heure(s)
#
# plusieurs possibilites au choix :
# date=19920901,time=00,
# date=19920901,time=00/to/18/by/6,
# date=19920901/to/19920906,time=00,
# date=19920901/to/19920907/by/3,time=00,
# date=19920901/to/19920906,time=00/to/18/by/6,
# date=19920901/to/19920907/by/3,time=00/to/18/by/6,
#################################################################
#  I - Extraction des donnees (code Mars)
# -------------------------------------------------------------
# DATE, HEURE,ZONE: le 26/12/1993 de 0h a 24h 
mars << *eor
ret,type=obs,obstype=t,repr=bufr,
date=20040101/to/20040131,time=0000/to/2400,target="Sond040101",
area=40/-45/-5/45
end
*eor
if [ $? != 0 ]
then
      echo " The MARS request failed."
      exit 1
fi
#################################################################
#  II - Donnees MARS dans input_datafile
# -------------------------------------------------------------
ln -sf Sond040101 input_datafile
################################################################ 
#  III - Execution du programme de decodage
#----------------------------------------------------------------------
$HOME/rdbufr.out
################################################################ 
#  IV - Sauvegarde des donnees dans le repertoire $SCRATCH
#---------------------------------------------------------------
rm Sond040101
compress fort.14
ls -al fort.14.Z
mv fort.14.Z Sond040101_don.Z
################################################################# 
#-------------------------------
# tidy up by deleting unwanted files
echo " Fin de la procedure cep_readsond.sh "
#-------------------------------
# This is done by the trap command.
 exit 0
#
# End of example script retpak  modifie le 2 juillet 1996
#################################################################