R1	acU*r1
U1	acU*r1t
Temp1	Tref*theta1
rhoLP1	rhoLcU*exp(qL1)
rhoVP1	rhoVcU*exp(qV1)
Vb1	4*pi*R1^3/3
pLP1	pfpp(Temp1/1[K],rhoLP1/(1[kg/m^3]))*1[Pa]
pVP1	pfpp(Temp1/1[K],rhoVP1/(1[kg/m^3]))*1[Pa]
uLP1	pfpu(Temp1/1[K],rhoLP1/(1[kg/m^3]))*1[J/kg]
uVP1	pfpu(Temp1/1[K],rhoVP1/(1[kg/m^3]))*1[J/kg]
gLP1	pfpg(Temp1/1[K],rhoLP1/(1[kg/m^3]))*1[J/kg]
gVP1	pfpg(Temp1/1[K],rhoVP1/(1[kg/m^3]))*1[J/kg]
fv1	(Vb1-mpcU/rhoLP1)/(mpcU*(1/rhoVP1-1/rhoLP1))
Rg1	(max(fv1*mpcU/rhoVP1,4*pi*(1[nm])^3/3)*3/(4*pi))^(1/3)
Sspec1	((1-fv1)*(uLP1+pLP1/rhoLP1-gLP1)+fv1*(uVP1+pVP1/rhoVP1-gVP1))/Temp1
qvol1	4*pi*R1^2*U1
pB1	pLP1-2*sigcU/R1-4*muL*U1/R1
En1	2*pi*rhoL*R1^3*U1^2*(1-R1/bcU)
Ei1	mpcU*((1-fv1)*uLP1+fv1*uVP1)+4*pi*sigIcU*Rg1^2-mpcU*uDcU-DP*mpcU/rhoDcU-4*pi*sigcU*rDcU^2
Ep1	DP*Vb1+4*pi*sigcU*R1^2
En	100*(En1)
Ei	100*(Ei1)
Ep	100*(Ep1)
Ein	E0cU
phiDt	acU^2/taucU*chit
pAc	-rhoL*phiDt
u2ac	(acU^2/taucU)^2*(chix^2+chiy^2+chiz^2)
Ea	ivcU(rhoL/2*(u2ac+phiDt^2/cL^2))