#!/bin/bash
#set -x
################################################################################################################################
#Script para contar a quantidade de observacoes GNSS assimiladas e rejeitadas
#
#
#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:
#  ./assim_sat.ksh 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:
#  10-03-2015 - Amarante, Lucas    - Codigo Inicial do Count_GNSS_MetopB_Assimilados.ksh
#  29-06-2017 - Banos, Ivette      - Criacao do script baseado no Count_GNSS_MetopB_Assimilados.ksh
################################################################################################################################

echo "Script para contar a quantidade de observacoes GNSS assimiladas e rejeitadas"

# 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}

data=${datai}
        
         dirOUT=${obsGnss_tabe}


       gnss_anl="EXP_Assimil_anl_${data}${dataf}.txt"
   rej_gnss_anl="EXP_Reject_anl_${data}${dataf}.txt"
#CNT_

   echo "Arquivos a gerar:"
   echo ${gnss_anl}
   echo ${rej_gnss_anl}
   echo

 echo "Ano    Mes    Dia    Hora  MetopB  MetopA  TerraSAR  GraceA  FM1  FM2  FM4  FM5  FM6  C/NOFS  TotalAssimilado" > ${dirOUT}/${gnss_anl}
 echo "Ano    Mes    Dia    Hora  MetopB  MetopA  TerraSAR  GraceA  FM1  FM2  FM4  FM5  FM6  C/NOFS  TotalReject"     > ${dirOUT}/${rej_gnss_anl}

while [ ${data} -le ${dataf} ]
do

 
  echo "Executando para a data: " ${data}
     
     ano=`echo ${data} | cut -c1-4`
     mes=`echo ${data} | cut -c5-6`
     dia=`echo ${data} | cut -c7-8`
    hora=`echo ${data} | cut -c9-10`

    dirIN=${pathDatain}/${ano}${mes}/${dia}${hora}/conv

########################################################################## 
#
#Separando os tipos de satelite GNSS
#Arquivo anl
#assimilados

 awk '{ if($2 == "'1.000000'") print  }' ${dirIN}/sso_gps.${data}.anl.dat > ${dirOUT}/use1_anl.dat
 awk '{ if($2 == "'-1.000000'") print  }' ${dirIN}/sso_gps.${data}.anl.dat > ${dirOUT}/usem1_anl.dat

     total=`grep "gps" ${dirIN}/soc.${data}.anl.dat | awk '{print $3}'`
 
    metopb=`awk '{ if($1 == 3.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
    metopa=`awk '{ if($1 == 4.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
  terrasar=`awk '{ if($1 == 42.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
    gracea=`awk '{ if($1 == 722.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
       fm1=`awk '{ if($1 == 740.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
       fm2=`awk '{ if($1 == 741.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
       fm4=`awk '{ if($1 == 743.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
       fm5=`awk '{ if($1 == 744.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
       fm6=`awk '{ if($1 == 745.000) print }' ${dirOUT}/use1_anl.dat | wc -l`
     cnofs=`awk '{ if($1 == 786.000) print }' ${dirOUT}/use1_anl.dat | wc -l`

   echo ${ano}'     '${mes}'     '${dia}'     '${hora}'    '${metopb}'     '${metopa}'     '${terrasar}'     '${gracea}'     '${fm1}'     '${fm2}'     '${fm4}'     '${fm5}'     '${fm6}'     '${cnofs}'     '${total} >> ${dirOUT}/${gnss_anl}

##############################
#Arquivo anl
#rejeitados

    reject=`grep "gps" ${dirIN}/soc.${data}.anl.dat | awk '{print $5}'` 

    metopb=`awk '{ if($1 == 3.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
    metopa=`awk '{ if($1 == 4.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
  terrasar=`awk '{ if($1 == 42.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
    gracea=`awk '{ if($1 == 722.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
       fm1=`awk '{ if($1 == 740.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
       fm2=`awk '{ if($1 == 741.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
       fm4=`awk '{ if($1 == 743.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
       fm5=`awk '{ if($1 == 744.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
       fm6=`awk '{ if($1 == 745.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`
     cnofs=`awk '{ if($1 == 786.000) print }' ${dirOUT}/usem1_anl.dat | wc -l`

   echo ${ano}'     '${mes}'     '${dia}'     '${hora}'     '${metopb}'     '${metopa}'     '${terrasar}'     '${gracea}'     '${fm1}'     '${fm2}'     '${fm4}'     '${fm5}'     '${fm6}'     '${cnofs}'     '${reject} >> ${dirOUT}/${rej_gnss_anl}

############################################################################################################################################################################################

   data=$(${inctime} ${data} +6hr %y4%m2%d2%h2)

done

cd ${dirOUT}
rm use*.dat

exit