38 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
47 int ill, ivstat, ierrth;
48 F77NAME(FPTHLR12)(&un, &T[ind * ncomp + id], &Hl, &ill, &ivstat, &ierrth);
52 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
56 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
57 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
60 double hvsp, hlsp, dhvsp1, dhlsp1;
61 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
65 F77NAME(FTLR12)(&un, &P[ind], &tsat, &rhog_sat, &hvsp , &Hl, &Tl, &rhol, &ill, &ivstat, &ierrth);
73 int ill, ivstat, ierrth;
77 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
80 double rhog, rhol, DP_rhoL, DP_rhoG;
82 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &T[ind * ncomp + id] , &rhog, &rhol , &DP_Tl, &DP_rhoL, &DP_rhoG, &ill, &ivstat, &ierrth);
88 for (
auto& val : res) val = 0;
94 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
103 int ill, ivstat, ierrth;
104 F77NAME(FPTHLR12)(&un, &T[ind * ncomp + id], &Hl, &ill, &ivstat, &ierrth);
108 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
112 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
113 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
116 double hvsp, hlsp, dhvsp1, dhlsp1;
117 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
121 F77NAME(FTLR12)(&un, &P[ind], &tsat, &rhog_sat, &hvsp , &Hl, &Tl, &rhol, &ill, &ivstat, &ierrth);
124 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
125 F77NAME(FCPLR12)(&un, &Tl, &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
127 val = dP_rhol - dHl_rhol* dP_Tl/dHl_Tl;
130 for (
auto& val : res)
133 int ill, ivstat, ierrth;
136 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
137 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
142 for (
auto& val : res) val = 0;
148 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
152 for (
auto& val : res)
158 int ill, ivstat, ierrth;
159 F77NAME(FPTHLR12)(&un, &T[ind * ncomp + id], &Hl, &ill, &ivstat, &ierrth);
163 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
167 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
168 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
171 double hvsp, hlsp, dhvsp1, dhlsp1;
172 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
176 F77NAME(FTLR12)(&un, &P[ind], &tsat, &rhog_sat, &hvsp , &Hl, &Tl, &rhol, &ill, &ivstat, &ierrth);
179 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
180 F77NAME(FCPLR12)(&un, &Tl, &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
182 val = dHl_rhol/dHl_Tl;
184 else for (
auto& val : res)
188 int ill, ivstat, ierrth;
189 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
190 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
192 val = dHl_rhol/dHl_Tl;
195 for (
auto& val : res) val = 0;
219 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
222 for (
auto& val : res)
228 int ill, ivstat, ierrth;
229 F77NAME(FPTHLR12)(&un, &T[ind * ncomp + id], &Hl, &ill, &ivstat, &ierrth);
233 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
237 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
238 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
241 double hvsp, hlsp, dhvsp1, dhlsp1;
242 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
246 F77NAME(FTLR12)(&un, &P[ind], &tsat, &rhog_sat, &hvsp , &Hl, &Tl, &rhol, &ill, &ivstat, &ierrth);
249 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
250 F77NAME(FCPLR12)(&un, &Tl, &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
254 else for (
auto& val : res)
257 int ill, ivstat, ierrth;
259 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
260 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
265 for (
auto& val : res) val = 0;
271 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
275 for (
auto& val : res)
281 int ill, ivstat, ierrth;
282 F77NAME(FPTHLR12)(&un, &T[ind * ncomp + id], &Hl, &ill, &ivstat, &ierrth);
286 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
290 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
291 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
294 double hvsp, hlsp, dhvsp1, dhlsp1;
295 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
299 F77NAME(FTLR12)(&un, &P[ind], &tsat, &rhog_sat, &hvsp , &Hl, &Tl, &rhol, &ill, &ivstat, &ierrth);
302 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
303 F77NAME(FCPLR12)(&un, &Tl, &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
309 double cpl2, dP_Tl2, dHl_Tl2, dP_rhol2, dHl_rhol2, dP_cpl2, dHl_cpl2;
310 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl2, &dP_Tl2, &dHl_Tl2, &dP_rhol2, &dHl_rhol2, &dP_cpl2, &dHl_cpl2, &ill, &ivstat, &ierrth);
311 if ((abs(1./dHl_Tl)>1.e-10) && (abs((1./dHl_Tl-1./dHl_Tl2)/(1./dHl_Tl)) > 1.e-4 )) Cout <<
"dT_h2 : " << 1./dHl_Tl2 <<
" dT_h : " << 1./dHl_Tl << finl ;
313 else for (
auto& val : res)
316 int ill, ivstat, ierrth;
317 double cpl, dP_Tl, dHl_Tl, dP_rhol, dHl_rhol, dP_cpl, dHl_cpl;
318 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl, &dP_Tl, &dHl_Tl, &dP_rhol, &dHl_rhol, &dP_cpl, &dHl_cpl, &ill, &ivstat, &ierrth);
322 for (
auto& val : res) val = 0;
328 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
331 for (
auto& val : res)
334 double cpl, dP_Tl, DHl_Tl, DP_Rhol, DHl_Rhol, DP_cpl, DHl_cpl ;
335 int ill, ivstat, ierrth;
336 F77NAME(FCPLR12)(&un, &T[ind * ncomp + id], &cpl, &dP_Tl, &DHl_Tl, &DP_Rhol, &DHl_Rhol, &DP_cpl, &DHl_cpl, &ill, &ivstat, &ierrth);
340 for (
auto& val : res) val = 0;