dispd = 'x'
id='20000716_0122'
st='2000-07-16 02:00:00'
ed='2000-07-16 02:10:00'
efile='efz20000716.021210'
coord = [-923, 182]
so_co = [251,261]
;id='20050825_0438
;st='25-Aug-05 04:00:00'
;ed='25-Aug-05 05:11:00'
;efile='efz20050825.044810'
;coord = [-943, 118]
;so_co = [249,263]

!p.thick=2.0 & !p.charsize=2.0 & !p.font=1

; 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, 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-5, 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 

;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'
   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'


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

zfile='ifz000716_020237'
;pt='2000-07-16 02:02:36'
;norh_trans,pt
;norh_synth,pt
;norh_trans,pt,freq=34
;norh_synth,pt,freq=34

norh_rd_img,zfile,zind,zdat
norh_index2map,zind,zdat,zzmap

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

;Make image
if 1 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

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)')
 
;gitfileh='020829_024827'
norh_rd_img, 'ifa'+gitfileh, indgia, datgia
norh_rd_img, 'ifz'+gitfileh, indgiz, datgiz
norh_index2map, indgia, datgia, giamap
norh_index2map, indgiz, datgiz, gizmap

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

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)') 

;gptfileh='020829_025333'
norh_rd_img, 'ifa'+gptfileh, indgpa, datgpa
norh_rd_img, 'ifz'+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)') 

;-kokokara-

;rptfileh='020829_025135'
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
;rdpix,datrpaR

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], format='(F10.3)')+' 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], format='(F10.3)')+' arcsec'

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

gi17pro = total(datgiaR[so_co[0]:so_co[1],*],1)
dummy = max(gi17pro, cnt)
print, ' 17GHz at GOES 30%  : '+string(x[cnt], format='(F10.3)')+' arcsec'
gi34pro = total(datgizR[so_co[0]*2:so_co[1]*2,*],1)
dummy = max(gi34pro, cnt)
print, ' 34GHz at GOES 30%  : '+string(xx[cnt], format='(F10.3)')+' 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=[800, 1100],yr=[1e2,1e10],$
;plot, x, rp17pro, /ylog, xr=[-800, -1100],yr=[1e2,1e10],$
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 Green: 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=7
oplot, x, rebin(gi34pro,512,/sample), color=7, linestyle=2

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


end