dnum_dtheta = teos(02)
& + ct*(2.*teos(03) + 3.*teos(04)*ct)
& + sa*(teos(06) + 2.*teos(07)*ct
& + sqrtsa*(teos(09) + ct*(2.*teos(10) + 3.*teos(11)*ct)))
& + p*( teos(13) + 2.*teos(14)*ct + sa*2.*teos(16)
& + p*(teos(18) + 2.*teos(19)*ct))
dden_dtheta = teos(22)
& + ct*(2.*teos(23) + ct*(3.*teos(24) + 4.*teos(25)*ct))
& + sa*(teos(27)
& + ct*(2.*teos(28) + ct*(3.*teos(29) + 4.*teos(30)*ct))
& + sqrtsa*(teos(32)
& + ct*(2.*teos(33) + ct*(3.*teos(34) + 4.*teos(35)*ct))))
& + p*(teos(38) + ct*(2.*teos(39) + 3.*teos(40)*ct)
& + teos(42)
& + p*(teos(44) + 2.*teos(45)*ct + teos(46)*sa
& + p*teos(48) ))
dnum_dsalt = teos(05) + ct*(teos(06) + teos(07)*ct)
& + 1.5*sqrtsa*(teos(08)
& + ct*(teos(09) + ct*(teos(10) + teos(11)*ct)))
& + p*(teos(15) + teos(16)*ct + p*teos(20))
dden_dsalt = teos(26)
& + ct*(teos(27) + ct*(teos(28) + ct*(teos(29) + teos(30)*ct)))
& + 2.*teos(36)*sa
& + 1.5*sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)
& + ct*(teos(34) + teos(35)*ct))))
& + p*(teos(41) + teos(42)*ct + p*teos(46))
rhoNum = teos(01)
& + ct*(teos(02) + ct*(teos(03) + teos(04)*ct))
& + sa*(teos(05) + ct*(teos(06) + teos(07)*ct)
& + sqrtsa*(teos(08) + ct*(teos(09)
& + ct*(teos(10) + teos(11)*ct))))
& + p*(teos(12) + ct*(teos(13) + teos(14)*ct)
& + sa*(teos(15) + teos(16)*ct)
& + p*(teos(17) + ct*(teos(18) + teos(19)*ct) + teos(20)*sa))
den = teos(21)
& + ct*(teos(22) + ct*(teos(23) + ct*(teos(24) + teos(25)*ct)))
& + sa*(teos(26) + ct*(teos(27) + ct*(teos(28)
& + ct*(teos(29) + teos(30)*ct)))
& + teos(36)*sa
% + sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)
& + ct*(teos(34) + teos(35)*ct)))))
% + p*(teos(37) + ct*(teos(38) + ct*(teos(39) + teos(40)*ct))
% + sa*(teos(41) + teos(42)*ct)
% + p*(teos(43) + ct*(teos(44) + teos(45)*ct + teos(46)*sa)
% + p*(teos(47) + teos(48)*ct)))
rhoNum(i,j) = teos(01)
& + sqrtsa*(teos(31) + ct*(teos(32) + ct*(teos(33)
& + p*(teos(37) + ct*(teos(38) + ct*(teos(39) + teos(40)*ct))
& + sa*(teos(41) + teos(42)*ct)
& + p*(teos(43) + ct*(teos(44) + teos(45)*ct + teos(46)*sa)
& + p*(teos(47) + teos(48)*ct)))
teos(01) = 9.998420897506056 _d +02
teos(02) = 2.839940833161907 _d 00
teos(03) = -3.147759265588511 _d -02
teos(04) = 1.181805545074306 _d -03
teos(05) = -6.698001071123802 _d 00
teos(06) = -2.986498947203215 _d -02
teos(07) = 2.327859407479162 _d -04
teos(08) = -3.988822378968490 _d -02
teos(09) = 5.095422573880500 _d -04
teos(10) = -1.426984671633621 _d -05
teos(11) = 1.645039373682922 _d -07
teos(12) = -2.233269627352527 _d -02
teos(13) = -3.436090079851880 _d -04
teos(14) = 3.726050720345733 _d -06
teos(15) = -1.806789763745328 _d -04
teos(16) = 6.876837219536232 _d -07
teos(17) = -3.087032500374211 _d -07
teos(18) = -1.988366587925593 _d -08
teos(19) = -1.061519070296458 _d -11
teos(20) = 1.550932729220080 _d -10
teos(21) = 1.000000000000000 _d 00
teos(22) = 2.775927747785646 _d -03
teos(23) = -2.349607444135925 _d -05
teos(24) = 1.119513357486743 _d -06
teos(25) = 6.743689325042773 _d -10
teos(26) = -7.521448093615448 _d -03
teos(27) = -2.764306979894411 _d -05
teos(28) = 1.262937315098546 _d -07
teos(29) = 9.527875081696435 _d -10
teos(30) = -1.811147201949891 _d -11
teos(31) = -3.303308871386421 _d -05
teos(32) = 3.801564588876298 _d -07
teos(33) = -7.672876869259043 _d -09
teos(34) = -4.634182341116144 _d -11
teos(35) = 2.681097235569143 _d -12
teos(36) = 5.419326551148740 _d -06
teos(37) = -2.742185394906099 _d -05
teos(38) = -3.212746477974189 _d -07
teos(39) = 3.191413910561627 _d -09
teos(40) = -1.931012931541776 _d -12
teos(41) = -1.105097577149576 _d -07
teos(42) = 6.211426728363857 _d -10
teos(43) = -1.119011592875110 _d -10
teos(44) = -1.941660213148725 _d -11
teos(45) = -1.864826425365600 _d -14
teos(46) = 1.119522344879478 _d -14
teos(47) = -1.200507748551599 _d -15
teos(48) = 6.057902487546866 _d -17
teos(k) = 0. _d 0
_RL teos(48)
& teos