c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc        parameters of fit2  for Fe                    ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 ) =  0.34888d0
        a( 1 ) = -0.11076d0
        a( 2 ) = -0.02349d0
        a( 3 ) = -0.02283d0
        a( 4 ) = -0.01250d0
        a( 5 ) = -0.00971d0
        b( 1 ) = -0.10661d0
        b( 2 ) = -0.02860d0
        b( 3 ) = -0.00785d0
        b( 4 ) = -0.00385d0
        b( 5 ) =  0.d0
        c      =  0.00766d0
        d      =  0.40991d0
        e( 0 ) =  0.23159d0
        e( 1 ) = -0.05891d0
        e( 2 ) = -0.02531d0
        e( 3 ) = -0.01747d0
        e( 4 ) = -0.01021d0
        e( 5 ) = -0.00778d0
        f( 1 ) = -0.07767d0
        f( 2 ) = -0.02516d0
        f( 3 ) = -0.00711d0
        f( 4 ) = -0.00260d0
        f( 5 ) =  0.d0
        g      = -0.03076d0
        h      =  0.36908d0
        i( 0 ) =  0.01892d0
        i( 1 ) = -0.01493d0
        i( 2 ) =  0.00125d0
        i( 3 ) =  0.00262d0
        i( 4 ) =  0.00055d0
        i( 5 ) =  0.00209d0
        j( 1 ) =  0.03034d0
        j( 2 ) = -0.01519d0
        j( 3 ) =  0.00204d0
        j( 4 ) = -0.00135d0
        j( 5 ) =  0.d0
        k      = -0.01494d0
        l      =  0.09395d0
        p( 0)  = -0.00011d0
        p( 1 ) = -0.00222d0
        p( 2 ) = -0.00104d0
        p( 3 ) =  0.00225d0
        p( 4 ) =  0.00037d0
        p( 5 ) =  0.00140d0
        q( 1 ) =  0.02013d0
        q( 2 ) = -0.00883d0
        q( 3 ) =  0.00090d0
        q( 4 ) = -0.00071d0
        q( 5 ) =  0.d0
        r      = -0.01249d0
        s      =  0.07850d0
        alpha0 = -0.07608d0
        alpha1 =  0.11559d0
        alpha2 =  1.75730d0
        alpha3 = -0.79677d0
        beta0  = -0.08034d0
        beta1  =  0.14368d0
        beta2  =  1.73140d0
        beta3  = -0.79467d0