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
