Merge pull request #731 from MetOffice/transect_fixes

Fix bug in transect operator

cset-workflow/meta/rose-meta.conf

@@ -419,8 +419,8 @@ help=This functionality extracts data for each variable defined in CS_VARS,
      such as transects through valleys, or through a front as it passes through.
 type=python_boolean
 trigger=template variables=CS_STARTCOORDS: True;
-trigger=template variables=CS_FINISHCOORDS: True;
-trigger=template variables=CS_VARS: True;
+        template variables=CS_FINISHCOORDS: True;
+        template variables=CS_VARS: True;
 compulsory=true
 sort-key=cs1
 

cset-workflow/rose-suite.conf.example

@@ -21,6 +21,9 @@ CSET_INPUT_FILE_PATH=""
 CSET_LOCAL_CSET_PATH="$HOME/CSET"
 CSET_RUNAHEAD_LIMIT=10
 !!CSET_TIMES_PER_FILE=0
+CS_FINISHCOORDS=""
+CS_STARTCOORDS=""
+CS_VARS=""
 DB_LONG_JOB=False
 DETERMINISTIC_PLOT_AGGREGATE_PRECIPITATION=False
 DETERMINISTIC_PLOT_CAPE_RATIO=False
@@ -34,6 +37,7 @@ DOMAIN_MEAN_VERTICAL_PROFILE_MODELLEVEL_SERIES=False
 DOMAIN_MEAN_VERTICAL_PROFILE_SERIES=False
 DOMAIN_SURFACE_HISTOGRAM_SERIES_FIELD=False
 ENSEMBLE_PLOT_SURFACE_AIR_TEMP=False
+EXTRACT_TRANSECT=False
 FETCH_FCST_OPT_CONF="filesystem"
 HOUSEKEEPING_MODE=2
 LFRIC_DOMAIN_MEAN_SURFACE_TIME_SERIES=False

src/CSET/operators/transect.py

@@ -94,8 +94,8 @@ def calc_transect(cube: iris.cube.Cube, startcoords: tuple, endcoords: tuple):
     )
 
     # Compute minimum gap between x/y spatial coords.
-    lon_min = np.min(lon_coord.points[1:] - lon_coord.points[:-1])
-    lat_min = np.min(lat_coord.points[1:] - lat_coord.points[:-1])
+    lon_min = np.abs(np.min(lon_coord.points[1:] - lon_coord.points[:-1]))
+    lat_min = np.abs(np.min(lat_coord.points[1:] - lat_coord.points[:-1]))
 
     # For scenarios where coord is at 90 degree to the grid (i.e. no
     # latitude/longitude change). Only xmin or ymin will be zero, not both