c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c       cccc       parameters of fit3 for S                       ccc
c       ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
        a( 0 )  =  0.03471d0
        a( 1 )  = -0.01341d0
        a( 2 )  = -0.00098d0
        a( 3 )  = -0.00204d0
        a( 4 )  = -0.00093d0
        b( 1 )  = -0.01107d0
        b( 2 )  = -0.00280d0
        b( 3 )  = -0.00100d0
        b( 4 )  =  0.d0
        c       =  0.00304d0
        d       =  0.13847d0
        e( 0 )  =  0.07593d0
        e( 1 )  = -0.02568d0
        e( 2 )  = -0.00575d0
        e( 3 )  = -0.00413d0
        e( 4 )  = -0.00240d0
        f( 1 )  = -0.01332d0
        f( 2 )  = -0.00495d0
        f( 3 )  = -0.00177d0
        f( 4 )  =  0.d0
        g       = -0.00913d0
        h       =  0.28946d0
        i( 0 )  =  0.00662d0
        i( 1 )  = -0.00432d0
        i( 2 )  =  0.00096d0
        i( 3 )  = -0.00004d0
        i( 4 )  =  0.00009d0
        j( 1 )  =  0.00839d0
        j( 2 )  = -0.00293d0
        j( 3 )  =  0.00068d0
        j( 4 )  =  0.d0
        k       = -0.00445d0
        l       =  0.02800d0
        p( 0 )  =  0.00848d0
        p( 1 )  = -0.00534d0
        p( 2 )  =  0.00091d0
        p( 3 )  =  0.00008d0
        p( 4 )  =  0.00011d0
        q( 1 )  =  0.01632d0
        q( 2 )  = -0.00506d0
        q( 3 )  =  0.00111d0
        q( 4 )  =  0.d0
        r       = -0.00928d0
        s       =  0.05837d0
        ap( 0 ) =  0.01927d0
        ap( 1 ) = -0.00567d0
        ap( 2 ) = -0.00143d0
        ap( 3 ) = -0.00091d0
        ap( 4 ) = -0.00041d0
        bp( 1 ) = -0.00365d0
        bp( 2 ) = -0.00162d0
        bp( 3 ) = -0.00045d0
        bp( 4 ) =  0.d0
        cp      = -0.00113d0
        dp      =  0.09407d0
        ip( 0 ) =  0.00277d0
        ip( 1 ) = -0.00137d0
        ip( 2 ) =  0.00008d0
        ip( 3 ) =  0.00009d0
        ip( 4 ) = -0.00003d0
        jp( 1 ) =  0.00535d0
        jp( 2 ) = -0.00159d0
        jp( 3 ) =  0.00034d0
        jp( 4 ) =  0.d0
        kp      = -0.00301d0
        lp      =  0.01896d0
        alpha0  =    0.5962d0
        alpha1  =   14.4670d0
        alpha2  = -144.0690d0
        alpha3  =  251.8580d0
        beta0   =    0.6814d0
        beta1   =   12.3005d0
        beta2   = -130.5920d0
        beta3   =  229.2470d0
        alphap0 =    0.2617d0
        alphap1 =  -12.4639d0
        alphap2 =  177.6570d0
        alphap3 = -475.7000d0
        betap0  =    0.3088d0
        betap1  =  -13.3769d0
        betap2  =  184.6140d0
        betap3  = -494.2780d0