pyhardisp 0.2.0__tar.gz → 0.2.2__tar.gz

{pyhardisp-0.2.0/pyhardisp.egg-info → pyhardisp-0.2.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyhardisp
-Version: 0.2.0
+Version: 0.2.2
 Summary: Ocean loading tidal displacement calculator - Python conversion of IERS HARDISP Fortran program
 Author: IERS Conventions, Python Conversion - Craig Miller
 License-Expression: LicenseRef-IERS
@@ -37,6 +37,7 @@ This document describes the Python conversion of the **HARDISP** (Harmonic Displ
 ## What is HARDISP?
 
 HARDISP is a scientific program developed by the International Earth Rotation and Reference Systems Service (IERS) for calculating tidal ocean loading effects at geodetic stations. Depending on the coefficients provided by the ocean loading provider, it computes either:
+
 - **Tidal displacements** (vertical, horizontal components) - when displacement coefficients are provided
 - **Tidal gravity effects** (gravity disturbances) - when gravity-type coefficients are provided
 
@@ -55,6 +56,7 @@ The program is part of the **IERS Conventions 2010** recommendations (Class 1 mo
 Given ocean loading coefficients in BLQ format (from the Bos-Scherneck loading service or equivalent providers), HARDISP computes a time series of tidal ocean loading effects in three components:
 
 **Output type depends on the ocean loading provider's coefficients:**
+
 - **Displacement coefficients**: Returns displacements
   - **dU**: Vertical (radial) displacement (meters)
   - **dS**: South (North-South) displacement (meters)
@@ -74,8 +76,8 @@ HARDISP processes the coefficients as-is, with no automatic conversions between
 The Python conversion consists of:
 **core.py** - Main Python module with optimized implementation
 
-
 Original Fortran files (for reference):
+
 - `HARDISP.F` - Main program
 - `ADMINT.F` - Admittance interpolation
 - `RECURS.F` - Recursive harmonic evaluation
@@ -221,6 +223,7 @@ for day_val in range(1, 8):  # Process days 1-7
 ### Coordinate System
 
 The output displacements are in a **local geodetic frame** at the station:
+
 - **dU**: Radial (positive upward)
 - **dS**: South component (positive southward, i.e., negative latitude direction)
 - **dW**: West component (positive westward, i.e., negative longitude direction)
@@ -228,6 +231,7 @@ The output displacements are in a **local geodetic frame** at the station:
 ### Tidal Constituents Used
 
 The 11 input harmonics (from BLQ format) represent:
+
 1. M₂ - Principal lunar semi-diurnal (12.42 hours)
 2. S₂ - Principal solar semi-diurnal (12.00 hours)
 3. N₂ - Lunar elliptic semi-diurnal (12.66 hours)
@@ -282,7 +286,6 @@ The Python code was converted from these 13 Fortran files:
 12. `ETUTC.F` (284 lines) - ET-UTC calculation
 13. `TOYS.F` - (not in workspace but referenced)
 
-
 ## Function Name Mapping (Original Fortran → New Python)
 
 | Original Fortran | New Python Name | Purpose |
@@ -306,6 +309,7 @@ The Python code was converted from these 13 Fortran files:
 ## Features Implemented
 
 ### Complete Conversion
+
 - [x] All 342 tidal constituents support
 - [x] Recursive harmonic evaluation (efficient)
 - [x] Cubic spline interpolation
@@ -317,6 +321,7 @@ The Python code was converted from these 13 Fortran files:
 - [x] Delaunay argument calculations
 
 ### Additional Features
+
 - [x] Object-oriented design with HardispComputer class
 - [x] NumPy integration for high performance
 - [x] Comprehensive error handling
@@ -327,6 +332,7 @@ The Python code was converted from these 13 Fortran files:
 ## Core Components Converted
 
 ### 1. Date/Time Functions (7 functions)
+
 ```
 ✓ is_leap_year(year) - Check leap year
 ✓ days_before_month(year, month) - Days before month start
@@ -336,6 +342,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 2. Spline Interpolation (3 functions)
+
 ```
 ✓ cublic_spline(x, u) - Compute cubic spline coefficients
 ✓ spline_eval(y, x, u, s) - Evaluate spline at point
@@ -343,6 +350,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 3. Tidal Frequency Calculations (3 functions)
+
 ```
 ✓ calculate_tidal_arguments(year, day, h, m, s) - Initialize tidal calculations
 ✓ tidal_frequency_and_phase(doodson_number) - Get frequency and phase from Doodson number
@@ -350,33 +358,34 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 4. Harmonic Recursion (1 function)
+
 ```
 ✓ recursion(n, hc, nf, om) - Efficient recursive harmonic evaluation
 ```
 
 ### 5. Utility Functions (2 functions)
+
 ```
 ✓ pyshells(x) - Sort array with indices (Shell sort)
 ✓ fortran_int_divide(a, b) - Fortran-style integer division
 ```
 
 ### 6. Main Class: HardispComputer
+
 ```
 ✓ read_blq_format(amp, phase) - Load ocean loading coefficients
 ✓ compute_ocean_loading(...) - Main computation engine
 ```
 
-
-
 ## References
 
 1. Petit, G. and Luzum, B. (eds.), **IERS Conventions (2010)**, IERS Technical Note No. 36, BKG (2010)
-   - Available at: https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php
+   - Available at: <https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php>
 
 2. Agnew, D. C., et al., **HARDISP**: Original algorithm and implementation
 
 3. Scherneck, H.-G., and M. S. Bos, **Ocean Loading Service**: BLQ format specification
-   - Available at: http://www.oso.chalmers.se/~loading/
+   - Available at: <http://www.oso.chalmers.se/~loading/>
 
 4. Cartwright, D. E., and A. C. Edden, **Tides of the Planet Earth**, Geophys. J. R. Astron. Soc. 65, 615-630, 1981
 
@@ -387,10 +396,11 @@ This Python conversion maintains the same IERS Conventions Software License as t
 ## Contact
 
 For questions or issues with the Fortran original:
-- IERS Conventions Center: https://www.iers.org/
-- Email: gpetit@bipm.org or brian.luzum@usno.navy.mil
+
+- IERS Conventions Center: <https://www.iers.org/>
+- Email: <gpetit@bipm.org> or <brian.luzum@usno.navy.mil>
 
 ---
 
-*Python conversion completed: 2024*  
+*Python conversion completed: 2026*  
 *Original Fortran: IERS Conventions 2010*

{pyhardisp-0.2.0 → pyhardisp-0.2.2}/README.md

@@ -1,10 +1,13 @@
 # pyhardisp Python Module - Quick Start Guide
+
 [![codecov](https://codecov.io/gh/craigmillernz/pyhardisp/branch/main/graph/badge.svg)](https://codecov.io/gh/craigmillernz/pyhardisp)
-[![Build wheels](https://github.com/craigmillernz/pyhardisp/actions/workflows/build.yml/badge.svg)](https://github.com/craigmillernz/pyhardisp/actions/workflows/build.yml)
+[![Build and Publish](https://github.com/craigmillernz/pyhardisp/actions/workflows/build.yml/badge.svg)](https://github.com/craigmillernz/pyhardisp/actions/workflows/build.yml)
+[![DOI](https://zenodo.org/badge/1170161431.svg)](https://doi.org/10.5281/zenodo.20189637)
 
 ## 📋 What is HARDISP vs pyhardisp?
 
 **HARDISP** computes ocean loading tidal effects for geodetic stations. It processes coefficients provided by ocean loading providers ([such as the Bos-Scherneck service](https://barre.oso.chalmers.se/loading/l.php)) to compute either:
+
 - **Displacements** (vertical and horizontal) when displacement-type coefficients are provided
 - **Gravity effects** (gravity (nm/s^2) and tilt(nrads)) when gravity-type coefficients are provided
 
@@ -12,33 +15,34 @@ It's part of the IERS Conventions 2010 recommended models for correcting space g
 
 **pyhardisp** is a python conversion of the original HARDISP fortran code from IERS (available in the src dir).
 
-
 ## Installation & Requirements
 
 ```bash
 # Requires:
-- Python 3.6+
-- NumPy (any recent version)
+- Python 3.8+
+- NumPy (> 1.19)
 
 # Installation:
-1. clone repo to your local machine
-2. pip install .
+1. pip install pyhardisp
+2. conda install -c conda-forge pyhardisp
 
 ```
 
 ## 📂 Files in This Directory
 
 ### Python Module
-- **`core.py`** 
+
+- **`core.py`**
   - Complete Python implementation of HARDISP
   - Vectorised routines for speed enhancement
 
 ### Documentation
+
 - **`README_PYTHON_CONVERSION.md`** - Complete user guide with examples
-- **`CONVERSION_SUMMARY.md`** - Technical overview and specifications
 - **`README.md`** - This file
 
 ### Original Fortran Code from IERS (Reference Only)
+
 - `HARDISP.F` - Main program
 - `ADMINT.F` - Admittance interpolation
 - `RECURS.F` - Recursive evaluation
@@ -52,16 +56,19 @@ It's part of the IERS Conventions 2010 recommended models for correcting space g
 ## 🚀 Quick Start (5 minutes)
 
 ### 1. Import the Module
+
 ```python
 import pyhardisp
 ```
 
 ### 2. Create a Computer
+
 ```python
 computer = pyhardisp.HardispComputer()
 ```
 
 ### 3. Load Your Data
+
 ```python
 # Ocean loading coefficients (from BLQ format)
 # 3 rows (vertical, east, north) × 11 columns (tidal constituents)
@@ -80,6 +87,7 @@ computer.read_blq_format(amplitudes, phases, units="nm/s2")
 ```
 
 ### 4. Compute Displacements
+
 ```python
 dz, ds, dw = computer.compute_ocean_loading(
     year=2018, month=4, day=5,
@@ -90,13 +98,15 @@ dz, ds, dw = computer.compute_ocean_loading(
 ```
 
 ### 5. Use the Results
+
 ```python
 # Results are in nm/s^2 for vertical or nrad for horizontal components
 for i in range(24):
     print(f"Hour {i}: Gravity={dz[i]:.6f} nm/s^2, South_tilt={ds[i]:.6f} nrad, West_tilt={dw[i]:.6f} nrad")
 ```
 
-## 📊 Example Output
+## Example Output
+
 ```
 Hour 0: Gravity=-46.218473 nm/s^2, South_tilt=31.320275 nrad, West_tilt=-62.009276 nrad
 Hour 1: Gravity=-46.166374 nm/s^2, South_tilt=31.286324 nrad, West_tilt=-62.111623 nrad
@@ -110,7 +120,7 @@ Hour 8: Gravity=-45.799106 nm/s^2, South_tilt=31.046809 nrad, West_tilt=-62.8244
 Hour 9: Gravity=-45.746273 nm/s^2, South_tilt=31.012328 nrad, West_tilt=-62.925713 nrad
 ```
 
-## 🔑 Key Functions Reference
+## Key Functions Reference
 
 | Function | Purpose |
 |----------|---------|
@@ -135,13 +145,7 @@ For detailed information, see:
    - Algorithm explanations
    - Validation results
 
-2. **`CONVERSION_SUMMARY.md`**
-   - Overview of conversion
-   - List of implemented features
-   - Quality metrics
-   - Performance information
-
-## 🔬 Technical Highlights
+## Technical Highlights
 
 - **342 tidal constituents** for precision of ~0.1%
 - **Recursive harmonic evaluation** for fast computation
@@ -158,7 +162,7 @@ Comparion with Quicktide Pro (QTP) and hardisp_x64.exe
 
 ![Conparison of ocean loading from Quicktide Pro, gtools (using hardisp_x64.exe) and pyhardisp](./ocean_loading_comparison.png)
 
-## 🎓 Understanding the Input Data
+## Understanding the Input Data
 
 The BLQ format contains 11 tidal constituents:
 
@@ -176,17 +180,18 @@ The BLQ format contains 11 tidal constituents:
 
 These are expanded to 342 constituents by the program for higher precision.
 
-## 🌍 Data Sources
+## Data Sources
 
 Ocean loading coefficients can be obtained from:
 
 **Bos & Scherneck Loading Service**
-- URL: http://https://barre.oso.chalmers.se/loading/l.php/
+
+- URL: <http://https://barre.oso.chalmers.se/loading/l.php/>
 - Format: BLQ (Binary Loading Queue)
 - Coverage: Global
 - Models: Various (CSR, FES, etc.)
 
-## 📚 References
+## References
 
 1. **IERS Conventions (2010)**
    - Petit, G. and Luzum, B. (eds.)
@@ -200,12 +205,12 @@ Ocean loading coefficients can be obtained from:
    - Cartwright & Edden (1981)
    - Doodson numbers and Darwin notation
 
-## 🔗 Related Tools
+## Related Tools
 
 - **SPOTL** - SPOTL ocean loading (original basis for HARDISP)
 - **IERS Conventions Software** - Complete IERS implementations
 
-## ❓ Common Questions
+## Common Questions
 
 **Q: What are BLQ files?**  
 A: Ocean loading coefficients in binary format from the Chalmers university loading service.
@@ -222,20 +227,12 @@ A: Yes, modify the IDT array in the HardispComputer class
 **Q: What about non-Gregorian calendars?**  
 A: The code assumes Gregorian calendar (since ~1582)
 
-## 📞 Support
+## Support
 
-- **For the original Fortran**: Contact IERS (gpetit@bipm.org)
-- **For this Python version**: Check the source code documentation
-- **For BLQ data**: Contact Hans-Georg Scherneck (scherneck@oso.chalmers.se)
+- **For the original Fortran**: Contact IERS (<gpetit@bipm.org>)
+- **For this Python version**: Make an [issue](https://github.com/craigmillernz/pyhardisp)
+- **For BLQ data**: Contact Hans-Georg Scherneck (<scherneck@oso.chalmers.se>)
 
-## 📜 License
+## License
 
 This Python version maintains the same IERS Conventions Software License as the original Fortran code. See copyright notices in the source files.
-
----
-
-**Happy computing! 🌐**
-
----
-
-*Based on: IERS Conventions 2010, Fortran HARDISP*  

{pyhardisp-0.2.0 → pyhardisp-0.2.2}/README_PYTHON_CONVERSION.md

@@ -7,6 +7,7 @@ This document describes the Python conversion of the **HARDISP** (Harmonic Displ
 ## What is HARDISP?
 
 HARDISP is a scientific program developed by the International Earth Rotation and Reference Systems Service (IERS) for calculating tidal ocean loading effects at geodetic stations. Depending on the coefficients provided by the ocean loading provider, it computes either:
+
 - **Tidal displacements** (vertical, horizontal components) - when displacement coefficients are provided
 - **Tidal gravity effects** (gravity disturbances) - when gravity-type coefficients are provided
 
@@ -25,6 +26,7 @@ The program is part of the **IERS Conventions 2010** recommendations (Class 1 mo
 Given ocean loading coefficients in BLQ format (from the Bos-Scherneck loading service or equivalent providers), HARDISP computes a time series of tidal ocean loading effects in three components:
 
 **Output type depends on the ocean loading provider's coefficients:**
+
 - **Displacement coefficients**: Returns displacements
   - **dU**: Vertical (radial) displacement (meters)
   - **dS**: South (North-South) displacement (meters)
@@ -44,8 +46,8 @@ HARDISP processes the coefficients as-is, with no automatic conversions between
 The Python conversion consists of:
 **core.py** - Main Python module with optimized implementation
 
-
 Original Fortran files (for reference):
+
 - `HARDISP.F` - Main program
 - `ADMINT.F` - Admittance interpolation
 - `RECURS.F` - Recursive harmonic evaluation
@@ -191,6 +193,7 @@ for day_val in range(1, 8):  # Process days 1-7
 ### Coordinate System
 
 The output displacements are in a **local geodetic frame** at the station:
+
 - **dU**: Radial (positive upward)
 - **dS**: South component (positive southward, i.e., negative latitude direction)
 - **dW**: West component (positive westward, i.e., negative longitude direction)
@@ -198,6 +201,7 @@ The output displacements are in a **local geodetic frame** at the station:
 ### Tidal Constituents Used
 
 The 11 input harmonics (from BLQ format) represent:
+
 1. M₂ - Principal lunar semi-diurnal (12.42 hours)
 2. S₂ - Principal solar semi-diurnal (12.00 hours)
 3. N₂ - Lunar elliptic semi-diurnal (12.66 hours)
@@ -252,7 +256,6 @@ The Python code was converted from these 13 Fortran files:
 12. `ETUTC.F` (284 lines) - ET-UTC calculation
 13. `TOYS.F` - (not in workspace but referenced)
 
-
 ## Function Name Mapping (Original Fortran → New Python)
 
 | Original Fortran | New Python Name | Purpose |
@@ -276,6 +279,7 @@ The Python code was converted from these 13 Fortran files:
 ## Features Implemented
 
 ### Complete Conversion
+
 - [x] All 342 tidal constituents support
 - [x] Recursive harmonic evaluation (efficient)
 - [x] Cubic spline interpolation
@@ -287,6 +291,7 @@ The Python code was converted from these 13 Fortran files:
 - [x] Delaunay argument calculations
 
 ### Additional Features
+
 - [x] Object-oriented design with HardispComputer class
 - [x] NumPy integration for high performance
 - [x] Comprehensive error handling
@@ -297,6 +302,7 @@ The Python code was converted from these 13 Fortran files:
 ## Core Components Converted
 
 ### 1. Date/Time Functions (7 functions)
+
 ```
 ✓ is_leap_year(year) - Check leap year
 ✓ days_before_month(year, month) - Days before month start
@@ -306,6 +312,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 2. Spline Interpolation (3 functions)
+
 ```
 ✓ cublic_spline(x, u) - Compute cubic spline coefficients
 ✓ spline_eval(y, x, u, s) - Evaluate spline at point
@@ -313,6 +320,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 3. Tidal Frequency Calculations (3 functions)
+
 ```
 ✓ calculate_tidal_arguments(year, day, h, m, s) - Initialize tidal calculations
 ✓ tidal_frequency_and_phase(doodson_number) - Get frequency and phase from Doodson number
@@ -320,33 +328,34 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 4. Harmonic Recursion (1 function)
+
 ```
 ✓ recursion(n, hc, nf, om) - Efficient recursive harmonic evaluation
 ```
 
 ### 5. Utility Functions (2 functions)
+
 ```
 ✓ pyshells(x) - Sort array with indices (Shell sort)
 ✓ fortran_int_divide(a, b) - Fortran-style integer division
 ```
 
 ### 6. Main Class: HardispComputer
+
 ```
 ✓ read_blq_format(amp, phase) - Load ocean loading coefficients
 ✓ compute_ocean_loading(...) - Main computation engine
 ```
 
-
-
 ## References
 
 1. Petit, G. and Luzum, B. (eds.), **IERS Conventions (2010)**, IERS Technical Note No. 36, BKG (2010)
-   - Available at: https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php
+   - Available at: <https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php>
 
 2. Agnew, D. C., et al., **HARDISP**: Original algorithm and implementation
 
 3. Scherneck, H.-G., and M. S. Bos, **Ocean Loading Service**: BLQ format specification
-   - Available at: http://www.oso.chalmers.se/~loading/
+   - Available at: <http://www.oso.chalmers.se/~loading/>
 
 4. Cartwright, D. E., and A. C. Edden, **Tides of the Planet Earth**, Geophys. J. R. Astron. Soc. 65, 615-630, 1981
 
@@ -357,10 +366,11 @@ This Python conversion maintains the same IERS Conventions Software License as t
 ## Contact
 
 For questions or issues with the Fortran original:
-- IERS Conventions Center: https://www.iers.org/
-- Email: gpetit@bipm.org or brian.luzum@usno.navy.mil
+
+- IERS Conventions Center: <https://www.iers.org/>
+- Email: <gpetit@bipm.org> or <brian.luzum@usno.navy.mil>
 
 ---
 
-*Python conversion completed: 2024*  
+*Python conversion completed: 2026*  
 *Original Fortran: IERS Conventions 2010*

{pyhardisp-0.2.0 → pyhardisp-0.2.2}/pyhardisp/__init__.py

@@ -54,4 +54,4 @@ __all__ = [
     "load_ocean_loading_coefficients",
 ]
 
-__version__ = "0.2.0"
+__version__ = "0.2.2"

{pyhardisp-0.2.0 → pyhardisp-0.2.2/pyhardisp.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyhardisp
-Version: 0.2.0
+Version: 0.2.2
 Summary: Ocean loading tidal displacement calculator - Python conversion of IERS HARDISP Fortran program
 Author: IERS Conventions, Python Conversion - Craig Miller
 License-Expression: LicenseRef-IERS
@@ -37,6 +37,7 @@ This document describes the Python conversion of the **HARDISP** (Harmonic Displ
 ## What is HARDISP?
 
 HARDISP is a scientific program developed by the International Earth Rotation and Reference Systems Service (IERS) for calculating tidal ocean loading effects at geodetic stations. Depending on the coefficients provided by the ocean loading provider, it computes either:
+
 - **Tidal displacements** (vertical, horizontal components) - when displacement coefficients are provided
 - **Tidal gravity effects** (gravity disturbances) - when gravity-type coefficients are provided
 
@@ -55,6 +56,7 @@ The program is part of the **IERS Conventions 2010** recommendations (Class 1 mo
 Given ocean loading coefficients in BLQ format (from the Bos-Scherneck loading service or equivalent providers), HARDISP computes a time series of tidal ocean loading effects in three components:
 
 **Output type depends on the ocean loading provider's coefficients:**
+
 - **Displacement coefficients**: Returns displacements
   - **dU**: Vertical (radial) displacement (meters)
   - **dS**: South (North-South) displacement (meters)
@@ -74,8 +76,8 @@ HARDISP processes the coefficients as-is, with no automatic conversions between
 The Python conversion consists of:
 **core.py** - Main Python module with optimized implementation
 
-
 Original Fortran files (for reference):
+
 - `HARDISP.F` - Main program
 - `ADMINT.F` - Admittance interpolation
 - `RECURS.F` - Recursive harmonic evaluation
@@ -221,6 +223,7 @@ for day_val in range(1, 8):  # Process days 1-7
 ### Coordinate System
 
 The output displacements are in a **local geodetic frame** at the station:
+
 - **dU**: Radial (positive upward)
 - **dS**: South component (positive southward, i.e., negative latitude direction)
 - **dW**: West component (positive westward, i.e., negative longitude direction)
@@ -228,6 +231,7 @@ The output displacements are in a **local geodetic frame** at the station:
 ### Tidal Constituents Used
 
 The 11 input harmonics (from BLQ format) represent:
+
 1. M₂ - Principal lunar semi-diurnal (12.42 hours)
 2. S₂ - Principal solar semi-diurnal (12.00 hours)
 3. N₂ - Lunar elliptic semi-diurnal (12.66 hours)
@@ -282,7 +286,6 @@ The Python code was converted from these 13 Fortran files:
 12. `ETUTC.F` (284 lines) - ET-UTC calculation
 13. `TOYS.F` - (not in workspace but referenced)
 
-
 ## Function Name Mapping (Original Fortran → New Python)
 
 | Original Fortran | New Python Name | Purpose |
@@ -306,6 +309,7 @@ The Python code was converted from these 13 Fortran files:
 ## Features Implemented
 
 ### Complete Conversion
+
 - [x] All 342 tidal constituents support
 - [x] Recursive harmonic evaluation (efficient)
 - [x] Cubic spline interpolation
@@ -317,6 +321,7 @@ The Python code was converted from these 13 Fortran files:
 - [x] Delaunay argument calculations
 
 ### Additional Features
+
 - [x] Object-oriented design with HardispComputer class
 - [x] NumPy integration for high performance
 - [x] Comprehensive error handling
@@ -327,6 +332,7 @@ The Python code was converted from these 13 Fortran files:
 ## Core Components Converted
 
 ### 1. Date/Time Functions (7 functions)
+
 ```
 ✓ is_leap_year(year) - Check leap year
 ✓ days_before_month(year, month) - Days before month start
@@ -336,6 +342,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 2. Spline Interpolation (3 functions)
+
 ```
 ✓ cublic_spline(x, u) - Compute cubic spline coefficients
 ✓ spline_eval(y, x, u, s) - Evaluate spline at point
@@ -343,6 +350,7 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 3. Tidal Frequency Calculations (3 functions)
+
 ```
 ✓ calculate_tidal_arguments(year, day, h, m, s) - Initialize tidal calculations
 ✓ tidal_frequency_and_phase(doodson_number) - Get frequency and phase from Doodson number
@@ -350,33 +358,34 @@ The Python code was converted from these 13 Fortran files:
 ```
 
 ### 4. Harmonic Recursion (1 function)
+
 ```
 ✓ recursion(n, hc, nf, om) - Efficient recursive harmonic evaluation
 ```
 
 ### 5. Utility Functions (2 functions)
+
 ```
 ✓ pyshells(x) - Sort array with indices (Shell sort)
 ✓ fortran_int_divide(a, b) - Fortran-style integer division
 ```
 
 ### 6. Main Class: HardispComputer
+
 ```
 ✓ read_blq_format(amp, phase) - Load ocean loading coefficients
 ✓ compute_ocean_loading(...) - Main computation engine
 ```
 
-
-
 ## References
 
 1. Petit, G. and Luzum, B. (eds.), **IERS Conventions (2010)**, IERS Technical Note No. 36, BKG (2010)
-   - Available at: https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php
+   - Available at: <https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.php>
 
 2. Agnew, D. C., et al., **HARDISP**: Original algorithm and implementation
 
 3. Scherneck, H.-G., and M. S. Bos, **Ocean Loading Service**: BLQ format specification
-   - Available at: http://www.oso.chalmers.se/~loading/
+   - Available at: <http://www.oso.chalmers.se/~loading/>
 
 4. Cartwright, D. E., and A. C. Edden, **Tides of the Planet Earth**, Geophys. J. R. Astron. Soc. 65, 615-630, 1981
 
@@ -387,10 +396,11 @@ This Python conversion maintains the same IERS Conventions Software License as t
 ## Contact
 
 For questions or issues with the Fortran original:
-- IERS Conventions Center: https://www.iers.org/
-- Email: gpetit@bipm.org or brian.luzum@usno.navy.mil
+
+- IERS Conventions Center: <https://www.iers.org/>
+- Email: <gpetit@bipm.org> or <brian.luzum@usno.navy.mil>
 
 ---
 
-*Python conversion completed: 2024*  
+*Python conversion completed: 2026*  
 *Original Fortran: IERS Conventions 2010*

{pyhardisp-0.2.0 → pyhardisp-0.2.2}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pyhardisp"
-version = "0.2.0"
+version = "0.2.2"
 description = "Ocean loading tidal displacement calculator - Python conversion of IERS HARDISP Fortran program"
 readme = "README_PYTHON_CONVERSION.md"
 requires-python = ">=3.8"