updating config files

WallyMaier · WallyMaier · commit a240d26f5f6a · 2022-07-04T17:17:40.000-07:00
diff --git a/TestCases/axisymmetric_rans/air_nozzle/air_nozzle.cfg b/TestCases/axisymmetric_rans/air_nozzle/air_nozzle.cfg
@@ -101,5 +101,5 @@ SOLUTION_FILENAME= solution_flow.dat
 OUTPUT_WRT_FREQ= 1000
 %
 % Note: This cfg is used for a primal and adjoint Testcase, therefore both residuals are present here.
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, RMS_TKE, RMS_DISSIPATION, \
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, RMS_TKE, RMS_DISSIPATION, TOTAL_HEATFLUX, \
                 RMS_ADJ_DENSITY, RMS_ADJ_ENERGY, RMS_ADJ_TKE, RMS_ADJ_DISSIPATION)
diff --git a/TestCases/ddes/flatplate/ddes_flatplate.cfg b/TestCases/ddes/flatplate/ddes_flatplate.cfg
@@ -226,4 +226,4 @@ SURFACE_ADJ_FILENAME= surface_adjoint
 OUTPUT_WRT_FREQ= 1000
 %
 % Screen output
-SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/incomp_navierstokes/bend/lam_bend.cfg b/TestCases/incomp_navierstokes/bend/lam_bend.cfg
@@ -255,5 +255,4 @@ SURFACE_ADJ_FILENAME= surface_adjoint
 % Writing solution file frequency
 OUTPUT_WRT_FREQ= 1000
 %
-%
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/incomp_navierstokes/cylinder/incomp_cylinder.cfg b/TestCases/incomp_navierstokes/cylinder/incomp_cylinder.cfg
@@ -181,7 +181,5 @@ SURFACE_FILENAME= surface_flow
 % Writing solution file frequency
 OUTPUT_WRT_FREQ= 100
 %
-%
 % Screen output fields
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)
-
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/incomp_navierstokes/cylinder/poly_cylinder.cfg b/TestCases/incomp_navierstokes/cylinder/poly_cylinder.cfg
@@ -287,4 +287,4 @@ OUTPUT_WRT_FREQ= 250
 %
 %
 % Screen output fields
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TEMPERATURE, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TEMPERATURE, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/incomp_navierstokes/streamwise_periodic/chtPinArray_3d/configMaster.cfg b/TestCases/incomp_navierstokes/streamwise_periodic/chtPinArray_3d/configMaster.cfg
@@ -22,7 +22,7 @@ OUTER_ITER = 31
 %
 CONV_RESIDUAL_MINVAL= -26
 %
-SCREEN_OUTPUT= (OUTER_ITER, BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], STREAMWISE_MASSFLOW[0], STREAMWISE_DP[0], AVG_TEMPERATURE[1] )
+SCREEN_OUTPUT= (OUTER_ITER, BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], STREAMWISE_MASSFLOW[0], STREAMWISE_DP[0], AVG_TEMPERATURE[1], TOTAL_HEATFLUX[0] )
 SCREEN_WRT_FREQ_OUTER= 100
 %
 HISTORY_OUTPUT= ( ITER, RMS_RES[0], RMS_RES[1], STREAMWISE_PERIODIC[0], FLOW_COEFF[0], HEAT[1] )
diff --git a/TestCases/incomp_rans/naca0012/naca0012.cfg b/TestCases/incomp_rans/naca0012/naca0012.cfg
@@ -242,5 +242,5 @@ OUTPUT_WRT_FREQ= 100
 %
 %
 % Screen output fields
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
 
diff --git a/TestCases/incomp_rans/naca0012/naca0012_SST_SUST.cfg b/TestCases/incomp_rans/naca0012/naca0012_SST_SUST.cfg
@@ -242,5 +242,5 @@ OUTPUT_WRT_FREQ= 100
 %
 %
 % Screen output fields
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, TOTAL_HEATFLUX)
 
diff --git a/TestCases/navierstokes/cylinder/cylinder_lowmach.cfg b/TestCases/navierstokes/cylinder/cylinder_lowmach.cfg
@@ -218,4 +218,4 @@ OUTPUT_WRT_FREQ= 250
 %
 %
 % Screen output fields
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/navierstokes/cylinder/lam_cylinder.cfg b/TestCases/navierstokes/cylinder/lam_cylinder.cfg
@@ -212,4 +212,4 @@ OUTPUT_WRT_FREQ= 250
 %
 % 
 % Screen output fields
-SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)
+SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/navierstokes/poiseuille/lam_poiseuille.cfg b/TestCases/navierstokes/poiseuille/lam_poiseuille.cfg
@@ -249,6 +249,6 @@ OUTPUT_WRT_FREQ= 500
 %
 %
 % Screen output 
-SCREEN_OUTPUT=(INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)
+SCREEN_OUTPUT=(INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG, TOTAL_HEATFLUX)
 %
 OUTPUT_FILES=(RESTART_ASCII)
diff --git a/TestCases/nonequilibrium/axi_visccone/axi_visccone.cfg b/TestCases/nonequilibrium/axi_visccone/axi_visccone.cfg
@@ -4,151 +4,76 @@
 % Case description: Mach 5 viscous flow over a 10deg wedge                     %
 % Author: C. Garbacz                                                           %
 % Institution: Strathclyde University                                          %
-% Date: 2020.11.01                                                             %
 % File Version 7.3.1 "Blackbird"                                               %
 %                                                                              %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 % ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
 %
-% Physical governing equations (EULER, NAVIER_STOKES,
-%                               WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
-%                               POISSON_EQUATION)
 SOLVER= NEMO_NAVIER_STOKES
-%
 GAS_MODEL= AIR-5
-%
 GAS_COMPOSITION= (0.77, 0.23, 0.0, 0.0, 0.0)
-%
-% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT)
 MATH_PROBLEM= DIRECT
-%
-% Restart solution (NO, YES)
 RESTART_SOL= NO
-%
-% Axisymmetric simulation, only compressible flows (NO, YES)
 AXISYMMETRIC= YES
 
 % ----------- COMPRESSIBLE AND INCOMPRESSIBLE FREE-STREAM DEFINITION ----------%
 %
-% Mach number (non-dimensional, based on the free-stream values)
 MACH_NUMBER= 5
-%
-% Angle of attack (degrees)
 AOA= 0.0
-%
-% Side-slip angle (degrees)
 SIDESLIP_ANGLE= 0.0
-%
-% Init option to choose between Reynolds (default) or thermodynamics quantities
-% for initializing the solution (REYNOLDS, TD_CONDITIONS)
 INIT_OPTION= TD_CONDITIONS
-%
-% Free-stream pressure (101325.0 N/m^2 by default, only Euler flows)  
 FREESTREAM_PRESSURE= 10.0
-%
-% Free-stream temperature (288.15 K by default)
 FREESTREAM_TEMPERATURE= 288.15
 FREESTREAM_TEMPERATURE_VE= 288.15
-%
-% Reynolds number (non-dimensional, based on the free-stream values)
 REYNOLDS_NUMBER= 10000
 KIND_TURB_MODEL= NONE
 
 % ---- NONEQUILIBRIUM GAS, IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
 %
-% Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
-%              CONSTANT_DENSITY, INC_IDEAL_GAS, INC_IDEAL_GAS_POLY, MUTATIONPP, SU2_NONEQ)
 FLUID_MODEL= SU2_NONEQ
 
 % -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
 %
-% Euler wall boundary marker(s) (NONE = no marker)
 MARKER_EULER= ( Euler)
 MARKER_ISOTHERMAL= (Wall, 300)
 MARKER_OUTLET= ( Exit, 5 )
 MARKER_FAR= ( Farfield, Inlet )
-%
-% Marker(s) of the surface to be plotted or designed
 MARKER_PLOTTING= (NONE )
-%
-% Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated
 MARKER_MONITORING= ( Wall )
 
 % ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
 %
-% Numerical method for spatial gradients (GREEN_GAUSS, LEAST_SQUARES, 
-%                                         WEIGHTED_LEAST_SQUARES)
 NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
-%
-% Courant-Friedrichs-Lewy condition of the finest grid
 CFL_NUMBER= 3.0
-%
-% Number of total iterations
 ITER= 11
-%
-% Linear solver for the implicit formulation (BCGSTAB, FGMRES)
 LINEAR_SOLVER= BCGSTAB
