load_image("center_box_0.25pix.fits");
image_data();

image_getregion();
notice2d("rotbox(88.16875,75.8625,70.416667,68.508334,0)");

load_psf("psf0", "psf_f1_norm_0.25pix.fits");

set_psf(psf0);
psf0.size=[32,32];
psf0.center=[128,129];
print(get_psf());
image_psf();

set_psf(psf0);
psf0.center=[128,129];
psf0.size=[72,72];
print(get_psf());
image_psf();

set_source(const2d.c1+gauss2d.g1);
c1.c0=1;
g1.fwhm = 6;
g1.xpos = 90;
g1.ypos = 80;
g1.ampl = 100;
show_source();

show_model();

set_stat("cash");
freeze(c1.c0);
set_method("neldermead");
set_method_opt("iquad",0);
set_method_opt("finalsimplex",0);
fit();


thaw(c1.c0);
fit();
projection();

contour_resid();
limits(X_AXIS,55,120);
limits(Y_AXIS,45,110);
get_contour_zrange();


image_crate = pack_image();
set_imagevals(image_crate,get_model_image().y);
write_file(image_crate,"2dpsf_model_image.fits");

set_imagevals(image_crate, get_resid_image.y);
write_file(image_crate, "2dpsf_resid_cnts.fits");