#!/bin/bash

# The job name
#SBATCH --job-name=mars_new_job
# Set the error and output files
#SBATCH --output=mars_new-%J.out
#SBATCH --error=mars_new-%J.out
# Set the initial working directory
#SBATCH --workdir=/scratch/ms/es/sp2b
# Choose the queue
#SBATCH --qos=long
# Wall clock time limit
##SBATCH --time=00:05:00
# Send an email on failure
##SBATCH -mail-type=FAIL

# Peticiones a MARS siendo un fichero por alcance.

# Fechas de los datos.
anno="2017"
mes="01"
dia="01"
pasada="00"

###############################################################################
###                                   OTROS                                 ###
###############################################################################

# TP en todos los perturbed.
############################

# Intervalos 3-horarios de 3 a 144.
for (( t=3; t<=144; t+=3 )) ; do
  if test ${t} -lt 10 ; then
    alcance="00${t}"
  elif test ${t} -lt 100 ; then
    alcance="0${t}"
  else
    alcance="${t}"
  fi
  mars>fcout_pf_ll_${anno}${mes}${dia}${pasada}_${alcance}<<EOF
retrieve,
  class=od,
  date=${anno}${mes}${dia},
  expver=1,
  levtype=sfc,
  number=1/2/3/4/5/6/7/8/9/10/11/12/13/14/15/16/17/18/19/20/21/22/23/24/25/26/27/28/29/30/31/32/33/34/35/36/37/38/39/40/41/42/43/44/45/46/47/48/49/50,
  param=228.128,
  step=${alcance},
  stream=enfo,
  time=${pasada},
  type=pf,
  grid=0.1/0.1,
  target="FC_PF_LL_${anno}${mes}${dia}${pasada}_${alcance}.grb"
EOF
done

# Intervalos 6-horarios de 150 a 360.
for (( t=150; t<=360; t+=6 )) ; do
  if test ${t} -lt 10 ; then
    alcance="00${t}"
  elif test ${t} -lt 100 ; then
    alcance="0${t}"
  else
    alcance="${t}"
  fi
  mars>fcout_pf_ll_${anno}${mes}${dia}${pasada}_${alcance}<<EOF
retrieve,
  class=od,
  date=${anno}${mes}${dia},
  expver=1,
  levtype=sfc,
  number=1/2/3/4/5/6/7/8/9/10/11/12/13/14/15/16/17/18/19/20/21/22/23/24/25/26/27/28/29/30/31/32/33/34/35/36/37/38/39/40/41/42/43/44/45/46/47/48/49/50,
  param=228.128,
  step=${alcance},
  stream=enfo,
  time=${pasada},
  type=pf,
  grid=0.1/0.1,
  target="FC_PF_LL_${anno}${mes}${dia}${pasada}_${alcance}.grb"
EOF
done

# Empirical distribution
########################

# Intervalos 6-horarios de 150 a 360.
for (( t=6; t<=360; t+=6 )) ; do
  if test ${t} -lt 10 ; then
    alcance="00${t}"
  elif test ${t} -lt 100 ; then
    alcance="0${t}"
  else
    alcance="${t}"
  fi
  mars>fcout_ed_ll_${anno}${mes}${dia}${pasada}_${alcance}<<EOF
retrieve,
  class=od,
  date=${anno}${mes}${dia},
  expver=1,
  levtype=sfc,
  number=51,
  param=164.128/167.128/207.128/228.128,
  step=${alcance},
  stream=enfo,
  time=${pasada},
  type=ed,
  grid=0.1/0.1,
  target="FC_ED_LL_${anno}${mes}${dia}${pasada}_${alcance}.grb"
EOF
done