_RL p0 , p25 , p5 , p2
parameter ( p0 = 0.0, p25 = 0.25, p5 = 0.5, p2 = 2.0 )
v(i,1) = p2 * v(i,1) +
zflag = p2 + KRi_range(1) * KRi_range(2)
v(i,k) = p2 * v(i,k) +
zflag = p2 + KRi_range(k)*(KRi_range(kp1)+KRi_range(km1))
p2 = p1*p1
& + p2*( eosJMDCKP(10)
& + p2*s1*( eosJMDCKP(13)
_RL t1, t2, t3, s1, s3o2, p1, p2, sp5, p1t1
& + p2 * ( eosP9 +
& + p2*s1* ( eosP12 +
& + p2*( eosJMDCKP(9)
& + p2*s1*( eosJMDCKP(12)
_RL t, t2, t3, t4, s, s3o2, p, p2
p2 = p*p
& + p2*s* ( eosP12 +
& + p2*s*( eosJMDCKP(12)
_RL t1, t2, t3, t4, s1, s3o2, p1, p2, sp5, p1t1