[zjj]

hello. i found when i deal with the unsteady shape optimization in naca0012 cases, the code can't run correctly with a thickness constraint. The message is shown as below.

Traceback (most recent call last):
File "/usr/local/bin/shape_optimization.py", line 183, in <module>
main()
File "/usr/local/bin/shape_optimization.py", line 92, in main
shape_optimization( options.filename ,
File "/usr/local/bin/shape_optimization.py", line 110, in shape_optimization
config = SU2.io.Config(filename)
File "/usr/local/bin/SU2/io/config.py", line 123, in __init__
tavg_keyGroup = "TAVG_" + historyOutFields[key_inner]["GROUP"]
KeyError: 'AIRFOIL_THICKNESS'


then i found the reason is that the 'if sentence' in below. thickness constraints is not in [historyOutFields]. who can tell me the reason ? thank you in advance!




if self.get("TIME_DOMAIN") == "YES":
objFuncsFields = self.get("OPT_OBJECTIVE")
histFields = self.get("HISTORY_OUTPUT")
diff_objective = self.get("OBJECTIVE_FUNCTION")
constrFuncFields = self.get("OPT_CONSTRAINT")

#OPT_OBJECTIVES
if bool (objFuncsFields):
for key in objFuncsFields:
tavg_keyGroup = "TAVG_" + historyOutFields[key]["GROUP"]
if not tavg_keyGroup in histFields:
histFields.append(tavg_keyGroup)

dtavg_keyGroup = "D_TAVG_" + historyOutFields[key]["GROUP"]
if not dtavg_keyGroup in histFields:
histFields.append(dtavg_keyGroup)

#OPT_CONSTRAINTS
if bool (constrFuncFields):
for key in constrFuncFields:
eqIneqConstrFunc = constrFuncFields.get(key)
for key_inner in eqIneqConstrFunc:
tavg_keyGroup = "TAVG_" + historyOutFields[key_inner]["GROUP"]
if not tavg_keyGroup in histFields:
histFields.append(tavg_keyGroup)