#! /usr/bin/env python3

#from ReadDiagTools import gsiDiagConv as gs
import gsiDiagEval2 as gs
import os

a=gs.open('/stornext/online6/das/gdad/SMG_eval/EX1ALL/gsi/dataout/2013010818/diag/pe%e.conv_02')

a.pgeomap('gps',3,'pbqc',2013010818,0,1000)
a.pgeomap('gps',4,'pbqc',2013010818,0,1000)
a.pgeomap('gps',740,'pbqc',2013010818,0,1000)
a.pgeomap('gps',745,'pbqc',2013010818,0,1000)
a.pgeomap('gps',744,'pbqc',2013010818,0,1000)
a.pgeomap('gps',741,'pbqc',2013010818,0,1000)
a.pgeomap('gps',722,'pbqc',2013010818,0,1000)

a.ptmap('gps', [3, 4, 740, 745, 744, 741, 722], 2013010818)
os.system("mv ptmap2013010818.png ptmap_ref_ROGNSS_2013010818.png")


dini = 2013010100 
dfin = 2013010818
vName = 'gps'
vType = 3

gsiTools = gs.gsiFunctions()
gsiTools.time_series(varName=vName, varType=vType, dateIni=dini, dateFin=dfin)