docs: def compute_center_of_mass(self, filename):
          """
          Center of mass is computed using the sensitivity data output from CRMod
          Data weights can be applied using command line options
          """
          sens = np.loadtxt(filename, skiprows=1)
  
          X = sens[:, 0]
          Z = sens[:, 1]
          # C = (np.abs(sens[:,2]))# ./ np.max(np.abs(sens[:,2]))
          C = sens[:, 2]
  
          x_center = 0
          z_center = 0
          sens_sum = 0
  
          for i in range(0, C.shape[0]):
              # unweighted
              if(self.weight == 0):
                  weight = (C[i])
              # abs
              if(self.weight == 1):
                  weight = np.abs(C[i])
              # log10
              if(self.weight == 2):
                  weight = np.log10(np.abs(C[i]))
              # sqrt
              if(self.weight == 3):
                  weight = np.sqrt(np.abs(C[i]))
  
              x_center += (X[i] * weight)
              z_center += (Z[i] * weight)
              sens_sum += weight
          x_center /= sens_sum
          z_center /= sens_sum
  
          return (x_center, z_center)