#!/bin/sh
####################################################################
#   Procédure cep_readsond.sh5
#  @(#)retobs     July 95 - ECMWF
#       This shell-script run the following commands on a workstation :
#         Retrieve observations from MARS, and unpack them.                     
#                                                                               
# modèle de programme pour la récuperation des données d'OBSERVATION
#      avec BUFREX au lieu de GRIBEX
#     RÉGION:  À CHOISIR
#  Guy Cautenet  Anne-Marie Lanquette     Y.P.    7/95  7/96   
#  Modele mis a jour le16 aout 2000 pour radoisondages   AML
####################################################################
#  setting the parameters for NQE
# ---------------------------
#QSUB -s /bin/sh         # Change to Bourne Shell
#QSUB -q ecgate          # Éxecution sur ecgate1 ou ecgate2
# 
#-------------------------------
# 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
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 données  (parmi d'autres)
#
# type=obs      données d'observation
# type=fb       données des fichiers "feedback"
####################################################################
#    2 - Type d'observations  (parmi d'autres)
#
# obstype=t               données de radiosondage
# obstype=s,ident=03955   données de surface ?
# obstype=vsns,           données de vertical sounding non satellite
####################################################################
#    3 - Sélection 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 négative vers le sud, positive vers le nord)
#  (longitude négative vers l'ouest, positive vers l'est)
####################################################################
#    4 - Date(s) et Heure(s)
#
# plusieurs possibilités 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 données (code Mars)
# -------------------------------------------------------------
# DATE, HEURE,ZONE: le 26/12/1993 de 0h à 24h 
mars << EOF
ret,type=obs,obstype=t,repr=bufr,
date=20010317/to/20010322,time=0000/to/2400,target="Sond010317",
area=45/3/42/7
EOF
if [ $? != 0 ]
then
      echo " The MARS request failed."
      exit 1
fi
#################################################################
#  II - Données MARS dans input_datafile
# -------------------------------------------------------------
ln -s Sond010317 input_datafile
################################################################ 
#  III - Exécution du programme de décodage
#----------------------------------------------------------------------
$HOME/rdbufr.out
################################################################ 
#  IV - Sauvegarde des données dans le répertoire $SCRATCH
#---------------------------------------------------------------
rm Sond010317
compress fort.14
ls -al fort.14.Z
mv fort.14.Z $SCRATCH/Sond010317_don.Z
################################################################# 
#-------------------------------
# tidy up by deleting unwanted files
echo " Fin de la procédure cep_readsond.sh5 - radiosondages"
#-------------------------------
# This is done by the trap command.
 exit 0
#
# End of example script cep_readsond.sh5 recopie dans pradmod
#################################################################