#!/bin/bash

#wget -nv ftp://ftp.hycom.org/datasets/GLBu0.08/expt_19.1/data/2011/hycom_GLBu0.08_191_2011111200_t003.nc
#ftps://ftp.hycom.org/datasets/GLBv0.08/expt_56.3
webdir=ftps://ftp.hycom.org/datasets/GLBu0.08/expt_56.3/
fileprefix=hycom_glb_563_

mn_d[1]=31
mn_d[2]=28
mn_d[3]=31
mn_d[4]=30
mn_d[5]=31
mn_d[6]=30
mn_d[7]=31
mn_d[8]=31
mn_d[9]=30
mn_d[10]=31
mn_d[11]=30
mn_d[12]=31

N_days=20
ini_yr=2016
ini_mn=1
ini_dd=1
ini_hh=00

for day_c in `seq 1 ${N_days}`
do
  yr=$ini_yr
  dd=$ini_dd
  mn=$ini_mn
  hh=$ini_hh
  if [ $ini_dd -lt 10 ]; then
    dd=0$ini_dd
  fi
  if [ $ini_mn -lt 10 ]; then
    mn=0$ini_mn
  fi

  for f_hr in t000 t003 t006 t009 t012 t015 t018 t021 
  do
    filename=${fileprefix}${yr}${mn}${dd}${hh}_${f_hr}.nc
    if [ ! -f $filename ]; then
      echo ${yr}${mn}${dd}${hh}_${f_hr}.nc downloading
      wget -nv ${webdir}${yr}/${filename}
    else 
      echo ${yr}${mn}${dd}${hh}_${f_hr}.nc exist
    fi
  done

  ini_dd=$((${ini_dd}+1))
  if [ $ini_dd -gt ${mn_d[$ini_mn]} ]; then
    ini_dd=1
    ini_mn=$((${ini_mn}+1))
    if [ $ini_mn -gt 12 ]; then
      ini_yr=$((${ini_yr}+1))
      ini_mn=1
    fi
  fi

done
#ncrcat hycom_GLBu0.08_191_20080* hycom_GLBu0.08_191_2008CS_event.nc hycom_GLBu0.08_191_2008JF.nc
echo done