;pro rat_aia
white=255&black=0&red=200&blue=100&green=40&pink=130&darkgreen=17&cyan=90&violet=125
cc=[white,red,green,blue,cyan,pink,violet,darkgreen,black,white];spisok cvetov dlja chastot
c2=red
rmag=0.998; sdo radius magnification factor
xshift=2; sdo shift
;xshift=-4; dlja sdo 1600 pjaten
;xshift=0
;fragment=[0,0,2400,2400]
;fragment=[0,0,1100,1100]
;fragment=[950,0,600,600]
fragment=[0,-150,1000,1000]; dlja 19/11/2014
;fragment=[-50,-250,400,400]; dlja 19/11/2014 jadro
pow=0.7;kontrastnost
ymin1=1000;0; naskolko pripodnjat skany
freqs=[4,7,10,14,16.4];spisok chastot
;freqs=[4,6,7,8,9,10,11,16.4];spisok chastot
cleanplot
device,decomp=0
window,xs=700,ys=700 ; razmer okna dlja fragmenta
;window,xs=800,ys=800 ; razmer okna dlja full disk
linethick=0; tolschina linij
chrsz=1.2 ; razmer schrifta



;f=pickfile()
;f='G:\2015-kaltman\sdo\saia_00171_fd_20141127_082959.fts.gz'
f='G:\2015-kaltman\sdo\saia_00171_fd_20141119_083059.fts.gz'
;fr='G:\2015-kaltman\20141127_120111_sun0_calib_iv_shift.fits'
;f='G:\2015-kaltman\sdo\saia_01600_fd_20141119_083328.fts.gz'
fbase=file_basename(f)
wavelength1=strmid(fbase,6,5)


;fr='G:\2015-kaltman\20141127sun0_17418.fits'
;fr='G:\2015-kaltman\20141127sun0_17418_sub.fits'
fr='G:\2015-kaltman\20141119sun0_22411_sub.fits'
fr='G:\2015-kaltman\20141119sun0_22411.fits'
;;fr='G:\2015-kaltman\20141127_120111_sun0_calib_iv_shift.fits'
;fr='G:\2015-kaltman\20141127sun0_17418_sub.fits'

if not file_test(fr) then goto,err
if not file_test(f) then goto,err

saia=readfits(f,h1); chtenie SDO AIA faila
n1=fxpar(h1,'NAXIS1')
n2=fxpar(h1,'NAXIS2')
n12=n1/2.0
n22=n2/2.0
crpix1= sxpar(h1, 'crpix1')
crpix2= sxpar(h1, 'crpix2')
cdelt1=sxpar(h1,'CDELT1')
cdelt2=sxpar(h1,'CDELT2')
solar_r=sxpar(h1, 'RSUN_OBS')
if solar_r eq 0 then solar_r=sxpar(h1, 'RSUN');

solar_rx=solar_r/cdelt1
saiashift=[n12-crpix1,n22-crpix2]
print, 'saiashift=',saiashift
if saiashift[0] gt 1 then saia=shift(saia,-saiashift)

rat=readfits(fr,h2); chtenie RATAN-600 faila
ee=mrdfits(fr,1,h3)
date_rat=sxpar(h2, 'DATE-OBS')
time_rat=sxpar(h2, 'TIME-OBS')
cdelt1r=sxpar(h2,'CDELT1')
n1r=fxpar(h2,'NAXIS1')
n3r=fxpar(h2,'NAXIS3')
n12r=n1r/2.0
crpix1r= sxpar(h2, 'crpix1')
solar_rr=sxpar(h2, 'SOLAR_R')
solar_rrx=solar_rr/cdelt1r
xshiftr=n12r-crpix1r
help,xshiftr
if abs(xshiftr) gt 1 then rat=shift(rat,[-xshiftr,0,0])
xmas=(findgen(n1r)-n12r)*cdelt1r
freq=10.0

ind=where(ee.freq ge freq-1e-4,count)
if count gt 0 then ifr0=ind[0] else ifr0=0
rbp0=get_rb0p(anytim(date_rat, /yy),/deg)
saia=shift(saia,xshift)
saia=rot(saia,-rbp0[2],rmag,n12,n22,/interp)


