;
pro getdata

restore,'mas-cr2068.sav'
restore,'pol_ch.sav'

save,filename='temp.sav'

end

;*********************************************
pro changedata

restore,'temp.sav'

chbound_mas=fltarr(360,180)

polch=congrid(pol_ch,360,180)

for i = 0,179 do begin
 chbound_mas(*,i) = polch(*,179-i)
endfor

test=where(chbound_mas eq 50.)
chbound_mas(test) = -1.
test=where(chbound_mas eq 250.)
chbound_mas(test) = 1.
test=where(chbound_mas eq 255.)
chbound_mas(test) = 0.


;brph=fltarr(100,220)
;brph(*,*)=br_corona(0,*,*)
;brph=transpose(brph)
;brph_mas=congrid(brph,360,180)
;sbr=brph_mas/abs(brph_mas)
;chbound_mas=sbr*chbound_mas

vr1au=fltarr(121,128)
for i = 0,120 do vr1au(i,*)=481.*vr_helio(ivr,120-i,*)
vr1au_mas=congrid(transpose(vr1au),360,180)

ibr=100
br1au=fltarr(120,128)
for i = 0,119 do br1au(i,*)=br_helio(ibr,119-i,*)
br1au_mas=congrid(transpose(br1au),360,180)

vr1au_mas_shift=fltarr(360,180)
br1au_mas_shift=fltarr(360,180)
for i = 45,359 do begin
  vr1au_mas_shift(i,*)=vr1au_mas(i-45,*)
  br1au_mas_shift(i,*)=br1au_mas(i-45,*)
endfor
for i = 0,44 do begin
  vr1au_mas_shift(i,*)=vr1au_mas(315+i,*)
  br1au_mas_shift(i,*)=br1au_mas(315+i,*)
endfor
 

ibr=118
brss=fltarr(100,220)
for i = 0,99 do brss(i,*)=br_corona(ibr,99-i,*)
brss_mas=transpose(brss)

temp=fltarr(100)
eventh=(findgen(100)*180./99.)*!dtor
for i = 0,219 do begin
 temp(*)=interpol(brss_mas(i,*),theta_brc,eventh,/lsquadratic)
 brss_mas(i,*)=temp(*)
endfor 
 
temp=fltarr(220)
evenph=(findgen(220)*360./219.)*!dtor
for i = 0,99 do begin
 temp(*)=interpol(brss_mas(*,i),phi_brc,evenph,/lsquadratic)
 brss_mas(*,i)=temp(*)
endfor 
 
 
brss_mas=congrid(brss_mas,360,180)
 

save,filename='sarahrevise_mas.sav',chbound_mas,vr1au_mas,br1au_mas,brss_mas,vr1au_mas_shift,br1au_mas_shift

stop
end