profile layer limits also included

pyaerocom/aeroval/obsentry.py

@@ -72,6 +72,7 @@ def __init__(self, **kwargs):
         self.is_superobs = False
         self.only_superobs = False
         self.colocation_layer_limts = None
+        self.profile_layer_limits = None
 
         self.read_opts_ungridded = {}
 

pyaerocom/colocation_3d.py

@@ -59,6 +59,7 @@ def colocate_vertical_profile_gridded(
     use_climatology_ref=False,
     resample_how=None,
     colocation_layer_limits=None,
+    profile_layer_limits=None,
     **kwargs,
 ):
     """
@@ -93,7 +94,11 @@ def colocate_vertical_profile_gridded(
 
     if not all(["start" and "end" in keys for keys in colocation_layer_limits]):
         raise KeyError(
-            "start and end must be provided for each vertical layer in colocate_vertical_profile_gridded"
+            "start and end must be provided for colocation in each vertical layer in colocate_vertical_profile_gridded"
+        )
+    if not all(["start" and "end" in keys for keys in profile_layer_limits]):
+        raise KeyError(
+            "start and end must be provided for displaying profiles in each vertical layer in colocate_vertical_profile_gridded"
         )
 
     dataset_ref = data_ref.contains_datasets[0]
@@ -259,6 +264,9 @@ def colocate_vertical_profile_gridded(
 
             # Make a copy of the station data and resample it to the mean based on hourly resolution. Needs testing!
             tmp_obs_stat = obs_stat.copy()
+            # add the station altitude to the altitudes so everything is against Above Sea Level (ASL)
+            tmp_obs_stat.altitude += tmp_obs_stat.station_coords["altitude"]
+
             tmp_obs_stat[var_ref] = (
                 tmp_obs_stat[var_ref][
                     (vertical_layer["start"] <= tmp_obs_stat.altitude)

pyaerocom/colocation_auto.py

@@ -352,6 +352,7 @@ def __init__(
         self.obs_filters = {}
         self._obs_is_vertical_profile = False
         self.colocation_layer_limits = None
+        self.profile_layer_limits = None
 
         self.read_opts_ungridded = {}
 
@@ -1417,7 +1418,8 @@ def _prepare_colocation_args(self, model_var: str, obs_var: str):
             ts_type = self._get_colocation_ts_type(model_data.ts_type, obs_data.ts_type)
             args.update(ts_type=ts_type)
         if self.obs_is_vertical_profile:
-            args.update(colocation_layer_limits=self.colocation_layer_limits)
+            args.update(colocation_layer_limits=self.colocation_layer_limits,
+                        profile_layer_limits=self.profile_layer_limits)
         return args
 
     def _check_dimensionality(self, args):