diff --git a/getDataFRF.py b/getDataFRF.py
index eb7fe77..cf04f15 100644
--- a/getDataFRF.py
+++ b/getDataFRF.py
@@ -2129,6 +2129,8 @@ def getLidarDEM(self, **kwargs):
                 'lat':       self.ncfile['latitude'][ys, xs],
                 'lon':       self.ncfile['longitude'][ys,xs]
                 }
+        DEMdata['elevation_mean'] = np.nanmean(DEMdata['elevation'], axis=0)
+        DEMdata['elevation_var'] = np.nanvar(DEMdata['elevation'], axis=0)
         return DEMdata
     
     def getBathyRegionalDEM(self, utmEmin, utmEmax, utmNmin, utmNmax):