# Gravity waves, 1s datasets (real2, imag2, ysig) August 10, 2000 #################################################################### model gravity; const M=701; var Yr[M],Yi[M], Ysig2[M],sr[M], si[M], psi[M], f[M], tau[M], N, phi0, tc, theta1, theta2, theta3, theta4, theta5, z, eta, m1, m2, mt; #data Yr in "real.txt", Yi in "imag.txt", Ysig2 in "ysig.txt"; # frequency range for this data set is 30-730 Hz # inits tc, m1, m2, phi0, N in "in1s.in"; { #distribution of Ys ################### eta <- m1*m2/((m1+m2)*(m1 + m2)); mt <- m1 + m2; z<- 0.000015481*(m1+m2); theta1<- (3/(128*eta)) *pow(z,-5/3); theta2<- (1/(384*eta)) * (3715/84 + 55*eta)/z; theta3<- (-1/(128*eta)) * 48 * 3.14159 * pow(z,-2/3); theta4<- (3/(128*eta)) *(15293365/508032 + (27145/504)*eta + (3085/72)*pow(eta,2)) *pow(z,-1/3); theta5<- (3.14159/(128*eta))*( 38645/252 + 5*eta); for (i in 1:M) { f[i]<- 29 + i; psi[i]<- theta1*pow(f[i],-5/3) + theta2*pow(f[i],-1.0) + theta3*pow(f[i],-2/3) + theta4*pow(f[i],-1/3) + theta5*log(f[i]); sr[i] <- N*pow(eta,1/2)*pow(m1 + m2,5/6)*pow(f[i],-7/6)*( cos(-phi0 + psi[i] + (2*3.14159*tc*f[i])) ); si[i] <- N*pow(eta,1/2)*pow(m1 + m2,5/6)*pow(f[i],-7/6)*( sin(-phi0 + psi[i] + (2*3.14159*tc*f[i])) ); tau[i] <- Ysig2[i]/10; Yr[i]~ dnorm(sr[i],tau[i]); Yi[i]~ dnorm(si[i],tau[i]); } #priors# ######## tc ~ dunif(-0.0021,0); # z ~ dunif(0,0.000150962); # eta ~ dunif(0,0.25); m1 ~ dunif(0.3,10); m2 ~ dunif(0.3,10); phi0 ~ dunif(-6.283185,6.283185); # N ~ dnorm(0,0.0001); N ~ dunif(-1000000,1000000); # N2 ~ dunif(-1000000,1000000); }