load_image("image2.fits")
image_data()

set_coord("physical")
set_stat("cash")
set_method("simplex")
show_all()

notice2d("circle(4072.46,4249.34,108)")

set_source(gauss2d.g1)
g1.ampl = 20
g1.fwhm = 20
g1.xpos = 4065.5
g1.ypos = 4250.5
g1.fwhm.max = 4300
g1.xpos.max = 4300
g1.ypos.max = 4300
g1.ampl.min = 1
g1.ampl.max = 1000
show_model()

set_source(g1+const2d.bgnd)
bgnd.c0 = 0.2
bgnd.c0.max = 100
show_model()

fit()
image_resid()

save_model("model.fits")
save_resid("resid.fits")

conf()

freeze(g1)
set_source(gauss2d.g2+g1+bgnd)
g2.fwhm = 10
g2.ampl = 100
g2.xpos = 4065.5
g2.ypos = 4250.5
g2.fwhm.max = 4300
g2.xpos.max = 4300
g2.ypos.max = 4300
g2.ampl.min = 1
g2.ampl.max = 1000

fit()
show_model()

thaw(g2.ellip)
thaw(g2.theta)
g2.ellip = 0.1
g2.theta = 1

fit()

conf(g2.ellip, g2.theta)

image_source_component(g2)

show_all()