diff --git a/src/ctapipe_io_lst/__init__.py b/src/ctapipe_io_lst/__init__.py
index 0b1e05d0..d763b0e2 100644
--- a/src/ctapipe_io_lst/__init__.py
+++ b/src/ctapipe_io_lst/__init__.py
@@ -45,7 +45,7 @@
     parse_tib_10MHz_counter,
 )
 from .constants import (
-    HIGH_GAIN, LST_LOCATIONS, N_GAINS, N_PIXELS, N_SAMPLES, LST1_LOCATION, REFERENCE_LOCATION,
+    HIGH_GAIN, LST_LOCATIONS, N_GAINS, N_PIXELS, N_SAMPLES, REFERENCE_LOCATION,
     PixelStatus, TriggerBits,
 )
 
@@ -367,7 +367,8 @@ def datalevels(self):
     @staticmethod
     def create_subarray(tel_id=1, reference_location=None):
         """
-        Obtain the subarray from the EventSource
+        Obtain a single-lst subarray description
+
         Returns
         -------
         ctapipe.instrument.SubarrayDescription
@@ -408,7 +409,7 @@ def create_subarray(tel_id=1, reference_location=None):
             name=f"LST-{tel_id} subarray",
             tel_descriptions=tel_descriptions,
             tel_positions=tel_positions,
-            reference_location=LST1_LOCATION,
+            reference_location=reference_location,
         )
 
         return subarray