38 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
48 int ill, ivstat, ierrth;
49 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
53 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
57 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
58 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
61 double hvsp, hlsp, dhvsp1, dhlsp1;
62 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
65 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
66 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
70 else for (
auto& val : res)
73 int ill, ivstat, ierrth;
77 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
81 double rhog, rhol, DP_rhoL, DP_rhoG;
82 F77NAME(FROVLR12)(&un, &P[ind], &T[ind * ncomp + id], &tsat , &rhog , &rhol, &DP_Tl, &DP_rhoL, &DP_rhoG, &ill, &ivstat, &ierrth);
87 for (
auto& val : res) val = 0;
93 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
102 int ill, ivstat, ierrth;
103 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
107 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
111 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
112 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
115 double hvsp, hlsp, dhvsp1, dhlsp1;
116 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
119 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
120 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
122 val = dP_rhog - dHg_rhog* dP_Tg/dHg_Tg;
125 for (
auto& val : res) val = 0;
131 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
134 for (
auto& val : res)
140 int ill, ivstat, ierrth;
141 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
145 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
149 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
150 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
153 double hvsp, hlsp, dhvsp1, dhlsp1;
154 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
157 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
158 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
160 val = dHg_rhog/dHg_Tg;
163 for (
auto& val : res) val = 0;
187 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
190 for (
auto& val : res)
196 int ill, ivstat, ierrth;
197 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
201 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
205 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
206 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
209 double hvsp, hlsp, dhvsp1, dhlsp1;
210 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
213 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
214 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
219 for (
auto& val : res) val = 0;
225 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
228 for (
auto& val : res)
234 int ill, ivstat, ierrth;
235 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
239 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
243 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
244 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
247 double hvsp, hlsp, dhvsp1, dhlsp1;
248 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
251 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
252 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
257 for (
auto& val : res) val = 0;
263 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
266 for (
auto& val : res)
272 int ill, ivstat, ierrth;
273 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
277 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
281 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
282 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
285 double hvsp, hlsp, dhvsp1, dhlsp1;
286 F77NAME(FHSATR12)(&un, &tsat, &rhog_sat, &rhol_sat, &P[ind], &hvsp, &hlsp, &dtsp1, &DP_rhoG_sat, &DP_rhoL_sat, &dhvsp1, &dhlsp1);
289 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
290 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg , &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog , &ill, &ivstat, &ierrth);
293 double cpg, dP_cpg, dhg_cpg;
294 F77NAME(FCPVR12)(&un, &Tg, &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog, &cpg, &dP_cpg, &dhg_cpg);
299 for (
auto& val : res) val = 0;
314 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
317 for (
auto& val : res)
322 int ill, ivstat, ierrth;
323 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
327 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
331 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
332 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
335 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
336 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg, &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog, &ill, &ivstat, &ierrth);
339 double mug, dP_mug, DHg_mug;
340 F77NAME(FMUVR12)(&un, &P[ind], &T[ind * ncomp + id], &dP_Tg, &dHg_Tg, &mug, &dP_mug, &DHg_mug);
345 for (
auto& val : res) val = 0;
351 assert((
int )T.size() == ncomp * (
int )P.size() && (
int )T.size() == ncomp * (
int )res.size());
354 for (
auto& val : res)
359 int ill, ivstat, ierrth;
360 F77NAME(FPTHGR12)(&un, &P[ind], &T[ind * ncomp + id], &Hg, &ill, &ivstat, &ierrth);
364 F77NAME(FPSATR12)(&un, &P[ind], &tsat, &dtsp1, &ill, &ivstat, &ierrth);
368 double rhog_sat, rhol_sat, DP_rhoL_sat, DP_rhoG_sat;
369 F77NAME(FROVLR12)(&un, &P[ind], &tsat, &tsat , &rhog_sat, &rhol_sat , &DP_Tl, &DP_rhoL_sat, &DP_rhoG_sat, &ill, &ivstat, &ierrth);
372 double Tg, dP_Tg, dHg_Tg, rhog, dP_rhog, dHg_rhog;
373 F77NAME(FTGR12)(&un, &P[ind], &Hg, &tsat, &rhog_sat, &Tg, &dP_Tg, &dHg_Tg, &rhog, &dP_rhog, &dHg_rhog, &ill, &ivstat, &ierrth);
376 double lambdag, dP_lambdag, DHg_lambdag;
377 F77NAME(FMUVR12)(&un, &P[ind], &T[ind * ncomp + id], &dP_Tg, &dHg_Tg, &lambdag, &dP_lambdag, &DHg_lambdag);
382 for (
auto& val : res) val = 0;