diff --git a/cset-workflow/includes/basic_qq_plot.cylc b/cset-workflow/includes/basic_qq_plot.cylc
index ae792ff8e..15d9d4f21 100644
--- a/cset-workflow/includes/basic_qq_plot.cylc
+++ b/cset-workflow/includes/basic_qq_plot.cylc
@@ -31,3 +31,20 @@
         """
 {% endfor %}
 {% endif %}
+
+{# Example code to use as basis for multi-model cycling:
+{% for equivalent_field in restructure_field_list(PRESSURE_LEVEL_MODEL_FIELDS) %}
+{% for model_number, field in equivalent_field.items() %}
+[runtime]
+    [[generic_basic_qq_plot_m{{model_number}}_{{sanitise_task_name(field)}}]]
+    inherit = PROCESS
+        [[[environment]]]
+        CSET_RECIPE_NAME = "generic_basic_qq_plot.yaml"
+        CSET_ADDOPTS = """
+            --VARNAME='{{field}}'
+            --MODEL_NAME='{{models[model_number-1]["name"]}}'
+        """
+        MODEL_NUMBER = {{model_number}}
+{% endfor %}
+{% endfor %}
+#}
diff --git a/cset-workflow/includes/mlevel_transect.cylc b/cset-workflow/includes/mlevel_transect.cylc
new file mode 100644
index 000000000..710074830
--- /dev/null
+++ b/cset-workflow/includes/mlevel_transect.cylc
@@ -0,0 +1,21 @@
+{% if EXTRACT_MLEVEL_TRANSECT %}
+{% for equivalent_field in restructure_field_list(MODEL_LEVEL_MODEL_FIELDS) %}
+{% for model_number, field in equivalent_field.items() %}
+[runtime]
+    [[mlevel_transect_m{{model_number}}_{{sanitise_task_name(field)}}]]
+    inherit = PROCESS
+    # Transect calculation can take a while, so increase the time limit.
+    execution time limit = PT60M
+        [[[environment]]]
+        CSET_RECIPE_NAME = "transect.yaml"
+        CSET_ADDOPTS = """
+            --VARNAME='{{field}}'
+            --VERTICAL_COORDINATE='model_level_number'
+            --MODEL_NAME='{{models[model_number-1]["name"]}}'
+            --START_COORDS='{{MLEVEL_TRANSECT_STARTCOORDS}}'
+            --FINISH_COORDS='{{MLEVEL_TRANSECT_FINISHCOORDS}}'
+        """
+        MODEL_NUMBER = {{model_number}}
+{% endfor %}
+{% endfor %}
+{% endif %}
diff --git a/cset-workflow/includes/plevel_transect.cylc b/cset-workflow/includes/plevel_transect.cylc
new file mode 100644
index 000000000..18b99f71f
--- /dev/null
+++ b/cset-workflow/includes/plevel_transect.cylc
@@ -0,0 +1,21 @@
+{% if EXTRACT_PLEVEL_TRANSECT %}
+{% for equivalent_field in restructure_field_list(PRESSURE_LEVEL_MODEL_FIELDS) %}
+{% for model_number, field in equivalent_field.items() %}
+[runtime]
+    [[plevel_transect_m{{model_number}}_{{sanitise_task_name(field)}}]]
+    inherit = PROCESS
+    # Transect calculation can take a while, so increase the time limit.
+    execution time limit = PT60M
+        [[[environment]]]
+        CSET_RECIPE_NAME = "transect.yaml"
+        CSET_ADDOPTS = """
+            --VARNAME='{{field}}'
+            --VERTICAL_COORDINATE='pressure'
+            --MODEL_NAME='{{models[model_number-1]["name"]}}'
+            --START_COORDS='{{PLEVEL_TRANSECT_STARTCOORDS}}'
+            --FINISH_COORDS='{{PLEVEL_TRANSECT_FINISHCOORDS}}'
+        """
+        MODEL_NUMBER = {{model_number}}
+{% endfor %}
+{% endfor %}
+{% endif %}
diff --git a/cset-workflow/includes/transect.cylc b/cset-workflow/includes/transect.cylc
deleted file mode 100644
index 69ab854fa..000000000
--- a/cset-workflow/includes/transect.cylc
+++ /dev/null
@@ -1,18 +0,0 @@
-{% if EXTRACT_TRANSECT %}
-{% for var in CS_VARS %}
-[runtime]
-    [[parallel_transect_{{var}}]]
-    inherit = PARALLEL
-    execution time limit = PT60M
-        [[[environment]]]
-        CSET_RECIPE_NAME = "transect.yaml"
-        CSET_ADDOPTS = "--CS_STARTCOORDS='{{CS_STARTCOORDS}}' --CS_FINISHCOORDS='{{CS_FINISHCOORDS}}' --CS_VAR='{{var}}' --CS_VERTLEV='{{CS_VERTLEV}}'"
-
-    [[collate_transect_{{var}}]]
-    inherit = COLLATE
-    execution time limit = PT60M
-        [[[environment]]]
-        CSET_RECIPE_NAME = "transect.yaml"
-        CSET_ADDOPTS = "--CS_STARTCOORDS='{{CS_STARTCOORDS}}' --CS_FINISHCOORDS='{{CS_FINISHCOORDS}}' --CS_VAR='{{var}}' --CS_VERTLEV='{{CS_VERTLEV}}'"
-{% endfor %}
-{% endif %}
diff --git a/cset-workflow/meta/diagnostics/rose-meta.conf b/cset-workflow/meta/diagnostics/rose-meta.conf
index cb2f69fa5..e3d672d1a 100644
--- a/cset-workflow/meta/diagnostics/rose-meta.conf
+++ b/cset-workflow/meta/diagnostics/rose-meta.conf
@@ -42,6 +42,42 @@ type=python_boolean
 compulsory=true
 sort-key=0surface4
 
+[template variables=SURFACE_SINGLE_POINT_TIME_SERIES]
+ns=Diagnostics/Quicklook
+description=Plot a time series at a single specified location in a surface field.
+help=Include values of desired longitude and latitude.
+type=python_boolean
+compulsory=true
+trigger=template variables=LATITUDE_POINT: True;
+        template variables=LONGITUDE_POINT: True;
+        template variables=SINGLE_POINT_METHOD: True;
+sort-key=0surface5
+
+[template variables=LATITUDE_POINT]
+ns=Diagnostics/Quicklook
+description=Latitude of selected point. Note that this could be rotated or not, depending on the data provided.
+help=The latitude must exist within the domain. Value should be a float: for example, -1.5.
+type=real
+compulsory=true
+sort-key=0surface6
+
+[template variables=LONGITUDE_POINT]
+ns=Diagnostics/Quicklook
+description=Longitude of selected point. Note that this could be rotated or not, depending on the data provided.
+help=The longitude must exist within the domain. Value should be a float: for example, 0.8.
+type=real
+compulsory=true
+sort-key=0surface6
+
+[template variables=SINGLE_POINT_METHOD]
+ns=Diagnostics/Quicklook
+description=Method used to map model data onto selected gridpoints.
+help=Method used to map model data onto selected gridpoints. These are regrid methods available in Iris.
+values="Nearest", "Linear"
+compulsory=true
+sort-key=0surface6
+
+
 # Pressure level fields.
 [template variables=PRESSURE_LEVEL_MODEL_FIELDS]
 ns=Diagnostics/Quicklook
@@ -97,6 +133,44 @@ values="step","barstacked","stepfilled"
 compulsory=true
 sort-key=1pressure4b
 
+[template variables=EXTRACT_PLEVEL_TRANSECT]
+ns=Diagnostics/Quicklook
+description=Plot a cross section of pressure for variables.
+help=This functionality extracts data for each variable and extracts a column
+     along each point of the transect. The number of points along the transect
+     is determined by calculating the length of the transect in degrees, and
+     dividing this by the minimum grid spacing within the domain. Additional
+     dimensions (time, ensemble) are accepted and returned as a cross section.
+     This allows the user to compute slices through features of meteorological
+     interest, such as transects through valleys, or through a front as it
+     passes through.
+type=python_boolean
+trigger=template variables=PLEVEL_TRANSECT_STARTCOORDS: True;
+        template variables=PLEVEL_TRANSECT_FINISHCOORDS: True;
+compulsory=true
+sort-key=1pressure5
+
+[template variables=PLEVEL_TRANSECT_STARTCOORDS]
+ns=Diagnostics/Quicklook
+description=Start latitude, longitude of the cross section.
+help=The latitude, longitude coordinate with respect to the model grid where the
+     cross section will start i.e. the furthest left hand point of the plot,
+     where the x axis is distance along transect, and y axis is pressure level).
+type=real,real
+compulsory=true
+sort-key=1pressure6
+
+[template variables=PLEVEL_TRANSECT_FINISHCOORDS]
+ns=Diagnostics/Quicklook
+description=Finish latitude, longitude of the cross section.
+help=The latitude, longitude coordinate with respect to the model grid where the
+     cross section will finish i.e. the furthest right hand point of the plot,
+     where the x axis is distance along transect, and y axis is pressure level).
+type=real,real
+compulsory=true
+sort-key=1pressure6
+
+
 # Model-level fields.
 [template variables=MODEL_LEVEL_MODEL_FIELDS]
 ns=Diagnostics/Quicklook
@@ -136,6 +210,43 @@ type=python_boolean
 compulsory=true
 sort-key=2modellevel4
 
+[template variables=EXTRACT_<LEVEL_TRANSECT]
+ns=Diagnostics/Quicklook
+description=Plot a cross section of model_level_number for variables.
+help=This functionality extracts data for each variable and extracts a column
+     along each point of the transect. The number of points along the transect
+     is determined by calculating the length of the transect in degrees, and
+     dividing this by the minimum grid spacing within the domain. Additional
+     dimensions (time, ensemble) are accepted and returned as a cross section.
+     This allows the user to compute slices through features of meteorological
+     interest, such as transects through valleys, or through a front as it
+     passes through.
+type=python_boolean
+trigger=template variables=MLEVEL_TRANSECT_STARTCOORDS: True;
+        template variables=MLEVEL_TRANSECT_FINISHCOORDS: True;
+compulsory=true
+sort-key=2modellevel5
+
+[template variables=MLEVEL_TRANSECT_STARTCOORDS]
+ns=Diagnostics/Quicklook
+description=Start latitude, longitude of the cross section.
+help=The latitude, longitude coordinate with respect to the model grid where the
+     cross section will start i.e. the furthest left hand point of the plot,
+     where the x axis is distance along transect, and y axis is pressure level).
+type=real,real
+compulsory=true
+sort-key=2modellevel6
+
+[template variables=MLEVEL_TRANSECT_FINISHCOORDS]
+ns=Diagnostics/Quicklook
+description=Finish latitude, longitude of the cross section.
+help=The latitude, longitude coordinate with respect to the model grid where the
+     cross section will finish i.e. the furthest right hand point of the plot,
+     where the x axis is distance along transect, and y axis is pressure level).
+type=real,real
+compulsory=true
+sort-key=2modellevel6
+
 
 # METplus.
 [template variables=METPLUS_POINT_STAT]
@@ -182,6 +293,8 @@ description=Path to directory containing analysis in MET netcdf format.
 type=quoted
 sort-key=met5
 
+
+# Convection diagnostics.
 [template variables=DETERMINISTIC_PLOT_CAPE_RATIO]
 ns=Diagnostics/Convection
 description=Extracts data required for, and calculates the CAPE ratio diagnostic, plotting on a map.
@@ -189,102 +302,63 @@ description=Extracts data required for, and calculates the CAPE ratio diagnostic
 help=See includes/deterministic_plot_cape_ratio.cylc
 type=python_boolean
 compulsory=true
-sort-key=cape1
 
-[template variables=EXTRACT_TRANSECT]
-ns=Diagnostics
-description=Plot a cross section for variables defined in CS_VARS.
-help=This functionality extracts data for each variable defined in CS_VARS,
-     and extracts a column along each point of the transect. The number of
-     points along the transect is determined by calculating the length of the
-     transect in degrees, and dividing this by the minimum grid spacing within the domain.
-     Additional dimensions (time, ensemble) are accepted and returned as a cross section.
-     This allows the user to compute slices through features of meteorological interest,
-     such as transects through valleys, or through a front as it passes through.
-type=python_boolean
-trigger=template variables=CS_STARTCOORDS: True;
-        template variables=CS_FINISHCOORDS: True;
-        template variables=CS_VARS: True;
-        template variables=CS_VERTLEV: True
-compulsory=true
-sort-key=cs1
-
-[template variables=CS_STARTCOORDS]
-ns=Diagnostics
-description=Start latitude, longitude of the cross section.
-help=The latitude, longitude coordinate with respect to the model grid where the
-     cross section will start i.e. the furthest left hand point of the plot, where the x
-     axis is distance along transect, and y axis is pressure level).
-type=real,real
-compulsory=true
-sort-key=cs2
-
-[template variables=CS_FINISHCOORDS]
-ns=Diagnostics
-description=Finish latitude, longitude of the cross section.
-help=The latitude, longitude coordinate with respect to the model grid where the
-     cross section will finish i.e. the furthest right hand point of the plot, where the x
-     axis is distance along transect, and y axis is pressure level).
-type=real,real
-compulsory=true
-sort-key=cs3
-
-[template variables=CS_VARS]
-ns=Diagnostics
-description=List of all variables to extract.
-help=A list of stash codes or variable names to extract from the data, and compute
-     the transect over. It must contain multi-level data (i.e. be at least 3D with pressure,
-     latitude and longitude coordinates).
-type=python_list
-compulsory=true
-sort-key=cs4
-
-[template variables=CS_VERTLEV]
-ns=Diagnostics
-description=Name of vertical coordinate to use for transect.
-help=A string, containing the vertical coordinate name to use for the transect. Typically
-     pressure or model_level_number, this is specified in case the source files contain
-     cubes where the variable name is the same, but they are on different vertical grids,
-     and it would be ambiguous which cube the user wants to compute the transect on.
-type=quoted
-compulsory=true
-sort-key=cs5
-
-[template variables=SURFACE_SINGLE_POINT_TIME_SERIES]
-ns=Diagnostics
-description=Plot a time series at a single specified location in a surface field.
-help=Include values of desired longitude and latitude.
+[template variables=DETERMINISTIC_PLOT_INFLOW_PROPERTIES]
+ns=Diagnostics/Convection
+description=Extracts data required for, and calculates the inflow properties diagnostic, plotting on a map.
+            Required STASH m01s20i119, m01s00i025, m01s00i033.
+help=See includes/deterministic_plot_inflow_properties.cylc
 type=python_boolean
 compulsory=true
-trigger=template variables=LATITUDE_POINT: True;
-        template variables=LONGITUDE_POINT: True;
-        template variables=SINGLE_POINT_METHOD: True;
-sort-key=point1
-
-[template variables=LATITUDE_POINT]
-ns=Diagnostics
-description=Latitude of selected point. Note that this could be rotated or not, depending on the data provided.
-help=The latitude must exist within the domain. Value should be a float: for example, -1.5.
-type=real
-compulsory=true
-sort-key=point2
-
-[template variables=LONGITUDE_POINT]
-ns=Diagnostics
-description=Longitude of selected point. Note that this could be rotated or not, depending on the data provided.
-help=The longitude must exist within the domain. Value should be a float: for example, 0.8.
-type=real
-compulsory=true
-sort-key=point2
 
-[template variables=SINGLE_POINT_METHOD]
-ns=Diagnostics
-description=Method used to map model data onto selected gridpoints.
-help=Method used to map model data onto selected gridpoints. These are regrid methods available in Iris.
-values="Nearest", "Linear"
-compulsory=true
-sort-key=point2
 
+# TODO: Reimplement in https://github.com/MetOffice/CSET/issues/807
+# [template variables=SELECT_SUBAREA]
+# ns=Diagnostics
+# description=Only analyse a subset of the region defined by the input data.
+# help=Select this option to run a recipe over a defined latitude-longitude range.
+# trigger=template variables=SUBAREA_LAT_BOUND_TOP: True;
+#        =template variables=SUBAREA_LAT_BOUND_BOTTOM: True;
+#        =template variables=SUBAREA_LON_BOUND_LEFT: True;
+#        =template variables=SUBAREA_LON_BOUND_RIGHT: True;
+# type=python_boolean
+# compulsory=true
+# sort-key=subsection1
+
+# [template variables=SUBAREA_LAT_BOUND_TOP]
+# ns=Diagnostics
+# description=Top edge coordinate of the sub-area, real.
+# help=Recommend looking at the input data to get these values. Uses the grid’s native units.
+# type=real
+# sort-key=subsection2
+# compulsory=true
+
+# [template variables=SUBAREA_LAT_BOUND_BOTTOM]
+# ns=Diagnostics
+# description=Bottom edge coordinate of the sub-area, real.
+# help=Recommend looking at the input data to get these values. Uses the grid’s native units.
+# type=real
+# sort-key=subsection2
+# compulsory=true
+
+# [template variables=SUBAREA_LON_BOUND_LEFT]
+# ns=Diagnostics
+# description=Left edge coordinate of the sub-area, real.
+# help=Recommend looking at the input data to get these values. Uses the grid’s native units.
+# type=real
+# sort-key=subsection2
+# compulsory=true
+
+# [template variables=SUBAREA_LON_BOUND_RIGHT]
+# ns=Diagnostics
+# description=Right edge coordinate of the sub-area, real.
+# help=Recommend looking at the input data to get these values. Uses the grid’s native units.
+# type=real
+# sort-key=subsection2
+# compulsory=true
+
+
+# TODO: Convert in https://github.com/MetOffice/CSET/issues/808
 [template variables=BASIC_QQ_PLOT]
 ns=Diagnostics
 description=Create a basic quantile-quantile plot for variables specified collapsing over specified coordinates.
@@ -352,9 +426,9 @@ sort-key=qq1d
 ns=Diagnostics
 description=Level coordinate name about which to constrain. The length of this
             list should be the same as MODEL_FIELDS B.
-help=Include a list of coordinate names in python list format["var1","var2"].
-     If single levels are required use any vertical coordinate in the cube and
-     make sure to put LEVELS_B as [].
+help=Include a list of coordinate names in python list format["var1","var2"]. If
+    single levels are required use any vertical coordinate in the cube and make
+    sure to put LEVELS_B as [].
 type=python_list
 compulsory=true
 sort-key=qq1e
@@ -364,17 +438,17 @@ ns=Diagnostics
 description=A list of levels to constrain over. The length of this
             list should be the same as MODEL_FIELDS B.
 help=Use "*" for all levels, INTEGERS for specific levels, and [] for no levels
-     (e.g., single level fields).
+    (e.g., single level fields).
 type=python_list
 compulsory=true
 sort-key=qq1f
 
 [template variables=COORDINATE_LIST]
 ns=Diagnostics
-description=List of coordinates to collapse the percentiles over. The coordinate list
-            must be of the form [\\"var1\\",\\"var2\\",...].The values should be
-            all of the coordinates in the cube to ensure the output is 1D with the
-            percentiles as the coordinate.
+description=List of coordinates to collapse the percentiles over. The coordinate
+    list must be of the form [\\"var1\\",\\"var2\\",...]. The values should be
+    all of the coordinates in the cube to ensure the output is 1D with the
+    percentiles as the coordinate.
 help=Include a list of coordinates in python list format[\\"var1\\",\\"var2\\"].
 type=quoted
 compulsory=true
@@ -387,55 +461,3 @@ help=Set to True if you want a one-to-one line; false if not.
 type=python_boolean
 compulsory=true
 sort-key=qq1h
-
-[template variables=SELECT_SUBAREA]
-ns=Diagnostics
-description=Only analyse a subset of the region defined by the input data.
-help=Select this option to run a recipe over a defined latitude-longitude range.
-trigger=template variables=SUBAREA_LAT_BOUND_TOP: True;
-       =template variables=SUBAREA_LAT_BOUND_BOTTOM: True;
-       =template variables=SUBAREA_LON_BOUND_LEFT: True;
-       =template variables=SUBAREA_LON_BOUND_RIGHT: True;
-type=python_boolean
-compulsory=true
-sort-key=subsection1
-
-[template variables=SUBAREA_LAT_BOUND_TOP]
-ns=Diagnostics
-description=Top edge coordinate of the sub-area, real.
-help=Recommend looking at the input data to get these values. Uses the grid’s native units.
-type=real
-sort-key=subsection2
-compulsory=true
-
-[template variables=SUBAREA_LAT_BOUND_BOTTOM]
-ns=Diagnostics
-description=Bottom edge coordinate of the sub-area, real.
-help=Recommend looking at the input data to get these values. Uses the grid’s native units.
-type=real
-sort-key=subsection2
-compulsory=true
-
-[template variables=SUBAREA_LON_BOUND_LEFT]
-ns=Diagnostics
-description=Left edge coordinate of the sub-area, real.
-help=Recommend looking at the input data to get these values. Uses the grid’s native units.
-type=real
-sort-key=subsection2
-compulsory=true
-
-[template variables=SUBAREA_LON_BOUND_RIGHT]
-ns=Diagnostics
-description=Right edge coordinate of the sub-area, real.
-help=Recommend looking at the input data to get these values. Uses the grid’s native units.
-type=real
-sort-key=subsection2
-compulsory=true
-
-[template variables=DETERMINISTIC_PLOT_INFLOW_PROPERTIES]
-ns=Diagnostics/Convection
-description=Extracts data required for, and calculates the inflow properties diagnostic, plotting on a map.
-            Required STASH m01s20i119, m01s00i025, m01s00i033.
-help=See includes/deterministic_plot_inflow_properties.cylc
-type=python_boolean
-compulsory=true
diff --git a/src/CSET/recipes/generic_basic_qq_plot.yaml b/src/CSET/recipes/generic_basic_qq_plot.yaml
index a36888e71..da0fff060 100644
--- a/src/CSET/recipes/generic_basic_qq_plot.yaml
+++ b/src/CSET/recipes/generic_basic_qq_plot.yaml
@@ -2,46 +2,9 @@ category: Quick Look
 title: basic q-q plot
 description: Extracts, quantiles, and plots two different variables ($VARNAME_A and $VARNAME_B) on a scatter plot.
 
-parallel:
+steps:
   - operator: read.read_cubes
-  - operator: filters.filter_multiple_cubes
-    constraint_A:
-        operator: constraints.combine_constraints
-        variable_constraint:
-            operator: constraints.generate_var_constraint
-            varname: $VARNAME_A
-        cell_method_constraint:
-            operator: constraints.generate_cell_methods_constraint
-            cell_methods: []
-        validity_time_constraint:
-            operator: constraints.generate_time_constraint
-            time_start: $VALIDITY_TIME
-        level_constraint:
-            operator: constraints.generate_level_constraint
-            coordinate: $VERTICAL_COORDINATE_A
-            levels: $LEVELS_A
-    constraint_B:
-        operator: constraints.combine_constraints
-        variable_constraint:
-            operator: constraints.generate_var_constraint
-            varname: $VARNAME_B
-        cell_method_constraint:
-            operator: constraints.generate_cell_methods_constraint
-            cell_methods: []
-        validity_time_constraint:
-            operator: constraints.generate_time_constraint
-            time_start: $VALIDITY_TIME
-        level_constraint:
-            operator: constraints.generate_level_constraint
-            coordinate: $VERTICAL_COORDINATE_B
-            levels: $LEVELS_B
 
-  - operator: write.write_cube_to_nc
-    filename: intermediate/qq_data
-
-collate:
-  - operator: read.read_cubes
-    filename_pattern: intermediate/*.nc
   - operator: plot.scatter_plot
     cube_y:
         operator: collapse.collapse
@@ -82,5 +45,6 @@ collate:
         method: 'PERCENTILE'
         additional_percent: [0,0.001,0.01,0.1,1,5,10,15,20,25,30,40,50,60,70,75,80,85,90,95,98,99,99.9,99.99,99.999,100]
     one-to-one: $ONE_TO_ONE
+
   - operator: write.write_cube_to_nc
     overwrite: True
diff --git a/src/CSET/recipes/transect.yaml b/src/CSET/recipes/transect.yaml
index 32dca1e63..0341a7f56 100644
--- a/src/CSET/recipes/transect.yaml
+++ b/src/CSET/recipes/transect.yaml
@@ -1,38 +1,31 @@
-title: Transect of $CS_VAR
+title: $MODEL_NAME Transect of $VARNAME
 description: |
-  Extracts a transect for a specified variable between two points and plots it.
+  Extracts a $VERTICAL_COORDINATE transect for $VARNAME between two points and
+  plots it.
 
-parallel:
+  Start coordinate: `$START_COORDS`
+  End coordinate: `$FINISH_COORDS`
+
+steps:
   - operator: read.read_cube
     constraint:
         operator: constraints.combine_constraints
-        validity_time_constraint:
-            operator: constraints.generate_time_constraint
-            time_start: $VALIDITY_TIME
+        cell_method_constraint:
+            operator: constraints.generate_cell_methods_constraint
+            cell_methods: []
         var_constraint:
             operator: constraints.generate_var_constraint
-            varname: $CS_VAR
+            varname: $VARNAME
         level_constraint:
             operator: constraints.generate_level_constraint
-            coordinate: $CS_VERTLEV
+            coordinate: $VERTICAL_COORDINATE
             levels: '*'
 
   - operator: transect.calc_transect
-    startcoords: $CS_STARTCOORDS
-    endcoords: $CS_FINISHCOORDS
-
-  - operator: write.write_cube_to_nc
-    filename: intermediate/transect
+    startcoords: $START_COORDS
+    endcoords: $FINISH_COORDS
 
-collate:
-  - operator: read.read_cube
-    constraint:
-        operator: constraints.generate_var_constraint
-        varname: $CS_VAR
-
-    filename_pattern: intermediate/*
+  - operator: plot.spatial_contour_plot
 
   - operator: write.write_cube_to_nc
     overwrite: True
-
-  - operator: plot.spatial_contour_plot