        real*8 FUNCTION F2_4(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  -1.415659669 ,
     2  -1.388699222 ,
     3  -1.321521607 ,
     4  -1.154824946 ,
     5  -0.7486276612 ,
     6  0.1168921239 ,
     7  1.433061347 ,
     8  1.606503951 ,
     9  -6.606066767 ,
     1  -50.75927175 ,
     2  -232.4794908 ,
     3  -895.6928536 ,
     4  -3175.193633 ,
     5  -10657.17457 ,
     6  -34869.92645 ,
     7  -111018.9687 /

        FT1 = 10.2343709
        FTN = -19.36216567
        FF = splint2(T,F,16,FT1,FTN,TT)
        F2_4 = FF
        RETURN
        END

        real*8 FUNCTION Fnn2_4(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  0.1603137234 ,
     2  0.2931852518 ,
     3  0.6426912722 ,
     4  1.469293043 ,
     5  3.576845911 ,
     6  10.99276461 ,
     7  43.31196946 ,
     8  152.9790238 ,
     9  452.2447282 ,
     1  1265.652764 ,
     2  3500.103361 ,
     3  9661.365475 ,
     4  26861.72288 ,
     5  56615.36657 ,
     6  227932.9797 ,
     7  560483.1322 /

        FT1 = 47.88926775
        FTN = 56.47796204
        FF = splint2(T,F,16,FT1,FTN,TT)
        Fnn2_4 = FF
        RETURN
        END

        real*8 FUNCTION F2(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  -1.413264692 ,
     2  -1.384679161 ,
     3  -1.313922475 ,
     4  -1.139016974 ,
     5  -0.7121623638 ,
     6  0.2235483621 ,
     7  1.834195213 ,
     8  3.038916337 ,
     9  -2.345274036 ,
     1  -38.79151123 ,
     2  -199.2823632 ,
     3  -804.2381425 ,
     4  -2921.96116 ,
     5  -10024.61436 ,
     6  -32834.52079 ,
     7  -105720.069 /

        FT1 = 10.89720317
        FTN = -18.53351751
        FF = splint2(T,F,16,FT1,FTN,TT)
        F2 = FF
        RETURN
        END

        real*8 FUNCTION F1_6(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  -1.404299206 ,
     2  -1.372861682 ,
     3  -1.296463277 ,
     4  -1.109861671 ,
     5  -0.6558917758 ,
     6  0.3653942173 ,
     7  2.315088157 ,
     8  4.750867433 ,
     9  2.836231281 ,
     1  -24.15309737 ,
     2  -158.6927255 ,
     3  -691.8698355 ,
     4  -2613.046201 ,
     5  -9160.55595 ,
     6  -30522.92861 ,
     7  -99225.23593 /

        FT1 = 12.15510022
        FTN = -17.56884471
        FF = splint2(T,F,16,FT1,FTN,TT)
        F1_6 = FF
        RETURN
        END

        real*8 FUNCTION F1_2(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  -1.359598487 ,
     2  -1.323541022 ,
     3  -1.237997904 ,
     4  -1.032851576 ,
     5  -0.5408238569 ,
     6  0.5921028145 ,
     7  2.948611673 ,
     8  6.908309762 ,
     9  9.449913979 ,
     1  -5.340088658 ,
     2  -106.3601915 ,
     3  -547.2504468 ,
     4  -2213.386607 ,
     5  -8052.606973 ,
     6  -27544.19479 ,
     7  -90833.20984 /

        FT1 = 14.13113252
        FTN = -16.24515961
        FF = splint2(T,F,16,FT1,FTN,TT)
        F1_2 = FF
        RETURN
        END

        real*8 FUNCTION F0_8(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  -1.071449432 ,
     2  -1.031952297 ,
     3  -0.9331049378 ,
     4  -0.6983659739 ,
     5  -0.1488085929 ,
     6  1.14131297 ,
     7  3.992316407 ,
     8  9.898560158 ,
     9  18.61522636 ,
     1  21.04149353 ,
     2  -32.67761918 ,
     3  -342.657355 ,
     4  -1649.781364 ,
     5  -6502.401136 ,
     6  -23257.16155 ,
     7  -78953.42378 /

        FT1 = 13.93093657
        FTN = -14.41716403
        FF = splint2(T,F,16,FT1,FTN,TT)
        F0_8 = FF
        RETURN
        END

        real*8 FUNCTION Fnn0_8(TT)

        implicit none
        real*8 TT,FF,FT1,FTN
        real*8 F(1:16)
        real*8 T(1:16)
        real*8 splint2

        data T /
     1  0.001583407607 ,
     2  0.004369348882 ,
     3  0.01205704051 ,
     4  0.03327091284 ,
     5  0.09180973061 ,
     6  0.2533452171 ,
     7  0.6990958211 ,
     8  1.929126481 ,
     9  5.323346054 ,
     1  14.6895569 ,
     2  40.53523473 ,
     3  111.8553314 ,
     4  308.660237 ,
     5  851.7353686 ,
     6  2350.329104 ,
     7  6485.637557 /
        data F /
     1  0.9462162342 ,
     2  0.959510795 ,
     3  0.9792052824 ,
     4  1.029070751 ,
     5  1.141488422 ,
     6  1.407566303 ,
     7  2.087434816 ,
     8  5.311714724 ,
     9  15.96996022 ,
     1  46.38385085 ,
     2  129.6894653 ,
     3  359.5069378 ,
     4  994.1065988 ,
     5  2770.065146 ,
     6  7561.872932 ,
     7  20959.65117 /

        FT1 = 5.078858229
        FTN = 3.269867863
        FF = splint2(T,F,16,FT1,FTN,TT)
        Fnn0_8 = FF
        RETURN
        END

        real*8 FUNCTION FALL(i,TT)

        implicit none
        integer i
        real*8 TT
        real*8 F2_4
        real*8 F2
        real*8 F1_6
        real*8 F1_2
        real*8 F0_8
        if (i.eq.1) FALL = F2_4(TT)
        if (i.eq.2) FALL = F2(TT)
        if (i.eq.3) FALL = F1_6(TT)
        if (i.eq.4) FALL = F1_2(TT)
        if (i.eq.5) FALL = F0_8(TT)

        RETURN
        END

        real*8 FUNCTION FnnALL(i,TT)

        implicit none
        integer i
        real*8 TT
        real*8 Fnn2_4
        real*8 Fnn0_8
        if (i.eq.1) FnnALL = Fnn2_4(TT)
        if (i.eq.5) FnnALL = Fnn0_8(TT)

        RETURN
        END

        real*8 FUNCTION nnALL(i)

        implicit none
        integer i
        real*8 nn
        if (i.eq.1) nn= 0.01726941657d0
        if (i.eq.2) nn= 0.02984155183d0
        if (i.eq.3) nn= 0.05828428092d0
        if (i.eq.4) nn= 0.1381553325d0
        if (i.eq.5) nn= 0.4662742473d0
        nnAll = nn

        RETURN
        END
