c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc       parameters of fit3 for Ca                      ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 )  =  0.03754d0
        a( 1 )  = -0.01521d0
        a( 2 )  = -0.00052d0
        a( 3 )  = -0.00218d0
        a( 4 )  = -0.00086d0
        b( 1 )  = -0.01349d0
        b( 2 )  = -0.00280d0
        b( 3 )  = -0.00115d0
        b( 4 )  =  0.d0
        c       =  0.00531d0
        d       =  0.13927d0
        e( 0 )  =  0.08106d0
        e( 1 )  = -0.02720d0
        e( 2 )  = -0.00613d0
        e( 3 )  = -0.00460d0
        e( 4 )  = -0.00260d0
        f( 1 )  = -0.01665d0
        f( 2 )  = -0.00559d0
        f( 3 )  = -0.00229d0
        f( 4 )  =  0.d0
        g       = -0.00698d0
        h       =  0.31671d0
        i( 0 )  =  0.00789d0
        i( 1 )  = -0.00520d0
        i( 2 )  =  0.00123d0
        i( 3 )  = -0.00010d0
        i( 4 )  =  0.00012d0
        j( 1 )  =  0.00865d0
        j( 2 )  = -0.00312d0
        j( 3 )  =  0.00073d0
        j( 4 )  =  0.d0
        k       = -0.00448d0
        l       =  0.02820d0
        p( 0 )  =  0.01025d0
        p( 1 )  = -0.00652d0
        p( 2 )  =  0.00123d0
        p( 3 )  =  0.00001d0
        p( 4 )  =  0.00016d0
        q( 1 )  =  0.01790d0
        q( 2 )  = -0.00556d0
        q( 3 )  =  0.00119d0
        q( 4 )  =  0.d0
        r       = -0.01015d0
        s       =  0.06388d0
        ap( 0 ) =  0.02137d0
        ap( 1 ) = -0.00661d0
        ap( 2 ) = -0.00134d0
        ap( 3 ) = -0.00103d0
        ap( 4 ) = -0.00043d0
        bp( 1 ) = -0.00477d0
        bp( 2 ) = -0.00175d0
        bp( 3 ) = -0.00054d0
        bp( 4 ) =  0.d0
        cp      = -0.00029d0
        dp      =  0.10004d0
        ip( 0 ) =  0.00343d0
        ip( 1 ) = -0.00178d0
        ip( 2 ) =  0.00018d0
        ip( 3 ) =  0.00007d0
        ip( 4 ) = -0.00003d0
        jp( 1 ) =  0.00579d0
        jp( 2 ) = -0.00176d0
        jp( 3 ) =  0.00038d0
        jp( 4 ) =  0.d0
        kp      = -0.00321d0
        lp      =  0.02018d0
        alpha0  =    0.5793d0
        alpha1  =   15.1152d0
        alpha2  = -151.6510d0
        alpha3  =  276.8910d0
        beta0   =    0.7331d0
        beta1   =   11.2258d0
        beta2   = -127.8900d0
        beta3   =  238.5240d0
        alphap0 =    0.2427d0
        alphap1 =  -12.2010d0
        alphap2 =  180.4920d0
        alphap3 = -496.3580d0
        betap0  =    0.3205d0
        betap1  =  -13.7451d0
        betap2  =  194.3840d0
        betap3  = -539.2810d0