Add mu omp fix

Author: AbbBallout

Solver/src/MultiphaseSolver/SpatialDiscretization.f90

@@ -1682,7 +1682,7 @@ subroutine compute_viscosity_at_faces(no_of_faces, no_of_sides, face_ids, mesh)
 
 
          if ( LESModel % Active ) then
-!$omp do schedule(runtime) private(i,j,delta,mu_smag)
+!$omp do schedule(runtime) private(i,j,delta,mu_smag,side)
             do iFace = 1, no_of_faces
                associate(f => mesh % faces(face_ids(iFace)))