Fix that finite_differences script pulls solution files for restarted run.

Author: TobiKattmann

SU2_PY/SU2/eval/gradients.py

@@ -838,6 +838,14 @@ def findiff( config, state=None ):
         name = su2io.expand_time(name,config)
         link.extend(name)
 
+    # files: restart solution for dual-time stepping first and second order
+    if 'RESTART_FILE_1' in files:
+        name = files['RESTART_FILE_1']
+        pull.append(name)
+    if 'RESTART_FILE_2' in files:
+        name = files['RESTART_FILE_2']
+        pull.append(name)
+
     # files: target equivarea distribution
     if 'EQUIV_AREA' in special_cases and 'TARGET_EA' in files:
         pull.append(files['TARGET_EA'])

SU2_PY/finite_differences.py

@@ -77,6 +77,21 @@ def finite_differences( filename           ,
     # State
     state = SU2.io.State()
     state.find_files(config)
+
+    # add restart files to state.FILES
+    if config.get('TIME_DOMAIN', 'NO') == 'YES' and config.get('RESTART_SOL', 'NO') == 'YES':
+        restart_name = config['RESTART_FILENAME'].split('.')[0]
+        restart_filename = restart_name + '_' + str(int(config['RESTART_ITER'])-1).zfill(5) + '.dat'
+        if not os.path.isfile(restart_filename): # throw, if restart files does not exist
+            sys.exit("Error: Restart file <" + restart_filename + "> not found.")
+        state['FILES']['RESTART_FILE_1'] = restart_filename
+
+        # use only, if time integration is second order
+        if config.get('TIME_MARCHING', 'NO') == 'DUAL_TIME_STEPPING-2ND_ORDER':
+            restart_filename = restart_name + '_' + str(int(config['RESTART_ITER'])-2).zfill(5) + '.dat'
+            if not os.path.isfile(restart_filename): # throw, if restart files does not exist
+                sys.exit("Error: Restart file <" + restart_filename + "> not found.")
+            state['FILES']['RESTART_FILE_2'] =restart_filename
 
     # Finite Difference Gradients
     SU2.eval.gradients.findiff(config,state)