c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc       parameters of fit3 for Mg                      ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 )  =  0.03191d0
        a( 1 )  = -0.01083d0
        a( 2 )  = -0.00211d0
        a( 3 )  = -0.00192d0
        a( 4 )  = -0.00109d0
        b( 1 )  = -0.00753d0
        b( 2 )  = -0.00281d0
        b( 3 )  = -0.00087d0
        b( 4 )  =  0.d0
        c       = -0.00110d0
        d       =  0.14075d0
        e( 0 )  =  0.06145d0
        e( 1 )  = -0.01751d0
        e( 2 )  = -0.00750d0
        e( 3 )  = -0.00340d0
        e( 4 )  = -0.00231d0
        f( 1 )  = -0.00322d0
        f( 2 )  = -0.00386d0
        f( 3 )  = -0.00148d0
        f( 4 )  =  0.d0
        g       = -0.01798d0
        h       =  0.29179d0
        i( 0 )  =  0.00468d0
        i( 1 )  = -0.00297d0
        i( 2 )  =  0.00053d0
        i( 3 )  =  0.00004d0
        i( 4 )  =  0.00005d0
        j( 1 )  =  0.00810d0
        j( 2 )  = -0.00261d0
        j( 3 )  =  0.00059d0
        j( 4 )  =  0.d0
        k       = -0.00451d0
        l       =  0.02840d0
        p( 0 )  =  0.00384d0
        p( 1 )  = -0.00219d0
        p( 2 )  =  0.00005d0
        p( 3 )  =  0.00017d0
        p( 4 )  =  0.00004d0
        q( 1 )  =  0.01526d0
        q( 2 )  = -0.00412d0
        q( 3 )  =  0.00086d0
        q( 4 )  =  0.d0
        r       = -0.00933d0
        s       =  0.05868d0
        ap( 0 ) =  0.01608d0
        ap( 1 ) = -0.00329d0
        ap( 2 ) = -0.00222d0
        ap( 3 ) = -0.00078d0
        ap( 4 ) = -0.00042d0
        bp( 1 ) = -0.00108d0
        bp( 2 ) = -0.00138d0
        bp( 3 ) = -0.00044d0 
        bp( 4 ) =  0.d0
        cp      = -0.00399d0
        dp      =  0.09525d0
        ip( 0 ) =  0.00137d0
        ip( 1 ) = -0.00044d0
        ip( 2 ) = -0.00015d0
        ip( 3 ) =  0.00011d0
        ip( 4 ) = -0.00003d0
        jp( 1 ) =  0.00500d0
        jp( 2 ) = -0.00129d0
        jp( 3 ) =  0.00025d0
        jp( 4 ) =  0.d0
        kp      = -0.00304d0
        lp      =  0.01915d0
        alpha0  =    0.5104d0
        alpha1  =   16.1851d0
        alpha2  = -147.1310d0
        alpha3  =  230.5050d0
        beta0   =    0.5855d0
        beta1   =   14.5626d0
        beta2   = -141.7690d0
        beta3   =  237.9640d0
        alphap0 =    0.2524d0
        alphap1 =  -12.4235d0
        alphap2 =  172.5400d0
        alphap3 = -446.9220d0
        betap0  =    0.2521d0
        betap1  =  -12.0769d0
        betap2  =  170.6080d0
        betap3  = -446.8310d0