#!/bin/ksh
####################################################################
#   Procedure cep_readsurf.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/98 12/98   
####################################################################
#  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
####################################################################
#    2 - Type d'observations  (parmi d'autres)
#
# obstype=s,    donnees de surface   
####################################################################
#    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
#  area=46./3./45./4.      # Clermont-Ferrand ou 
#  ident=07460             # CLERMONT-FERRAND
#
#  (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=19980101,time=00,
# date=19980101,time=00/to/18/by/6,
# date=19980101/to/19980106,time=00,
# date=19980101/to/19980107/by/3,time=00,
# date=19980101/to/19980106,time=00/to/18/by/6,
# date=19980101/to/19980107/by/3,time=00/to/18/by/6,
#################################################################
#  I - Extraction des donnees (code Mars)
# -------------------------------------------------------------
mars <<EOF
retrieve,date=20030802/to/20030816,time=00/to/24,target = "Surf030802",
type=obs,obstype=s,repr=bufr,
area=56./-15./35./25.      # Clermont-Ferrand
EOF
if [ $? != 0 ]
then
      echo " The MARS request failed."
      exit 1
fi
#################################################################
#  II - Donnees MARS dans input_datafile
# -------------------------------------------------------------
rm input_datafile
ln -sf  Surf030802 input_datafile
################################################################ 
#  III - Execution du programme de decodage
#----------------------------------------------------------------------
#  Differents parametres a demander
#
#PRESSURE          : Pressure
#MSLPRESSU         : Mean Sea Level Pressure
#3HPRESSURE        : 3 Hours Pressure change
#TENDPRESSU        : Pressure Tendency (characteristic)
#WINDDIR           : Wind Direction at 10 m
#WINDSPEED         : Wind Speed at 10 m
#DRYTEMPERA        : Dry bulb Temperature at 2 m
#DEWTEMPERA        : Dew point Temperature at 2 m
#RELHUMIDIT        : Relative Humidity
#HORIZVISIB        : Horizontal Visibility
#PRESWEATHE        : Present Weather
#PASTWEATH1        : Past Weather 1
#PASTWEATH2        : Past Weather 2
#CLOUDCOVER        : Cloud Cover
#VERT.SIGN.        : 
#CLOUDAMOUN        : Cloud Amount
#CLOUDBASE         : Cloud Base height
#CLOUDTYPEA        : Cloud Type 
#CLOUDTYPEB        : 
#CLOUDTYPEC        : 
#VERT.SIGN1        : 
#CLOUDAMOU1        : 
#CLOUDTYPE1        : 
#CLOUDBAS1         : 
#VERT.SIGN2        : 
#CLOUDAMOU2        : 
#CLOUDTYPE2        : 
#CLOUDBAS2         : 
#VERT.SIGN3        : 
#CLOUDAMOU3        : 
#CLOUDTYPE3        : 
#CLOUDBAS3         : 
#VERT.SIGN4        : 
#CLOUDAMOU4        : 
#CLOUDTYPE4        : 
#CLOUDBAS4         : 
#TOTPRECIP.        : Total Precipitation past 24 h
#TOTSNOW           : Total snow depth
#QUL.INFO.         : 
#DATAPRESE1        : 
#DATAPRESE2        : 
#----------------------------------------------------------------------
$HOME/rdbufrsrf.out << EOF
PRESSURE          : Pressure
WINDDIR           : Wind Direction at 10 m
WINDSPEED         : Wind Speed at 10 m
DRYTEMPERA        : Dry bulb Temperature at 2 m
DEWTEMPERA        : Dew point Temperature at 2 m
RELHUMIDIT        : Relative Humidity
TOTPRECIP.        : Total Precipitation past 24 h
FIN               : Fin du choix des parametres
EOF
################################################################ 
#  IV - Sauvegarde des donnees dans le repertoire $SCRATCH
#---------------------------------------------------------------
rm Surf030802
compress fort.12
ls -al fort.12.Z
mv fort.12.Z Surf030802_don.Z
################################################################# 
#-------------------------------
# tidy up by deleting unwanted files
echo " Fin de la procedure cep_readsurf.sh "
#-------------------------------
# This is done by the trap command.
 exit 0
#
#################################################################