hmi_x = (findgen(n1)-n1/2.0)*cdelt1
hmi_y = (findgen(n2)-n2/2.0)*cdelt2
a=fragment
xr1=[a[0]-a[2]*0.5,a[0]+a[2]*0.5]
yr1=[a[1]-a[3]*0.5,a[1]+a[3]*0.5]
x1=xr1[0]/cdelt1r+crpix1r
x2=xr1[1]/cdelt1r+crpix1r

;risovanie SDO AIA
aspect=float(a[2])/a[3]
aia_lct, rgbr, rgbg, rgbb, wavelnth=wavelength1, /load
if chrsz le 1 then xmar1=[8,8] else xmar1=[10,6]
gvi_multiplot, [1,1], asp = aspect, xmargin=xmar1,ymargin=xmar1
tvplot, saia^pow,charsize=chrsz, $
hmi_x,hmi_y,xra=xr1,yra=yr1, miss = 0,title='SUN '+date_rat+' '+strmid(time_rat,0,8)
xyouts,20,!d.y_size-20,file_basename(f,'.fts.gz'),/device,charsize=chrsz
rads=[-solar_r,0,solar_r]
for jj=0,2 do plots,rads[jj],!y.crange, color=colorr,linest=1,thick=linethick


loadct,12 ; nadrisovat cvetnye skany  RATAN-600
n4=n_elements(freqs)
ifreqs=findgen(n4)
for i=0,n4-1 do begin
ifr=where(ee.freq ge freqs[i]-1e-4,count)
if count gt 0 then ifreqs[i]=ifr[0] else ifreqs[i]=0
endfor

y1=min(rat(x1:x2,0,ifreqs),max=y2)-ymin1
plot,xmas,rat[*,0,ifr0],xra=xr1,/noerase,xst=5,yst=5, $
yra=[y1,y2],/nodata;,pos=[0., 0.3, 1., 1.0]
for i=0,n4-1 do begin
oplot, xmas,rat(*,0, ifreqs[i]), color=cc[i],thick=linethick ;,linest=i
endfor
ss=string(ee[ifreqs].freq,format='(f5.1)')+'GHz'

rads=[-solar_rr,0,solar_rr]
for jj=0,2 do plots,rads[jj],!y.crange, color=white, $
linest=1,thick=linethick

legend,ss,colors=cc[0:n4-1],linest=0, $
charsize= chrsz,thick=2,charthick=chrsz,/right,number=0.1


;if part eq 0 then begin
;lats = [-80, -45, -30, -20, 0, 15, 27, 35, 45, 55, 75]
;latnames = strtrim(lats, 2)
;map_set, bb,0,-pp, /ortho, /noerase, pos=[0,0,1,1], /nobor
;MAP_GRID, LABEL=2, LATS=lats, LATNAMES=latnames,  $
;LATLAB=7,latdel=10, londel=10,LONLAB=-2.5,  LONS=-15, $
;GLINESTYLE=0, color=c0;12
;endif


; posmotret otdelno, pravilno li otcentrovan disk solnca na 2D
window,1,xs=800,ys=800
a=[0,0,2400,2400]
hmi_x = (findgen(n1)-n1/2.0)*cdelt1
hmi_y = (findgen(n2)-n2/2.0)*cdelt2
xr1=[a[0]-a[2]*0.5,a[0]+a[2]*0.5]
yr1=[a[1]-a[3]*0.5,a[1]+a[3]*0.5]
aspect=float(a[2])/a[3]
aia_lct, rgbr, rgbg, rgbb, wavelnth=wavelength1, /load
gvi_multiplot, [1,1], asp = aspect
tvplot, saia^pow,charsize=chrsz, $
hmi_x,hmi_y,xra=xr1,yra=yr1, miss = 0,title='SUN '+date_rat+' '+strmid(time_rat,0,8)
xyouts,20,!d.y_size-20,file_basename(f,'.fts.gz'),/device,charsize=chrsz
rads=[-solar_r,0,solar_r]
for jj=0,2 do plots,rads[jj],!y.crange, color=colorr,linest=1,thick=linethick

;posmotret' otdelno, pravilno li otcentrovan skan- sovpadajut li limby s raschetnymi
window,2,xs=1200,ys=800
loadct,12
plot,xmas,rat[*,0,n3r-2],xra=xr1,/noerase,xst=1,yst=1,title=file_basename(fr)
for jj=0,2 do plots,rads[jj],!y.crange, color=red, $
linest=1,thick=linethick

goto,kon
err:
print, 'file not found'
kon:
end