(FEATURE) adding argument to skip variogram

If there aren't enough points 1) variogram is pointless
2) the variogram will error. So adding ability to skip it
probably can add a check as well

Author: lachlangrose

LoopStructural/modelling/core/geological_model.py

@@ -666,7 +666,7 @@ def create_and_add_folded_foliation(self, foliation_data, fold_frame=None,
             # allow for predefined functions to be used
             fold_limb_rotation.set_function(kwargs['limb_function'])
         else:
-            fold_limb_rotation.fit_fourier_series(wl=l_wl)
+            fold_limb_rotation.fit_fourier_series(wl=l_wl,**kwargs)
         fold.fold_limb_rotation = fold_limb_rotation
         # fold_limb_fitter = kwargs.get("fold_limb_function",
         # _interpolate_fold_limb_rotation_angle)

LoopStructural/modelling/fold/fold_rotation_angle.py

@@ -31,7 +31,7 @@ def __init__(self, rotation_angle, fold_frame_coordinate, svario=False):
             self.svario = SVariogram(self.fold_frame_coordinate, self.rotation_angle)
         self.fitted_params = None
 
-    def fit_fourier_series(self, wl = None, lags = None, nlag = None, lag = None):
+    def fit_fourier_series(self, wl = None, lags = None, nlag = None, lag = None, skip_variogram=False,**kwargs):
         """
 
         Parameters
@@ -48,8 +48,8 @@ def fit_fourier_series(self, wl = None, lags = None, nlag = None, lag = None):
         if self.svario is None:
             self.svario = SVariogram(self.fold_frame_coordinate,
                                      self.rotation_angle)
-
-        self.svario.calc_semivariogram(lags=lags, nlag=nlag, lag=lag)
+        if skip_variogram != False:
+            self.svario.calc_semivariogram(lags=lags, nlag=nlag, lag=lag)
         if wl is None:
             wl = self.svario.find_wavelengths(lags=lags, nlag=nlag, lag=lag)
             # for now only consider single fold wavelength