acstools 3.7.0__tar.gz → 3.7.2__tar.gz

acstools-3.7.2/.github/dependabot.yml

@@ -0,0 +1,28 @@
+version: 2
+updates:
+
+  # Maintain dependencies for GitHub Actions
+  - package-ecosystem: "github-actions"
+    directory: ".github/workflows"
+    target-branch: "master"
+    schedule:
+      interval: monthly
+      time: "03:10"
+    groups:
+      actions:
+        patterns:
+          - "*"
+    open-pull-requests-limit: 10
+
+  # Maintain dependencies for pip
+  - package-ecosystem: "pip"
+    directory: "/"
+    target-branch: "master"
+    schedule:
+      interval: monthly
+      time: "03:00"
+    groups:
+      actions:
+        patterns:
+          - "*"
+    open-pull-requests-limit: 10

acstools-3.7.2/.github/labeler.yml

@@ -0,0 +1,61 @@
+docs:
+- changed-files:
+  - any-glob-to-any-file:
+    - .github/**/*.md
+    - .readthedocs.yaml
+    - doc/**/*
+    - '*.md'
+  - all-globs-to-any-file:
+    - '*.rst'
+    - '!CHANGES.rst'
+
+testing:
+- changed-files:
+  - any-glob-to-any-file:
+    - .bandit.yaml
+    - .github/**/*.yml
+    - conftest.py
+    - acstools/tests/**/*
+
+CALACS:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/acs2d.py
+    - acstools/acsccd.py
+    - acstools/acscteforwardmodel.py
+    - acstools/acscte.py
+    - acstools/acsrej.py
+    - acstools/acssum.py
+    - acstools/calacs.py
+    - acstools/utils_calib.py
+
+CTE:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/acscteforwardmodel.py
+    - acstools/acscte.py
+    - acstools/acsphotcte.py
+
+Destripe:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/acs_destripe_plus.py
+    - acstools/acs_destripe.py
+
+Polarization:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/polarization_tools.py
+    - acstools/data/polarizer_tables.yaml
+
+Satellite Detection:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/satdet.py
+    - acstools/findsat_mrt.py
+    - acstools/utils_findsat_mrt.py
+
+Zeropoints:
+- changed-files:
+  - any-glob-to-any-file:
+    - acstools/acszpt.py

{acstools-3.7.0 → acstools-3.7.2}/.github/workflows/ci_workflows.yml

@@ -21,7 +21,7 @@ jobs:
     runs-on: ubuntu-latest
     steps:
     - name: Check base branch
-      uses: actions/github-script@v6
+      uses: actions/github-script@v7
       if: github.event_name == 'pull_request'
       with:
         script: |
@@ -59,6 +59,6 @@ jobs:
           linux: py39-test-alldeps
           posargs: --remote-data -v
 
-        - name: Python 3.12 with remote data
+        - name: Python 3.12 with remote data and dev dependencies
           linux: py312-test-devdeps
           posargs: --remote-data -v

{acstools-3.7.0 → acstools-3.7.2}/.github/workflows/codeql-analysis.yml

@@ -41,7 +41,7 @@ jobs:
 
     # Initializes the CodeQL tools for scanning.
     - name: Initialize CodeQL
-      uses: github/codeql-action/init@v2
+      uses: github/codeql-action/init@v3
       with:
         languages: ${{ matrix.language }}
         # If you wish to specify custom queries, you can do so here or in a config file.
@@ -52,7 +52,7 @@ jobs:
     # Autobuild attempts to build any compiled languages  (C/C++, C#, or Java).
     # If this step fails, then you should remove it and run the build manually (see below)
     - name: Autobuild
-      uses: github/codeql-action/autobuild@v2
+      uses: github/codeql-action/autobuild@v3
 
     # ℹ️ Command-line programs to run using the OS shell.
     # 📚 https://git.io/JvXDl
@@ -66,4 +66,4 @@ jobs:
     #   make release
 
     - name: Perform CodeQL Analysis
-      uses: github/codeql-action/analyze@v2
+      uses: github/codeql-action/analyze@v3

{acstools-3.7.0 → acstools-3.7.2}/.github/workflows/open_actions.yml

