#! /usr/bin/env python3

#from ReadDiagTools import gsiDiagConv as gs
import gsidiag as gd
import matplotlib.pyplot as plt
import os

file=gd.read_diag('/pesq/share/das/dist/joao.gerd/EXP18/GSI/dataout/2020022918/diag_conv_03.2020022918')

#Dados disponiveis espacial
file.ptmap('ps', [120])
plt.title('Distribuicao Espacial dos Dados Assimilados - 2020022918', fontsize='10')
plt.savefig('ptmap_ps_adpupa_2020022918.png')
file.ptmap('t', [120])
plt.title('Distribuicao Espacial dos Dados Assimilados - 2020022918', fontsize='10')
plt.savefig('ptmap_t_adpupa_2020022918.png')
file.ptmap('q', [120])
plt.title('Distribuicao Espacial dos Dados Assimilados - 2020022918', fontsize='10')
plt.savefig('ptmap_q_adpupa_2020022918.png')
file.ptmap('uv', [220, 221])
plt.title('Distribuicao Espacial dos Dados Assimilados - 2020022918', fontsize='10')
plt.savefig('ptmap_uv_adpupa_2020022918.png')


#file.ptmap('ps', [120, 132 ])
#os.system("mv ptmap_ps_2020022918.png ptmap_ps_adpupa_2020022918.png >& /dev/null")
#file.ptmap('t', [120, 132 ])
#file.ptmap('uv', [220, 221, 232])

#os.system("mv ptmap_uv_2020022918.png ptmap_uv_adpupa_2020022918.png >& /dev/null")

# Espacial usado nao usado por variavel
#file.pvmap(['t', 'ps', 'q','uv'], mask='kx==120')
#file.plot('ps', [120, 132 ], 'obs', mask='iuse==1',legend=True)