function f = fluidPVT(Ppr, BG, MUG, Bwi, prw, cw, vwi, cvw, SW, KRG, KRW, PCGL, SWF, KRGF, KRWF, PCGLF, density_g_sc, ifpcgl, rpt, cs_data, csa_data, cb_data, cba_data)
% f = fluidPVT(Bopb, pb, co, Bwi, prw, cw, vwi, cvw, visopb, cvo, SW, KRO, KRW, PCOW, ifpcow, SWF, KROF, KRWF, PCOWF, rpt,Dosi,Dwsi)
f.Bw = @(p) Bw(p, Bwi, prw, cw);
f.Bg = @(p) Bg(p, BG, Ppr);
f.muw = @(p) muw(p, vwi, prw, cvw);
f.mug = @(p) mug(p, MUG, Ppr);
f.krrg = @(sw) krrg(sw, SW, KRG, SWF, KRGF, rpt);
f.krw = @(sw) krw(sw, SW, KRW, SWF, KRWF, rpt);
f.pcgl = @(sw) pcgl(sw, SW, PCGL, SWF, PCGLF, rpt);
f.cs_absorb = @(cs) cs_absorb_f(cs, cs_data, csa_data);
f.cb_absorb = @(cb) cb_absorb_f(cb, cb_data, cba_data);
f.density_g_sc = density_g_sc;
f.ifpcgl = ifpcgl;
% f.Dosi=Dosi;
% f.Dwsi=Dwsi;