Merge branch 'develop' into chore_catch_ortho_fail

Author: pcarruscag

SU2_CFD/src/numerics/flow/flow_diffusion.cpp

@@ -163,7 +163,7 @@ void CAvgGrad_Base::AddTauWall(const su2double *UnitNormal,
   GeometryToolbox::TangentProjection(nDim, tau, UnitNormal, TauTangent);
 
   su2double WallShearStress = GeometryToolbox::Norm(nDim, TauTangent);
-  su2double Scale = TauWall / WallShearStress;
+  su2double Scale = TauWall / fmax(WallShearStress, EPS);
 
   /*--- Scale the stress tensor by the ratio of the wall shear stress
    (from wall functions) to the one computed above. ---*/

SU2_CFD/src/solvers/CRadP1Solver.cpp

@@ -575,7 +575,7 @@ void CRadP1Solver::SetTime_Step(CGeometry *geometry, CSolver **solver_container,
   su2double Area, Vol, Lambda;
   su2double Global_Delta_Time = 1E6, Local_Delta_Time = 0.0, K_v = 0.25;
   su2double CFL = config->GetCFL_Rad();
-  su2double GammaP1 = 1.0 / (3.0*(Absorption_Coeff + Scattering_Coeff));
+  su2double GammaP1 = 1.0 / (3.0*fmax(Absorption_Coeff + Scattering_Coeff, EPS));
   const su2double* Normal;
 
   /*--- Compute spectral radius based on thermal conductivity ---*/

SU2_PY/SU2/eval/gradients.py

@@ -259,7 +259,7 @@ def adjoint(func_name, config, state=None):
             name = files["RESTART_FILE_1"]
             name = su2io.expand_part(name, config)
             link.extend(name)
-        if "RESTART_FILE_1" in files:  # not the case for 1st order time stepping
+        if "RESTART_FILE_2" in files:  # not the case for 1st order time stepping
             name = files["RESTART_FILE_2"]
             name = su2io.expand_part(name, config)
             link.extend(name)

SU2_PY/SU2/io/tools.py

@@ -1064,9 +1064,11 @@ def expand_time(name, config):
             name_pat = add_suffix(name, "%05d")
             names = [name_pat % i for i in range(n_start_time, n_time)]
         else:
+            names = []  # Initialize empty list before loop
             for n in range(len(name)):
                 name_pat = add_suffix(name[n], "%05d")
-                names = [name_pat % i for i in range(n_start_time, n_time)]
+                # Use extend() to accumulate all filenames (consistent with expand_zones)
+                names.extend([name_pat % i for i in range(n_start_time, n_time)])
     else:
         if not isinstance(name, list):
             names = [name]

SU2_PY/SU2/run/interface.py

@@ -38,7 +38,13 @@
 #  Setup
 # ------------------------------------------------------------
 
-SU2_RUN = os.environ["SU2_RUN"]
+try:
+    SU2_RUN = os.environ["SU2_RUN"]
+except KeyError as exc:
+    raise RuntimeError(
+        'Environment variable "SU2_RUN" is not set. Please set SU2_RUN to the SU2 installation bin directory.'
+    ) from exc
+
 sys.path.append(SU2_RUN)
 quote = '"' if sys.platform == "win32" else ""
 
@@ -263,14 +269,17 @@ def run_command(Command):
     """
     sys.stdout.flush()
 
+    # Avoid potential deadlocks when a child process writes heavily to stderr:
+    # read stderr via communicate() while streaming stdout to the console.
     # Use communicate() to continuously drain stderr and avoid deadlocks if the
     # subprocess writes a lot of output to stderr.
     proc = subprocess.Popen(
         Command, shell=True, stdout=sys.stdout, stderr=subprocess.PIPE
     )
     _, stderr = proc.communicate()
+
     return_code = proc.returncode
-    message = stderr.decode(errors="replace")
+    message = (stderr or b"").decode(errors="replace")
 
     if return_code < 0:
         message = "SU2 process was terminated by signal '%s'\n%s" % (