dispd = 'x'
id='20031024_0246
st='24-Oct-03 02:00:00'
ed='24-Oct-03 03:15:00'
efile='efz20031024.030011'
coord = [-864, -349]
so_co = [250,260]

; Make Lightcurve, fixed the peak & initial time
gs = ogoes()
gs -> set, tstart=st, tend =ed, mode=0
gs_str = gs -> getdata(/str)
gs_tim = anytim(gs_str.utbase,/sec) + gs_str.tarray
gs_low = gs_str.yclean


max_goes = max(gs_low, cnt1)
print, 'GOES Peak Time:'+anytim(gs_tim(cnt1), /yohkoh)
goes_pt = anytim(gs_tim(cnt1), /yohkoh)

qv = max_goes*0.3
qv_ss = min(abs(gs_low[0:cnt1]-qv), cnt2)
print, 'Initial Time (GOES 30%):'+anytim(gs_tim(cnt2[0]),/yohkoh)
goes_it = anytim(gs_tim(cnt2[0]),/yohkoh)

norh_rd_tcx,st, ed, index,data
max_tcx =max(data[0:18000], cnt3)
print, 'Radio Peak Time:'+anytim(index[cnt3],/yohkoh)
radio_pt = anytim(index[cnt3],/yohkoh)

loadct,0
if dispd eq 'ps' then begin
   set_plot,'ps'
   device, filename=id+'_lc.ps'
   !p.font=0
endif else begin
   window,1
endelse

utplot, gs_tim, gs_low, /ylog, $
        title=id+' GOES_Low, NoRH_tcx, Lightcurve', $
        subtitle='SXR Peak:'+goes_pt+' / SXR 30%:'+goes_it+' / 17GHz Peak:'+radio_pt
outplot, index, data*1e-4, linestyle=3
outplot, [goes_it, goes_it],[10e-10, 1]
outplot, [goes_pt, goes_pt],[10e-10, 1]
outplot, [radio_pt, radio_pt],[10e-10, 1]

if dispd eq 'ps' then begin
   device,/close
   set_plot,'x'
   !p.font=-1
endif 

;stop

;Make movie
int = 30
if 0 then begin
   norh_trans, st, ed, int
   norh_trans, st, ed, int, freq=34
   norh_synth, st, ed, int, cenfnl=coord, size=128
   norh_synth, st, ed, int, cenfnl=coord, size=512, freq=34
endif

   files17 = findfile('ipa*')
   files34 = findfile('ipz*')
   norh_rd_img, files17, ind17, dat17
   norh_rd_img, files34, ind34, dat34
   norh_index2map, ind17, dat17, amap
   norh_index2map, ind34, dat34, zmap
   
   loadct,0
   set_plot,'z'
   device, set_resolution=[512,512]
   for i=0, n_elements(zmap)-1 do begin
      plot_map, zmap[i], fov=3, /log, dmin=8e3, grid_spacing=5
      plot_map, amap[i], /over
      img = tvrd() 
      if i eq 0 then vtr=img else vtr=[[[vtr]],[[img]]] 
   endfor
   set_plot,'x'

;stop

mreadfits, efile, index, data
index2map, index, data, emap
eemap = map2earth(emap)
dummy=min(abs(anytim(zmap.time, /sec)-anytim(radio_pt, /sec)), cnt)

wdef,512
loadct,1
plot_map, eemap, fov=2, cen=zmap[cnt],/log
plot_map, zmap[cnt], /over

;Make image
if 0 then begin
   norh_trans, goes_it, goes_pt, 10
   norh_trans, goes_it, goes_pt, 10, freq=34
   norh_synth, goes_it, /usetsx
   norh_synth, goes_it, freq=34, /usetsx
   norh_synth, goes_pt, /usetsx
   norh_synth, goes_pt, freq=34, /usetsx
   norh_synth, radio_pt, /usetsx
   norh_synth, radio_pt, freq=34, /usetsx
endif

; make goes initial time 
git = anytim(goes_it,/ex)
gitfileh = string(git[6]-2000, format='(I02)') + $
           string(git[5], format='(I02)') + $
           string(git[4], format='(I02)') + '_' +$
           string(git[0], format='(I02)') + $
           string(git[1], format='(I02)') + $
           string(git[2]+round(git[3]/1000.), format='(I02)') 

;make map at goes initial time
norh_rd_img, 'ifa'+gitfileh, indgia, datgia
norh_rd_img, 'ifz'+gitfileh, indgiz, datgiz
norh_index2map, indgia, datgia, giamap
norh_index2map, indgiz, datgiz, gizmap

;plot map
loadct,10
wdef,512
plot_map, gizmap, fov=2, cen=zmap[0], dmin=8e3, grid_spacing=5
plot_map, giamap, /over

; make goes peak time 
;gpt = anytim(goes_pt,/ex)
;gptfileh = string(gpt[6]-2000, format='(I02)') + $
;           string(gpt[5], format='(I02)') + $
;           string(gpt[4], format='(I02)') + '_' +$
;           string(gpt[0], format='(I02)') + $
;           string(gpt[1], format='(I02)') + $
;           string(gpt[2]+round(gpt[3]/1000.), format='(I02)') 

