# IMAGE PARAMETERS
A) galaxy.fits              # Input data image (FITS file)
B) result.fits       # Name for the output image
C) var.fits               # Sigma image name (made from data if blank or "none") 
D) psf.fits     # Input PSF image and (optional) diffusion kernel
E) 1                   # PSF oversampling factor relative to data
F) none           # Pixel mask (ASCII file or FITS file with non-0 values)
G) none         # Parameter constraint file (ASCII)
H) 1 79 1 79     # Image region to fit (xmin xmax ymin ymax)
I) 80 80             # Size of convolution box (x y)
J) 22.5              # Magnitude photometric zeropoint 
K) 0.05   0.05         # Plate scale (dx dy)  [arcsec/pix]
O) regular                # Display type (regular, curses, both)
P) 0                   # Create output image only? (1=yes; 0=optimize) 
S) 0                   # Modify/create objects interactively?

# Start a list of objects below here......

# Component number: 1
 0) sersic                 #  Component type
 1) 40.6   40.7   1 1  #  Position x, y
 3) 19.0     1          #  Integrated magnitude 
 4) 1         1          #  R_e (effective radius)   [pix]
 5) 4         1          #  Sersic index n (de Vaucouleurs n=4) 
 6) 0.0000      0          #     ----- 
 7) 0.0000      0          #     ----- 
 8) 0.0000      0          #     ----- 
 9) 1.0000      0          #  Axis ratio (b/a)  
10) 71      1          #  Position angle (PA) [deg: Up=0, Left=90]
 Z) 0                      #  Skip this model in output image?  (yes=1, no=0)

# Component number: 1
 0) sersic                 #  Component type
 1) 40.6   40.7   1 1  #  Position x, y
 3) 18.0     1          #  Integrated magnitude 
 4) 8          1          #  R_e (effective radius)   [pix]
 5) 1         1          #  Sersic index n (de Vaucouleurs n=4) 
 6) 0.0000      0          #     ----- 
 7) 0.0000      0          #     ----- 
 8) 0.0000      0          #     ----- 
 9) 1      1          #  Axis ratio (b/a)  
10) 75      1          #  Position angle (PA) [deg: Up=0, Left=90]
 Z) 0                      #  Skip this model in output image?  (yes=1, no=0)

# Component number: 1
# 0) sersic                 #  Component type
# 1) 40.6   40.7   1 1  #  Position x, y
# 3) 16.0     1          #  Integrated magnitude 
# 4) 10          1          #  R_e (effective radius)   [pix]
# 5) 2         1          #  Sersic index n (de Vaucouleurs n=4) 
# 6) 0.0000      0          #     ----- 
# 7) 0.0000      0          #     ----- 
# 8) 0.0000      0          #     ----- 
# 9) 0.4      1          #  Axis ratio (b/a)  
#10) 75      1          #  Position angle (PA) [deg: Up=0, Left=90]
# Z) 0                      #  Skip this model in output image?  (yes=1, no=0)


# Component number: 2
 0) expdisk                #  Component type
 1) 40.6   40.7   1 1  #  Position x, y
 3) 16      1          #  Integrated magnitude 
 4) 10      1          #  R_s (disk scale-length)   [pix]
 5) 0.0000      0          #     ----- 
6) 0.0000      0          #     ----- 
 7) 0.0000      0          #     ----- 
8) 0.0000      0          #     ----- 
 9) 0.4      1          #  Axis ratio (b/a)  
10) 75      1          #  Position angle (PA) [deg: Up=0, Left=90]
 Z) 0                      #  Skip this model in output image?  (yes=1, no=0)

 0) sky
 1) 0.000       1       # sky background       [ADU counts]
 2) 0.000      0       # dsky/dx (sky gradient in x) 
 3) 0.000      0       # dsky/dy (sky gradient in y) 
 Z) 0                  #  Skip this model in output image?  (yes=1, no=0)