def idlfitstocasaim(imname, reffreq, freqincrement):
    """One of IDL's many crimes is claiming to write FITS files when it doesn't.
    This function (based on what Kumar Golap suggested) tries to get an IDL
    'fits' file into casa.

    imname - this function copies imname.fits in the current directory to
             imname.im
    reffreq - Reference frequency in Hz.
    freqincrement - Frequency increment in Hz.
    """

    fitsfile = imname + '.fits'
    casaname = imname + '.im'
    
    ia.fromfits(outfile = 'temp.im', infile = fitsfile)

    ### add stokes axis
    ib = ia.adddegaxes(stokes='I')

    #### add spectral axis
    ic = ib.adddegaxes(outfile = casaname, spectral=True)
    ic.summary()

    ###Now change the reference frequency and increments
    csys = ic.coordsys()
    a = csys.referencevalue()
    a['numeric'][3] = reffreq
    csys.setreferencevalue(a)
    a = csys.increment()
    a['numeric'][3] = freqincrement
    csys.setincrement(a)
    ic.setcoordsys(csys=csys.torecord())
    ic.summary()