-%
-% Min error of the linear solver for the implicit formulation
 LINEAR_SOLVER_ERROR= 1E-6
-%
-% Max number of iterations of the linear solver for the implicit formulation
 LINEAR_SOLVER_ITER= 5
 
 % -----------------------------------------------------------------------%
 %
-% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
-%                              TURKEL_PREC, MSW)
 CONV_NUM_METHOD_FLOW= AUSM
-%
-% Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER)
-%
 MUSCL_FLOW= YES
-%
-% Slope limiter (VENKATAKRISHNAN, MINMOD)
 SLOPE_LIMITER_FLOW= VENKATAKRISHNAN
-%
-% Coefficient for the Venkat's limiter (upwind scheme). A larger values decrease
-%             the extent of limiting, values approaching zero cause
-%             lower-order approximation to the solution (0.05 by default)
 VENKAT_LIMITER_COEFF= 0.05
-%
-% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
+
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
-% Min value of the residual (log10 of the residual)
 CONV_RESIDUAL_MINVAL= -15
-%
-% Start convergence criteria at iteration number
 CONV_STARTITER= 10
-%
+
 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
-% Mesh input file
 MESH_FILENAME= viscwedge.su2
-%
-% Mesh input file format (SU2, CGNS)
 MESH_FORMAT= SU2
-%
-% Restart flow input file
 SOLUTION_FILENAME= restart_flow.dat  
-%
 TABULAR_FORMAT= TECPLOT
-%
-% Output file convergence history (w/o extension) 
 CONV_FILENAME= convergence
-%
-% Output file restart flow
 RESTART_FILENAME= restart_flow.dat
-%
-% Output file flow (w/o extension) variables
 VOLUME_FILENAME= soln_volume
-%
-% Output file surface flow coefficient (w/o extension)
 SURFACE_FILENAME= soln_surface     
-%
 OUTPUT_WRT_FREQ= 100
-%
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY_0, RMS_DENSITY_1, RMS_DENSITY_2, RMS_DENSITY_3, RMS_DENSITY_4, RMS_ENERGY, RMS_ENERGY_VE, LIFT, DRAG, TOTAL_HEATFLUX)
 OUTPUT_FILES= (RESTART_ASCII, PARAVIEW_ASCII)
diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py
@@ -97,8 +97,8 @@ def main():
     visc_cone.cfg_dir   = "nonequilibrium/axi_visccone"
     visc_cone.cfg_file  = "axi_visccone.cfg"
     visc_cone.test_iter = 10
-    visc_cone.test_vals         = [-5.222278, -5.746529, -20.569425, -20.633787, -20.547644, -1.928717, -2.247306, 1.255759, -3.208374, 0.0]
-    visc_cone.test_vals_aarch64 = [-5.222267, -5.746522, -20.569408, -20.633783, -20.546393, -1.928714, -2.247317, 1.255761, -3.208360, 0.0]
+    visc_cone.test_vals         = [-5.222278, -5.746529, -20.569425, -20.633787, -20.547644, -1.928717, -2.247306, 1.255759, -3.208374, 32633.000000]
+    visc_cone.test_vals_aarch64 = [-5.222267, -5.746522, -20.569408, -20.633783, -20.546393, -1.928714, -2.247317, 1.255761, -3.208360, 32633.000000]
     visc_cone.su2_exec  = "mpirun -n 2 SU2_CFD"
     visc_cone.timeout   = 1600
     visc_cone.new_output = True
@@ -1425,8 +1425,8 @@ def main():
     solid_periodic_pins.cfg_dir   = "solid_heat_conduction/periodic_pins"
     solid_periodic_pins.cfg_file  = "configSolid.cfg"
     solid_periodic_pins.test_iter = 750
-    solid_periodic_pins.test_vals         = [-15.878958, -14.569206, 300.900000, 425.320000, 0.000000, 5.000000, -1.672714, 0.0] #last 7 lines
-    solid_periodic_pins.test_vals_aarch64 = [-15.879010, -14.569206, 300.900000, 425.320000, 0.000000, 5.000000, -1.672630, 0.0] #last 7 lines
+    solid_periodic_pins.test_vals         = [-15.878958, -14.569206, 300.900000, 425.320000, 0.000000, 5.000000, -1.672714] #last 7 lines
+    solid_periodic_pins.test_vals_aarch64 = [-15.879010, -14.569206, 300.900000, 425.320000, 0.000000, 5.000000, -1.672630] #last 7 lines
     solid_periodic_pins.su2_exec  = "mpirun -n 2 SU2_CFD"
     solid_periodic_pins.timeout   = 1600
     solid_periodic_pins.tol       = 0.00001
