1 if (runTime.outputTime())
19 totalRhoN.internalField() = totalRhoN_sum;
38 totalRhoM.internalField() = totalRhoM_sum;
57 forAll(singleSpeciesVelocity, sSV)
61 singleSpeciesVelocity[sSV] =
65 totalMomentum_sum[sSV] +=
74 singleSpeciesVelocity[sSV] = vector::zero;
79 forAll(totalVelocity.internalField(), tV)
81 if (totalMass_sum[tV] > VSMALL)
83 totalVelocity.internalField()[tV] =
89 totalVelocity.internalField()[tV] =
111 forAll(singleSpeciesTemp, sST)
115 singleSpeciesTemp[sST] =
130 totalTemperatureVTerms_sum[sST] +=
148 singleSpeciesTemp[sST] = 0.0;
153 forAll(totalTemperature.internalField(), tT)
155 if(totalN_sum[tT] > 0)
157 totalTemperature.internalField()[tT] =
158 totalTemperatureVTerms_sum[tT]
159 /(3.0 * moleculeCloud::kb * totalN_sum[tT]);
163 totalTemperature.internalField()[tT] = 0.0;
182 forAll(singleSpeciesMeanKE, sSMKE)
186 singleSpeciesMeanKE[sSMKE] =
194 totalKE_sum[sSMKE] +=
204 singleSpeciesMeanKE[sSMKE] = 0.0;
209 forAll(totalMeanKE.internalField(), tMKE)
211 if(totalN_sum[tMKE] > 0)
213 totalMeanKE.internalField()[tMKE] =
219 totalMeanKE.internalField()[tMKE] = 0.0;