#!/bin/sh

sd0=$1
nhours=$2
show=$3

year=`echo $sd0 | cut -c1-4`
month=`echo $sd0 | cut -c5-6`

case ${year} in
  1998) section=19980101-19981231; lev=10 ;;
  1999) section=19990101-19991231; lev=10 ;;
  2000) section=20000101-20001231; lev=18 ;;
  2001)   if [ ${month} -le 04 ] ; then section=20010101-20010430 ; lev=18 ;
        elif [ ${month} -le 08 ] ; then section=20010501-20010831 ; lev=18 ;
        else                            section=20010901-20011231 ; lev=18 ;
          fi ;;
  2002)   if [ ${month} -le 04 ] ; then section=20020101-20020430 ; lev=18 ;
        elif [ ${month} -le 08 ] ; then section=20020501-20020831 ; lev=18 ;
        else                            section=20020901-20021231 ; lev=18 ;
          fi ;;
     *) echo "*** DATE OF SIMULATION OUT OF ECMWF DATA BASE" ; exit ;;
esac

# Meteo interpolation

echo "Meteo Horizontal Interpolation..."
tcompile exdom-ecm
echo "Interpolating Zonal Wind..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/u_bin.${year}"' '"$TMP/u-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Meridional Wind..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/v_bin.${year}"' '"$TMP/v-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Temperature..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/t_bin.${year}"' '"$TMP/t-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Specific Humidity..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/q_bin.${year}"' '"$TMP/q-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Liquid water..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/l_bin.${year}"' '"$TMP/l-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Ice water..."
echo "$lev $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/i_bin.${year}"' '"$TMP/i-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Surface pressure..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/surp_bin.${year}"' '"$TMP/surp-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating 2m Temperature..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/temp_bin.${year}"' '"$TMP/temp-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Medium cloudiness..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/medc_bin.${year}"' '"$TMP/medc-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating High cloudiness..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/higc_bin.${year}"' '"$TMP/higc-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Surface sensible heat flux..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/sshf_bin.${year}"' '"$TMP/sshf-ecm"'" | $TMP/exdom-ecm.e
echo "Interpolating Latent heat flux..."
echo "1 $sd0 $nhours '"$FNCOO"' '"$METEODIR/${section}/slhf_bin.${year}"' '"$TMP/slhf-ecm"'" | $TMP/exdom-ecm.e

# Meteo processing

echo "Meteo ECMMET diagnostic interface..."
echo ${sd0} ${lev} ${show}  > $TMP/ecmmet.p
echo "'"$TMP/u-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/v-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/t-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/q-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/l-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/i-ecm"'"      >> $TMP/ecmmet.p
echo "'"$TMP/surp-ecm"'"   >> $TMP/ecmmet.p
echo "'"$TMP/medc-ecm"'"   >> $TMP/ecmmet.p
echo "'"$TMP/higc-ecm"'"   >> $TMP/ecmmet.p
echo "'"$TMP/temp-ecm"'"   >> $TMP/ecmmet.p
echo "'"$TMP/sshf-ecm"'"   >> $TMP/ecmmet.p
echo "'"$TMP/slhf-ecm"'"   >> $TMP/ecmmet.p
echo "'"$FNCOO"'"          >> $TMP/ecmmet.p
echo "'"$FNLDU"'"          >> $TMP/ecmmet.p
echo "'"$FNLAP"'"          >> $TMP/ecmmet.p
echo "'"$FNBIO"'"          >> $TMP/ecmmet.p
echo "'"$FNVCO"'"          >> $TMP/ecmmet.p
echo "'"$TMP/METEO"'"      >> $TMP/ecmmet.p
echo "'"$TMP/BEMISSIONS"'" >> $TMP/ecmmet.p
tcompile ecmmet
$TMP/ecmmet.e < $TMP/ecmmet.p
