#!/usr/bin/env python3

import pyfits
import sys


if __name__ == "__main__":

    in_fn = sys.argv[1]

    filtername = sys.argv[2]
    
    out_fn = sys.argv[3]


    hdu = pyfits.open(in_fn)

    hdr = hdu[0].header
    hdr['AUTHOR'] = "Ralf Kotulla"
    hdr['ORIGIN'] = "Dept of Astronomy, U Wisconsin-Madison"
    hdr['EMAIL'] = 'kotulla@wisc.edu'
    hdr['FILTER'] = (sys.argv[2], "filter name")
    hdr['MAGZERO'] = (26.0, 'photometric zeropoint')

    
    for key in [
            'DATE',
            'SOFTNAME',
            'SOFTVERS',
            'SOFTDATE',
            'SOFTAUTH',
            'SOFTINST',
            'COMBINET',
            'RESAMPT1',
            'CENTERT1',
            'PSCALET1',
            'RESAMPT2',
            'CENTERT2',
            'PSCALET2',
            'EXPTIME', 'GAIN',
            ]:
        del hdr[key]

    try:
        while (True):
            del hdr['COMMENT']
    except KeyError:
        pass
    
    hdu.writeto(out_fn, clobber=True)
    
