My Code
for WS in WSlistCVT_Pa: TW = q_cruise_Pa*(CDmin/WSlistCVT_Pa[i] + WSlistCVT_Pa[i]\\ *k*(n_cvt_/q_cruise_Pa)**2)