#!/bin/sh

# General Script for ECMWF data extraction ofr a given period

year=$1
sdate=$2
ndays=$3
index=$4
diri=$5
diro=$6

echo "Extracting data from the Summer $year ECMWF archive..."
echo "Starting date = ${sdate}"

# Running the meteo interface

echo ${sdate} ${ndays} > eurmet_${year}.p
echo "'"${diri}/u_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/v_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/t_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/q_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/surp_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/lowc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/medc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/higc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/temp_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/totc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/copc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diri}/lspc_bin.summer${year}"'" >> eurmet_${year}.p
echo "'"${diro}/ZW_eur"'" >> eurmet_${year}.p
echo "'"${diro}/MW_eur"'" >> eurmet_${year}.p
echo "'"${diro}/TE_eur"'" >> eurmet_${year}.p
echo "'"${diro}/SH_eur"'" >> eurmet_${year}.p
echo "'"${diro}/DE_eur"'" >> eurmet_${year}.p
echo "'"${diro}/MH_eur"'" >> eurmet_${year}.p
echo "'"${diro}/LH_eur"'" >> eurmet_${year}.p
echo "'"${diro}/RI_eur"'" >> eurmet_${year}.p
echo "'"${diro}/T2_eur"'" >> eurmet_${year}.p
echo "'"${diro}/CL_eur"'" >> eurmet_${year}.p

compile eurmet_${year}
eurmet_${year}.e < eurmet_${year}.p

#  Biogenic emissions

echo "'"${diro}/CL_eur"'" > calbio.p
echo "'"${diro}/T2_eur"'" >> calbio.p
echo "'"${diro}/BE_eur"'" >> calbio.p

compile calbio
calbio.e < calbio.p

#  Extraction of the model domain

compile exCONT

echo "5 '${diro}/DE_eur' '${diro}/DE.${index}'" | exCONT.e
echo "1 '${diro}/MH_eur' '${diro}/MH.${index}'" | exCONT.e
echo "5 '${diro}/MW_eur' '${diro}/MW.${index}'" | exCONT.e
echo "5 '${diro}/SH_eur' '${diro}/SH.${index}'" | exCONT.e
echo "5 '${diro}/TE_eur' '${diro}/TE.${index}'" | exCONT.e
echo "5 '${diro}/ZW_eur' '${diro}/ZW.${index}'" | exCONT.e
echo "5 '${diro}/LH_eur' '${diro}/LH.${index}'" | exCONT.e
echo "2 '${diro}/RI_eur' '${diro}/RI.${index}'" | exCONT.e
echo "1 '${diro}/T2_eur' '${diro}/T2.${index}'" | exCONT.e
echo "4 '${diro}/CL_eur' '${diro}/CL.${index}'" | exCONT.e
echo "3 '${diro}/BE_eur' '${diro}/BE.${index}'" | exCONT.e

# Housecleaning

rm ${diro}/ZW_eur
rm ${diro}/MW_eur
rm ${diro}/TE_eur
rm ${diro}/SH_eur
rm ${diro}/DE_eur
rm ${diro}/MH_eur
rm ${diro}/LH_eur
rm ${diro}/RI_eur
rm ${diro}/T2_eur
rm ${diro}/CL_eur
rm ${diro}/BE_eur
