Merge pull request #46 from MTgeophysics/simpeg

For now only 1D and 2D are "working" but need to include in next release.

Author: kujaku11

.github/workflows/testing_in_conda.yml

@@ -23,21 +23,26 @@ jobs:
       uses: s-weigand/setup-conda@v1
       with:
         update-conda: false
+        conda-channels: conda-forge
         python-version: ${{ matrix.python-version }}
 
     - name: Install Env
       shell: bash
       run: |
         python --version
-        conda create -n mtpy-v2-test python=${{ matrix.python-version }}
+        conda env create -f environment.yml
         source activate mtpy-v2-test
+        pip install git+https://github.com/simpeg/pydiso.git
         conda install pytest
         conda install pytest-subtests
         conda install pytest-cov
         pip install git+https://github.com/kujaku11/mt_metadata.git@main
-        pip install git+https://github.com/kujaku11/mth5.git@add_aurora_tools
-        pip install git+https://github.com/MTgeophysics/mtpy_data.git@main
+        pip install git+https://github.com/kujaku11/mth5.git@patches
         pip install git+https://github.com/simpeg/aurora@move_to_mth5
+        git clone https://github.com/MTgeophysics/mtpy_data.git
+        cd mtpy_data
+        pip install -e .
+        cd ..
         
     - name: Install Our Package
       shell: bash

environment.yml

@@ -0,0 +1,18 @@
+name: mtpy-v2-test
+channels:
+  - conda-forge
+  - defaults
+dependencies:
+  - matplotlib
+  - numpy=1.26
+  - git
+  - pip
+  - simpeg
+  - dill
+  - jupyterlab
+  - ipywidgets
+  - mkl-devel
+  - pkg-config
+  - cmake
+  - cython 
+  - ninja

mtpy/core/mt.py

@@ -94,6 +94,7 @@ def clone_empty(self):
         new_mt_obj.model_north = self.model_north
         new_mt_obj.model_elevation = self.model_elevation
         new_mt_obj._rotation_angle = self._rotation_angle
+        new_mt_obj.profile_offset = self.profile_offset
 
         return new_mt_obj
 
@@ -958,9 +959,7 @@ def add_white_noise(self, value, inplace=True):
             ] = self._transfer_function.transfer_function.real * (
                 noise_real
             ) + (
-                1j
-                * self._transfer_function.transfer_function.imag
-                * noise_imag
+                1j * self._transfer_function.transfer_function.imag * noise_imag
             )
 
             self._transfer_function["transfer_function_error"] = (
@@ -974,16 +973,20 @@ def add_white_noise(self, value, inplace=True):
             ] = self._transfer_function.transfer_function.real * (
                 noise_real
             ) + (
-                1j
-                * self._transfer_function.transfer_function.imag
-                * noise_imag
+                1j * self._transfer_function.transfer_function.imag * noise_imag
             )
 
             self._transfer_function["transfer_function_error"] = (
                 self._transfer_function.transfer_function_error + value
             )
             return new_mt_obj
 
+    def edit_curve(self, method="default", tolerance=0.05):
+        """
+        try to remove bad points in a scientific way.
+
+        """
+
     def to_occam1d(self, data_filename=None, mode="det"):
         """Write an Occam1DData data file.
         :param data_filename: Path to write file, if None returns Occam1DData