imcombine beta*_spectra_?.fts  medflat combine=median

blkavg medflat[*,1:32] avcol 1 32
blkrep avcol avcol_out 1 32
imarith medflat / avcol_out masterflat calctype=real pixtype=real