From 2f4d3c81f8a8e38e01026e8d91b8eef8f44be75c Mon Sep 17 00:00:00 2001
From: Enrico Ferreguti <enricofer@gmail.com>
Date: Thu, 21 Oct 2021 18:03:56 +0200
Subject: [PATCH] detect grib wind  dateset metadata capabilities

---
 wind_forecast_routing_algorithm.py |  2 ++
 wind_forecast_routing_launch.py    | 17 +++++++++--------
 2 files changed, 11 insertions(+), 8 deletions(-)

diff --git a/wind_forecast_routing_algorithm.py b/wind_forecast_routing_algorithm.py
index 0376be5..2e39b18 100644
--- a/wind_forecast_routing_algorithm.py
+++ b/wind_forecast_routing_algorithm.py
@@ -149,6 +149,8 @@ def __init__(self, gribLayer, wind_idx, destinationCrs=None):
             destinationCrs = QgsCoordinateReferenceSystem(4326)
         transform = QgsCoordinateTransform(self.grib.crs(), destinationCrs, QgsProject.instance().transformContext())
         self.grib.updateTriangularMesh(transform)
+        metadata = self.grib.datasetGroupMetadata(self.grib.datasetIndexAtRelativeTime(QgsInterval(1),wind_idx))
+        print ("GRIB-METADATA",metadata.dataType(),metadata.isScalar(),metadata.isTemporal(),metadata.isVector())
         self.wind_idx = wind_idx
 
     def getWindAt(self, t, lat, lon):
diff --git a/wind_forecast_routing_launch.py b/wind_forecast_routing_launch.py
index 0df163a..b40a43c 100644
--- a/wind_forecast_routing_launch.py
+++ b/wind_forecast_routing_launch.py
@@ -96,14 +96,6 @@
         "limits": None,
         "resolution": 0.5
     },
-    {
-        "service": 'icon_eu_p06_',
-        "context": "Europe - ICON_EU",
-        "interval": '1',
-        "days": '5',
-        "limits": [-23.5, 29.5, 45.0, 70.5],
-        "resolution": 0.06
-    },
     {
         "service": 'arpege_eu_p10_',
         "context": "Europe - ARPEGE_EU",
@@ -112,6 +104,14 @@
         "limits": [-32.0, 20.0, 42.0, 72.0],
         "resolution": 0.1
     },
+    {
+        "service": 'icon_eu_p06_',
+        "context": "Europe - ICON_EU",
+        "interval": '1',
+        "days": '5',
+        "limits": [-23.5, 29.5, 45.0, 70.5],
+        "resolution": 0.06
+    },
     {
         "service": 'nam_conus_12km_',
         "context": "North America - NAM_CONUS",
@@ -306,6 +306,7 @@ def processAlgorithm(self, parameters, context, feedback):
         rawReplyObject = manager.blockingGet(request)
         j = QJsonDocument.fromJson(rawReplyObject.content())
         replyObject = j.toVariant()
+        print("replyObject", replyObject)
         if replyObject["status"]:
 
             download_params = {