21 static double a[15]={-26.6446988,-26.9066506,-27.0619832,-27.1826903,
22 -27.2783527,-27.3595949,-27.569406,-27.611159,-27.654748,-27.70479,
23 -27.745398,-27.776126,-27.807261,-27.833093,-27.866134};
24 static double b[15]={-0.40511045,-0.41644707,-0.45834359,-0.49137946,
25 -0.51931762,-0.54971231,-0.18555807,-0.29204736,-0.36741085,
26 -0.45843009,-0.49753713,-0.51418739,-0.52287028,-0.52445456,
28 static double c[15]={11.29232731,11.71035693,12.89309608,13.85569087,
29 14.67354775,15.56090323,6.147461,9.0304953,11.094731,13.630431,
30 14.721959,15.172335,15.413946,15.458123,15.428761};
31 static double d[15]={.067257375,.07638384,.089925637,.102252192,
32 .111016071,.119518918,0.0093832482,0.028119606,0.039357697,0.050378417,
33 0.051644049,0.051367182,0.04938724,0.050139066,0.043085968};
34 static double e[15]={-1.99108378,-2.26898352,-2.65163846,-3.02333001,
35 -3.29462338,-3.56633674,-1.0019228,-1.5128672,-1.8327058,-2.1866371,
36 -2.2286257,-2.1932699,-2.1205891,-2.1317169,-1.9175186};
37 static double f[15]={-0.0050802618,-0.005849291,-0.0074854405,-0.0085677543,
38 -0.0093067267,-0.0098455637,0.040903604,0.037491802,0.035618861,
39 0.034132954,0.032418252,0.02947883,0.027393564,0.027607009,0.02433868};
40 static double g[15]={.166267838,.196780541,.242675042,.282237211,
41 .310204623,.335160025,-0.81087376,-0.73435108,-0.69164333,-0.64907209,
42 -0.61216299,-0.55239109,-0.51048669,-0.51963194,-0.4504203};
43 static double h[15]={.00020528663,.00027588492,.00033980652,.00041445515,
44 .00046423276,.0005121808,-0.0011986559,-0.0011333973,-0.0010992935,
45 -0.0010878727,-0.0010412678,-0.00095539899,-0.00089141547,-0.00089294364,
47 static double i[15]={-0.0071357493,-0.0099630604,-0.01178647,-0.014696455,
48 -0.01670318,-0.01887373,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.};
63 if( n > 15 || x < 0.2 )
82 if( x >
phycon.TEMP_LIMIT_HIGH_LOG )
84 fprintf(
ioQQQ,
" HydroRecCool called with invalid temperature=%e nelem=%li\n",
100 (1. + b[nm1]*
phycon.telogn[0] +
112 i[nm1]*
powi( x,4) ) /
116 h[nm1]*
powi( x ,4) );
119 hclf_v = pow(10.,fac)*
POW3(nelem+1.);