function sza,log,mtime,oldlog

; Test a given test log and mission time, calculate the current rotation period
; and solar zenith angle (using Marty Tomasko's formulas)
; For DISR#2.

; v3 modified to correct Roof 9 data by Lyn 
; v4 modified to work on PC 3/9/4

common already_done,doneit,t,t2,in_lock_time,out_of_lock_time,lock_group

lost_lock=strarr(100)
using=strarr(100)
if n_elements(oldlog) eq 0 then goto,initialize
if log ne oldlog then goto,initialize
if n_elements(doneit) ne 0 then goto,process

initialize:
; Does the sun pulse file already exist?
fn=findfile(log+'/post/anal_results/sun_pulses')
if fn(0) eq '' then list_sun_pulses,log

str=''
t=dblarr(3,2000)
ttemp=dblarr(3)
openr,unit,log+'/post/anal_results/sun_pulses',/get_lun
readf,unit,str
i=0
while not eof(unit) do begin
  readf,unit,ttemp
  t(*,i)=ttemp
  i=i+1
endwhile
free_lun,unit

t=t(*,0:i-1)
t2=reform(t(1,*))

; Read the message_log to find in-lock and out-of-lock times
in_lock_time=dblarr(100)
out_of_lock_time=dblarr(100)
results='string'

openr,ol_in,log+'/message_log',/get_lun
halt=0
n_usings=0
while not eof(ol_in) do begin
  readf,ol_in,results
  junk=strpos(results,'using_sun_pulses')
  if junk ne -1 then begin
    using(n_usings)=results   
    n_usings=n_usings+1
  endif
endwhile
using=using(0:n_usings-1)

for i=0,n_usings-1 do begin
  pieces=strsplit(using(i),'Time: ',/extract,/regex)
  time_extract,pieces(1),hr,min,sec,ilt
  in_lock_time(i)=double(ilt)
endfor
in_lock_time=in_lock_time(0:i-1)
openr,ol_in2,log+'/post/anal_results/sun_pulses',/get_lun
while not (eof(ol_in2)) do begin
  readf,ol_in2,results
endwhile
reads,results,o_time,format='(13x,f12.0)'
last_time=o_time+0.01

openr,ol_in3,log+'/message_log',/get_lun
halt=0
ll_flag=0
n_lost=0
while not eof(ol_in3) do begin
  readf,ol_in3,results
  junk=strpos(results,'SunSensor_sig_lost ')
  if junk ne -1 then begin
    ll_flag=1
    lost_lock(n_lost)=results
    n_lost=n_lost+1
  endif
endwhile
lost_lock=lost_lock(0:n_lost-1)
if results ne 'code:  55 - SunSensor_sig_lost         Id:    99  Time: 00:47:26.3645' and $
  n_lost ne 0 then begin
  k=0
  for j=0,(n_usings<n_lost)-1 do begin
    pieces=strsplit(lost_lock(j),'Time: ',/extract,/regex)
    time_extract,pieces(1),hr,min,sec,olt
    olt=double(olt)
    if j eq 0 then begin 
      out_of_lock_time(k)=olt
      k=k+1
    endif
    out_of_lock_time(k)=last_time
    if j gt 0 and olt gt in_lock_time(k) then begin
      out_of_lock_time(k)=olt
      k=k+1
    endif
  endfor
  out_of_lock_time=out_of_lock_time(0:i-1)
endif else begin ; if never lost lock!
  openr,ol_in4,log+'/post/anal_results/sun_pulses',/get_lun
  while not (eof(ol_in4)) do begin
    readf,ol_in4,results
  endwhile
  reads,results,o_time,format='(13x,f12.0)'
  out_of_lock_time(0)=o_time+0.01D0
endelse
if out_of_lock_time(i-1) eq 0.0 then begin
  openr,ol_in5,log+'/post/anal_results/sun_pulses',/get_lun
  while not (eof(ol_in5)) do begin
    readf,ol_in5,results
  endwhile
  reads,results,o_time,format='(13x,f12.0)'
  out_of_lock_time(i-1)=o_time+0.01D0
endif
; Assign a "lock group" number to each triplet
lock_group=intarr(n_elements(t))-1
for j=0,i-1 do begin
  index=where(t2 ge in_lock_time(j) and t2 le out_of_lock_time(j))
  if index(0) ne -1 then lock_group(index)=j
endfor
doneit=1
oldlog=log

process:
for l=0,n_elements(in_lock_time)-1 do begin
  if mtime ge in_lock_time(l) and mtime le out_of_lock_time(l) then goto,locked
endfor
results=[-1.0,-1.0]
return,results

locked:
index=where(t2 ge mtime)
index=index(0)
if index eq -1 then return,[-1.0,-1.0]
; Case 1 -- three consecutive pulses in lock
if (lock_group(index-1) eq lock_group(index)) and (lock_group(index) eq $
  lock_group(index+1)) then begin
  p1=t2(index)-t2(index-1)
  p2=t2(index+1)-t2(index)
  tt=dblarr(3)
  p=dblarr(3)
  tt(0)=t2(index-1)
  tt(1)=t2(index)
  tt(2)=t2(index+1)
  p(0)=(t2(index-1)-0.5D0*(tt(0)+tt(1)))/(0.5D0*(tt(2)-tt(0)))*(p2-p1)+p1
  p(1)=(t2(index)-0.5D0*(tt(0)+tt(1)))/(0.5D0*(tt(2)-tt(0)))*(p2-p1)+p1
  p(2)=(t2(index+1)-0.5D0*(tt(0)+tt(1)))/(0.5D0*(tt(2)-tt(0)))*(p2-p1)+p1
  deltaz=dblarr(3)
  deltaz(0)=(t(0,index-1)-t(2,index-1))/p(0)*360.D0
  deltaz(1)=(t(0,index)-t(2,index))/p(1)*360.D0
  deltaz(2)=(t(0,index+1)-t(2,index+1))/p(2)*360.D0
  sza=100.91D0+4.2575D0*deltaz+0.11029D0*deltaz^2+0.0014921D0*deltaz^3+ $
    8.0614D-06*deltaz^4
  coefs=poly_fit(tt-tt(0),sza,2)
  x=coefs(0)+coefs(1)*(mtime-tt(0))+coefs(2)*(mtime-tt(0))^2
  p=(mtime-0.5D0*(tt(0)+tt(1)))/(0.5D0*(tt(2)-tt(0)))*(p2-p1)+p1
  results=[x,p]
  return,results
endif

; Case 2 only first two pulses in lock
if (lock_group(index-1) eq lock_group(index)) and (lock_group(index) lt $
  lock_group(index+1)) then begin
  if mtime gt t2(index) then begin
    results=[-1.0,-1.0] & STOP
    return,results
  endif
  p=t2(index)-t2(index-1)
  deltaz=0.5D0*(t(0,index-1)-t(2,index-1)+t(0,index)-t(2,index))
  ; deltaz should now be converted to degrees
  sza=100.91D0+4.2575D0*deltaz+0.11029D0*deltaz^2+0.0014921D0*deltaz^3+ $
    8.0614D-06*deltaz^4
  results=[sza,p]
  return,results
endif

; Case 3 only last two pulses in lock
if (lock_group(index-1) lt lock_group(index)) and (lock_group(index) eq $
  lock_group(index+1)) then begin
  if mtime lt t2(index) then begin
    results=[-1.0,-1.0]
    return,results
  endif
  p=t2(index+1)-t2(index)
  deltaz=0.5D0*(t(0,index)-t(2,index)+t(0,index+1)-t(2,index+1))
  ; deltaz should now be converted to degrees
  sza=100.91D0+4.2575D0*deltaz+0.11029D0*deltaz^2+0.0014921D0*deltaz^3+ $
    8.0614D-06*deltaz^4
  results=[sza,p]
  return,results
endif
close,/all
return,[-1.0,-1.0]
end

; cmd=string(log,format='("tail -n1 ",a,"/post/anal_results/sun_pulses")')
; spawn,cmd,results
; cmd=string(log,format='("tail -n1 ",a,"/post/anal_results/sun_pulses")')
; spawn,cmd,results
;cmd=string(log+'/message_log',format='("grep SunSensor_sig_lost ",a)')
;spawn,cmd,results
;n_lost=n_elements(results)
;if results(0) ne 'code:  55 - SunSensor_sig_lost         Id:    99  Time: 00:47:26.3645' and $
;  results(0) ne 0 then begin
;   reads,pieces(1),hr,min,sec,format='(56x,f2.0,1x,f2.0,1x,f7.4)'
;   olt=3600.D0*double(hr)+60.D0*double(min)+double(sec)
; reads,results(i),hr,min,sec,format='(56x,f2.0,1x,f2.0,1x,f7.4)'
; in_lock_time(i)=3600.D0*double(hr)+60.D0*double(min)+double(sec)