temp

Author: henryleberre

src/simulation/include/inline_riemann.fpp

@@ -59,12 +59,9 @@
             end if
         end do
         gamma_avg = Cp_avg/Cv_avg
-        c_avggg = 0.0d0
 
-        do i = 1, num_species
-            Phi_avg(i) = (gamma_avg - 1.d0)*(vel_avg_rms/2.0d0 - h_avg_2(i)) + gamma_avg*gas_constant/mol_weights(i)*T_avg
-            c_avggg = c_avggg + Yi_avg(i)*Phi_avg(i)
-        end do
+        Phi_avg(:) = (gamma_avg - 1.d0)*(vel_avg_rms/2.0d0 - h_avg_2(:)) + gamma_avg*gas_constant/mol_weights(:)*T_avg
+        c_avggg = sum(Yi_avg(:)*Phi_avg(:))
     #:endif
 
 #:enddef roe_avg