debug

Author: avsm

scripts/debug_tile_lookup.py

@@ -2,6 +2,7 @@
 """Debug script to investigate tile lookup failures."""
 
 import sys
+import platform
 import numpy as np
 
 def test_tile(registry, year, lon, lat):
@@ -50,6 +51,9 @@ def main():
     print("=" * 60)
     print("TILE LOOKUP DEBUG SCRIPT")
     print("=" * 60)
+    print(f"Platform: {platform.system()} {platform.release()}")
+    print(f"Python: {sys.version}")
+    print(f"NumPy: {np.__version__}")
 
     # Target tile that's failing
     target_year = 2024
@@ -70,6 +74,17 @@ def main():
     print(f"Embeddings registry: {len(gdf)} rows")
     print(f"Landmasks registry: {len(ldf) if ldf is not None else 'None'} rows")
 
+    # Check for negative lon_i values (potential Windows issue)
+    neg_lon_tiles = gdf[gdf['lon_i'] < 0]
+    print(f"Tiles with negative lon_i: {len(neg_lon_tiles)}")
+    if len(neg_lon_tiles) > 0:
+        print(f"  Sample negative lon_i values: {sorted(neg_lon_tiles['lon_i'].unique())[:10]}")
+        # Check if -15 exists
+        lon_minus_15 = gdf[gdf['lon_i'] == -15]
+        print(f"  Tiles with lon_i == -15: {len(lon_minus_15)}")
+        lon_minus_5 = gdf[gdf['lon_i'] == -5]
+        print(f"  Tiles with lon_i == -5: {len(lon_minus_5)}")
+
     # Check columns
     print("\n--- Registry Columns ---")
     print(f"Embeddings columns: {list(gdf.columns)}")