# Run parasite drag solver
PD_analysis = "ParasiteDrag"
vsp.SetAnalysisInputDefaults(PD_analysis)
vsp.SetIntAnalysisInput(PD_analysis, "RecomputeGeom", (0,)) # do not recompute
vsp.SetIntAnalysisInput(PD_analysis, "LengthUnit", (2,)) # meters
vsp.SetDoubleAnalysisInput(PD_analysis, "Sref", (Sref,))
vsp.SetIntAnalysisInput(PD_analysis, "FreestreamPropChoice", (5,)) # Re/L + Mach control
vsp.SetDoubleAnalysisInput(PD_analysis, "SpecificHeatRatio", (gamma_PD,))
vsp.SetDoubleAnalysisInput(PD_analysis, "Re_L", (Re_l_PD,))
vsp.SetDoubleAnalysisInput(PD_analysis, "Mach", (Machinf_PD,))
vsp.PrintAnalysisInputs(PD_analysis)
vsp.Update()
PD_results = vsp.ExecAnalysis(PD_analysis)
#vsp.PrintResults(PD_results)
vsp.WriteResultsCSVFile(PD_results, fname_csvdata_PD)
cD_tot = vsp.GetDoubleResults(PD_results, "Total_CD_Total", 0)
print(f"\nCD_tot = {cD_tot[0]}\n")