pro rf, cme_name

print, 'cme_name: ', cme_name

cme_img = rfits (cme_name, header=cme_hdu, /scale)

print, cme_hdu

END