data_ini=$1
data_end=$2

dir_out=/pesq/share/oceanwave/dist/Valdir/era5_waves
data=$data_ini

while [[ $data -le $data_end ]] ; do

echo ${data}
year=${data:0:4}
month=${data:4:2}

\rm -f wave.py

cat > wave.py << EOF

import cdsapi

dataset = "reanalysis-era5-single-levels"
request = {
    "product_type": ["reanalysis"],
    "variable": [
        "mean_wave_direction",
        "mean_wave_period",
        "significant_height_of_combined_wind_waves_and_swell",
        "peak_wave_period"
    ],
    "year": ["${year}"],
    "month": ["${month}"],
    "day": [
        "01", "02", "03",
        "04", "05", "06",
        "07", "08", "09",
        "10", "11", "12",
        "13", "14", "15",
        "16", "17", "18",
        "19", "20", "21",
        "22", "23", "24",
        "25", "26", "27",
        "28", "29", "30",
        "31"
    ],
    "time": [
        "00:00", "01:00", "02:00",
        "03:00", "04:00", "05:00",
        "06:00", "07:00", "08:00",
        "09:00", "10:00", "11:00",
        "12:00", "13:00", "14:00",
        "15:00", "16:00", "17:00",
        "18:00", "19:00", "20:00",
        "21:00", "22:00", "23:00"
    ],
    "data_format": "netcdf",
    "download_format": "unarchived",
    "area": [20, -80, -8, -30]
}

client = cdsapi.Client()
client.retrieve(dataset, request).download(
        "${dir_out}/era5_wave_hourly_${month}${year}.nc"
)
EOF
chmod u+wx wave.py
python wave.py
\rm -f wave.py

data=$(date -u "+%Y%m" -d "${year}${month}01 00:00 1 month")

done