diff --git a/properties.py b/properties.py
index be523b5..2273afc 100644
--- a/properties.py
+++ b/properties.py
@@ -60,6 +60,13 @@ def plot_ylabel(self):
         return r"Mass Flux $(M_\odot yr^{-1})$"
 
     def _get_profile(self, halo, maxrad):
+        halo.g['metals']
+        halo.s['metals']
+        halo.g['FeMassFrac']
+        halo.s['FeMassFrac']
+        halo.g['OxMassFrac']
+        halo.s['OxMassFrac']
+
         delta = self.plot_xdelta()
         nbins = int(maxrad / delta)
         maxrad = delta * nbins
@@ -186,6 +193,12 @@ def _get_profile(self, halo, maxrad):
         ox_pro_g = None
         ox_pro_s = None
 
+        halo.g['metals']
+        halo.s['metals']
+        halo.g['FeMassFrac']
+        halo.s['FeMassFrac']
+        halo.g['OxMassFrac']
+        halo.s['OxMassFrac']
         if len(halo.g)>10:
             pro_g = pynbody.analysis.profile.Profile(halo.g, type='lin', ndim=3,
                                                min=0, max=maxrad, nbins=nbins, weight='mass')
@@ -227,6 +240,13 @@ def _get_profile(self, halo, maxrad):
         fe_pro = None
         ox_pro = None
 
+        halo.g['metals']
+        halo.s['metals']
+        halo.g['FeMassFrac']
+        halo.s['FeMassFrac']
+        halo.g['OxMassFrac']
+        halo.s['OxMassFrac']
+
         if len(halo)>10:
             if 'massHot' in halo.loadable_keys():
                 twophase = pynbody.filt.HighPass('massHot', 0)