@@ -1441,7 +1441,7 @@ def main():
     cht_incompressible.cfg_dir   = "coupled_cht/incomp_2d"
     cht_incompressible.cfg_file  = "cht_2d_3cylinders.cfg"
     cht_incompressible.test_iter = 10
-    cht_incompressible.test_vals = [-2.128826, -0.588813, -0.588813, -0.588813, 0.0] #last 4 columns
+    cht_incompressible.test_vals = [-2.128826, -0.588813, -0.588813, -0.588813] #last 4 columns
     cht_incompressible.su2_exec  = "SU2_CFD"
     cht_incompressible.timeout   = 1600
     cht_incompressible.multizone = True
@@ -1453,7 +1453,7 @@ def main():
     cht_compressible.cfg_dir   = "coupled_cht/comp_2d"
     cht_compressible.cfg_file  = "cht_2d_3cylinders.cfg"
     cht_compressible.test_iter = 10
-    cht_compressible.test_vals = [-4.256032, -0.532728, -0.532729, -0.532728]
+    cht_compressible.test_vals = [-4.256032, -0.532728, -0.532729, -0.532728, 0.0]
     cht_compressible.su2_exec  = "SU2_CFD"
     cht_compressible.timeout   = 1600
     cht_compressible.multizone = True
diff --git a/TestCases/radiation/p1model/configp1.cfg b/TestCases/radiation/p1model/configp1.cfg
@@ -96,7 +96,7 @@ CONV_CAUCHY_EPS= 1E-6
 % OUTPUT
 %%%%%%%%%%%%%%%%%%%%%%%
 
-SCREEN_OUTPUT = (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, RMS_TEMPERATURE, RMS_RAD_ENERGY)
+SCREEN_OUTPUT = (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, RMS_TEMPERATURE, RMS_RAD_ENERGY, TOTAL_HEATFLUX)
 
 MESH_FORMAT= SU2
 MESH_FILENAME= meshp1.su2
diff --git a/TestCases/rans/naca0012/turb_NACA0012_sst_2003m.cfg b/TestCases/rans/naca0012/turb_NACA0012_sst_2003m.cfg
@@ -92,4 +92,4 @@ VOLUME_FILENAME= flow
 SURFACE_FILENAME= surface_flow
 OUTPUT_WRT_FREQ= 10000
 %
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, LINSOL_RESIDUAL)
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, LINSOL_RESIDUAL, TOTAL_HEATFLUX)
diff --git a/TestCases/rans/oneram6/turb_ONERAM6.cfg b/TestCases/rans/oneram6/turb_ONERAM6.cfg
@@ -283,4 +283,4 @@ OUTPUT_WRT_FREQ= 250
 %
 %
 % Screen output fields
-SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/rans/rae2822/turb_SA_RAE2822.cfg b/TestCases/rans/rae2822/turb_SA_RAE2822.cfg
@@ -251,7 +251,7 @@ OUTPUT_WRT_FREQ= 250
 %
 %
 % Screen output fields 
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
 
 % ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
 %
diff --git a/TestCases/rans/rae2822/turb_SST_RAE2822.cfg b/TestCases/rans/rae2822/turb_SST_RAE2822.cfg
@@ -227,5 +227,5 @@ SURFACE_ADJ_FILENAME= surface_adjoint
 OUTPUT_WRT_FREQ= 250
 %
 % Screen output fields
-SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG)
+SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, TOTAL_HEATFLUX)
 
diff --git a/TestCases/serial_regression.py b/TestCases/serial_regression.py
@@ -1457,7 +1457,7 @@ def main():
     p1rad.cfg_file  = "configp1.cfg"
     p1rad.new_output= True
     p1rad.test_iter = 100
-    p1rad.test_vals = [-7.751309, -7.923059, -2.119084, 0.091733] #last 4 columns
+    p1rad.test_vals = [-7.751309, -7.923059, -2.119084, 0.091733, 0.0] #last 4 columns
     p1rad.su2_exec  = "SU2_CFD"
     p1rad.timeout   = 1600
     p1rad.tol       = 0.00001
