Fix for missing fgd values
The fgd data file I'm using has some missing plants.
@@ -606,6 +606,9 @@ def atb_fixed_var_om_existing(
fgd = coal_fgd_df.query(
"plant_id_eia == @plant_id & generator_id in @gen_ids"
)["fgd"].values
+ if not np.any(fgd):
+ # If FGD isn't found, use average of with/without FGD
+ fgd = np.ones_like(age) * 0.5
# https://www.eia.gov/analysis/studies/powerplants/generationcost/pdf/full_report.pdf
annual_capex = (16.53 + (0.126 * age) + (5.68 * fgd)) * 1000
