minor update

Author: Hyeoksu Lee

examples/3D_turb_mixing/turbulence_stat/run_turbulence.m

@@ -305,7 +305,7 @@ function f_compute_Reynolds_stress(rho, vel_fluc, y_norm_vth, timestep)
     S23 = squeeze(0.5*(dvel_ds(2,3,:,:,:) + dvel_ds(3,2,:,:,:)));
     S33 = squeeze(dvel_ds(3,3,:,:,:));
 
-    SS = S11.^2 + S12.^2 + S13.^2 + S12.^2 + S22.^2 + S23.^2 + S13.^2 + S23.^2 + S33.^2;
+    SS = S11.^2 + 2*S12.^2 + 2*S13.^2 + S22.^2 + 2*S23.^2 + S33.^2;
 end
 
 % Compute TKE budget