# dados = ws wg wd at rh dw pr st bh cl tu sr hs hm tp dp sp # numero de parametros processados npa = 17 # ================================================================================== # # cria vetores de flags das series processadas (depende das qtdade de variaveis a serem consistidas + data) flagp = np.zeros((len(dados), npa + 1), dtype="|S32") flagp[:, 0] = data.astype(str) # ================================================================================== # # Testes de consistencia dos dados processados # #Teste 1 - faixa flagp[:, 1] = consiste_proc.faixa(dados[:, 0], 0.25, 30, 0.5, 15, flagp[:, 1]) # ws flagp[:, 2] = consiste_proc.faixa(dados[:, 1], 0.25, 40, 0.5, 25, flagp[:, 2]) # wg flagp[:, 3] = consiste_proc.faixa(dados[:, 2], 0, 360, 0, 360, flagp[:, 3]) # wd flagp[:, 4] = consiste_proc.faixa(dados[:, 3], 0, 40, 10, 30, flagp[:, 4]) # at flagp[:, 5] = consiste_proc.faixa(dados[:, 4], 25, 102, 40, 100, flagp[:, 5]) # rh # flagp[:,6] = consiste_proc.faixa(dados[:,5],0,40,10,30,flagp[:,6]) #dw flagp[:, 7] = consiste_proc.faixa(dados[:, 6], 900, 1100, 950, 1050, flagp[:, 7]) # pr flagp[:, 8] = consiste_proc.faixa(dados[:, 7], 12, 38, 15, 35, flagp[:, 8]) # st flagp[:, 9] = consiste_proc.faixa(dados[:, 8], 0, 360, 0, 360, flagp[:, 9]) # bh # flagp[:,10] = consiste_proc.faixa(dados[:,9],???,flagp[:,10]) #cl # flagp[:,11] = consiste_proc.faixa(dados[:,10],???,flagp[:,11]) #tu # flagp[:,12] = consiste_proc.faixa(dados[:,11],0,1500,10,1450,flagp[:,12]) #sr flagp[:, 13] = consiste_proc.faixa(dados[:, 12], 0, 20, 0.25, 8, flagp[:, 13]) # hs # flagp[:,14] = consiste_proc.faixa(dados[:,13],0,35,0.5,20,flagp[:,14]) #hm flagp[:, 15] = consiste_proc.faixa(dados[:, 14], 3, 30, 4, 18, flagp[:, 15]) # tp flagp[:, 16] = consiste_proc.faixa(dados[:, 15], 0, 360, 0, 360, flagp[:, 16]) # dp
hm02 = matondab[:,14] tp2 = matondab[:,15] dp2 = matondab[:,16] gam = matondab[:,17] gam1 = matondab[:,18] gam2 = matondab[:,19] #cria vetores de flags das series processadas (depende das qtdade de variaveis a serem consistidas + data) flagp = np.zeros((len(listap),npa+1),dtype='|S32') flagp[:,0] = listap[:] # ================================================================================== # # Testes de consistencia dos dados processados #Teste 1 - faixa flagp[:,1] = consiste_proc.faixa(hs,0,20,0.25,8,flagp[:,1]) flagp[:,2] = consiste_proc.faixa(h10,0,20,0.25,8,flagp[:,2]) flagp[:,3] = consiste_proc.faixa(hmax,0,35,0.5,20,flagp[:,3]) flagp[:,4] = consiste_proc.faixa(tmed,3,30,4,18,flagp[:,4]) flagp[:,5] = consiste_proc.faixa(thmax,3,30,4,18,flagp[:,5]) flagp[:,6] = consiste_proc.faixa(hm0,0,20,0.25,8,flagp[:,6]) flagp[:,7] = consiste_proc.faixa(tp,3,30,4,18,flagp[:,7]) flagp[:,8] = consiste_proc.faixa(dp,0,360,0,360,flagp[:,8]) flagp[:,9] = consiste_proc.faixa(sigma1p,0,360,0,360,flagp[:,9]) flagp[:,10] = consiste_proc.faixa(sigma2p,0,360,0,360,flagp[:,10]) flagp[:,11] = consiste_proc.faixa(hm01,0,20,0.25,8,flagp[:,11]) flagp[:,12] = consiste_proc.faixa(tp1,3,30,4,18,flagp[:,12]) flagp[:,13] = consiste_proc.faixa(dp1,0,360,0,360,flagp[:,13]) flagp[:,14] = consiste_proc.faixa(hm02,0,20,0.25,8,flagp[:,14]) flagp[:,15] = consiste_proc.faixa(tp2,3,30,4,18,flagp[:,15]) flagp[:,16] = consiste_proc.faixa(dp2,0,360,0,360,flagp[:,16])
print 'Iniciando processamento em ... ' + local # ================================================================================== # #cria vetores de flags das series processadas (depende das qtdade de variaveis a serem consistidas + data) flagp = np.zeros((len(dados),npa+1),dtype='|S32') flagp[:,0] = dados[:,0].astype(int).astype(str) # ================================================================================== # # Testes de consistencia dos dados processados #Teste 1 - faixa #sintaxe: faixa(serie,linf_inst,lmax_inst,lmin_reg,lmax_reg,flag) print 'Realizando teste de Faixa' flagp[:,1] = consiste_proc.faixa(dados[:,1],0.1,35,0.2,25,flagp[:,1]) #ws flagp[:,2] = consiste_proc.faixa(dados[:,2],0.1,35,0.3,25,flagp[:,2]) #wg flagp[:,3] = consiste_proc.faixa(dados[:,3],0,360,0,360,flagp[:,3]) #wd flagp[:,4] = consiste_proc.faixa(dados[:,4],0,40,6,30,flagp[:,4]) #at flagp[:,5] = consiste_proc.faixa(dados[:,5],980,1040,985,1035,flagp[:,5]) #pr flagp[:,6] = consiste_proc.faixa(dados[:,6],10,35,13,30,flagp[:,6]) #wt flagp[:,7] = consiste_proc.faixa(dados[:,7],0,20,0.15,7,flagp[:,7]) #hs flagp[:,8] = consiste_proc.faixa(dados[:,8],1.95,26,3,22,flagp[:,8]) #tp - ndbc09 flagp[:,9] = consiste_proc.faixa(dados[:,9],0,360,0,360,flagp[:,9]) #dp #Teste 2 - Variabilidade #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
hmax = sens[:, 21] tp = sens[:, 22] dp = sens[:, 23] # numero de parametros a serem processados npa = 4 # cria vetores de flags das series processadas (depende das qtdade de variaveis a serem consistidas + data) flagp = np.zeros((len(data), npa + 1), dtype="|S32") flagp[:, 0] = data # listap[:] # ================================================================================== # # 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])