@@ -1480,7 +1480,7 @@ def main():
     cht_incompressible.cfg_dir   = "coupled_cht/incomp_2d"
     cht_incompressible.cfg_file  = "cht_2d_3cylinders.cfg"
     cht_incompressible.test_iter = 10
-    cht_incompressible.test_vals = [-2.128826, -0.588813, -0.588813, -0.588813] #last 4 columns
+    cht_incompressible.test_vals = [-2.128826, -0.588813, -0.588813, -0.588813, 0.0] #last 4 columns
     cht_incompressible.su2_exec  = "SU2_CFD"
     cht_incompressible.timeout   = 1600
     cht_incompressible.multizone = True
@@ -1492,7 +1492,7 @@ def main():
     cht_incompressible.cfg_dir   = "coupled_cht/comp_2d"
     cht_incompressible.cfg_file  = "cht_2d_3cylinders.cfg"
     cht_incompressible.test_iter = 10
-    cht_incompressible.test_vals = [-4.256032, -0.532728, -0.532729, -0.532728]
+    cht_incompressible.test_vals = [-4.256032, -0.532728, -0.532729, -0.532728, 0.0]
     cht_incompressible.su2_exec  = "SU2_CFD"
     cht_incompressible.timeout   = 1600
     cht_incompressible.multizone = True
diff --git a/TestCases/transition/Schubauer_Klebanoff/transitional_BC_model_ConfigFile.cfg b/TestCases/transition/Schubauer_Klebanoff/transitional_BC_model_ConfigFile.cfg
@@ -266,6 +266,5 @@ SURFACE_ADJ_FILENAME= surface_adjoint
 % Writing solution file frequency
 OUTPUT_WRT_FREQ= 500
 %
-%
 % Screen output 
-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/turbulence_models/sa/rae2822/turb_SA_NEG_RAE2822.cfg b/TestCases/turbulence_models/sa/rae2822/turb_SA_NEG_RAE2822.cfg
@@ -185,4 +185,4 @@ SURFACE_FILENAME= surface_flow
 OUTPUT_WRT_FREQ= 250
 %
 % Screen output fields 
-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_NU_TILDE, LIFT, DRAG)
+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)
diff --git a/TestCases/unsteady/square_cylinder/turb_square.cfg b/TestCases/unsteady/square_cylinder/turb_square.cfg
@@ -254,6 +254,6 @@ OUTPUT_WRT_FREQ= 1
 %
 % 
 % Screen output 
-SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, TAVG_LIFT, TAVG_DRAG)
+SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_TKE, RMS_DISSIPATION, LIFT, DRAG, TAVG_LIFT, TAVG_DRAG, TOTAL_HEATFLUX)
 %
 OUTPUT_FILES= (RESTART_ASCII)
diff --git a/TestCases/wallfunctions/flatplate/incompressible_SST/turb_SST_flatplate.cfg b/TestCases/wallfunctions/flatplate/incompressible_SST/turb_SST_flatplate.cfg
@@ -99,7 +99,7 @@ MESH_FILENAME= ../mesh_flatplate_140x100.su2
 MESH_FORMAT= SU2
 %
 SCREEN_OUTPUT= (INNER_ITER, WALL_TIME, RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TEMPERATURE, RMS_TKE, RMS_DISSIPATION, \
-                LINSOL_ITER, LINSOL_RESIDUAL, LIFT, DRAG)
+                LINSOL_ITER, LINSOL_RESIDUAL, LIFT, DRAG, TOTAL_HEATFLUX)
 SCREEN_WRT_FREQ_INNER= 10
 %
 HISTORY_OUTPUT= (ITER, RMS_RES, AERO_COEFF, FLOW_COEFF)

