#!/bin/bash

data1=$1 #Data inicial do experimento.
data2=$2 #Data Final do experimento.

Perio=${data1}${data2}
echo 
echo ">> Usando o GNUplot para gerar a figura com as series "
echo ">> para os diferentes dominios no periodo"
echo 
echo "Data inicial: "${data1}
echo "Data final:   "${data2}
echo
echo "Chamando o GNUplot >>>>"
echo

hh=`echo ${data1} |cut -c 9-10`
YYYYMMANL=`echo ${6} | cut -c 1-8`

#Usando o GNUplot para gerar as figuras das series do somatório dos valores de umidade Não física
gnuplot << EOF 
reset
set xtics format ""
set xdata time
set timefmt "%Y%m%d%H"
set format x "%d"

set terminal png size 1360,900 enhanced 
set output 'humNphy.series_${Perio}.png'

set multiplot layout 5, 2 title "Sum of negative humidity values in the    profile                                          Sum of supersaturate humidity values in the profile"
set tmargin 1
set bmargin 4

set object 2 rect from graph 0, graph 0 to graph 1, graph 1 behind
set object 2 rect fc rgb "#FFF9D9" fillstyle solid 1.0
set style fill pattern 2 bo 1
set grid lc rgb "black" 

set xrange ["${dat}":"${YYYYMMANL}${hh}"]
set xtics "${dat}",259200,"${YYYYMMANL}${hh}"

# Definicao dos estilos padrão dos protocolos
set style line 1 lc 1 lt 9 lw 3 pi -1 ps 1.0  # FisrtGuess Experimento
set style line 2 lc 2 lt 17 lw 3 pi -1 ps 1.0 # Analise Experimento
set style line 3 lc -1 lt -1 lw 2.5           # FisrtGuess Controle
set style line 4 lc 5  lt -1 lw 2             # Analise Controle
set style line 5 lc 14 lt 6 lw 3 pi -1 ps 1.0 # Diferenca FisrtGuess - Analise
set style line 6 lc 4 lt 12 lw 3 pi -1 ps 1.0 # Diferenca Experimento - controle

set style data linespoints
set key reverse bottom outside horizontal
#set key reverse Left outside
set ylabel "Artic Region " 
set lmargin 13
set rmargin 5
set xzeroaxis linetype 3 linewidth 3 lc rgb "black"
#set yrange [-0 : -20000]
plot 'URnegV_BAM.fct_${Perio}.txt' using 1:6  ls 1 title 'Guess',\
     'URnegV_BAM.anl_${Perio}.txt' using 1:6  ls 2 title 'Anali',\
     'CNT_URnegV_BAM.fct_${Perio}.txt' using 1:6  ls 3  title 'Guess CNT',\
     'CNT_URnegV_BAM.anl_${Perio}.txt' using 1:6  ls 4  title 'Anali CNT'

set style data linespoints
set lmargin 13
set rmargin 5
#set yrange [-0 : -20000]
plot 'URsatV_BAM.fct_${Perio}.txt' using 1:6  ls 1 title 'Guess',\
     'URsatV_BAM.anl_${Perio}.txt' using 1:6  ls 2 title 'Anali',\
     'CNT_URsatV_BAM.fct_${Perio}.txt' using 1:6  ls 3 title 'Guess CNT',\
     'CNT_URsatV_BAM.anl_${Perio}.txt' using 1:6  ls 4  title 'Anali CNT'

set style data linespoints
set key off
#set key reverse Left outside
set ylabel "Hemisphere North " 
set lmargin 13
set rmargin 5
#set yrange [-0 : -20000]
plot 'URnegV_BAM.fct_${Perio}.txt' using 1:5  ls 1 title 'Guess',\
     'URnegV_BAM.anl_${Perio}.txt' using 1:5  ls 2  title 'Anali',\
     'CNT_URnegV_BAM.fct_${Perio}.txt' using 1:5  ls 3 title 'Guess CNT',\
     'CNT_URnegV_BAM.anl_${Perio}.txt' using 1:5  ls 4  title 'Anali CNT'

