diff --git a/slcl1butils/scripts/do_IW_L1C_SAFE_from_L1B_SAFE.py b/slcl1butils/scripts/do_IW_L1C_SAFE_from_L1B_SAFE.py
index 1fa8bfb..722d1c7 100644
--- a/slcl1butils/scripts/do_IW_L1C_SAFE_from_L1B_SAFE.py
+++ b/slcl1butils/scripts/do_IW_L1C_SAFE_from_L1B_SAFE.py
@@ -271,8 +271,9 @@ def append_ancillary_field(ancillary, ds_intra, ds_inter):
     else:
         raise ValueError("%s ancillary name not handled" % ancillary["name"])
     # Get the polygons of the swath data
+    first_pola_available = ds_intra.coords['pol'].data[0]
     polygons, coordinates, variables = get_swath_tiles_polygons_from_l1bgroup(
-        ds_intra, swath_only=True
+        ds_intra,polarization=first_pola_available, swath_only=True
     )
     # Crop the raster to the swath bounding box limit
     raster_bb_ds = raster_cropping_in_polygon_bounding_box(
diff --git a/slcl1butils/scripts/do_WV_L1C_SAFE_from_L1B_SAFE.py b/slcl1butils/scripts/do_WV_L1C_SAFE_from_L1B_SAFE.py
index 623230d..c4bd30b 100644
--- a/slcl1butils/scripts/do_WV_L1C_SAFE_from_L1B_SAFE.py
+++ b/slcl1butils/scripts/do_WV_L1C_SAFE_from_L1B_SAFE.py
@@ -224,8 +224,9 @@ def append_ancillary_field(ancillary, ds_intra):
             raster_ds = ww3_global_yearly_3h(filename, closest_date)
 
         # Get the polygons of the swath data
+        first_pola_available = ds_intra.coords['pol'].data[0]
         polygons, coordinates, variables = get_swath_tiles_polygons_from_l1bgroup(
-            ds_intra, swath_only=True
+            ds_intra,polarization=first_pola_available, swath_only=True
         )
         # Crop the raster to the swath bounding box limit