# coding: utf-8

# In[3]:


import math


# In[32]:


#t1 = 3.82715371E-07
t1 =  -9.68575478E-07 / (1000) # to convert msec # sac data thought this is sec data
t1 = t1IN / (1000)
print("# t1 = "+str(t1))


# In[33]:


sampring_dt = 1/48000
#sampring_dt = 1/(48*1000*1000) # interp
sampring_dt = 1/(48*1000) # 48 Hz in sac to conver 48kHz
sampring_dt = 1/(SAMPIN*1000) # 48 Hz in sac to conver 48kHz

print("# sampling_dt = "+str(sampring_dt))


# In[18]:


y0=0.94772073784559008
y1=0.99999962063224257
y2=0.94775696178821955


# In[22]:


y0=0.999956012
y1=0.999996126 
y2=0.999971807


# In[35]:


y0= 0.947327495
y1=  0.999996185
y2= 0.947768807


# In[ ]:


y0=y0IN
y1=y1IN
y2=y2IN

print("# y0 = "+str(y0))
print("# y1 = "+str(y1))
print("# y2 = "+str(y2))

# In[36]:


temp1 = y0+y2
temp2 = 2*y1
temp3 = temp1 / temp2
alpha = math.acos(temp3)
#print(alpha)
print("# alpha = "+str(alpha))


# In[37]:


temp4 = temp2 * math.sin(alpha)
temp5 = y0 - y2
temp6 = temp5  / temp4
temp7 = math.atan(temp6)
# dt2 is correct
# dt1 is incrrect
#dt = alpha / temp7
#print(dt)
dt2 = -1 * (temp7 / alpha )
#print(dt2)
print("# dt2 = "+str(dt2))


# In[38]:


#tau = dt2 * sampring_dt;
#tau_micro = tau * 1000 * 1000
#print(tau_micro)


# In[39]:


tau = dt2 * sampring_dt
#print(tau)
print("# tau = "+str(tau))


# In[40]:


tauOUT=t1+tau
#print(tauOUT)
print("# tauOUT = "+str(tauOUT))


# In[41]:


tau_microOUT = tauOUT * 1000 * 1000
#print(tau_microOUT)
print("# tau_microOUT = "+str(tau_microOUT))