dnum_dtheta = eosMDJWFnum(1)
& + t1*(2.*eosMDJWFnum(2) + 3.*eosMDJWFnum(3)*t1)
& + eosMDJWFnum(5)*s1
& + p1t1*(2.*eosMDJWFnum(8) + 2.*eosMDJWFnum(11)*p1)
dnum_dsalt = eosMDJWFnum(4)
& + eosMDJWFnum(5)*t1
& + 2.*eosMDJWFnum(6)*s1 + eosMDJWFnum(9)*p1
rhoNum = eosMDJWFnum(0)
& + t1*(eosMDJWFnum(1)
& + t1*(eosMDJWFnum(2) + eosMDJWFnum(3)*t1) )
& + s1*(eosMDJWFnum(4)
& + eosMDJWFnum(5)*t1 + eosMDJWFnum(6)*s1)
& + p1*(eosMDJWFnum(7) + eosMDJWFnum(8)*t2
& + eosMDJWFnum(9)*s1
& + p1*(eosMDJWFnum(10) + eosMDJWFnum(11)*t2) )
rhoNum(i,j) = eosMDJWFnum(0)
eosMDJWFnum( 0) = 9.99843699 _d +02
eosMDJWFnum( 1) = 7.35212840 _d +00
eosMDJWFnum( 2) = -5.45928211 _d -02
eosMDJWFnum( 3) = 3.98476704 _d -04
eosMDJWFnum( 4) = 2.96938239 _d +00
eosMDJWFnum( 5) = -7.23268813 _d -03
eosMDJWFnum( 6) = 2.12382341 _d -03
eosMDJWFnum( 7) = 1.04004591 _d -02
eosMDJWFnum( 8) = 1.03970529 _d -07
eosMDJWFnum( 9) = 5.18761880 _d -06
eosMDJWFnum(10) = -3.24041825 _d -08
eosMDJWFnum(11) = -1.23869360 _d -11
eosMDJWFnum(k) = 0. _d 0
_RL eosMDJWFnum(0:11), eosMDJWFden(0:12)
& eosMDJWFnum, eosMDJWFden