C
C
SUBROUTINE GENCKWCFS(T, C, WDOT)
C
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER (I-N)
C
DIMENSION C(*), WDOT(*)
C
DIMENSION RKF( 34)
C
ALOGT = LOG(T)
C
CTOT = C(1)
& + C( 2)
& + C( 3)
& + C( 4)
& + C( 5)
& + C( 6)
& + C( 7)
& + C( 8)
& + C( 9)
C
RKF( 1) = 0.1860000000000000D+15*EXP(
& - ( 0.8449924820065498D+04 / T ))
& * C( 3)
& * C( 1)
RKF( 2) = 0.1480000000000000D+14*EXP(
& - ( 0.3422244715690613D+03 / T ))
& * C( 4)
& * C( 5)
RKF( 3) = 0.1820000000000000D+11*EXP(
& ALOGT
& - ( 0.4479114407300949D+04 / T ))
& * C( 6)
& * C( 4)
RKF( 4) = 0.8320000000000000D+10*EXP(
& ALOGT
& - ( 0.3497735407948494D+04 / T ))
& * C( 3)
& * C( 5)
RKF( 5) = 0.3390000000000000D+14*EXP(
& - ( 0.9235028019547463D+04 / T ))
& * C( 2)
& * C( 4)
RKF( 6) = 0.3160000000000000D+13*EXP(
& - ( 0.5535984098911285D+03 / T ))
& * C( 5)** 2
RKF( 7) = 0.9550000000000000D+14*EXP(
& - ( 0.1021640701889992D+05 / T ))
& * C( 2)
& * C( 3)
RKF( 8) = 0.2190000000000000D+14*EXP(
& - ( 0.2591847100853920D+04 / T ))
& * C( 6)
& * C( 5)
RKF( 9) = 0.1000000000000000D+14*EXP(
& - ( 0.9058883070945740D+03 / T ))
& * C( 8)
& * C( 5)
RKF( 10) = 0.2820000000000000D+14*EXP(
& - ( 0.1650226532757282D+05 / T ))
& * C( 2)
& * C( 7)
RKF( 11) = 0.5010000000000000D+14*EXP(
& - ( 0.5032712817192078D+03 / T ))
& * C( 7)
& * C( 4)
RKF( 12) = 0.6460000000000000D+14*EXP(
& - ( 0.2826371518135071D+05 / T ))
& * C( 5)
& * C( 1)
RKF( 13) = 0.2510000000000000D+15*EXP(
& - ( 0.9562154352664948D+03 / T ))
& * C( 7)
& * C( 3)
RKF( 14) = 0.1200000000000000D+14*EXP(
& - ( 0.2018117839694023D+05 / T ))
& * C( 5)** 2
RKF( 15) = 0.2510000000000000D+14*EXP(
& - ( 0.3522898972034454D+03 / T ))
& * C( 7)
& * C( 3)
RKF( 16) = 0.5500000000000000D+14*EXP(
& - ( 0.2908908008337021D+05 / T ))
& * C( 6)
& * C( 1)
RKF( 17) = 0.5010000000000000D+14*EXP(
& - ( 0.5032712817192078D+03 / T ))
& * C( 7)
& * C( 5)
RKF( 18) = 0.6310000000000000D+15*EXP(
& - ( 0.3717161686778069D+05 / T ))
& * C( 2)
& * C( 1)
RKF( 19) = 0.3980000000000000D+14*EXP(
& - ( 0.2145948745250702D+05 / T ))
& * C( 8)
& * C( 1)
RKF( 20) = 0.1000000000000000D+14*EXP(
& - ( 0.5032712817192078D+03 / T ))
& * C( 7)** 2
RKF( 21) = 0.1700000000000000D+13*EXP(
& - ( 0.1887267306447029D+04 / T ))
& * C( 8)
& * C( 3)
RKF( 22) = 0.7240000000000000D+12*EXP(
& - ( 0.9411172968149185D+04 / T ))
& * C( 7)
& * C( 6)
RKF( 23) = 0.2190000000000000D+17*EXP(
& - ( 0.5284348458051682D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 2)
RKF( 24) = 0.1410000000000000D+24*EXP(
& -0.2000000000000000D+01* ALOGT
& - ( 0.0000000000000000D+00 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 3)
& * C( 5)
RKF( 25) = 0.1660000000000000D+16*EXP(
& - ( -0.5032712817192078D+03 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 3)
& * C( 1)
RKF( 26) = 0.2290000000000000D+16*EXP(
& - ( 0.2310015183091164D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 7)
RKF( 27) = 0.1200000000000000D+18*EXP(
& - ( 0.2289884331822395D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 8)
RKF( 28) = 0.9120000000000000D+15*EXP(
& - ( -0.2551585398316383D+04 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 5)** 2
RKF( 29) = 0.1000000000000000D+17*EXP(
& - ( 0.0000000000000000D+00 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 4)
& * C( 3)
RKF( 30) = 0.7940000000000000D+20*EXP(
& -ALOGT
& - ( 0.5219929733991623D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 5)
RKF( 31) = 0.5130000000000000D+16*EXP(
& - ( 0.5787619739770889D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 1)
RKF( 32) = 0.4680000000000000D+16*EXP(
& -0.2800000000000000D+00* ALOGT
& - ( 0.0000000000000000D+00 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 4)** 2
RKF( 33) = 0.2190000000000000D+15*EXP(
& - ( 0.4831404304504395D+05 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 6)
RKF( 34) = 0.3020000000000000D+16*EXP(
& - ( 0.0000000000000000D+00 / T ))
& * ( CTOT
& + ( 0.4000000000000000D+00-1.D0)*C( 1)
& + ( 0.6500000000000000D+01-1.D0)*C( 2)
& )
& * C( 3)** 2
C
WDOT( 1)=0.0
WDOT( 2)=0.0
WDOT( 3)=0.0
WDOT( 4)=0.0
WDOT( 5)=0.0
WDOT( 6)=0.0
WDOT( 7)=0.0
WDOT( 8)=0.0
WDOT( 9)=0.0
C
WDOT( 3) = WDOT( 3) - RKF( 1)
WDOT( 4) = WDOT( 4) - RKF( 2)
WDOT( 6) = WDOT( 6) - RKF( 3)
WDOT( 3) = WDOT( 3) - RKF( 4)
WDOT( 2) = WDOT( 2) - RKF( 5)
WDOT( 5) = WDOT( 5) +
& ( -0.2000000000000000D+01 * RKF( 6) )
WDOT( 2) = WDOT( 2) - RKF( 7)
WDOT( 6) = WDOT( 6) - RKF( 8)
WDOT( 8) = WDOT( 8) - RKF( 9)
WDOT( 2) = WDOT( 2) - RKF( 10)
WDOT( 7) = WDOT( 7) - RKF( 11)
WDOT( 5) = WDOT( 5) - RKF( 12)
WDOT( 7) = WDOT( 7) - RKF( 13)
WDOT( 5) = WDOT( 5) +
& ( -0.2000000000000000D+01 * RKF( 14) )
WDOT( 7) = WDOT( 7) - RKF( 15)
WDOT( 6) = WDOT( 6) - RKF( 16)
WDOT( 7) = WDOT( 7) - RKF( 17)
WDOT( 2) = WDOT( 2) - RKF( 18)
WDOT( 8) = WDOT( 8) - RKF( 19)
WDOT( 7) = WDOT( 7) +
& ( -0.2000000000000000D+01 * RKF( 20) )
WDOT( 8) = WDOT( 8) - RKF( 21)
WDOT( 7) = WDOT( 7) - RKF( 22)
WDOT( 2) = WDOT( 2) - RKF( 23)
WDOT( 3) = WDOT( 3) - RKF( 24)
WDOT( 3) = WDOT( 3) - RKF( 25)
WDOT( 7) = WDOT( 7) - RKF( 26)
WDOT( 8) = WDOT( 8) - RKF( 27)
WDOT( 5) = WDOT( 5) +
& ( -0.2000000000000000D+01 * RKF( 28) )
WDOT( 4) = WDOT( 4) - RKF( 29)
WDOT( 5) = WDOT( 5) - RKF( 30)
WDOT( 1) = WDOT( 1) - RKF( 31)
WDOT( 4) = WDOT( 4) +
& ( -0.2000000000000000D+01 * RKF( 32) )
WDOT( 6) = WDOT( 6) - RKF( 33)
WDOT( 3) = WDOT( 3) +
& ( -0.2000000000000000D+01 * RKF( 34) )
WDOT( 1) = WDOT( 1) - RKF( 1)
WDOT( 5) = WDOT( 5) - RKF( 2)
WDOT( 4) = WDOT( 4) - RKF( 3)
WDOT( 5) = WDOT( 5) - RKF( 4)
WDOT( 4) = WDOT( 4) - RKF( 5)
WDOT( 3) = WDOT( 3) - RKF( 7)
WDOT( 5) = WDOT( 5) - RKF( 8)
WDOT( 5) = WDOT( 5) - RKF( 9)
WDOT( 7) = WDOT( 7) - RKF( 10)
WDOT( 4) = WDOT( 4) - RKF( 11)
WDOT( 1) = WDOT( 1) - RKF( 12)
WDOT( 3) = WDOT( 3) - RKF( 13)
WDOT( 3) = WDOT( 3) - RKF( 15)
WDOT( 1) = WDOT( 1) - RKF( 16)
WDOT( 5) = WDOT( 5) - RKF( 17)
WDOT( 1) = WDOT( 1) - RKF( 18)
WDOT( 1) = WDOT( 1) - RKF( 19)
WDOT( 3) = WDOT( 3) - RKF( 21)
WDOT( 6) = WDOT( 6) - RKF( 22)
WDOT( 5) = WDOT( 5) - RKF( 24)
WDOT( 1) = WDOT( 1) - RKF( 25)
WDOT( 3) = WDOT( 3) - RKF( 29)
WDOT( 4) = WDOT( 4) + RKF( 1)
WDOT( 3) = WDOT( 3) + RKF( 2)
WDOT( 3) = WDOT( 3) + RKF( 3)
WDOT( 6) = WDOT( 6) + RKF( 4)
WDOT( 5) = WDOT( 5) +
& ( 0.2000000000000000D+01 * RKF( 5) )
WDOT( 2) = WDOT( 2) + RKF( 6)
WDOT( 6) = WDOT( 6) + RKF( 7)
WDOT( 2) = WDOT( 2) + RKF( 8)
WDOT( 2) = WDOT( 2) + RKF( 9)
WDOT( 8) = WDOT( 8) + RKF( 10)
WDOT( 5) = WDOT( 5) + RKF( 11)
WDOT( 7) = WDOT( 7) + RKF( 12)
WDOT( 5) = WDOT( 5) +
& ( 0.2000000000000000D+01 * RKF( 13) )
WDOT( 7) = WDOT( 7) + RKF( 14)
WDOT( 6) = WDOT( 6) + RKF( 15)
WDOT( 7) = WDOT( 7) + RKF( 16)
WDOT( 2) = WDOT( 2) + RKF( 17)
WDOT( 7) = WDOT( 7) + RKF( 18)
WDOT( 7) = WDOT( 7) +
& ( 0.2000000000000000D+01 * RKF( 19) )
WDOT( 8) = WDOT( 8) + RKF( 20)
WDOT( 7) = WDOT( 7) + RKF( 21)
WDOT( 8) = WDOT( 8) + RKF( 22)
WDOT( 3) = WDOT( 3) + RKF( 23)
WDOT( 2) = WDOT( 2) + RKF( 24)
WDOT( 7) = WDOT( 7) + RKF( 25)
WDOT( 3) = WDOT( 3) + RKF( 26)
WDOT( 5) = WDOT( 5) +
& ( 0.2000000000000000D+01 * RKF( 27) )
WDOT( 8) = WDOT( 8) + RKF( 28)
WDOT( 5) = WDOT( 5) + RKF( 29)
WDOT( 4) = WDOT( 4) + RKF( 30)
WDOT( 4) = WDOT( 4) +
& ( 0.2000000000000000D+01 * RKF( 31) )
WDOT( 1) = WDOT( 1) + RKF( 32)
WDOT( 3) = WDOT( 3) +
& ( 0.2000000000000000D+01 * RKF( 33) )
WDOT( 6) = WDOT( 6) + RKF( 34)
WDOT( 5) = WDOT( 5) + RKF( 1)
WDOT( 1) = WDOT( 1) + RKF( 2)
WDOT( 5) = WDOT( 5) + RKF( 3)
WDOT( 4) = WDOT( 4) + RKF( 4)
WDOT( 4) = WDOT( 4) + RKF( 6)
WDOT( 5) = WDOT( 5) + RKF( 7)
WDOT( 3) = WDOT( 3) + RKF( 8)
WDOT( 7) = WDOT( 7) + RKF( 9)
WDOT( 5) = WDOT( 5) + RKF( 10)
WDOT( 1) = WDOT( 1) + RKF( 11)
WDOT( 4) = WDOT( 4) + RKF( 12)
WDOT( 3) = WDOT( 3) + RKF( 14)
WDOT( 1) = WDOT( 1) + RKF( 15)
WDOT( 3) = WDOT( 3) + RKF( 16)
WDOT( 1) = WDOT( 1) + RKF( 17)
WDOT( 5) = WDOT( 5) + RKF( 18)
WDOT( 1) = WDOT( 1) + RKF( 20)
WDOT( 6) = WDOT( 6) + RKF( 21)
WDOT( 3) = WDOT( 3) + RKF( 22)
WDOT( 5) = WDOT( 5) + RKF( 23)
WDOT( 1) = WDOT( 1) + RKF( 26)
WDOT( 3) = WDOT( 3) + RKF( 30)
C
RETURN
END
C