;make map at goes pt
;norh_rd_img, 'ipa'+gptfileh, indgpa, datgpa
;norh_rd_img, 'ipz'+gptfileh, indgpz, datgpz
;norh_index2map, indgpa, datgpa, gpamap
;norh_index2map, indgpz, datgpz, gpzmap

;wdef,512
;plot_map, gpzmap, fov=2, cen=zmap[0], dmin=8e3, grid_spacing=5
;plot_map, gpamap, /over


rpt = anytim(radio_pt,/ex)
rptfileh = string(rpt[6]-2000, format='(I02)') + $
          string(rpt[5], format='(I02)') + $
          string(rpt[4], format='(I02)') + '_' +$
          string(rpt[0], format='(I02)') + $
          string(rpt[1], format='(I02)') + $
          string(rpt[2]+round(rpt[3]/1000.), format='(I02)') 

norh_rd_img, 'ifa'+rptfileh, indrpa, datrpa
norh_rd_img, 'ifz'+rptfileh, indrpz, datrpz
norh_index2map, indrpa, datrpa, rpamap
norh_index2map, indrpz, datrpz, rpzmap

wdef,512
plot_map, rpzmap, fov=2, cen=zmap[0],dmin=8e3, grid_spacing=5
plot_map, rpamap, /over


;

dummy = max(datrpa,cnt)
xcen = (cnt mod 512) - 256
ycen = fix(cnt / 512) - 256

ang_rot = atan(float(xcen)/float(ycen))*180/!pi*(-1)

datrpaR = rot(datrpa, ang_rot, /interp)
datrpzR = rot(datrpz, ang_rot, /interp)
;datgpaR = rot(datgpa, ang_rot, /interp)
;datgpzR = rot(datgpz, ang_rot, /interp)
datgiaR = rot(datgia, ang_rot, /interp)
datgizR = rot(datgiz, ang_rot, /interp)

wdef,512
tvscl, bytscl(datrpaR, max=5e5)
plots,[256,256],[0,512],/dev
plots,[so_co[0],so_co[0]],[0,512], /dev, linestyle=2
plots,[so_co[1],so_co[1]],[0,512], /dev, linestyle=2

wdef,512
tvscl, rebin(bytscl(datrpzR, max=5e5),512,512)
plots,[256,256],[0,512],/dev
plots,[so_co[0],so_co[0]],[0,512], /dev, linestyle=2
plots,[so_co[1],so_co[1]],[0,512], /dev, linestyle=2

x = (float(findgen(512))-256.)*amap[0].dx
xx = (float(findgen(1024))-512.)*zmap[0].dx

print, '---Distance from Sun Center---'
rp17pro = total(datrpaR[so_co[0]:so_co[1],*],1)
dummy = max(rp17pro, cnt)
print, ' 17GHz at Radio Peak: '+string(x[cnt])+' arcsec'
rp34pro = total(datrpzR[so_co[0]*2:so_co[1]*2,*],1)
dummy = max(rp34pro, cnt)
print, ' 34GHz at Radio Peak: '+string(xx[cnt])+' arcsec'

;gp17pro = total(datgpaR[so_co[0]:so_co[1],*],1)
;dummy = max(gp17pro[0:60], cnt)
;print, ' 17GHz at GOES Peak : '+string(x[cnt])+' arcsec'
;gp34pro = total(datgpzR[so_co[0]*2:so_co[1]*2,*],1)
;dummy = max(gp34pro[0:120], cnt)
;print, ' 34GHz at GOES Peak : '+string(xx[cnt])+' arcsec'

gi17pro = total(datgiaR[so_co[0]:so_co[1],*],1)
dummy = max(gi17pro[0:60], cnt)
print, ' 17GHz at GOES 30%  : '+string(x[cnt])+' arcsec'
gi34pro = total(datgizR[so_co[0]*2:so_co[1]*2,*],1)
dummy = max(gi34pro[0:120], cnt)
print, ' 34GHz at GOES 30%  : '+string(xx[cnt])+' arcsec'

if dispd eq 'ps' then begin
   set_plot,'ps'
   device, filename=id+'_sp.ps',/color
   !p.font=0
endif else begin
   window,2
endelse

loadct,0 
linecolors


plot, x, rp17pro, /ylog, xr=[-1100, 800],$
title='Radio Source Profiles',$
ytitle='Brightness Temp. (K)', $
xtitle='Distance from Sun Center (arcsec)', $
subtitle='Solid: 17GHz, Dashed: 34GHz, Black:Radio Peak, Red: GOES Peak Blue: 30% GOES'
oplot, x, rebin(rp34pro,512,/sample), linestyle=2
;oplot, x, gp17pro, color=2
;oplot, x, rebin(gp34pro,512,/sample), color=2, linestyle=2
oplot, x, gi17pro, color=10
oplot, x, rebin(gi34pro,512,/sample), color=10, linestyle=2

if dispd eq 'ps' then begin
   device,/close
   set_plot,'x'
   !p.font=-1
endif 


end