#!/bin/bash
#set -x
################################################################################################################################
#Script para realizar a média e stddev de OmA e OmF no perfil veritcal para cada satélite usado
#
#
#gps      003    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  METOP-B
#gps      004    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  METOP-A
#gps      042    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  TerraSAR
#gps      722    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  GRACE-A
#gps      740    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  COSMIC
#gps      741    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  COSMIC
#gps      743    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  COSMIC
#gps      744    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  COSMIC
#gps      745    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.000000     0    0.     0.      0  COSMIC
#gps      786    0    1     3.0      0      0      0  10.0  10.0   1.0  10.0  0.00000      0    0.     0.      0  C/NOFS
#
#
# !INTERFACE:
#  ./fit_vert.sh datai dataf pathDatain
#
#  datai = Data inicial YYYYMMDDHH 
#  dataf = Data final   YYYYMMDDHH  
#  pathDatain = Enderećo dos dados de saida do STAT rodado do GSI 
#
# !REVISION HISTORY:
#  04-07-2017 - Banos, Ivette      - Criacao do script para a criacao de tabelas com todos os dados de todos os horarios
#                                     e realizacao da media para cada nivel de pressao
################################################################################################################################
# Exportando variaveis do eval
source  ../../../config_eval.sh vars_export ${nome_sys} ${EXPE}

# Exportando variaveis do obsGnss
source ${obsGnss_home}/config_obsGnss.ksh vars_export ${nome_sys} ${EXPE}

datai=${1}
dataf=${2}
pathDatain=${3}

#dirIN=/scratchin/grupos/assim_dados/home/eduardo.khamis/eval/Ex02ALL/stat/binary/temporal/conv
dirIN=${pathDatain}/temporal/conv

dirOUT=${obsGnss_tabe}
#dirOUT=${obsGnss_CNT}

   echo "Dir input:"
   echo ${dirIN}
   echo "Dir output:"   
   echo ${dirOUT} 

echo ""
echo "Lista de arquivos gerados:" 

cd ${dirIN}

#Criando tabelas anl
#METOPB
 awk '{  print  }' *anl*03* > ${dirOUT}/metb_anl.dat
#METOPA
 awk '{  print  }' *anl*0004* > ${dirOUT}/meta_anl.dat
#TERRASAR-X
 awk '{  print  }' *anl*042* > ${dirOUT}/terra_anl.dat
#GRACEA
 awk '{  print  }' *anl*722* > ${dirOUT}/gracea_anl.dat
#COSMIC
 awk '{  print  }' *anl*074*  > ${dirOUT}/cosm_anl.dat
#CNOFS
 awk '{  print  }' *anl*0786* > ${dirOUT}/cnof_anl.dat

#Criando tabelas fgs
#METOPB
 awk '{  print  }' *fgs*03* > ${dirOUT}/metb_fgs.dat
#METOPA
 awk '{  print  }' *fgs*0004* > ${dirOUT}/meta_fgs.dat
#TERRASAR-X
 awk '{  print  }' *fgs*042* > ${dirOUT}/terra_fgs.dat
#GRACEA
 awk '{  print  }' *fgs*722* > ${dirOUT}/gracea_fgs.dat
#COSMIC
 awk '{  print  }' *fgs*074*  > ${dirOUT}/cosm_fgs.dat
#CNOFS
 awk '{  print  }' *fgs*0786* > ${dirOUT}/cnof_fgs.dat

cd ${dirOUT}

#Calculando as medias e desvio padrão de cada nivel por cada satelite

#TERRASAR-X
terr=("terra_anl.dat" "terra_fgs.dat" )
for i in $(seq 0 1)
do
arq=${terr[i]}
echo " Satelites=   ${arq}"
echo "NaN                NaN                " > ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
/scratchin/grupos/assim_dados/home/gdad/public/gnuplot5/bin/gnuplot << EOF
set print "${dirOUT}/stat_${arq}" append
do for [i=4:15] {
stats  '${arq}' u i nooutput ; 
print STATS_mean, STATS_stddev
} 
exit
EOF
done

