t = t[0:nlin] eta = eta[0:nlin] etax = etax[0:nlin] etay = etay[0:nlin] #cria matriz com serie de heave (para o calculo do multiplicador do desvio padrao) eta_mat.append(eta) eta_med.append(np.mean(eta)) eta_dp.append(np.std(eta)) #lista nome dos arquivos consistentes listac.append(listap[cont]) #processamento no dominio do tempo hs,h10,hmax,tmed,thmax = proconda.ondat(t,eta,h) #processamento no dominio da frequencia hm0, tp, dp, sigma1, sigma2, sigma1p, sigma2p, freq, df, k, sn, snx,sny, snn, snnx, snny, snxny, snxnx, snyny, a1, b1, a2, b2, dire1, dire2 = proconda.ondaf( eta,etax,etay,h,nfft,fs) #calcula o espectro de fase (fase e coerencia) #acha o indice da fp indfp = pl.find(sn[:,0]==sn[sn[:,1]==max(sn[:,1]),0]) fase_nnx.append(np.real(snnx[indfp,4])[0]) #fase de heave e dspx fase_nny.append(np.real(snny[indfp,4])[0]) #fase de heave e dspx fase_nxny.append(np.real(snxny[indfp,4])[0]) #fase de heave e dspx coer_nnx.append(np.real(snnx[indfp,5])[0]) #coerencia de heave e dspx coer_nny.append(np.real(snny[indfp,5])[0]) #coerencia de heave e dspx
dc = -1 #loop para processar os dados for arq in lista: dc += 1 print arq dd = np.loadtxt(pathname + arq) pr, vx, vy, vz = dd[:, [14, 2, 3, 4]].T eta = np.concatenate((eta, vz), axis=1) #processamento no dominio do tempo hs, h10, hmax, tmed, thmax = proconda.ondat(t, vz, h) #processamento no dominio da frequencia #hm0_pr, tp_pr, dp_pr, sigma1, sigma2, sigma1p, sigma2p, freq, df, k, sn_pr, snx,sny, snn, snnx, snny, snxny, snxnx, snyny, a1, b1, a2, b2, dire1_pr, dire2 = proconda.ondaf( #pr,vx,vy,h,nfft,fs) ##pressao hm0, tp, dp, sigma1, sigma2, sigma1p, sigma2p, freq, df, k, sn, snx, sny, snn, snnx, snny, snxny, snxnx, snyny, a1, b1, a2, b2, dire1, dire2 = proconda.ondaf( vz, vx, vy, h, nfft, fs) #calcula os espectros de pr, vz, vx e vy apr = espec.espec1(pr, nfft, fs) avz = espec.espec1(vz, nfft, fs) avx = espec.espec1(vx, nfft, fs) avy = espec.espec1(vy, nfft, fs) #corrige a declinacao magnetica
pl.close('all') #Hs=6,2, Tp= , Dp=130 - Onda com maior Hs dados = np.loadtxt(os.environ['HOME'] + '/Dropbox/lioc/dados/pnboia/triaxys/rio_grande/HNE/200912130600.HNE', skiprows = 11) t = dados[:,0] - dados[0,0] eta = dados[:,1] # ------------------------------------------------- # #calcula parametros de onda para a serie normal #param = [Hs,H10,Hmax,Tmed,THmax] #parametros de onda no tempo etat = proconda.ondat(t,eta,200) #parametros de onda na frequencia etaf = espec.espec1(eta,328,1.28) hm0 = 4.01 * sp.sqrt(sum(etaf[:,1]) * (etaf[0,0])) # ------------------------------------------------- # #cria serie perturbada etap = cp.copy(eta) #etap1 = cp.copy(eta) #numero de perturbacoes npt = 50 #duracao da perturbacao