PyStormTracker 0.2.1__tar.gz

pystormtracker-0.2.1/.gitattributes

@@ -0,0 +1 @@
+data/test/*.nc filter=lfs diff=lfs merge=lfs -text export-ignore

pystormtracker-0.2.1/.github/workflows/python-publish.yml

@@ -0,0 +1,65 @@
+# This workflow will upload a Python Package to PyPI when a release is created
+# For more information see: https://docs.github.com/en/actions/automating-builds-and-tests/building-and-testing-python#publishing-to-package-registries
+
+name: Upload Python Package
+
+on:
+  release:
+    types: [published]
+
+permissions:
+  contents: read
+
+jobs:
+  release-build:
+    runs-on: ubuntu-latest
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.x"
+
+      - name: Build release distributions
+        run: |
+          python -m pip install build
+          python -m build
+
+      - name: Upload distributions
+        uses: actions/upload-artifact@v4
+        with:
+          name: release-dists
+          path: dist/
+
+  pypi-publish:
+    runs-on: ubuntu-latest
+    needs:
+      - release-build
+    permissions:
+      # IMPORTANT: this permission is mandatory for trusted publishing
+      id-token: write
+
+    # Dedicated environments with protections for publishing are strongly recommended.
+    environment:
+      name: pypi
+      url: https://pypi.org/project/PyStormTracker/${{ steps.get_version.outputs.version }}
+
+    steps:
+      - name: Get version from tag
+        id: get_version
+        run: |
+          # Strips 'v' prefix from tag_name (e.g. v0.2.1 -> 0.2.1)
+          VERSION=${{ github.event.release.tag_name }}
+          echo "version=${VERSION#v}" >> $GITHUB_OUTPUT
+
+      - name: Retrieve release distributions
+        uses: actions/download-artifact@v4
+        with:
+          name: release-dists
+          path: dist/
+
+      - name: Publish release distributions to PyPI
+        uses: pypa/gh-action-pypi-publish@release/v1
+        with:
+          packages-dir: dist/

pystormtracker-0.2.1/.gitignore

@@ -0,0 +1,70 @@
+# Byte-compiled / optimized / DLL files
+__pycache__/
+*.py[cod]
+
+# C extensions
+*.so
+
+# Distribution / packaging
+.Python
+env/
+build/
+develop-eggs/
+dist/
+downloads/
+eggs/
+.eggs/
+lib/
+lib64/
+parts/
+sdist/
+var/
+*.egg-info/
+.installed.cfg
+*.egg
+
+# PyInstaller
+#  Usually these files are written by a python script from a template
+#  before PyInstaller builds the exe, so as to inject date/other infos into it.
+*.manifest
+*.spec
+
+# Installer logs
+pip-log.txt
+pip-delete-this-directory.txt
+
+# Unit test / coverage reports
+htmlcov/
+.tox/
+.coverage
+.coverage.*
+.cache
+nosetests.xml
+coverage.xml
+*,cover
+
+# Translations
+*.mo
+*.pot
+
+# Django stuff:
+*.log
+
+# Sphinx documentation
+docs/_build/
+
+# PyBuilder
+target/
+
+# IPython intermediate checkpoints
+.ipynb_checkpoints
+
+# Python Pickle files
+*.pickle
+
+# netCDF files
+*.nc
+!data/test/*.nc
+
+# CSV files
+*.csv

pystormtracker-0.2.1/.pre-commit-config.yaml

@@ -0,0 +1,24 @@
+repos:
+  - repo: https://github.com/pre-commit/pre-commit-hooks
+    rev: v4.5.0
+    hooks:
+      - id: trailing-whitespace
+      - id: end-of-file-fixer
+      - id: check-yaml
+      - id: check-added-large-files
+      - id: debug-statements
+
+  - repo: https://github.com/astral-sh/ruff-pre-commit
+    rev: v0.3.4
+    hooks:
+      # Run the linter.
+      - id: ruff
+        args: [ --fix ]
+      # Run the formatter.
+      - id: ruff-format
+
+  - repo: https://github.com/pre-commit/mirrors-mypy
+    rev: v1.9.0
+    hooks:
+      - id: mypy
+        additional_dependencies: [ types-all ]

pystormtracker-0.2.1/.zenodo.json

@@ -0,0 +1,34 @@
+{
+    "title": "PyStormTracker",
+    "description": "A Parallel Object-Oriented Cyclone Tracker in Python",
+    "creators": [
+        {
+            "name": "Yau, Albert M. W.",
+            "orcid": "00009-0000-7825-9627"
+        }
+    ],
+    "keywords": [
+        "cyclone tracking",
+        "climate variability",
+        "dask",
+        "mpi"
+    ],
+    "license": "BSD-3-Clause",
+    "upload_type": "software",
+    "version": "0.2.1",
+    "language": "python",
+    "code_repository": "https://github.com/mwyau/PyStormTracker",
+    "development_status": "active",
+    "related_identifiers": [
+        {
+            "identifier": "10.1175/JCLI-D-20-0393.1",
+            "relation": "isSupplementTo",
+            "scheme": "doi"
+        },
+        {
+            "identifier": "10.5281/zenodo.18868625",
+            "relation": "isDerivedFrom",
+            "scheme": "doi"
+        }
+    ]
+}

pystormtracker-0.2.1/CITATION.cff

@@ -0,0 +1,46 @@
+cff-version: 1.2.0
+message: "If you use this software, please cite it as below."
+authors:
+  - family-names: Yau
+    given-names: Albert M. W.
+    email: albert@mwyau.com
+title: "PyStormTracker"
+version: 0.2.1
+doi: 10.5281/zenodo.18764813
+date-released: 2026-03-04
+url: "https://github.com/mwyau/PyStormTracker"
+repository-code: "https://github.com/mwyau/PyStormTracker"
+license: BSD-3-Clause
+keywords:
+  - cyclone tracking
+  - climate variability
+  - dask
+  - mpi
+preferred-citation:
+  type: article
+  authors:
+    - family-names: Yau
+      given-names: Albert M. W.
+    - family-names: Chang
+      given-names: Edmund K. M.
+  title: "Finding Storm Track Activity Metrics That Are Highly Correlated with Weather Impacts. Part I: Frameworks for Evaluation and Accumulated Track Activity"
+  journal: "Journal of Climate"
+  volume: 33
+  issue: 23
+  start: 10169
+  end: 10186
+  year: 2020
+  doi: 10.1175/JCLI-D-20-0393.1
+references:
+  - type: poster
+    authors:
+      - family-names: Yau
+        given-names: Albert M. W.
+      - family-names: Paul
+        given-names: Kevin
+      - family-names: Dennis
+        given-names: John
+    title: "PyStormTracker: A Parallel Object-Oriented Cyclone Tracker in Python"
+    collection-title: "96th American Meteorological Society Annual Meeting"
+    year: 2016
+    doi: 10.5281/zenodo.18868625

pystormtracker-0.2.1/LICENSE

@@ -0,0 +1,40 @@
+Copyright (c) 2015, University Corporation for Atmospheric Research (UCAR)
+Copyright (c) 2026, Albert M. W. Yau
+All rights reserved.
+
+Developed by: Application Scalability and Performance Group
+              University Corporation for Atmospheric Research -
+              National Center for Atmospheric Research
+              https://www.cisl.ucar.edu/
+
+Developed by: Albert M. W. Yau
+              https://github.com/mwyau/PyStormTracker
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the "Software"),
+to deal with the Software without restriction, including without limitation
+the rights to use, copy, modify, merge, publish, distribute, sublicense,
+and/or sell copies of the Software, and to permit persons to whom
+the Software is furnished to do so, subject to the following conditions:
+
+- Redistributions of source code must retain the above copyright notice,
+  this list of conditions and the following disclaimers.
+- Redistributions in binary form must reproduce the above copyright notice,
+  this list of conditions and the following disclaimers in the documentation
+  and/or other materials provided with the distribution.
+- Neither the names of the Application Scalability and Performance (ASAP)
+  Group, UCAR, nor the names of its contributors may be used to endorse or
+  promote products derived from this Software without specific prior written
+  permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
+LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+POSSIBILITY OF SUCH DAMAGE.

pystormtracker-0.2.1/PKG-INFO

@@ -0,0 +1,146 @@
+Metadata-Version: 2.4
+Name: PyStormTracker
+Version: 0.2.1
+Summary: A Parallel Object-Oriented Cyclone Tracker in Python
+Project-URL: Homepage, https://github.com/mwyau/PyStormTracker
+Project-URL: Repository, https://github.com/mwyau/PyStormTracker.git
+Project-URL: Issues, https://github.com/mwyau/PyStormTracker/issues
+Author-email: "Albert M. W. Yau" <albert@mwyau.com>
+Maintainer-email: "Albert M. W. Yau" <albert@mwyau.com>
+License: BSD-3-Clause
+License-File: LICENSE
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: BSD License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: Atmospheric Science
+Requires-Python: >=3.10
+Requires-Dist: dask>=2024.2.0
+Requires-Dist: distributed>=2024.2.0
+Requires-Dist: netcdf4>=1.6.5
+Requires-Dist: numpy>=1.24.0
+Requires-Dist: scipy>=1.10.0
+Provides-Extra: dev
+Requires-Dist: mpi4py>=3.1.5; extra == 'dev'
+Requires-Dist: mypy>=1.9.0; extra == 'dev'
+Requires-Dist: netcdf4>=1.6.5; extra == 'dev'
+Requires-Dist: pandas>=2.2.0; extra == 'dev'
+Requires-Dist: pre-commit>=3.7.0; extra == 'dev'
+Requires-Dist: pytest>=8.1.1; extra == 'dev'
+Requires-Dist: ruff>=0.3.4; extra == 'dev'
+Description-Content-Type: text/markdown
+
+# PyStormTracker
+
+[![DOI](https://zenodo.org/badge/36328800.svg)](https://doi.org/10.5281/zenodo.18764813)
+
+PyStormTracker provides the implementation of the "Simple Tracker" algorithm used for cyclone trajectory analysis in **Yau and Chang (2020)**. It was originally developed at the **National Center for Atmospheric Research (NCAR)** as part of the **2015 Summer Internships in Parallel Computational Science (SIParCS)** program, utilizing a task-parallel strategy with temporal decomposition and a tree reduction algorithm to process large climate datasets.
+
+## Features
+
+- **Modern Python 3 Support**: Fully migrated from Python 2 with comprehensive type hints.
+- **Flexible Data Support**: Works with `netCDF4` and handles various coordinate naming conventions (`lat`/`lon` vs `latitude`/`longitude`).
+- **Parallel Backends**:
+  - **Dask (Default)**: Automatically scales to all available CPU cores on local machines.
+  - **MPI**: Supports distributed execution via `mpi4py`.
+  - **Serial**: Standard sequential execution for smaller datasets or debugging.
+- **Robust Detection**: Handles masked/missing data correctly and includes automated unit/integration tests.
+- **User-Friendly Output**: Results are exported directly to CSV with readable datetime strings and formatted numeric values.
+
+## Technical Methodology
+
+PyStormTracker treats meteorological fields as 2D images and leverages `scipy.ndimage` for robust feature detection:
+
+- **Local Extrema Detection**: Uses `generic_filter` with a sliding window (default 5x5) to identify local minima (cyclones) or maxima (anticyclones/vorticity).
+- **Intensity & Refinement**: Applies the **Laplacian operator** (`laplace`) to measure the "sharpness" of the field at each detected center point. This is used to resolve duplicates and ensure only the most physically intense point is kept when multiple adjacent pixels are flagged.
+- **Spherical Continuity**: Utilizes `mode='wrap'` for all filters to correctly handle periodic boundaries across the Prime Meridian, enabling seamless tracking across the entire globe.
+- **Heuristic Linking**: Implements a nearest-neighbor linking strategy to connect detected centers into trajectories across successive time steps.
+
+## Installation
+
+### Prerequisites
+- Python 3.10+
+- (Optional) MS-MPI or OpenMPI for MPI support.
+
+### From PyPI (Recommended)
+You can install the latest stable version of PyStormTracker directly from PyPI:
+```bash
+pip install PyStormTracker
+```
+
+### From Source
+1. Clone the repository:
+   ```bash
+   git clone https://github.com/mwyau/PyStormTracker.git
+   cd PyStormTracker
+   ```
+
+2. Install the package in editable mode:
+   ```bash
+   pip install -e .
+   ```
+
+## Usage
+
+Once installed, you can use the `stormtracker` command directly:
+
+```bash
+stormtracker -i data/test/slp.2012.nc -v slp -o my_tracks
+```
+
+### Command Line Arguments
+
+| Argument | Short | Description |
+| :--- | :--- | :--- |
+| `--input` | `-i` | **Required.** Path to the input NetCDF file. |
+| `--var` | `-v` | **Required.** Variable name to track (e.g., `slp`, `vo`). |
+| `--output` | `-o` | **Required.** Path to the output CSV file (appends `.csv` if missing). |
+| `--num` | `-n` | Number of time steps to process. |
+| `--mode` | `-m` | `min` (default) for low pressure, `max` for vorticity/high pressure. |
+| `--backend` | `-b` | `dask` (default), `serial`, or `mpi`. |
+| `--workers` | `-w` | Number of Dask workers (defaults to CPU core count). |
+
+### Examples
+
+**Run with Dask (Auto-detected cores):**
+```bash
+stormtracker -i input.nc -v slp -o tracks
+```
+
+**Run with MPI (Distributed):**
+```bash
+mpiexec -n 4 stormtracker -i input.nc -v slp -o tracks -b mpi
+```
+
+## Project Structure
+
+- `src/pystormtracker/models/`: Data structures (`Center`, `Grid`, `Tracks`).
+- `src/pystormtracker/simple/`: Implementation of the Simple Tracker logic (`SimpleDetector`, `SimpleLinker`).
+- `src/pystormtracker/stormtracker.py`: CLI orchestration and parallel backends.
+
+## Testing
+
+Run the full test suite (unit and integration tests) using `pytest`:
+
+```bash
+pytest
+```
+
+## Citations
+
+If you use this software in your research, please cite the following:
+
+- **Yau, A. M. W., and E. K. M. Chang**, 2020: Finding Storm Track Activity Metrics That Are Highly Correlated with Weather Impacts. Part I: Frameworks for Evaluation and Accumulated Track Activity. *J. Climate*, **33**, 10169–10186, https://doi.org/10.1175/JCLI-D-20-0393.1.   
+
+- **Yau, A. M. W.**, 2026: mwyau/PyStormTracker. *Zenodo*, https://doi.org/10.5281/zenodo.18764813.
+
+## References
+
+ - **Yau, A. M. W., K. Paul and J. Dennis**, 2016: PyStormTracker: A Parallel Object-Oriented Cyclone Tracker in Python. *96th American Meteorological Society Annual Meeting*, New Orleans, LA. *Zenodo*, https://doi.org/10.5281/zenodo.18868625.
+
+## License
+
+This project is licensed under the terms found in the `LICENSE` file.

pystormtracker-0.2.1/README.md

@@ -0,0 +1,111 @@
+# PyStormTracker
+
+[![DOI](https://zenodo.org/badge/36328800.svg)](https://doi.org/10.5281/zenodo.18764813)
+
+PyStormTracker provides the implementation of the "Simple Tracker" algorithm used for cyclone trajectory analysis in **Yau and Chang (2020)**. It was originally developed at the **National Center for Atmospheric Research (NCAR)** as part of the **2015 Summer Internships in Parallel Computational Science (SIParCS)** program, utilizing a task-parallel strategy with temporal decomposition and a tree reduction algorithm to process large climate datasets.
+
+## Features
+
+- **Modern Python 3 Support**: Fully migrated from Python 2 with comprehensive type hints.
+- **Flexible Data Support**: Works with `netCDF4` and handles various coordinate naming conventions (`lat`/`lon` vs `latitude`/`longitude`).
+- **Parallel Backends**:
+  - **Dask (Default)**: Automatically scales to all available CPU cores on local machines.
+  - **MPI**: Supports distributed execution via `mpi4py`.
+  - **Serial**: Standard sequential execution for smaller datasets or debugging.
+- **Robust Detection**: Handles masked/missing data correctly and includes automated unit/integration tests.
+- **User-Friendly Output**: Results are exported directly to CSV with readable datetime strings and formatted numeric values.
+
+## Technical Methodology
+
+PyStormTracker treats meteorological fields as 2D images and leverages `scipy.ndimage` for robust feature detection:
+
+- **Local Extrema Detection**: Uses `generic_filter` with a sliding window (default 5x5) to identify local minima (cyclones) or maxima (anticyclones/vorticity).
+- **Intensity & Refinement**: Applies the **Laplacian operator** (`laplace`) to measure the "sharpness" of the field at each detected center point. This is used to resolve duplicates and ensure only the most physically intense point is kept when multiple adjacent pixels are flagged.
+- **Spherical Continuity**: Utilizes `mode='wrap'` for all filters to correctly handle periodic boundaries across the Prime Meridian, enabling seamless tracking across the entire globe.
+- **Heuristic Linking**: Implements a nearest-neighbor linking strategy to connect detected centers into trajectories across successive time steps.
+
+## Installation
+
+### Prerequisites
+- Python 3.10+
+- (Optional) MS-MPI or OpenMPI for MPI support.
+
+### From PyPI (Recommended)
+You can install the latest stable version of PyStormTracker directly from PyPI:
+```bash
+pip install PyStormTracker
+```
+
+### From Source
+1. Clone the repository:
+   ```bash
+   git clone https://github.com/mwyau/PyStormTracker.git
+   cd PyStormTracker
+   ```
+
+2. Install the package in editable mode:
+   ```bash
+   pip install -e .
+   ```
+
+## Usage
+
+Once installed, you can use the `stormtracker` command directly:
+
+```bash
+stormtracker -i data/test/slp.2012.nc -v slp -o my_tracks
+```
+
+### Command Line Arguments
+
+| Argument | Short | Description |
+| :--- | :--- | :--- |
+| `--input` | `-i` | **Required.** Path to the input NetCDF file. |
+| `--var` | `-v` | **Required.** Variable name to track (e.g., `slp`, `vo`). |
+| `--output` | `-o` | **Required.** Path to the output CSV file (appends `.csv` if missing). |
+| `--num` | `-n` | Number of time steps to process. |
+| `--mode` | `-m` | `min` (default) for low pressure, `max` for vorticity/high pressure. |
+| `--backend` | `-b` | `dask` (default), `serial`, or `mpi`. |
+| `--workers` | `-w` | Number of Dask workers (defaults to CPU core count). |
+
+### Examples
+
+**Run with Dask (Auto-detected cores):**
+```bash
+stormtracker -i input.nc -v slp -o tracks
+```
+
+**Run with MPI (Distributed):**
+```bash
+mpiexec -n 4 stormtracker -i input.nc -v slp -o tracks -b mpi
+```
+
+## Project Structure
+
+- `src/pystormtracker/models/`: Data structures (`Center`, `Grid`, `Tracks`).
+- `src/pystormtracker/simple/`: Implementation of the Simple Tracker logic (`SimpleDetector`, `SimpleLinker`).
+- `src/pystormtracker/stormtracker.py`: CLI orchestration and parallel backends.
+
+## Testing
+
+Run the full test suite (unit and integration tests) using `pytest`:
+
+```bash
+pytest
+```
+
+## Citations
+
+If you use this software in your research, please cite the following:
+
+- **Yau, A. M. W., and E. K. M. Chang**, 2020: Finding Storm Track Activity Metrics That Are Highly Correlated with Weather Impacts. Part I: Frameworks for Evaluation and Accumulated Track Activity. *J. Climate*, **33**, 10169–10186, https://doi.org/10.1175/JCLI-D-20-0393.1.   
+
+- **Yau, A. M. W.**, 2026: mwyau/PyStormTracker. *Zenodo*, https://doi.org/10.5281/zenodo.18764813.
+
+## References
+
+ - **Yau, A. M. W., K. Paul and J. Dennis**, 2016: PyStormTracker: A Parallel Object-Oriented Cyclone Tracker in Python. *96th American Meteorological Society Annual Meeting*, New Orleans, LA. *Zenodo*, https://doi.org/10.5281/zenodo.18868625.
+
+## License
+
+This project is licensed under the terms found in the `LICENSE` file.

pystormtracker-0.2.1/notebooks/detect_center_latlon.ipynb

@@ -0,0 +1,149 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "%pylab inline\n",
+    "\n",
+    "import sys\n",
+    "import numpy as np\n",
+    "import netCDF4\n",
+    "from mpl_toolkits.basemap import Basemap\n",
+    "import matplotlib.pyplot as plt\n",
+    "from ipywidgets import interact\n",
+    "\n",
+    "sys.path.append('../src/')\n",
+    "from pystormtracker import *"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "# Plot chart 0 only\n",
+    "chart = 0\n",
+    "size = 11\n",
+    "\n",
+    "var, time, lat, lon = get_var_handle(filename=\"../data/test/slp.2012.nc\")\n",
+    "lons, lats = np.meshgrid(lon[:], lat[:])\n",
+    "\n",
+    "m = Basemap(projection='cyl',llcrnrlat=-90,urcrnrlat=90,\\\n",
+    "            llcrnrlon=0,urcrnrlon=360,resolution='c')\n",
+    "x, y = m(lons, lats)\n",
+    "\n",
+    "data = var[chart,:,:]\n",
+    "centers = detect_center_latlon(data, size=size)\n",
+    "center_indices = centers.nonzero()\n",
+    "# print(np.transpose(center_indices))\n",
+    "\n",
+    "xcenter = x[center_indices]\n",
+    "ycenter = y[center_indices]\n",
+    "\n",
+    "fig = plt.figure(figsize=(16,10))\n",
+    "ax = fig.add_axes([0.1,0.1,0.9,0.9])\n",
+    "m.drawcoastlines()\n",
+    "\n",
+    "clevs = np.arange(950,1080.,5.)\n",
+    "cs = m.contourf(x, y, data/100, clevs)\n",
+    "cb = m.colorbar(cs)\n",
+    "for xi,yi in zip(xcenter, ycenter):\n",
+    "    plt.text(xi,yi,'L',fontsize=18,fontweight='bold',\n",
+    "                    ha='center',va='center',color='k')\n",
+    "plt.title('All Detected Centers, size=' + str(size) + ', chart=' + str(chart))\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "size = 5\n",
+    "\n",
+    "var, time, lat, lon = get_var_handle(filename=\"../data/test/slp.2012.nc\")\n",
+    "lons, lats = np.meshgrid(lon[:], lat[:])\n",
+    "\n",
+    "m = Basemap(projection='cyl',llcrnrlat=-90,urcrnrlat=90,\\\n",
+    "            llcrnrlon=0,urcrnrlon=360,resolution='c')\n",
+    "x, y = m(lons, lats)\n",
+    "\n",
+    "def pltchart(chart):\n",
+    "    data = var[chart,:,:]\n",
+    "    centers = detect_center_latlon(data, size=size, threshold=10.)\n",
+    "    center_indices = centers.nonzero()\n",
+    "    # print(np.transpose(center_indices))\n",
+    "    \n",
+    "    xcenter = x[center_indices]\n",
+    "    ycenter = y[center_indices]\n",
+    "    \n",
+    "    fig = plt.figure(figsize=(16,10))\n",
+    "    ax = fig.add_axes([0.1,0.1,0.9,0.9])\n",
+    "    m.drawcoastlines()\n",
+    "    \n",
+    "    clevs = np.arange(975,1100.,1.)\n",
+    "    cs = m.contourf(x, y, data/100, clevs)\n",
+    "    cb = m.colorbar(cs)\n",
+    "    for xi,yi in zip(xcenter, ycenter):\n",
+    "        plt.text(xi,yi,'L',fontsize=18,fontweight='bold',\n",
+    "                        ha='center',va='center',color='k')\n",
+    "    plt.title('All Detected Centers, size=' + str(size) + ', chart=' + str(chart))\n",
+    "    plt.show()\n",
+    "\n",
+    "chart = (0, len(time))\n",
+    "    \n",
+    "interact(pltchart, chart=chart)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": []
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.0"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}