There's still some hardcoding of coordinates

README.md

@@ -1 +1,20 @@
 # Canada datastore
+
+An example crontab call for the script is:
+
+```bash
+SHELL=/bin/bash
+BASH_ENV=~/.bashrc_conda
+LOC=/home/jose/data/Canada_datastore_2024
+EUMETSAT_KEY=wZbhgx5sf1qMYKsCpJ9Bf8lE7Aga
+EUMETSAT_SECRET=DrBYHdM3AQwhxFiMaR0MIRbbX48a
+QT_QPA_PLATFORM=offscreen
+PATH=/bin:/usr/bin:/home/jose/mamba/bin/
+MAILTO=jose.gomez-dans@kcl.ac.uk
+
+0       3       *       *       *       conda activate py312 ; get_arable -f /home/jose/data/KCL_Arable/raw_data/ >> /home/jose/data/KCL_Arable/logs/logfile.$(date +\%Y-\%m-\%d).log 2>&1
+0 6 * * * conda activate py312 ; get_sentinel3 -f1 ${LOC}/Sentinel3 -ff ${LOC}/FIRMS/ >> ${LOC}/logs/logfile.$(date +\%Y-\%m-\%d-\%H-\%M).log 2>&1
+0 12 * * * conda activate py312 ; get_sentinel3 -f1 ${LOC}/Sentinel3 -ff ${LOC}/FIRMS/ >> ${LOC}/logs/logfile.$(date +\%Y-\%m-\%d-\%H-\%M).log 2>&1
+45 17 * * * conda activate py312 ; get_sentinel3 -f1 ${LOC}/Sentinel3 -ff ${LOC}/FIRMS/ >> ${LOC}/logs/logfile.$(date +\%Y-\%m-\%d-\%H-\%M).log 2>&1
+
+```

src/canada_datastore/utils.py

@@ -17,7 +17,7 @@
 logger = logging.getLogger("canada_datastore")
 
 # URL = "http://10.81.205.13:3333/Canada_datastore/"
-URL = "https://gws-access.jasmin.ac.uk/public/global_fire_models/"
+URL = "https://gws-access.jasmin.ac.uk/public/global_fire_models/Canada_datastore_2024/"
 # URL = "/home/jose/data/"
 
 
@@ -113,8 +113,7 @@ def add_layers_to_project(
             layer = QgsVectorLayer(the_url, f"FIRMS {label}", "ogr")
             layerz.append(layer)
             if layer.isValid():
-                print(f"Adding {label}")
-                print(the_url)
+                print(f"Adding {label} <> {the_url}")
                 project.addMapLayer(
                     layer, False
                 )  # Add the layer without adding it to the layer tree
@@ -126,6 +125,7 @@ def add_layers_to_project(
         # Add layers to the group based on their type (raster or vector)
         for group_label, [label, the_url] in layer_files.items():
             if label.find("F1") >= 0:
+                print(f"Adding {label} <> {the_url}")
                 layer = add_qgis_raster_layer(the_url, label)
                 layerz.append(layer)
                 project.addMapLayer(layer, False)
@@ -139,7 +139,7 @@ def add_layers_to_project(
         print(the_url, label)
         layerz.append(layer)
         if layer.isValid():
-            print(f"Adding {label}")
+            print(f"Adding {label} <> {the_url}")
             project.addMapLayer(
                 layer, False
             )  # Add the layer without adding it to the layer tree
@@ -233,9 +233,6 @@ def create_all_project_files(
         l1_today = {k: v for k, v in l1.items() if k.date() == today}
         l2_today = {k: v for k, v in l2.items() if k.date() == today}
         firms_today = {k: v for k, v in f1.items() if k.date() == today}
-        import pdb
-
-        pdb.set_trace()
         add_layers_to_project(l1_today, l2_today, firms_today, output_fname)
 
     qgs.exit()