;      dBz_detrend_tmp = dBz
;      weights=1./dbz(inan)
;      A = [max(dBz(inan)), dbz(inan(0))]
;      ycfit=curvefit(time(inan),dbz(inan),weight,A,Function_Name='sinfunct')
;        dBz_detrend_tmp(inan) = dbz(inan)-ycfit
;        result = poly_fit(time(inan),dBz_detrend_tmp(inan),jfitBz,yfit=yfit)
;        dBz_detrend(inan) = dBz_detrend_tmp(inan)-yfit

;jfitBx = 5
;jfitBz = 5
;dBz_detrend = dbz
;inan = where(finite(dbz),cnan)
;  result = poly_fit(time(inan),dbz(inan),jfitBz,yfit=yfit)
;  dBz_detrend(inan) = dbz(inan)-yfit
;
;dBx_detrend=dbx
;inan = where(finite(dbx),cnan)
;  result = poly_fit(time(inan),dbx(inan),jfitBx,yfit=yfit)
;  dBx_detrend(inan) = dbx(inan)-yfit
;
;ilat2 = where(abs(glat) lt 45. and finite(dBz_detrend), clat2)
;meandBz = total(dBz_detrend(ilat2))/clat2
;dBz_detrend = dBz_detrend-meandBz
;ilat2 = where(abs(glat) lt 45. and finite(dBx_detrend), clat2)
;meandBx = total(dBx_detrend(ilat2))/clat2
;dBx_detrend = dBx_detrend-meandBx

      jfitBy = 5 
      dBy_detrend_tmp = dBy
      inan = where(finite(dBy),cnan)
      weights=1./dby(inan)
      A = [max(dBy(inan)), dby(inan(0))]
      ycfit=curvefit(time(inan),dby(inan),weight,A,Function_Name='msinfunct')
;      ycfit=curvefit(time(inan),dby(inan),weight,A,Function_Name='sinfunct')
        dBy_detrend_tmp(inan) = dby(inan)-ycfit
        result = poly_fit(time(inan),dBy_detrend_tmp(inan),jfitBy,yfit=yfit)
        dBy_detrend(inan) = dBy_detrend_tmp(inan)-yfit

;window, 0
;plot, dbZ
;oplot, ycfit
;oplot, dBz_detrend_tmp, color=254
;oplot, dBz_detrend, color=60
;
;      ilat2 = where(abs(glat) lt 45. and finite(dBz_detrend), clat2)
;      meandBz = total(dBz_detrend(ilat2))/clat2
;      dBz_detrend = dBz_detrend-meandBz

      ilat2 = where(abs(glat) lt 45. and finite(dBy_detrend), clat2)
      meandBy = total(dBy_detrend(ilat2))/clat2
      dBy_detrend = dBy_detrend-meandBy

      inanx = where(finite(dbx_detrend), cnanx)
      if cnanx gt 1 then dbx_tmp = interpol(dbx_detrend(inanx),time(inanx),time_utd)
      inany = where(finite(dby_detrend), cnany)
      if cnany gt 1 then dby_tmp = interpol(dby_detrend(inany),time(inany),time_utd)
      inanz = where(finite(dbz_detrend), cnanz)
      if cnanz gt 1 then dbz_tmp = interpol(dbz_detrend(inanz),time(inanz),time_utd)

    Sx = 1e-6*(Ey*dBz_tmp-Ez*dBy_tmp)/mu ; erg/(sec-cm^2)
    Sx_ey = 1e-6*Ey*dBz_tmp/mu ; erg/(sec-cm^2)
    Sx_ez = -1e-6*Ez*dBy_tmp/mu ; erg/(sec-cm^2)
    Sx0 = 1e-6*(Ey0*dBz_tmp-Ez1*dBy_tmp)/mu ; erg/(sec-cm^2)

  outfn = 'poynt_F'+string(dat.sat,format='(i2.2)')+'-'+string(outyear,format='(i4.4)')+$
          string(outdoy,format='(i3.3)')+string(outhour,format='(i2.2)')+$
          string(outmin,format='(i2.2)')+string(outsec,format='(i2.2)')

  f = { name: 'Poynting Flux'}
  add_str_element, f, 'fn', outfn
  add_str_element, f, 'sat', dat.sat
  add_str_element, f, 'year', outyear
  add_str_element, f, 'doy', outdoy
  add_str_element, f, 'hour', outhour
  add_str_element, f, 'min', outmin
  add_str_element, f, 'sec', outsec
  add_str_element, f, 'time', time
  add_str_element, f, 'NS', NS
  add_str_element, f, 'glat', glat
  add_str_element, f, 'glon', glon
  add_str_element, f, 'alt', alt
  add_str_element, f, 'eflx', eflx
  add_str_element, f, 'iflx', iflx
  add_str_element, f, 'dbx_o', dbx
  add_str_element, f, 'dby_o', dby
  add_str_element, f, 'dbz_o', dbz
  add_str_element, f, 'dbx_d', dbx_detrend
  add_str_element, f, 'dby_d', dby_detrend
  add_str_element, f, 'dbz_d', dbz_detrend
  add_str_element, f, 'fitdBx', jfitBx
  add_str_element, f, 'fitdBy', jfitBy
  add_str_element, f, 'fitdBz', jfitBz 
  add_str_element, f, 'time_utd', time_utd
  add_str_element, f, 'lat_utd', glat_utd
  add_str_element, f, 'lon_utd', glon_utd
  add_str_element, f, 'alt_utd', alt_utd
  add_str_element, f, 'mlt_utd', mlt_utd
  add_str_element, f, 'mlat_utd', mlat_utd
  add_str_element, f, 'vx', vx
  add_str_element, f, 'vy', vy
  add_str_element, f, 'vz', vz
  add_str_element, f, 'Ni', ni
  add_str_element, f, 'denO', denO
  add_str_element, f, 'QuaR', quaR
  add_str_element, f, 'QuaI', QuaI
  add_str_element, f, 'QuaE', QuaE
  add_str_element, f, 'Ex', Ex
  add_str_element, f, 'Ey', Ey
  add_str_element, f, 'Ez', Ez
  add_str_element, f, 'Ey0', Ey0
  add_str_element, f, 'dBx', dBx_tmp
  add_str_element, f, 'dBy', dBy_tmp
  add_str_element, f, 'dBz', dBz_tmp
  add_str_element, f, 'Sx', Sx
  add_str_element, f, 'Sx0', Sx0
  add_str_element, f, 'Sx_ey', Sx_ey
  add_str_element, f, 'Sx_ez', Sx_ez
  add_str_element, f, 't60', time60
  add_str_element, f, 'mlt60', mlt60
  add_str_element, f, 'mlat60', mlat60
  add_str_element, f, 'mlon60', mlon60

plot_poynt_utd, f, dir=dir, fileType='x'

;    save, filename=dir+outfn+'.sav', f
;    plot_poynt_utd, f, dir=dir, fileType=fileType