debug

Author: avsm

scripts/debug_tile_lookup.py

@@ -246,7 +246,22 @@ def main():
     print("\nFirst few raw tiles (showing lon, lon_i, lon_i/100):")
     for idx, row in raw_unique.head(5).iterrows():
         lon_back = row['lon_i'] / 100.0
-        print(f"  lon={row['lon']!r}, lon_i={row['lon_i']}, lon_i/100={lon_back!r} (hex: {lon_back.hex()})")
+        print(f"  lon={row['lon']!r}, lon_i={row['lon_i']} (type: {type(row['lon_i']).__name__}), lon_i/100={lon_back!r}")
+
+    # Also check what .values gives us
+    print("\nChecking .values array types:")
+    values_array = raw_unique[["lon_i", "lat_i"]].values
+    print(f"  values_array.dtype: {values_array.dtype}")
+    print(f"  values_array[0]: {values_array[0]} (types: {type(values_array[0][0]).__name__}, {type(values_array[0][1]).__name__})")
+
+    # Check if there's float conversion happening
+    print("\nDirect iteration test:")
+    for year_val, lon_val, lat_val, lon_i_val, lat_i_val in raw_unique.values[:3]:
+        print(f"  lon_i_val={lon_i_val!r} (type: {type(lon_i_val).__name__})")
+        divided = lon_i_val / 100.0
+        print(f"    lon_i_val / 100.0 = {divided!r} (type: {type(divided).__name__})")
+        back_to_int = coord_to_grid_int(divided)
+        print(f"    coord_to_grid_int({divided!r}) = {back_to_int!r} (type: {type(back_to_int).__name__})")
 
     tiles = list(registry.iter_tiles_in_region(bounds, year))
     print(f"\nFound {len(tiles)} tiles in region (from iter_tiles_in_region)")
@@ -258,8 +273,8 @@ def main():
         lon_i_computed = coord_to_grid_int(tile_lon)
         lat_i_computed = coord_to_grid_int(tile_lat)
         print(f"\n[{i+1}/{len(tiles)}] year={tile_year}")
-        print(f"  lon={tile_lon!r} (hex: {tile_lon.hex()}) → lon_i={lon_i_computed}")
-        print(f"  lat={tile_lat!r} (hex: {tile_lat.hex()}) → lat_i={lat_i_computed}")
+        print(f"  lon={tile_lon!r} (hex: {tile_lon.hex()}) -> lon_i={lon_i_computed}")
+        print(f"  lat={tile_lat!r} (hex: {tile_lat.hex()}) -> lat_i={lat_i_computed}")
         print(f"  lon*100={tile_lon * 100!r}, round={np.round(tile_lon * 100)}")
 
         # Test get_tile_file_size