pro niris_inv_script

;read the sample.fts file
img=readfits(dp(),h)
;run the ME inversion
out=main_test(img,h)
help, out
window, xs=800, ys=600
;B_total
tv, congrid(bytscl(out[*,*,0], 0, 3000), 200, 200), 0
loadct, 13
;Inclination
tv, congrid(bytscl(out[*,*,1], 0, !pi), 200, 200), 1
;Azimuth
tv, congrid(bytscl(out[*,*,2], 0, !pi), 200, 200), 2
;Doppler
tv, congrid(bytscl(out[*,*,6], -0.2,0.2), 200, 200), 3
;read original.fts
o=readfits(dp())
;read fit_array.fts
f=readfits(dp())
!p.multi=[0,2,2]
x=20 & y=20
for i=0, 3 do begin & plot, o[x,y,*,i] & oplot, f[x,y,*,i], linest=2 & endfor
x=(findgen(50)-25.)/10.
x=[x,x,x,x]
help, x
;X              FLOAT     = Array[200]
;create a model
par=[1000, !pi/6, !pi/3, 1.3, 0.15,1., 0., 0.6, 0.2]
niris_MEsinglet, x, par, f
help, x, f
;X              FLOAT     = Array[200]
;F               DOUBLE    = Array[200]
for i=0, 3 do begin & plot, x[0:49], f[i*50:i*50+49] & endfor
stop
;type .cont
par=[1000, !pi/6, !pi/3, 5, 0.15,1., 0., 0.6, 0.2]
niris_MEsinglet, xx, par, f
for i=0, 3 do begin & plot, x[0:49], f[i*50:i*50+49] & endfor
stop
;type .cont
par=[1000, !pi/6, !pi/3, 5, 0.15,1., 0., 0.6, 0.4]
niris_MEsinglet, xx, par, f
for i=0, 3 do begin & plot, xx[0:49], f[i*50:i*50+49] & endfor


end