#GRACE-A
grac=("gracea_anl.dat" "gracea_fgs.dat" )
for i in $(seq 0 1)
do
arq=${grac[i]}
echo " Satelites=   ${arq}"
echo "NaN                NaN                " > ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
/scratchin/grupos/assim_dados/home/gdad/public/gnuplot5/bin/gnuplot << EOF
set print "${dirOUT}/stat_${arq}" append
do for [i=5:15] {
stats  '${arq}' u i nooutput ; 
print STATS_mean, STATS_stddev
} 
exit
EOF
done

#METOPB/A
met=("metb_anl.dat" "metb_fgs.dat" "meta_anl.dat" "meta_fgs.dat" )
for i in $(seq 0 3)
do
arq=${met[i]}
echo " Satelites=   ${arq}"
echo "NaN                NaN                " > ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
/scratchin/grupos/assim_dados/home/gdad/public/gnuplot5/bin/gnuplot << EOF 
set print "${dirOUT}/stat_${arq}" append
do for [i=9:15] {
stats  '${arq}' u i nooutput ; 
print STATS_mean, STATS_stddev
}
exit
EOF
done

#COSMIC
cos=("cosm_anl.dat" "cosm_fgs.dat" )
for i in $(seq 0 1)
do
arq=${cos[i]}
echo " Satelites=   ${arq}"
echo "NaN                NaN                " > ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
/scratchin/grupos/assim_dados/home/gdad/public/gnuplot5/bin/gnuplot << EOF 
set print "${dirOUT}/stat_${arq}" append
do for [i=7:15] {
stats  '${arq}' u i nooutput ; 
print STATS_mean, STATS_stddev
}
exit
EOF
done

#CNOFS
cnf=("cnof_anl.dat" "cnof_fgs.dat" )
for i in $(seq 0 1)
do
arq=${cnf[i]}
echo " Satelites=   ${arq}"
echo "NaN                NaN                " > ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
echo "NaN                NaN                " >> ${dirOUT}/stat_${arq}
/scratchin/grupos/assim_dados/home/gdad/public/gnuplot5/bin/gnuplot << EOF
set print "${dirOUT}/stat_${arq}" append
do for [i=8:15] {
stats  '${arq}' u i nooutput ; 
print STATS_mean, STATS_stddev
}
exit
EOF
done

#Gerando arquivo com os niveis
vetor=("1000" "900" "800" "700" "600" "500" "400" "300" "250" "200" "150" "100" "50" "0" )
for j in $(seq 0 13)
do
var=${vetor[j]}
echo ${var} >> ${dirOUT}/niveis.dat
done

medias_stdev_anl="EXP_medias_stdev_anl_${datai}${dataf}.txt"
medias_stdev_fgs="EXP_medias_stdev_fgs_${datai}${dataf}.txt"
#CNT_

#Colocando todo em uma mesma tabela
echo "Niv    	metbM                metbS                metaM                metaS                terrM                terrS                graceaM                graceaS                cosmM                cosmS                cnofM                cnofS" > ${dirOUT}/${medias_stdev_anl}
echo "Niv       metbM                metbS                metaM                metaS                terrM                terrS                graceaM                graceaS                cosmM                cosmS                cnofM                cnofS" > ${dirOUT}/${medias_stdev_fgs}

pr -m -t -J niveis.dat stat_metb_anl.dat stat_meta_anl.dat stat_terra_anl.dat stat_gracea_anl.dat stat_cosm_anl.dat stat_cnof_anl.dat >> ${dirOUT}/${medias_stdev_anl}

pr -m -t -J niveis.dat stat_metb_fgs.dat stat_meta_fgs.dat stat_terra_fgs.dat stat_gracea_fgs.dat stat_cosm_fgs.dat stat_cnof_fgs.dat >> ${dirOUT}/${medias_stdev_fgs}

rm cnof_* cosm_* me* niveis.dat terra_* gracea* stat*

exit