dden_dtheta = eosMDJWFden(1)
& + t1*(2.*eosMDJWFden(2)
& + t1*(3.*eosMDJWFden(3)
& + 4.*eosMDJWFden(4)*t1 ) )
& + s1*(eosMDJWFden(6)
& + t1*(3.*eosMDJWFden(7)*t1
& + 2.*eosMDJWFden(9)*sp5 ) )
& + p1*p1*(3.*eosMDJWFden(11)*t2 + eosMDJWFden(12)*p1)
dden_dsalt = eosMDJWFden(5)
& + t1*( eosMDJWFden(6) + eosMDJWFden(7)*t2 )
& + 1.5*sp5*(eosMDJWFden(8) + eosMDJWFden(9)*t2)
den = eosMDJWFden(0)
& + t1*(eosMDJWFden(1)
& + t1*(eosMDJWFden(2)
& + t1*(eosMDJWFden(3) + t1*eosMDJWFden(4) ) ) )
& + s1*(eosMDJWFden(5)
& + t1*(eosMDJWFden(6)
& + eosMDJWFden(7)*t2)
& + sp5*(eosMDJWFden(8) + eosMDJWFden(9)*t2) )
& + p1*(eosMDJWFden(10)
& + p1t1*(eosMDJWFden(11)*t2 + eosMDJWFden(12)*p1) )
eosMDJWFden( 0) = 1.00000000 _d +00
eosMDJWFden( 1) = 7.28606739 _d -03
eosMDJWFden( 2) = -4.60835542 _d -05
eosMDJWFden( 3) = 3.68390573 _d -07
eosMDJWFden( 4) = 1.80809186 _d -10
eosMDJWFden( 5) = 2.14691708 _d -03
eosMDJWFden( 6) = -9.27062484 _d -06
eosMDJWFden( 7) = -1.78343643 _d -10
eosMDJWFden( 8) = 4.76534122 _d -06
eosMDJWFden( 9) = 1.63410736 _d -09
eosMDJWFden(10) = 5.30848875 _d -06
eosMDJWFden(11) = -3.03175128 _d -16
eosMDJWFden(12) = -1.27934137 _d -17
eosMDJWFden(k) = 0. _d 0
_RL eosMDJWFnum(0:11), eosMDJWFden(0:12)
& eosMDJWFnum, eosMDJWFden