Read params Open input SPECTRAL file, extract EXPOSURE, BACKSCAL. Open input template IMAGE file, extract keywords and 2D array. Configure lines and gaussian normalization (gnorm) values. Open and read RMF. Open and read ARF. Loop through selected lines chosen from OVII OVIII OVIIIB NEIX MGXI. call CAL_swcxmap to get selected swcxMAP corresponding to that line. rnorm = gnorms(line) * area (from ARF) * exposure (from SPEC) * backscal calculate fractional area (frac) from RMF MATRIX rnorm = frac * rnorm / total_swcx_exposure do i=1,780 do j=1,780 outSWCXmap(i,j) = outSWCXmap(i,j) + rnorm*swcxMAP(i,j) enddo enddo endLoop Mask chips if desired (ccd not selected) Create output fits image based on template (copying Attributes).