pro baseline_phase,filename,dbx,dby,dbz
   restore,filename
   ha = fltarr(n_elements(tavg))
   for i = 0, n_elements(tavg)-1 do begin
      ; Calculate HA in degrees
      ha[i] = (ovsa_lst(cvdoy(geometry.year,geometry.doy[0],/JUL) + (tavg[i]+geometry.dut)/8.64D7)*15/3.6e6 - calinfo.ra/10000.)
   endfor
   ha = ha*!dtor
   dec = (calinfo.dec/10000.)*!dtor
   dphi = 2*!pi*(dbx*cos(dec)*cos(ha) - dby*cos(dec)*sin(ha) + dbz*sin(dec))/!dtor
   window,0,xsiz=500,ysiz=800
   !p.multi=[0,1,2]
   j = 5
   plot,ha/!dtor,lobe(atan(avg[5,j,*,0],avg[6,j,*,0])/!dtor- dphi*f[j],/mid),yran=[-180,180],psym=1
   plot,ha/!dtor,lobe(atan(avg[7,j,*,0],avg[8,j,*,0])/!dtor- dphi*f[j],/mid),yran=[-180,180],psym=1
   ;for i = 1, 5 do oplot,ha/!dtor,atan(avg[5,j,*,0],avg[6,j,*,0])/!dtor+360*i,psym=1
   ;oplot,ha/!dtor,lobe(dphi*f[j])
;   plot,ha/!dtor,atan(avg[7,j,*,0],avg[8,j,*,0])/!dtor-180,yran=[0,2000],psym=1
;   for i = 1, 5 do oplot,ha/!dtor,atan(avg[7,j,*,0],avg[8,j,*,0])/!dtor+360*i-180,psym=1
;   oplot,ha/!dtor,dphi*f[j]
return
end