PARAMETER ( B1 = 16.6 )
PARAMETER ( B3 = 1. / B1 )
QE(I,1) = B1*QE(I,1)
_RL A1, A2, B1, B2, C1
PARAMETER( B1 = 16.6D0 )
gam1 = 1. _d 0/3. _d 0 - 2. _d 0 *A1/B1
gam2 = (B2+6. _d 0*A1)/B1
beta1 = A1*B1*(gam1-C1)
beta2 = A1*(B1*(gam1-C1) + 6. _d 0 *A1 + 3. _d 0 *A2)
beta3 = A2*B1*gam1
beta4 = A2*(B1*(gam1+gam2) - 3. _d 0 *A1)
_RL B1
B1 = -1636.75 _d 0 + 12.0408 _d 0*tk - 0.0327957 _d 0*tk*tk
B = B1 + 3.16528 _d 0*tk*tk*tk*(0.00001 _d 0)