cd "$(dirname "$0")"
curl "http://grbsrv.opengribs.org/getmygribs2.php?osys=osx&ver=1.2.6&model=gfs_p25_&la1=0&la2=35&lo1=-120&lo2=-50&intv=3&days=7&cyc=last&par=W;P;R;C;T;H;c;G;r;8;&wmdl=none&wpar=" -A "XyGrib_mac/1.2.6" > gfs.json
jq .message.url gfs.json > gfs.url
cat gfs.url 
sed 's/"//g' gfs.url > gfs2.url
curl $(cat gfs2.url) -A "XyGrib_mac/1.2.6" > gfs.grib2

##distintas variables
echo 'ref_time'
ref_time=$(gdalinfo gfs.grib2 | grep -i grib_ref_time | sort | uniq)
echo $ref_time
ref_time=$(echo $ref_time | cut -d' ' -f 2)
echo $ref_time,$(date -d @$ref_time)
#echo 'valid_time'
valid_time=$(gdalinfo gfs.grib2 | grep -i grib_valid_time | sort | uniq)
echo $valid_time
valid_time=$(echo $valid_time | cut -d' ' -f 2)
echo $valid_time,$(date -d @$valid_time)
echo "Variables"
echo $(gdalinfo gfs.grib2 | grep -i grib_comment | sort | uniq)
##extraer las capas

##lluvia 
j=0
tifs=""
for i in {14..70}
do
    if [ $i != 15 ]; then
        echo "Banda: $i pcpn_$j"
        gdal_translate gfs.grib2 -B $i -a_srs EPSG:4326 -r bilinear -tr 0.05 0.05 /var/www/nube1.on.gt/gfs/pcpn_$j.tif
        newfile="/var/www/nube1.on.gt/gfs/pcpn_$j.tif"
        tifs+=" ${newfile}"
        j=$((j+3))
    fi    
done

# agrupaciones
bands="ABCDEFGH"
j=1
tifs=""
    for j in {0..4}
    do
    tifs=""
    for i in {0..7}
    do  
        tif=$((($i+1)*3))
        tif2=$(($tif+$j*24))
        newfile="-${bands:i:1} /var/www/nube1.on.gt/gfs/pcpn_$tif2.tif"
        tifs+=" ${newfile}"
    done 
    gdal_calc.py $tifs --outfile=/var/www/nube1.on.gt/gfs/a_pcpn_$tif2.tif --calc='A+B+C+D+E+F+G+H'
done
# j=0
# for i in {0..55}
# do
#         echo "Banda: $i pcpn_$j"
#         gdal_translate gfs.grib2 -B $i -a_srs EPSG:4326 -r bilinear -tr 0.05 0.05 /var/www/nube1.on.gt/gfs/pcpn_$j.tif
#         newfile="/var/www/nube1.on.gt/gfs/pcpn_$j.tif"
#         tifs+=" ${newfile}"
#         j=$((j+3))
# done

#echo " tifs: $tifs"
### demasiado grande 450mb gdalbuildvrt -separate -o pcpn.vrt /var/www/nube1.on.gt/gfs/pcpn_*.tif
### demasiado grande gdal_translate pcpn.vrt /var/www/nube1.on.gt/gfs/pcpn.tif
#gdal_translate gfs.grib2 -B 24 -a_srs EPSG:4326 -r bilinear -tr 0.05 0.05 /var/www/nube1.on.gt/gfs/b24.tif
#gdal_translate gfs.grib2 -B 166 -a_srs EPSG:4326 -r bilinear -tr 0.05 0.05 /var/www/nube1.on.gt/gfs/temp_12z.tif