CCDPROC the dome flats (cross-talk, overscan, trim, fixpix, bias frame, dark frame)
 (You might have already done some of these steps in your initial round of processing, but would
still need to apply the bias frame correction and the dark frame correction if you need to apply it -- i.e.
if there is a measurable dark current. We are currently not applying a dark correction to your images).

      images = "@dflat.list"         List of Mosaic CCD images to process
     (output = "")                   List of output processed images
    (ccdtype = "flat")               CCD image type to process
     (noproc = no)                   List processing steps only?\n
   (xtalkcor = yes)                  Apply crosstalk correction?
   (overscan = yes)                  Apply overscan strip correction?
       (trim = yes)                  Trim the image?
     (fixpix = yes)                   Apply bad pixel mask correction?
    (zerocor = yes)                  Apply zero level correction?
    (darkcor = no)                   Apply dark count correction?
    (flatcor = no)                   Apply flat field correction?
   (sflatcor = no)                   Apply sky flat field correction?
      (merge = no)                   Merge amplifiers from same CCD?\n
  (xtalkfile = "mscdb$noao/CCDMosaThin1/xtalk.dat") Crosstalk file
    (biassec = "!biassec")           Overscan strip image section
    (trimsec = "!trimsec")           Trim data section
    (fixfile = "BPM")                List of bad pixel masks
       (zero = "Zero990413")         List of zero level calibration images
       (dark = "Dark*")              List of dark count calibration images
       (flat = "Dflat*")             List of flat field images
      (sflat = "Sflat*")             List of secondary flat field images
 (minreplace = 1.)                   Minimum flat field value\n
(interactive = no)                   Fit overscan interactively?
   (function = "legendre")           Fitting function
      (order = 1)                    Number of polynomial terms or spline pieces
     (sample = "*")                  Sample points to fit
   (naverage = 1)                    Number of sample points to combine
   (niterate = 1)                    Number of rejection iterations
 (low_reject = 3.)                   Low sigma rejection factor
(high_reject = 3.)                   High sigma rejection factor
       (grow = 0.)                   Rejection growing radius
         (fd = "")
        (fd2 = "")
       (mode = "ql")
 

Combine the dome flats using flatcombine. (For MOSAIC-I-E, Chip4 should not go above 7000 ADU---non-linear region.)   If you had not already processed the images, flatcombine would check
and if necessary call ccdproc to process the images before combining them into your new domeflat.

Parameters for mscred task flatcombine

     input = "@DflatR.list"      List of flat field images to combine
   (output = "Dflat981011")      Output flat field root name
  (combine = "average")          Type of combine operation
   (reject = "ccdclip")          Type of rejection
  (ccdtype = "flat")             CCD image type to combine
  (process = yes)                Process images before combining?
  (subsets = yes)                Combine images by subset parameter?
   (delete = no)                 Delete input images after combining?
    (scale = "mode")             Image scaling
  (statsec = "")                 Image section for computing statistics
     (nlow = 1)                  minmax: Number of low pixels to reject
    (nhigh = 1)                  minmax: Number of high pixels to reject
    (nkeep = 1)                  Minimum to keep (pos) or maximum to reject (neg
    (mclip = no)                 Use median in sigma clipping algorithms?
   (lsigma = 3.)                 Lower sigma clipping factor
   (hsigma = 3.)                 Upper sigma clipping factor
  (rdnoise = "rdnoise")          ccdclip: CCD readout noise (electrons)
     (gain = "gain")             ccdclip: CCD gain (electrons/DN)
   (snoise = "0.")               ccdclip: Sensitivity noise (fraction)
    (pclip = -0.5)               pclip: Percentile clipping parameter
    (blank = 1.)                 Value if there are no pixels
     (mode = "ql")