Fix particle tracking bug, add version to cli output

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "tc1d"
-version = "0.3.0"
+version = "0.3.1"
 authors = [
     {name = "Dave Whipp", email = "david.whipp@helsinki.fi"},
     {name= "Benjamin Gérard", email = "benjamin.gerard@helsinki.fi"},
@@ -126,3 +126,4 @@ ruff = ">=0.12.12,<0.13"
 hatchling = ">=1.27.0,<2"
 mypy = ">=1.17.1,<2"
 pip = ">=25.2,<26"
+pre-commit = ">=4.3.0,<5"

src/tc1d/tc1d.py

@@ -3209,6 +3209,9 @@ def run_model(params):
             # Loop over all times and use -dt to run erosion model backwards
             # TODO: Make this a function???
             while curtime > 0.0:
+                # Increment current time
+                curtime -= dt
+
                 # Find particle velocities and move incrementally to starting depths
                 vx_pts, vx, vx_max, fault_depth = calculate_erosion_rate(
                     params,
@@ -3228,8 +3231,7 @@ def run_model(params):
                 if params["plot_fault_depth_history"]:
                     fault_depth_history[idx] = fault_depth
 
-                # Increment current time and idx
-                curtime -= dt
+                # Increment idx
                 idx += 1
 
             # Reverse exhumation history for plotting

src/tc1d/tc1d_cli.py

@@ -2,6 +2,7 @@
 
 # Import libraries we need
 import argparse
+from importlib.metadata import version
 import sys
 import tc1d
 
@@ -15,6 +16,7 @@ def main():
         description="Calculates transient 1D temperatures and thermochronometer ages",
         formatter_class=argparse.ArgumentDefaultsHelpFormatter,
     )
+    parser.add_argument("-v", "--version", action="version", version=version("tc1d"))
     general = parser.add_argument_group(
         "General options", "Options for various general features"
     )