set style data linespoints
set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
#set yrange [-0 : -20000]
plot 'URsatV_BAM.fct_${Perio}.txt' using 1:5  ls 1 title 'Guess',\
     'URsatV_BAM.anl_${Perio}.txt' using 1:5  ls 2  title 'Anali',\
     'CNT_URsatV_BAM.fct_${Perio}.txt' using 1:5  ls 3 title 'Guess CNT',\
     'CNT_URsatV_BAM.anl_${Perio}.txt' using 1:5  ls 4  title 'Anali CNT'

set ylabel "Tropical Region"
set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
set tics out
set ytics nomirror
#set yrange [-0 : -20000]
plot 'URnegV_BAM.fct_${Perio}.txt' using 1:4  ls 1 title 'Guess',\
     'URnegV_BAM.anl_${Perio}.txt' using 1:4  ls 2 title 'Anali',\
     'CNT_URnegV_BAM.fct_${Perio}.txt' using 1:4  ls 3 title 'Guess CNT',\
     'CNT_URnegV_BAM.anl_${Perio}.txt' using 1:4  ls 4 title 'Anali CNT'

set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
set tics out
set ytics nomirror
#set yrange [-0 : -20000]
set xzeroaxis linetype 3 linewidth 3 lc rgb "black"
plot 'URsatV_BAM.fct_${Perio}.txt' using 1:4  ls 1 title 'Guess',\
     'URsatV_BAM.anl_${Perio}.txt' using 1:4  ls 2 title 'Anali',\
     'CNT_URsatV_BAM.fct_${Perio}.txt' using 1:4  ls 3 title 'Guess CNT',\
     'CNT_URsatV_BAM.anl_${Perio}.txt' using 1:4  ls 4 title 'Anali CNT'

set style data linespoints
set key off
#set key reverse Left outside
set ylabel "Hemisphere South" 
set lmargin 13
set rmargin 5
#set yrange [-0 : -20000]
plot 'URnegV_BAM.fct_${Perio}.txt' using 1:3  ls 1 title 'Guess',\
     'URnegV_BAM.anl_${Perio}.txt' using 1:3  ls 2  title 'Anali',\
     'CNT_URnegV_BAM.fct_${Perio}.txt' using 1:3  ls 3 title 'Guess CNT',\
     'CNT_URnegV_BAM.anl_${Perio}.txt' using 1:3  ls 4  title 'Anali CNT'


set style data linespoints
set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
#set yrange [-0 : -20000]
plot 'URsatV_BAM.fct_${Perio}.txt' using 1:3  ls 1 title 'Guess',\
     'URsatV_BAM.anl_${Perio}.txt' using 1:3  ls 2  title 'Anali',\
     'CNT_URsatV_BAM.fct_${Perio}.txt' using 1:3  ls 3 title 'Guess CNT',\
     'CNT_URsatV_BAM.anl_${Perio}.txt' using 1:3  ls 4  title 'Anali CNT'


set ylabel "Antartic Region "
set xlabel "Tempo (Dia)"
set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
set tics out
set ytics nomirror
#set yrange [-0 : -20000]
plot 'URnegV_BAM.fct_${Perio}.txt' using 1:2  ls 1 title 'Guess',\
     'URnegV_BAM.anl_${Perio}.txt' using 1:2  ls 2 title 'Anali',\
     'CNT_URnegV_BAM.fct_${Perio}.txt' using 1:2  ls 3 title 'Guess CNT',\
     'CNT_URnegV_BAM.anl_${Perio}.txt' using 1:2  ls 4 title 'Anali CNT'

set key off
#set key reverse Left outside
set lmargin 13
set rmargin 5
set tics out
set ytics nomirror
#set yrange [-0 : -20000]
plot 'URsatV_BAM.fct_${Perio}.txt' using 1:2  ls 1 title 'Guess',\
     'URsatV_BAM.anl_${Perio}.txt' using 1:2  ls 2 title 'Anali',\
     'CNT_URsatV_BAM.fct_${Perio}.txt' using 1:2  ls 3 title 'Guess CNT',\
     'CNT_URsatV_BAM.anl_${Perio}.txt' using 1:2  ls 4 title 'Anali CNT'

exit
EOF

# Verificando se o arquivo foi criado
if [ -e humNphy.series_${Perio}.png  ]; then 
    echo "Arquivo gerado: "humNphy.series_${Perio}.png
    echo 
  else 
    echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
    echo "! Problemas na geracao do arquivo no GNUplot     !"
    echo "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"
    echo 
  fi