R1	acP*r1
U1	acP*r1t
Temp1	Tref*theta1
rhoLP1	rhoLcP*exp(qL1)
rhoVP1	rhoVcP*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-mpcP/rhoLP1)/(mpcP*(1/rhoVP1-1/rhoLP1))
Rg1	(max(fv1*mpcP/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*sigcP/R1-4*muL*U1/R1
En1	2*pi*rhoL*R1^3*U1^2*(1-R1/bcP)
Ei1	mpcP*((1-fv1)*uLP1+fv1*uVP1)+4*pi*sigIcP*Rg1^2-mpcP*uDcP-DP*mpcP/rhoDcP-4*pi*sigcP*rDcP^2
Ep1	DP*Vb1+4*pi*sigcP*R1^2
En	1*(En1)
Ei	1*(Ei1)
Ep	1*(Ep1)
Ein	E0cP
phiDt	acP^2/taucP*chit
pAc	-rhoL*phiDt
u2ac	(acP^2/taucP)^2*(chix^2+chiy^2+chiz^2)
Ea	ivcP(rhoL/2*(u2ac+phiDt^2/cL^2))