@@ -17,13 +17,12 @@ jobs:
     runs-on: ubuntu-latest
     steps:
     - name: Label PR
-      uses: actions/labeler@v4
+      uses: actions/labeler@v5
       if: github.event_name == 'pull_request_target' || github.event_name == 'pull_request'
       with:
         repo-token: "${{ secrets.GITHUB_TOKEN }}"
-        sync-labels: ''
     - name: 'Comment Draft PR'
-      uses: actions/github-script@v6
+      uses: actions/github-script@v7
       if: github.event.pull_request.draft == true
       with:
         github-token: ${{ secrets.GITHUB_TOKEN }}

{acstools-3.7.0 → acstools-3.7.2}/.github/workflows/predeps_workflows.yml

@@ -18,8 +18,8 @@ jobs:
         jref: "https://ssb.stsci.edu/trds_open/jref"
       submodules: false
       coverage: ''
-    # TODO: Enable more OSes when possible.
-    envs: |
-      - name: Python 3.11 with pre-release
-        linux: py311-test-predeps
-        posargs: --remote-data -v
+      # TODO: Enable more OSes when possible.
+      envs: |
+        - name: Python 3.11 with pre-release
+          linux: py311-test-predeps
+          posargs: --remote-data -v

{acstools-3.7.0 → acstools-3.7.2}/.github/workflows/publish-to-pypi.yml

@@ -36,7 +36,7 @@ jobs:
       with:
         fetch-depth: 0
 
-    - uses: actions/setup-python@v4
+    - uses: actions/setup-python@v5
       with:
         python-version: '3.11'
 

{acstools-3.7.0 → acstools-3.7.2}/CITATION.md

@@ -1,7 +1,9 @@
 If you use acstools, please cite as:
 
