Merge pull request #704 from akrherz/artcc

add manual ARTCC assignment for PIREPs
akrherz · Feb 13, 2023 · 954aa45 · 954aa45
2 parents 34d4e67 + 4b06d36
commit 954aa45
Showing 1 changed file with 11 additions and 3 deletions.
diff --git a/src/pyiem/nws/products/pirep.py b/src/pyiem/nws/products/pirep.py
@@ -315,12 +315,20 @@ def sql(self, txn):
                 continue
             if report.longitude is None:
                 geom = "POINT EMPTY"
+                artcc = "null"
             else:
                 geom = f"SRID=4326;POINT({report.longitude} {report.latitude})"
+                artcc = (
+                    "(select ident from airspaces where st_dwithin(geom, "
+                    f"ST_MakePoint({report.longitude}, {report.latitude}), 0) "
+                    "and type_code = 'ARTCC' LIMIT 1)"
+                )
             txn.execute(
-                "INSERT into pireps(valid, geom, is_urgent, "
-                "aircraft_type, report) VALUES (%s, "
-                "ST_GeographyFromText(%s), %s, %s, %s)",
+                f"""
+                INSERT into pireps(valid, geom, is_urgent, aircraft_type,
+                report, artcc) VALUES (%s,
+                ST_GeographyFromText(%s), %s, %s, %s, {artcc})
+                """,
                 (
                     report.valid,
                     geom,