R1	ac3*r1
U1	ac3*r1t
Temp1	Tref*theta1
rhoLP1	rhoLc3*exp(qL1)
rhoVP1	rhoVc3*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-mpc3/rhoLP1)/(mpc3*(1/rhoVP1-1/rhoLP1))
Rg1	(max(fv1*mpc3/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*sigc3/R1-4*muL*U1/R1
En1	2*pi*rhoL*R1^3*U1^2*(1-R1/bc3)
Ei1	mpc3*((1-fv1)*uLP1+fv1*uVP1)+4*pi*sigIc3*Rg1^2-mpc3*uDc3-DP*mpc3/rhoDc3-4*pi*sigc3*rDc3^2
Ep1	DP*Vb1+4*pi*sigc3*R1^2
R2	ac3*r2
U2	ac3*r2t
Temp2	Tref*theta2
rhoLP2	rhoLc3*exp(qL2)
rhoVP2	rhoVc3*exp(qV2)
Vb2	4*pi*R2^3/3
pLP2	pfpp(Temp2/1[K],rhoLP2/(1[kg/m^3]))*1[Pa]
pVP2	pfpp(Temp2/1[K],rhoVP2/(1[kg/m^3]))*1[Pa]
uLP2	pfpu(Temp2/1[K],rhoLP2/(1[kg/m^3]))*1[J/kg]
uVP2	pfpu(Temp2/1[K],rhoVP2/(1[kg/m^3]))*1[J/kg]
gLP2	pfpg(Temp2/1[K],rhoLP2/(1[kg/m^3]))*1[J/kg]
gVP2	pfpg(Temp2/1[K],rhoVP2/(1[kg/m^3]))*1[J/kg]
fv2	(Vb2-mpc3/rhoLP2)/(mpc3*(1/rhoVP2-1/rhoLP2))
Rg2	(max(fv2*mpc3/rhoVP2,4*pi*(1[nm])^3/3)*3/(4*pi))^(1/3)
Sspec2	((1-fv2)*(uLP2+pLP2/rhoLP2-gLP2)+fv2*(uVP2+pVP2/rhoVP2-gVP2))/Temp2
qvol2	4*pi*R2^2*U2
pB2	pLP2-2*sigc3/R2-4*muL*U2/R2
En2	2*pi*rhoL*R2^3*U2^2*(1-R2/bc3)
Ei2	mpc3*((1-fv2)*uLP2+fv2*uVP2)+4*pi*sigIc3*Rg2^2-mpc3*uDc3-DP*mpc3/rhoDc3-4*pi*sigc3*rDc3^2
Ep2	DP*Vb2+4*pi*sigc3*R2^2
R3	ac3*r3
U3	ac3*r3t
Temp3	Tref*theta3
rhoLP3	rhoLc3*exp(qL3)
rhoVP3	rhoVc3*exp(qV3)
Vb3	4*pi*R3^3/3
pLP3	pfpp(Temp3/1[K],rhoLP3/(1[kg/m^3]))*1[Pa]
pVP3	pfpp(Temp3/1[K],rhoVP3/(1[kg/m^3]))*1[Pa]
uLP3	pfpu(Temp3/1[K],rhoLP3/(1[kg/m^3]))*1[J/kg]
uVP3	pfpu(Temp3/1[K],rhoVP3/(1[kg/m^3]))*1[J/kg]
gLP3	pfpg(Temp3/1[K],rhoLP3/(1[kg/m^3]))*1[J/kg]
gVP3	pfpg(Temp3/1[K],rhoVP3/(1[kg/m^3]))*1[J/kg]
fv3	(Vb3-mpc3/rhoLP3)/(mpc3*(1/rhoVP3-1/rhoLP3))
Rg3	(max(fv3*mpc3/rhoVP3,4*pi*(1[nm])^3/3)*3/(4*pi))^(1/3)
Sspec3	((1-fv3)*(uLP3+pLP3/rhoLP3-gLP3)+fv3*(uVP3+pVP3/rhoVP3-gVP3))/Temp3
qvol3	4*pi*R3^2*U3
pB3	pLP3-2*sigc3/R3-4*muL*U3/R3
En3	2*pi*rhoL*R3^3*U3^2*(1-R3/bc3)
Ei3	mpc3*((1-fv3)*uLP3+fv3*uVP3)+4*pi*sigIc3*Rg3^2-mpc3*uDc3-DP*mpc3/rhoDc3-4*pi*sigc3*rDc3^2
Ep3	DP*Vb3+4*pi*sigc3*R3^2
En	1*(En1+En2+En3)
Ei	1*(Ei1+Ei2+Ei3)
Ep	1*(Ep1+Ep2+Ep3)
Ein	E0c3
phiDt	ac3^2/tauc3*chit
pAc	-rhoL*phiDt
u2ac	(ac3^2/tauc3)^2*(chix^2+chiy^2+chiz^2)
Ea	ivc3(rhoL/2*(u2ac+phiDt^2/cL^2))