-[![ascl:2011.024](https://img.shields.io/badge/ascl-2011.024-blue.svg?colorB=262255)](https://ascl.net/2011.024)
+[![NASA/ADS-ACStools](https://img.shields.io/badge/NASA%2FADS-ACStools-blue)](https://ui.adsabs.harvard.edu/abs/2020ascl.soft11024L/abstract)
 
 #### Other forms of citation
 
+[![ascl:2011.024](https://img.shields.io/badge/ascl-2011.024-blue.svg?colorB=262255)](https://ascl.net/2011.024)
+
 [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.7406933.svg)](https://doi.org/10.5281/zenodo.7406933)

{acstools-3.7.0 → acstools-3.7.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: acstools
-Version: 3.7.0
+Version: 3.7.2
 Summary: Python Tools for HST ACS
 Home-page: https://github.com/spacetelescope/acstools
 Author: Matt Davis, Warren Hack, Norman Grogin, Pey Lian Lim, Sara Ogaz, Leonardo Ubeda, Mihai Cara, David Borncamp, Nathan Miles, Tyler Desjardins, Jenna Ryon, David Stark, Gagandeep Anand, Yotam Cohen
@@ -39,7 +39,7 @@ Requires-Dist: pytest-astropy-header; extra == "test"
 Requires-Dist: pytest-remotedata; extra == "test"
 Requires-Dist: ci-watson; extra == "test"
 Provides-Extra: docs
-Requires-Dist: sphinx<7; extra == "docs"
+Requires-Dist: sphinx; extra == "docs"
 Requires-Dist: sphinx-automodapi; extra == "docs"
 Requires-Dist: sphinx-rtd-theme; extra == "docs"
 

{acstools-3.7.0 → acstools-3.7.2}/acstools/findsat_mrt.py

@@ -1290,8 +1290,8 @@ class WfcWrapper(TrailFinder):
 
     @binsize.setter
     def binsize(self, value):
-        if value is not None and not isinstance(value, int):
-            raise ValueError(f"binsize must be None or int but got: {value}")
+        if value is not None and not (isinstance(value, int) and value >= 1):
+            raise ValueError(f"binsize must be None or int >= 1 but got: {value}")  # noqa
         self._binsize = value
 
     @property
@@ -1360,6 +1360,12 @@ class WfcWrapper(TrailFinder):
 
         LOG.info('Rebinning the data by {}'.format(self.binsize))
 
+        # trigger a warning if the image dimensions are not even
+        if (((self.image.shape[0] % self.binsize != 0)) |
+           ((self.image.shape[1] % self.binsize != 0))):
+            LOG.warning('Image dimensions do not evenly divide by bin size.'
+                        '\nExtra rows/columns will be trimmed first.')
+
         # setting a warning filter for the nansum calculations. These
         # warnings are inconsequential and already accounted for in the code
         with warnings.catch_warnings():
@@ -1406,4 +1412,4 @@ class WfcWrapper(TrailFinder):
             raise ValueError(f'DQ array can only be updated for FLC/FLT images, not {self.image_type}')  # noqa
 
         update_dq(self.image_file, self.extension + 2, self.mask, dqval=dqval,
-                  verbose=verbose)
+                  verbose=verbose, expand_mask=True)

{acstools-3.7.0 → acstools-3.7.2}/acstools/satdet.py

@@ -285,7 +285,7 @@ def _detsat_one(filename, ext, sigma=2.0, low_thresh=0.1, h_thresh=0.5,
         trail_angle = trail_angle[mask]
         result = result[mask]
 
-        ang, num = stats.mode(round_angle)
+        ang, num = stats.mode(round_angle, keepdims=True)
 
         # do the filtering
         truth = round_angle == ang[0]

{acstools-3.7.0 → acstools-3.7.2}/acstools/tests/test_wfc_destripe.py

@@ -26,7 +26,9 @@ class TestDestripe(BaseACSTOOLS):
                                mask1=masks[0], mask2=masks[1])
 
         # Compare results
-        self.compare_outputs([(outputfile, truthfile)], rtol=1e-6)
+        # NOTE: atol instead of rtol was set due to numpy 2.0,
+        # see https://github.com/spacetelescope/acstools/issues/184
+        self.compare_outputs([(outputfile, truthfile)], atol=0.01)
 
     # jb5g05ubq = calibrated post-SM4 WFC full-frame exposures without masks
     # ja0x03ojq = calibrated polarizer WFC subarray exposures

{acstools-3.7.0 → acstools-3.7.2}/acstools/tests/test_wfc_findsat_mrt.py

@@ -1,11 +1,14 @@
 '''Test satellite trail detection and masking using findsat_mrt.'''
 import os
 
+import pytest
 from astropy.io.fits import FITSDiff
 from astropy.utils.data import get_pkg_data_filename
 
 from acstools.findsat_mrt import WfcWrapper
 
+pytest.importorskip("skimage")
+
 
 def test_WfcWrapper(tmp_path):
     """Identify and mask trails in WFC extension 4."""

{acstools-3.7.0 → acstools-3.7.2}/acstools/tests/test_wfc_satdet.py

@@ -3,10 +3,14 @@
 .. note:: Cannot test ``detsat()`` because PHT results change from run to run!
 
 """
+import pytest
+
 from acstools import satdet
 from acstools.utils_findsat_mrt import update_dq
 from acstools.tests.helpers import BaseACSTOOLS
 
+pytest.importorskip("skimage")
+
 
 class TestSatDet(BaseACSTOOLS):
     detector = 'wfc'

{acstools-3.7.0 → acstools-3.7.2}/acstools/utils_findsat_mrt.py

@@ -10,7 +10,7 @@ import numpy as np
 from astropy.convolution import convolve, Gaussian2DKernel
 from astropy.io import fits
 from astropy.modeling import models, fitting
-from astropy.nddata import Cutout2D
+from astropy.nddata import Cutout2D, block_replicate
 from astropy.stats import sigma_clip
 from astropy.table import Table, vstack
 from astropy.utils.exceptions import AstropyUserWarning
@@ -20,8 +20,8 @@ try:
     from skimage import transform
     from skimage.transform._warps import warp
     from skimage._shared.utils import convert_to_float
-except ImportError:
-    warnings.warn('skimage not installed. MRT calculation will not work.')
+except ImportError as e:
+    warnings.warn(f'skimage not installed. MRT calculation will not work: {repr(e)}')  # noqa
 
 # check for scipy
 try:
@@ -108,8 +108,8 @@ def merge_tables(tbls, theta_sep=10, rho_sep=10):
             drho = np.abs(t['ycentroid'] - s['ycentroid'])
             sel = (dtheta < theta_sep) & (drho < rho_sep)
             keep[sel] = False
-        # Silence MergeConflictWarning for "date" metadata because
-        # the tables were created with different timestamps. The last entry is kept.
+        # Silence MergeConflictWarning for "date" metadata because the tables
+        # were created with different timestamps. The last entry is kept.
         src = vstack([src, t[keep]], metadata_conflicts="silent")
     return src
 
@@ -1179,7 +1179,8 @@ def radon(image, theta=None, circle=False, *, preserve_range=False,
         return radon_image
 
 
-def update_dq(filename, ext, mask, dqval=16384, verbose=True):
+def update_dq(filename, ext, mask, dqval=16384, verbose=True,
+              expand_mask=False):
     """Update the given image and DQ extension with the given
     satellite trails mask and flag.
 
@@ -1199,10 +1200,85 @@ def update_dq(filename, ext, mask, dqval=16384, verbose=True):
         tailored for ACS/WFC.
     verbose : bool, optional
         Print extra information to the terminal.
+    expand_mask : bool, optional
+        Allows the mask to be expanded to match the size of the original mask.
+        This is relevant for masks that were generated used binned versions of
+        the original image.
 
     """
     with fits.open(filename, mode='update') as pf:
         dqarr = pf[ext].data
+
+        # check that the old mask and new mask are the same size
+        same_size = mask.shape == dqarr.shape
+        mask_larger = np.any([mask.shape[i] > dqarr.shape[i]
+                              for i in range(mask.ndim)])
+
+        # provide informative error messages depending on mask/dqarr size
+        if not same_size:
+            LOG.info('Inconsistent mask sizes: \n'
+                     'Input mask: {} \n'
+                     'DQ array: {}'.format(mask.shape, dqarr.shape))
+
+            if mask_larger:
+                LOG.error('Mask is larger than the DQ array in one or more '
+                          'dimensions.\nThis routine cannot currently rescale '
+                          'masks that are larger than the origin input image')
+                return
+
+            elif not mask_larger and not expand_mask:
+                LOG.error('Cannot proceed due to size mismatch.\n'
+                          'Set expand_mask=True if the mask was generated '
+                          'with binned data and needs to be enlarged to the '
+                          'original size.')
+                return
+            else:
+                LOG.info('Enlarging mask to original size.')
+
+                # Determine factor by which to scale-up image. Scales should
+                # be the same for the x and y dimensions, after we round down
+                # to the nearest integer in case there are extra rows/columns
+                # that do not divide evenly. These extra rows/columns will be
+                # added back at the end. This approach is the reverse of
+                # astropy.nddata.block_reduce, which trims extra rows/cols
+                # first.
+
+                factor_x = np.floor(dqarr.shape[1]/mask.shape[1]).astype(int)
+                factor_y = np.floor(dqarr.shape[0]/mask.shape[0]).astype(int)
+                if factor_x != factor_y:
+
+                    mask_aspect = mask.shape[1]/mask.shape[0]
+                    dq_aspect = dqarr.shape[1]/dqarr.shape[0]
+
+                    raise ValueError('Mask and DQ array axis aspect ratios are'
+                                     ' too different\n'
+                                     'Mask dimensions (x,y) = {}\n'
+                                     'Mask aspect ratio (x/y) = {}\n'
+                                     'DQ array dimensions (y,x) = {}\n'
+                                     'DQ array aspect ratio (x/y) = {}'
+                                     .format(mask.shape, mask_aspect,
+                                             dqarr.shape, dq_aspect))
+
+                else:
+                    factor = factor_x
+
+                # assess whether there will be extra rows/columns after scaling
+                extra_cols = dqarr.shape[1] - mask.shape[1]*factor
+                extra_rows = dqarr.shape[0] - mask.shape[0]*factor
+
+                mask = block_replicate(mask, factor,
+                                       conserve_sum=False)
+
+                if extra_rows > 0:
+                    LOG.info('Duplicating end rows to match orig size')
+                    while mask.shape[0] != dqarr.shape[0]:
+                        mask = np.vstack([mask, mask[-1, :]])
+
+                if extra_cols > 0:
+                    LOG.info('Duplicating end columns to match orig size')
+                    while mask.shape[1] != dqarr.shape[1]:
+                        mask = np.hstack([mask, mask[:, -1].reshape(mask.shape[0], 1)])  # noqa
+
         old_mask = (dqval & dqarr) != 0  # Existing flagged trails
         new_mask = mask & ~old_mask  # Only flag previously unflagged trails
         npix_updated = np.count_nonzero(new_mask)
@@ -1356,16 +1432,24 @@ def create_mrt_line_kernel(width, sigma, outfile=None, shape=(1024, 2048),
         LOG.info('Inteprolating onto new grid to center kernel')
         theta_arr = np.arange(cutout.shape[1])
         rho_arr = np.arange(cutout.shape[0])
-        theta_grid, rho_grid = np.meshgrid(theta_arr, rho_arr)
 
         new_theta_arr = theta_arr + theta_shift
         new_rho_arr = rho_arr + rho_shift
         new_theta_grid, new_rho_grid = np.meshgrid(new_theta_arr, new_rho_arr)
 
-        # inteprolate onto new grid
-        f = interpolate.interp2d(theta_grid, rho_grid, cutout.data,
-                                 kind='cubic')
-        cutout = f(new_theta_arr, new_rho_arr)  # overwrite old cutout
+        f = interpolate.RegularGridInterpolator((rho_arr, theta_arr),
+                                                cutout.data,
+                                                bounds_error=False,
+                                                fill_value=0,
+                                                method='cubic')
+
+        # unable to provide 2D grid of points for interpolation. Switching to
+        # 1D, and then will switch back
+        points = [[r, t] for r, t in zip(np.ravel(new_rho_grid),
+                                         np.ravel(new_theta_grid))]
+
+        # return to original shape and replace the old cutout
+        cutout = np.reshape(f(points), cutout.data.shape)
 
     if ax:
         fig4, ax4 = plt.subplots()

{acstools-3.7.0 → acstools-3.7.2}/acstools/version.py

@@ -12,5 +12,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '3.7.0'
-__version_tuple__ = version_tuple = (3, 7, 0)
+__version__ = version = '3.7.2'
+__version_tuple__ = version_tuple = (3, 7, 2)

{acstools-3.7.0 → acstools-3.7.2}/acstools.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: acstools
-Version: 3.7.0
+Version: 3.7.2
 Summary: Python Tools for HST ACS
 Home-page: https://github.com/spacetelescope/acstools
 Author: Matt Davis, Warren Hack, Norman Grogin, Pey Lian Lim, Sara Ogaz, Leonardo Ubeda, Mihai Cara, David Borncamp, Nathan Miles, Tyler Desjardins, Jenna Ryon, David Stark, Gagandeep Anand, Yotam Cohen
@@ -39,7 +39,7 @@ Requires-Dist: pytest-astropy-header; extra == "test"
 Requires-Dist: pytest-remotedata; extra == "test"
 Requires-Dist: ci-watson; extra == "test"
 Provides-Extra: docs
-Requires-Dist: sphinx<7; extra == "docs"
+Requires-Dist: sphinx; extra == "docs"
 Requires-Dist: sphinx-automodapi; extra == "docs"
 Requires-Dist: sphinx-rtd-theme; extra == "docs"
 

{acstools-3.7.0 → acstools-3.7.2}/acstools.egg-info/requires.txt

@@ -13,7 +13,7 @@ photutils
 dask
 
 [docs]
-sphinx<7
+sphinx
 sphinx-automodapi
 sphinx-rtd-theme
 

{acstools-3.7.0 → acstools-3.7.2}/doc/source/index.rst

@@ -5,7 +5,7 @@ ACSTOOLS
 ########
 
 Software tools for
-`Advanced Camera for Surveys (ACS) <http://www.stsci.edu/hst/acs/>`_.
+`Advanced Camera for Surveys (ACS) <https://www.stsci.edu/hst/acs/>`_.
 
 Different ways to install the latest release of this package::
 
@@ -52,3 +52,10 @@ Indices and tables
 * :ref:`genindex`
 * :ref:`modindex`
 * :ref:`search`
+
+******************
+Citation
+******************
+
+* Please cite `acstools <https://ui.adsabs.harvard.edu/abs/2020ascl.soft11024L/abstract>`_ if you have used it in your work. 
+

{acstools-3.7.0 → acstools-3.7.2}/setup.cfg

@@ -52,7 +52,7 @@ test =
 	pytest-remotedata
 	ci-watson
 docs = 
-	sphinx<7
+	sphinx
 	sphinx-automodapi
 	sphinx-rtd-theme
 
@@ -83,6 +83,7 @@ filterwarnings =
 	ignore:pkg_resources is deprecated as an API:DeprecationWarning
 	ignore:.*utcfromtimestamp:DeprecationWarning
 	ignore:.*numpy\.core.*:DeprecationWarning
+	ignore:skimage not installed
 
 [flake8]
 ignore = E221,E226,E262,E265,E501,E704,W504

{acstools-3.7.0 → acstools-3.7.2}/tox.ini

@@ -11,7 +11,7 @@ envlist =
 passenv = HOME,CI,jref
 
 setenv =
-    devdeps: PIP_EXTRA_INDEX_URL = https://pypi.anaconda.org/astropy/simple https://pypi.anaconda.org/scientific-python-nightly-wheels/simple
+    devdeps: PIP_EXTRA_INDEX_URL = https://pypi.anaconda.org/astropy/simple https://pypi.anaconda.org/liberfa/simple https://pypi.anaconda.org/scientific-python-nightly-wheels/simple
 
 # Run the tests in a temporary directory to make sure that we don't import
 # package from the source tree
@@ -38,8 +38,9 @@ deps =
     devdeps: scipy>=0.0.dev0
     devdeps: matplotlib>=0.0.dev0
     devdeps: scikit-image>=0.0.dev0
+    devdeps: pyerfa>=0.0.dev0
     devdeps: astropy>=0.0.dev0
-    devdeps: git+https://github.com/astropy/photutils.git
+    devdeps: photutils>=0.0.dev0
     devdeps: git+https://github.com/spacetelescope/ci_watson.git
     devdeps: stsci.tools
 
@@ -49,8 +50,6 @@ extras =
     predeps: all
 
 commands =
-    # Force numpy-dev after matplotlib downgrades it (https://github.com/matplotlib/matplotlib/issues/26847)
-    devdeps: python -m pip install --pre --upgrade --extra-index-url https://pypi.anaconda.org/scientific-python-nightly-wheels/simple numpy
     pip freeze
     pytest --pyargs acstools {toxinidir}/doc {posargs}
 

acstools-3.7.0/.github/dependabot.yml

@@ -1,8 +0,0 @@
-version: 2
-updates:
-- package-ecosystem: pip
-  directory: "/"
-  schedule:
-    interval: monthly
-    time: "03:00"
-  open-pull-requests-limit: 10

acstools-3.7.0/.github/labeler.yml

@@ -1,43 +0,0 @@
-docs:
-  - .github/**/*.md
-  - .readthedocs.yaml
-  - doc/**/*
-  - '*.md'
-  - any: ['*.rst', '!CHANGES.rst']
-
-testing:
-  - .bandit.yaml
-  - .github/**/*.yml
-  - conftest.py
-  - acstools/tests/**/*
-  - Jenkinsfile
-  - JenkinsfileRT
-
-CALACS:
-  - acstools/acs2d.py
-  - acstools/acsccd.py
-  - acstools/acscteforwardmodel.py
-  - acstools/acscte.py
-  - acstools/acsrej.py
-  - acstools/acssum.py
-  - acstools/calacs.py
-  - acstools/utils_calib.py
-
-CTE:
-  - acstools/acscteforwardmodel.py
-  - acstools/acscte.py
-  - acstools/acsphotcte.py
-
-Destripe:
-  - acstools/acs_destripe_plus.py
-  - acstools/acs_destripe.py
-
-Polarization:
-  - acstools/polarization_tools.py
-  - acstools/data/polarizer_tables.yaml
-
-Satellite Detection:
-  - acstools/satdet.py
-
-Zeropoints:
-  - acstools/acszpt.py