fix: adding logger to intrusion body

added bool flag for intrusion import

Author: Lachlan Grose

LoopStructural/modelling/core/geological_model.py

@@ -29,11 +29,14 @@
 from LoopStructural.utils.helper import (all_heading, gradient_vec_names,
                                          strike_dip_vector)
 
-from LoopStructural.modelling.intrusions import IntrusionNetwork
-from LoopStructural.modelling.intrusions import IntrusionBuilder
-from LoopStructural.modelling.intrusions import IntrusionBody
-from LoopStructural.modelling.intrusions import IntrusionFeature 
-
+intrusions = True
+try:
+    from LoopStructural.modelling.intrusions import IntrusionNetwork
+    from LoopStructural.modelling.intrusions import IntrusionBuilder
+    from LoopStructural.modelling.intrusions import IntrusionBody
+    from LoopStructural.modelling.intrusions import IntrusionFeature 
+except ImportError:
+     intrusions = False
 from LoopStructural.utils import getLogger, log_to_file
 logger = getLogger(__name__)
 
@@ -1011,6 +1014,9 @@ def create_and_add_intrusion(self, intrusion_name, intrusion_frame_name = None,
         intrusion feature
         
         """
+        if intrusions == False:
+            logger.error("Libraries not installed")
+            raise Exception("Libraries not installed")
         feature_data = self.data[self.data['feature_name'] == intrusion_name].copy()
 
         # Create and build Intrusion Network

LoopStructural/modelling/intrusions/intrusion_body.py

@@ -1,6 +1,8 @@
 import numpy as np
 import pandas as pd
 
+from LoopStructural.utils import getLogger
+logger = getLogger(__name__)
 #GSLIB library        
 try:
     import geostatspy.GSLIB as GSLIB                          # GSLIB utilities, viz and wrapped functions
@@ -516,6 +518,8 @@ def simulate_s_thresholds(self):
         ndmin = self.simulation_s_parameters.get('ndmin'); ndmax = self.simulation_s_parameters.get('ndmax')         
         radius = self.simulation_s_parameters.get('radius') 
 
+        if geostats is None:
+            raise Exception('geostats is not installed')    
         s_min_simulation = geostats.sgsim(inputsimdata_minS,'coord1','ref_coord','s_residual',wcol=-1,scol=-1,
                                           tmin=tmin,tmax=tmax,itrans=itrans,ismooth=0,dftrans=0,tcol=0,twtcol=0,
                                           zmin=zmin,zmax=zmax,ltail=1,ltpar=0.0,utail=1,utpar=0.3,nsim=1,