R1	ac5*r1
U1	ac5*r1t
Temp1	Tref*theta1
rhoLP1	rhoLc5*exp(qL1)
rhoVP1	rhoVc5*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-mpc5/rhoLP1)/(mpc5*(1/rhoVP1-1/rhoLP1))
Rg1	(max(fv1*mpc5/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*sigc5/R1-4*muL*U1/R1
En1	2*pi*rhoL*R1^3*U1^2*(1-R1/bc5)
Ei1	mpc5*((1-fv1)*uLP1+fv1*uVP1)+4*pi*sigIc5*Rg1^2-mpc5*uDc5-DP*mpc5/rhoDc5-4*pi*sigc5*rDc5^2
Ep1	DP*Vb1+4*pi*sigc5*R1^2
R2	ac5*r2
U2	ac5*r2t
Temp2	Tref*theta2
rhoLP2	rhoLc5*exp(qL2)
rhoVP2	rhoVc5*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-mpc5/rhoLP2)/(mpc5*(1/rhoVP2-1/rhoLP2))
Rg2	(max(fv2*mpc5/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*sigc5/R2-4*muL*U2/R2
En2	2*pi*rhoL*R2^3*U2^2*(1-R2/bc5)
Ei2	mpc5*((1-fv2)*uLP2+fv2*uVP2)+4*pi*sigIc5*Rg2^2-mpc5*uDc5-DP*mpc5/rhoDc5-4*pi*sigc5*rDc5^2
Ep2	DP*Vb2+4*pi*sigc5*R2^2
R3	ac5*r3
U3	ac5*r3t
Temp3	Tref*theta3
rhoLP3	rhoLc5*exp(qL3)
rhoVP3	rhoVc5*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-mpc5/rhoLP3)/(mpc5*(1/rhoVP3-1/rhoLP3))
Rg3	(max(fv3*mpc5/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*sigc5/R3-4*muL*U3/R3
En3	2*pi*rhoL*R3^3*U3^2*(1-R3/bc5)
Ei3	mpc5*((1-fv3)*uLP3+fv3*uVP3)+4*pi*sigIc5*Rg3^2-mpc5*uDc5-DP*mpc5/rhoDc5-4*pi*sigc5*rDc5^2
Ep3	DP*Vb3+4*pi*sigc5*R3^2
R4	ac5*r4
U4	ac5*r4t
Temp4	Tref*theta4
rhoLP4	rhoLc5*exp(qL4)
rhoVP4	rhoVc5*exp(qV4)
Vb4	4*pi*R4^3/3
pLP4	pfpp(Temp4/1[K],rhoLP4/(1[kg/m^3]))*1[Pa]
pVP4	pfpp(Temp4/1[K],rhoVP4/(1[kg/m^3]))*1[Pa]
uLP4	pfpu(Temp4/1[K],rhoLP4/(1[kg/m^3]))*1[J/kg]
uVP4	pfpu(Temp4/1[K],rhoVP4/(1[kg/m^3]))*1[J/kg]
gLP4	pfpg(Temp4/1[K],rhoLP4/(1[kg/m^3]))*1[J/kg]
gVP4	pfpg(Temp4/1[K],rhoVP4/(1[kg/m^3]))*1[J/kg]
fv4	(Vb4-mpc5/rhoLP4)/(mpc5*(1/rhoVP4-1/rhoLP4))
Rg4	(max(fv4*mpc5/rhoVP4,4*pi*(1[nm])^3/3)*3/(4*pi))^(1/3)
Sspec4	((1-fv4)*(uLP4+pLP4/rhoLP4-gLP4)+fv4*(uVP4+pVP4/rhoVP4-gVP4))/Temp4
qvol4	4*pi*R4^2*U4
pB4	pLP4-2*sigc5/R4-4*muL*U4/R4
En4	2*pi*rhoL*R4^3*U4^2*(1-R4/bc5)
Ei4	mpc5*((1-fv4)*uLP4+fv4*uVP4)+4*pi*sigIc5*Rg4^2-mpc5*uDc5-DP*mpc5/rhoDc5-4*pi*sigc5*rDc5^2
Ep4	DP*Vb4+4*pi*sigc5*R4^2
R5	ac5*r5
U5	ac5*r5t
Temp5	Tref*theta5
rhoLP5	rhoLc5*exp(qL5)
rhoVP5	rhoVc5*exp(qV5)
Vb5	4*pi*R5^3/3
pLP5	pfpp(Temp5/1[K],rhoLP5/(1[kg/m^3]))*1[Pa]
pVP5	pfpp(Temp5/1[K],rhoVP5/(1[kg/m^3]))*1[Pa]
uLP5	pfpu(Temp5/1[K],rhoLP5/(1[kg/m^3]))*1[J/kg]
uVP5	pfpu(Temp5/1[K],rhoVP5/(1[kg/m^3]))*1[J/kg]
gLP5	pfpg(Temp5/1[K],rhoLP5/(1[kg/m^3]))*1[J/kg]
gVP5	pfpg(Temp5/1[K],rhoVP5/(1[kg/m^3]))*1[J/kg]
fv5	(Vb5-mpc5/rhoLP5)/(mpc5*(1/rhoVP5-1/rhoLP5))
Rg5	(max(fv5*mpc5/rhoVP5,4*pi*(1[nm])^3/3)*3/(4*pi))^(1/3)
Sspec5	((1-fv5)*(uLP5+pLP5/rhoLP5-gLP5)+fv5*(uVP5+pVP5/rhoVP5-gVP5))/Temp5
qvol5	4*pi*R5^2*U5
pB5	pLP5-2*sigc5/R5-4*muL*U5/R5
En5	2*pi*rhoL*R5^3*U5^2*(1-R5/bc5)
Ei5	mpc5*((1-fv5)*uLP5+fv5*uVP5)+4*pi*sigIc5*Rg5^2-mpc5*uDc5-DP*mpc5/rhoDc5-4*pi*sigc5*rDc5^2
Ep5	DP*Vb5+4*pi*sigc5*R5^2
En	1*(En1+En2+En3+En4+En5)
Ei	1*(Ei1+Ei2+Ei3+Ei4+Ei5)
Ep	1*(Ep1+Ep2+Ep3+Ep4+Ep5)
Ein	E0c5
phiDt	ac5^2/tauc5*chit
pAc	-rhoL*phiDt
u2ac	(ac5^2/tauc5)^2*(chix^2+chiy^2+chiz^2)
Ea	ivc5(rhoL/2*(u2ac+phiDt^2/cL^2))