#sintaxe: variab(serie,lag(hr),lim,flag) print 'Realizando teste de Variabilidade' flagp[:,1] = consiste_proc.variab(dados[:,1],1,10,flagp[:,1]) #ws flagp[:,2] = consiste_proc.variab(dados[:,2],1,10,flagp[:,2]) #wg flagp[:,3] = consiste_proc.variab(dados[:,3],1,360,flagp[:,3]) #wd flagp[:,4] = consiste_proc.variab(dados[:,4],1,5,flagp[:,4]) #at flagp[:,5] = consiste_proc.variab(dados[:,5],1,10,flagp[:,5]) #pr flagp[:,6] = consiste_proc.variab(dados[:,6],1,5,flagp[:,6]) #wt flagp[:,7] = consiste_proc.variab(dados[:,7],1,3,flagp[:,7]) #hs flagp[:,8] = consiste_proc.variab(dados[:,8],1,15,flagp[:,8]) #tp flagp[:,9] = consiste_proc.variab(dados[:,9],1,360,flagp[:,9]) #dp #Teste 3 - Consec. iguais #sintaxe = iguais(var,nci,flag) print 'Realizando teste de Consec. Iguais' flagp[:,1] = consiste_proc.iguais(dados[:,1],3,flagp[:,1]) #ws flagp[:,2] = consiste_proc.iguais(dados[:,2],3,flagp[:,2]) #wg flagp[:,3] = consiste_proc.iguais(dados[:,3],3,flagp[:,3]) #wd flagp[:,4] = consiste_proc.iguais(dados[:,4],3,flagp[:,4]) #at flagp[:,5] = consiste_proc.iguais(dados[:,5],3,flagp[:,5]) #pr flagp[:,6] = consiste_proc.iguais(dados[:,6],3,flagp[:,6]) #wt flagp[:,7] = consiste_proc.iguais(dados[:,7],3,flagp[:,7]) #hs flagp[:,8] = consiste_proc.iguais(dados[:,8],3,flagp[:,8]) #tp flagp[:,9] = consiste_proc.iguais(dados[:,9],3,flagp[:,9]) #dp #Teste 4 - Media e Desvio Padrao #sintaxe = iguais(var,per_meddp,mult_dp,flag) print 'Realizando teste de Media e Desvio Padrao' flagp[:,1] = consiste_proc.meddp(dados[:,1],100,3,flagp[:,1]) #ws flagp[:,2] = consiste_proc.meddp(dados[:,2],100,3,flagp[:,2]) #wg flagp[:,3] = consiste_proc.meddp(dados[:,3],100,3,flagp[:,3]) #wd
flagp[:, 5] = consiste_proc.variab(dados[:, 4], 1, 5, flagp[:, 5]) # rh # flagp[:,6] = consiste_proc.variab(dados[:,5],1,5,flagp[:,6]) #dw flagp[:, 7] = consiste_proc.variab(dados[:, 6], 1, 5, flagp[:, 7]) # pr flagp[:, 8] = consiste_proc.variab(dados[:, 7], 1, 5, flagp[:, 8]) # st flagp[:, 9] = consiste_proc.variab(dados[:, 8], 1, 5, flagp[:, 9]) # bh # flagp[:,10] = consiste_proc.variab(dados[:,9],???,flagp[:,10]) #cl # flagp[:,11] = consiste_proc.variab(dados[:,10],???,flagp[:,11]) #tu # flagp[:,12] = consiste_proc.variab(dados[:,11],1,5,flagp[:,12]) #sr flagp[:, 13] = consiste_proc.variab(dados[:, 12], 1, 5, flagp[:, 13]) # hs # flagp[:,14] = consiste_proc.variab(dados[:,13],1,5,flagp[:,14]) #hm flagp[:, 15] = consiste_proc.variab(dados[:, 14], 1, 5, flagp[:, 15]) # tp flagp[:, 16] = consiste_proc.variab(dados[:, 15], 1, 5, flagp[:, 16]) # dp # flagp[:,17] = consiste_proc.variab(dados[:,16],1,5,flagp[:,17]) #sp # #Teste 1 - Consec. iguais flagp[:, 1] = consiste_proc.iguais(dados[:, 0], 5, flagp[:, 1]) # ws flagp[:, 2] = consiste_proc.iguais(dados[:, 1], 5, flagp[:, 2]) # wg flagp[:, 3] = consiste_proc.iguais(dados[:, 2], 5, flagp[:, 3]) # wd flagp[:, 4] = consiste_proc.iguais(dados[:, 3], 5, flagp[:, 4]) # at flagp[:, 5] = consiste_proc.iguais(dados[:, 4], 5, flagp[:, 5]) # rh # flagp[:,6] = consiste_proc.iguais(dados[:,5],5,flagp[:,6]) #dw flagp[:, 7] = consiste_proc.iguais(dados[:, 6], 5, flagp[:, 7]) # pr flagp[:, 8] = consiste_proc.iguais(dados[:, 7], 5, flagp[:, 8]) # st flagp[:, 9] = consiste_proc.iguais(dados[:, 8], 5, flagp[:, 9]) # bh # flagp[:,10] = consiste_proc.iguais(dados[:,9],???,flagp[:,10]) #cl # flagp[:,11] = consiste_proc.iguais(dados[:,10],???,flagp[:,11]) #tu # flagp[:,12] = consiste_proc.iguais(dados[:,11],5,flagp[:,12]) #sr flagp[:, 13] = consiste_proc.iguais(dados[:, 12], 5, flagp[:, 13]) # hs # flagp[:,14] = consiste_proc.iguais(dados[:,13],5,flagp[:,14]) #hm flagp[:, 15] = consiste_proc.iguais(dados[:, 14], 5, flagp[:, 15]) # tp flagp[:, 16] = consiste_proc.iguais(dados[:, 15], 5, flagp[:, 16]) # dp
# Testes de consistencia dos dados processados
# Teste 1 - faixa
flagp[:, 1] = consiste_proc.faixa(hs, t1a, t1b, t1c, t1d, flagp[:, 1])
# flagp[:,2] = consiste_proc.faixa(hmax,0,35,0.5,20,flagp[:,2])
# flagp[:,3] = consiste_proc.faixa(tp,3,30,4,18,flagp[:,3])
# flagp[:,4] = consiste_proc.faixa(dp,0,360,0,360,flagp[:,4])
# Teste 2 - Variabilidade temporal
flagp[:, 1] = consiste_proc.variab(hs, t2a, t2b, flagp[:, 1])
# flagp[:,2] = consiste_proc.variab(hmax,1,5,flagp[:,2])
# flagp[:,3] = consiste_proc.variab(tp,1,20,flagp[:,3])
# flagp[:,4] = consiste_proc.variab(dp,1,360,flagp[:,4])
# Teste 3 - Valores consecutivos iguais (*verificar num de arquivos em 'listac')
flagp[:, 1] = consiste_proc.iguais(hs, t3a, flagp[:, 1])
# flagp[:,2] = consiste_proc.iguais(hmax,15,flagp[:,2])
# flagp[:,3] = consiste_proc.iguais(tp,20,flagp[:,3])
# flagp[:,4] = consiste_proc.iguais(dp,20,flagp[:,4])
# ================================================================================== #
# Coloca nan nos dados reprovados
# data em julianos
# matondap = np.array([pl.date2num(data),hs,hmax,tp,dp]).T
matondap = np.array([data, hs, hmax, tp, dp]).T
for c in range(1, flagp.shape[1]):
for i in range(len(flagp)):