sqrtsa = 0. _d 0
sqrtsa = SQRT(sa)
& + sqrtsa*(teos(09) + ct*(2.*teos(10) + 3.*teos(11)*ct)))
& + sqrtsa*(teos(32)
_RL ct, sa, sqrtsa, p
& + 1.5*sqrtsa*(teos(08)
& + 1.5*sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)
& + sqrtsa*(teos(08) + ct*(teos(09)
% + sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)
& + sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)