Original file line number	Diff line number	Diff line change
`@@ -101,5 +101,5 @@ SOLUTION_FILENAME= solution_flow.dat`
`101`	`101`	`OUTPUT_WRT_FREQ= 1000`
`102`	`102`	`%`
`103`	`103`	`% Note: This cfg is used for a primal and adjoint Testcase, therefore both residuals are present here.`
`104`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, RMS_TKE, RMS_DISSIPATION, \`
	`104`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, RMS_TKE, RMS_DISSIPATION, TOTAL_HEATFLUX, \`
`105`	`105`	`RMS_ADJ_DENSITY, RMS_ADJ_ENERGY, RMS_ADJ_TKE, RMS_ADJ_DISSIPATION)`
Original file line number	Diff line number	Diff line change
`@@ -226,4 +226,4 @@ SURFACE_ADJ_FILENAME= surface_adjoint`
`226`	`226`	`OUTPUT_WRT_FREQ= 1000`
`227`	`227`	`%`
`228`	`228`	`% Screen output`
`229`		`-SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG)`
	`229`	`+SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)`
Original file line number	Diff line number	Diff line change
`@@ -255,5 +255,4 @@ SURFACE_ADJ_FILENAME= surface_adjoint`
`255`	`255`	`% Writing solution file frequency`
`256`	`256`	`OUTPUT_WRT_FREQ= 1000`
`257`	`257`	`%`
`258`		`-%`
`259`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)`
	`258`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG, TOTAL_HEATFLUX)`
Original file line number	Diff line number	Diff line change
`@@ -181,7 +181,5 @@ SURFACE_FILENAME= surface_flow`
`181`	`181`	`% Writing solution file frequency`
`182`	`182`	`OUTPUT_WRT_FREQ= 100`
`183`	`183`	`%`
`184`		`-%`
`185`	`184`	`% Screen output fields`
`186`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG)`
`187`		`-`
	`185`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_VELOCITY-X, LIFT, DRAG, TOTAL_HEATFLUX)`
Original file line number	Diff line number	Diff line change
`@@ -287,4 +287,4 @@ OUTPUT_WRT_FREQ= 250`
`287`	`287`	`%`
`288`	`288`	`%`
`289`	`289`	`% Screen output fields`
`290`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TEMPERATURE, LIFT, DRAG)`
	`290`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_TEMPERATURE, LIFT, DRAG, TOTAL_HEATFLUX)`
Original file line number	Diff line number	Diff line change
`@@ -22,7 +22,7 @@ OUTER_ITER = 31`
`22`	`22`	`%`
`23`	`23`	`CONV_RESIDUAL_MINVAL= -26`
`24`	`24`	`%`
`25`		`-SCREEN_OUTPUT= (OUTER_ITER, BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], STREAMWISE_MASSFLOW[0], STREAMWISE_DP[0], AVG_TEMPERATURE[1] )`
	`25`	`+SCREEN_OUTPUT= (OUTER_ITER, BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], STREAMWISE_MASSFLOW[0], STREAMWISE_DP[0], AVG_TEMPERATURE[1], TOTAL_HEATFLUX[0] )`
`26`	`26`	`SCREEN_WRT_FREQ_OUTER= 100`
`27`	`27`	`%`
`28`	`28`	`HISTORY_OUTPUT= ( ITER, RMS_RES[0], RMS_RES[1], STREAMWISE_PERIODIC[0], FLOW_COEFF[0], HEAT[1] )`
Original file line number	Diff line number	Diff line change
`@@ -242,5 +242,5 @@ OUTPUT_WRT_FREQ= 100`
`242`	`242`	`%`
`243`	`243`	`%`
`244`	`244`	`% Screen output fields`
`245`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG)`
	`245`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_PRESSURE, RMS_NU_TILDE, LIFT, DRAG, TOTAL_HEATFLUX)`
`246`	`246`
Original file line number	Diff line number	Diff line change
`@@ -218,4 +218,4 @@ OUTPUT_WRT_FREQ= 250`
`218`	`218`	`%`
`219`	`219`	`%`
`220`	`220`	`% Screen output fields`
`221`		`-SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)`
	`221`	`+SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG, TOTAL_HEATFLUX)`
Original file line number	Diff line number	Diff line change
`@@ -212,4 +212,4 @@ OUTPUT_WRT_FREQ= 250`
`212`	`212`	`%`
`213`	`213`	`%`
`214`	`214`	`% Screen output fields`
`215`		`-SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG)`
	`215`	`+SCREEN_OUTPUT = (INNER_ITER, RMS_DENSITY, RMS_ENERGY, LIFT, DRAG, TOTAL_HEATFLUX)`