diff --git a/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstGFS_HofX_obsIODAv2_PyEmbed/read_iodav2_mpr.py b/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstGFS_HofX_obsIODAv2_PyEmbed/read_iodav2_mpr.py
index de79edf64a..84a0140fb3 100644
--- a/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstGFS_HofX_obsIODAv2_PyEmbed/read_iodav2_mpr.py
+++ b/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstGFS_HofX_obsIODAv2_PyEmbed/read_iodav2_mpr.py
@@ -83,7 +83,7 @@
                     var_name+'@ObsType','na','na','lead','na','na','na','na','MPR',
                     'nobs','index','station_id','latitude','longitude',
                     'air_pressure','na',var_name+'@hofx',var_name+'@ObsValue',
-                    var_name+'@EffectiveQC','na','na']
+                    var_name+'@EffectiveQC','na','na','na','na','na']
             
             ioda_df_var = ioda_df_var[cols]
 
diff --git a/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstHAFS_obsPrepBufr_JEDI_IODA_interface/read_ioda_mpr.py b/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstHAFS_obsPrepBufr_JEDI_IODA_interface/read_ioda_mpr.py
index 538c364cca..b062697bf7 100644
--- a/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstHAFS_obsPrepBufr_JEDI_IODA_interface/read_ioda_mpr.py
+++ b/parm/use_cases/model_applications/data_assimilation/StatAnalysis_fcstHAFS_obsPrepBufr_JEDI_IODA_interface/read_ioda_mpr.py
@@ -83,7 +83,7 @@ def read_netcdfs(files, dim):
                     var_name+'@ObsType','na','na','lead','na','na','na','na','MPR',
                     'nobs','index','station_id@MetaData','latitude@MetaData','longitude@MetaData',
                     'air_pressure@MetaData','na',var_name+'@hofx',var_name+'@ObsValue',
-                    var_name+'@PreQC','na','na']
+                    var_name+'@PreQC','na','na','na','na','na']
             
             ioda_df_var = ioda_df_var[cols]
 
diff --git a/parm/use_cases/model_applications/unstructured_grids/StatAnalysis_fcstLFRIC_UGRID_obsASCII_PyEmbed/ugrid_lfric_mpr.py b/parm/use_cases/model_applications/unstructured_grids/StatAnalysis_fcstLFRIC_UGRID_obsASCII_PyEmbed/ugrid_lfric_mpr.py
index 78d1994208..7e9c2c78fa 100644
--- a/parm/use_cases/model_applications/unstructured_grids/StatAnalysis_fcstLFRIC_UGRID_obsASCII_PyEmbed/ugrid_lfric_mpr.py
+++ b/parm/use_cases/model_applications/unstructured_grids/StatAnalysis_fcstLFRIC_UGRID_obsASCII_PyEmbed/ugrid_lfric_mpr.py
@@ -185,7 +185,7 @@ def to_centers(mesh: pv.PolyData) -> pv.PolyData:
                 'var_name','na','na','lead','na','na','na','na','MPR',
                 'nobs','index','station_id','obs_lat','obs_lon',
                 'level','na','fcst_value','obs_value',
-                'QC','na','na']
+                'QC','na','na','na','na','na']
         
         match_df = match_df[cols]