c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc        parameters of fit2  for  Mg                   ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 ) =  0.25567d0
        a( 1 ) = -0.08158d0
        a( 2 ) = -0.01632d0
        a( 3 ) = -0.01667d0
        a( 4 ) = -0.00952d0
        a( 5 ) = -0.00748d0
        b( 1 ) = -0.07491d0
        b( 2 ) = -0.02220d0
        b( 3 ) = -0.00653d0
        b( 4 ) = -0.00345d0
        b( 5 ) =  0.d0
        c      =  0.00899d0
        d      =  0.38006d0
        e( 0 ) =  0.13455d0
        e( 1 ) = -0.02828d0
        e( 2 ) = -0.01846d0
        e( 3 ) = -0.01087d0
        e( 4 ) = -0.00693d0
        e( 5 ) = -0.00547d0
        f( 1 ) = -0.04508d0
        f( 2 ) = -0.01782d0
        f( 3 ) = -0.00607d0
        f( 4 ) = -0.00200d0
        f( 5 ) =  0.d0
        g      = -0.02832d0
        h      =  0.31541d0
        i( 0 ) =  0.01231d0
        i( 1 ) = -0.01030d0
        i( 2 ) =  0.00016d0
        i( 3 ) =  0.00259d0
        i( 4 ) =  0.00048d0
        i( 5 ) =  0.00183d0
        j( 1 ) =  0.02669d0
        j( 2 ) = -0.01291d0
        j( 3 ) =  0.00154d0
        j( 4 ) = -0.00108d0
        j( 5 ) =  0.d0
        k      = -0.01386d0
        l      =  0.08711d0
        p( 0)  = -0.00583d0
        p( 1 ) =  0.00192d0
        p( 2 ) = -0.00179d0
        p( 3 ) =  0.00193d0
        p( 4 ) =  0.00034d0
        p( 5 ) =  0.00102d0
        q( 1 ) =  0.01501d0
        q( 2 ) = -0.00582d0
        q( 3 ) =  0.00024d0
        q( 4 ) = -0.00038d0
        q( 5 ) =  0.d0
        r      = -0.01059d0
        s      =  0.06657d0
        alpha0 = -0.07003d0
        alpha1 =  0.07808d0
        alpha2 =  1.75870d0
        alpha3 = -0.76675d0
        beta0  = -0.06960d0
        beta1  =  0.06577d0
        beta2  =  1.81180d0
        beta3  = -0.80797d0