def fill(self, fase, estado):
fase.M = self.M
fase.rho = unidades.Density(estado.rho)
fase.v = unidades.SpecificVolume(estado.v)
fase.Z = unidades.Dimensionless(self.P*estado.v/R/1000*self.M/self.T)
fase.h = unidades.Enthalpy(estado.h)
fase.s = unidades.SpecificHeat(estado.s)
fase.u = unidades.Enthalpy(estado.u)
fase.a = unidades.Enthalpy(fase.u-self.T*fase.s)
fase.g = unidades.Enthalpy(fase.h-self.T*fase.s)
fase.cv = unidades.SpecificHeat(estado.cv)
fase.cp = unidades.SpecificHeat(estado.cp)
fase.cp_cv = unidades.Dimensionless(fase.cp/fase.cv)
fase.w = unidades.Speed(estado.w)
fase.mu = unidades.Viscosity(estado.mu)
fase.k = unidades.ThermalConductivity(estado.k)
fase.nu = unidades.Diffusivity(fase.mu/fase.rho)
fase.dielec = unidades.Dimensionless(iapws._Dielectric(estado.rho, self.T))
fase.Prandt = unidades.Dimensionless(estado.mu*estado.cp/estado.k)
# fase.joule=unidades.TemperaturePressure(self.Liquido["hjt"], "KkPa")
# fase.xkappa=unidades.InvPressure(DerivTerms("IsothermalCompressibility", self.T, fase.rho, self.name), "kPa")
# fase.alfav=unidades.InvTemperature(-estado.PFC.drhodT_constp()/estado.rho)
cp0 = iapws.prop0(self.T, self.P)
fase.v0 = unidades.SpecificVolume(cp0.v)
fase.h0 = unidades.Enthalpy(cp0.h)
fase.u0 = unidades.Enthalpy(fase.h0-self.P*1000*fase.v0)
fase.s0 = unidades.SpecificHeat(cp0.s)
fase.a0 = unidades.Enthalpy(fase.u0-self.T*fase.s0)
fase.g0 = unidades.Enthalpy(fase.h0-self.T*fase.s0)
fase.cp0 = unidades.SpecificHeat(cp0.cp)
fase.cv0 = unidades.SpecificHeat(cp0.cv)
fase.cp0_cv = unidades.Dimensionless(fase.cp0/fase.cv0)
fase.w0 = cp0.w
# fase.gamma0 = cp0.gamma
fase.f = unidades.Pressure(self.P*exp((fase.g-fase.g0)/R/self.T))