@@ -121,42 +121,42 @@ namespace aspect
121121 in.requested_properties
122122 = MaterialModel::MaterialProperties::uninitialized;
123123
124- for (unsigned int i= 0 ; i<property_names. size (); ++i )
124+ for (const auto &property_name : property_names )
125125 {
126- if (property_names[i] == " viscosity" )
126+ if (property_name == " viscosity" )
127127 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::viscosity;
128128
129- else if (property_names[i] == " density" )
129+ else if (property_name == " density" )
130130 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::density;
131131
132- else if (property_names[i] == " thermal expansivity" )
132+ else if (property_name == " thermal expansivity" )
133133 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::thermal_expansion_coefficient;
134134
135- else if (property_names[i] == " specific heat" )
135+ else if (property_name == " specific heat" )
136136 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::specific_heat;
137137
138- else if (property_names[i] == " thermal conductivity" )
138+ else if (property_name == " thermal conductivity" )
139139 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::thermal_conductivity;
140140
141- else if (property_names[i] == " thermal diffusivity" )
141+ else if (property_name == " thermal diffusivity" )
142142 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::density
143143 | MaterialModel::MaterialProperties::specific_heat
144144 | MaterialModel::MaterialProperties::thermal_conductivity;
145145
146- else if (property_names[i] == " compressibility" )
146+ else if (property_name == " compressibility" )
147147 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::compressibility;
148148
149- else if (property_names[i] == " entropy derivative pressure" )
149+ else if (property_name == " entropy derivative pressure" )
150150 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::entropy_derivative_pressure;
151151
152- else if (property_names[i] == " entropy derivative temperature" )
152+ else if (property_name == " entropy derivative temperature" )
153153 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::entropy_derivative_temperature;
154154
155- else if (property_names[i] == " reaction terms" )
155+ else if (property_name == " reaction terms" )
156156 in.requested_properties = in.requested_properties | MaterialModel::MaterialProperties::reaction_terms;
157157
158158 // We don't know what a material model needs to compute the melt fraction.
159- else if (property_names[i] == " melt fraction" )
159+ else if (property_name == " melt fraction" )
160160 in.requested_properties = MaterialModel::MaterialProperties::all_properties;
161161 }
162162
@@ -201,48 +201,50 @@ namespace aspect
201201 for (unsigned int q=0 ; q<n_quadrature_points; ++q)
202202 {
203203 unsigned output_index = 0 ;
204- for (unsigned int i= 0 ; i<property_names. size (); ++i, ++output_index )
204+ for (const auto &property_name : property_names )
205205 {
206- if (property_names[i] == " viscosity" )
206+ if (property_name == " viscosity" )
207207 computed_quantities[q][output_index] = out.viscosities [q];
208208
209- else if (property_names[i] == " density" )
209+ else if (property_name == " density" )
210210 computed_quantities[q][output_index] = out.densities [q];
211211
212- else if (property_names[i] == " thermal expansivity" )
212+ else if (property_name == " thermal expansivity" )
213213 computed_quantities[q][output_index] = out.thermal_expansion_coefficients [q];
214214
215- else if (property_names[i] == " specific heat" )
215+ else if (property_name == " specific heat" )
216216 computed_quantities[q][output_index] = out.specific_heat [q];
217217
218- else if (property_names[i] == " thermal conductivity" )
218+ else if (property_name == " thermal conductivity" )
219219 computed_quantities[q][output_index] = out.thermal_conductivities [q];
220220
221- else if (property_names[i] == " thermal diffusivity" )
221+ else if (property_name == " thermal diffusivity" )
222222 computed_quantities[q][output_index] = out.thermal_conductivities [q]/(out.densities [q]*out.specific_heat [q]);
223223
224- else if (property_names[i] == " compressibility" )
224+ else if (property_name == " compressibility" )
225225 computed_quantities[q][output_index] = out.compressibilities [q];
226226
227- else if (property_names[i] == " entropy derivative pressure" )
227+ else if (property_name == " entropy derivative pressure" )
228228 computed_quantities[q][output_index] = out.entropy_derivative_pressure [q];
229229
230- else if (property_names[i] == " entropy derivative temperature" )
230+ else if (property_name == " entropy derivative temperature" )
231231 computed_quantities[q][output_index] = out.entropy_derivative_temperature [q];
232232
233- else if (property_names[i] == " reaction terms" )
233+ else if (property_name == " reaction terms" )
234234 {
235235 for (unsigned int k=0 ; k<this ->n_compositional_fields (); ++k, ++output_index)
236236 {
237237 computed_quantities[q][output_index] = out.reaction_terms [q][k];
238238 }
239239 --output_index;
240240 }
241- else if (property_names[i] == " melt fraction" )
241+ else if (property_name == " melt fraction" )
242242 computed_quantities[q][output_index] = melt_fractions[q];
243243 else
244244 AssertThrow (false ,
245245 ExcMessage (" Material property not implemented for this postprocessor." ));
246+
247+ ++output_index;
246248 }
247249 }
248250 }
0 commit comments